-
0 votesanswersviews
Rollapply用于滚动回归，具有不同的窗口大小，保存矩阵中的系数
我正在尝试执行滚动线性回归，重叠窗口以及不同的窗口大小 . 我想将结果保存在两个表中，一个是窗口大小和斜率，一个是窗口大小a和截取 . 我的数据集现在保存为xts对象 . 我正在尝试使用for循环生成不同的窗口大小，然后在两个变量之间进行滚动回归，使用rollapply重叠窗口 . 然后我希望每个回归的系数被保存为矩阵中的一行，可以构建，以便每个窗口大小对应一行 . 我看过很多以前的问题，但我无法... -
0 votesanswersviews
在SQL Server R服务中使用lm-model进行预测时，“因素具有新级别”错误
我是R的新手，但试图将它与SQL Server R服务一起使用 . 作为R脚本的输入，我选择几个数字列和一个带有分类数据的文本列 . 让它成为“company_name” . 然后我执行： declare @script_r nvarchar(max) = ' dat <- InputDataSet; model <- lm(first_commission_received... -
2 votesanswersviews
plot.lm错误：$运算符对原子向量无效
我有以下回归模型与转换： fit <- lm( I(NewValue ^ (1 / 3)) ~ I(CurrentValue ^ (1 / 3)) + Age + Type - 1, data = dataReg) plot(fit) 但是 plot 给了我以下错误： Error: $ operator is invalid for atomic vectors 关... -
0 votesanswersviews
解释R中的回归系数[关闭]
我正在尝试将x * log（x）模型拟合到数据中 . 拟合成功完成但我在解释结果系数时遇到困难 . 这是我的代码的快照 . x <- c(6, 11, 16, 21, 26, 31, 36, 41, 46, 51) y <- c(5.485, 6.992, 7.447, 8.134, 8.524, 8.985, 9.271, 9.647, 10.561, 9.971) fit <... -
4 votesanswersviews
获取`lm（）`返回的“mlm”对象的回归系数的标准误差
我想对同一个回归量运行10次回归，然后拉出所有标准错误 without using a loop . depVars <- as.matrix(data[,1:10]) # multiple dependent variables regressor <- as.matrix([,11]) # independent variable allModels <- lm(depVa... -
27 votesanswersviews
从R中的lm中提取标准化系数
我为这个愚蠢的问题道歉......但我似乎无法找到一个简单的解决方案 我想从拟合的线性模型中提取标准化系数（在R中）必须有一个简单的方法或函数来做到这一点 . 你能告诉我它是什么吗？ 编辑（以下一些评论）：我应该提供有关我的问题的更多上下文信息 . 我正在为一群心理学家教授一个介绍性的R工作室 . 对于他们来说，没有能够获得标准化系数的线性模型就好像你根本没有运行模型一样（好吧，这有点夸张，但你明... -
2 votesanswersviews
最小二乘回归系数非线性函数的标准误差和置信区间
我在R中运行OLS回归，从中得到几个系数 . 这是代码的一部分： Attacks <- Treat.Terr.Dataset$Attacks[2:30] Attackslag <- Treat.Terr.Dataset$Attacks[1:29] TreatmentEffect <- Treat.Terr.Dataset$TreatmentEffect[2:30] Treatme... -
0 votesanswersviews
R：在函数内部将变量传递给lm
我想编写一个调用 lm 的函数，并使用 ggplot2 绘制带有回归线的散点图 . 来自here，这是我的代码： fun <- function(m, n, o, p) { library(ggplot2) data <- as.data.frame(read.table(file = m, header = T, dec = ".", sep = &quo... -
1 votesanswersviews
线性模型中的分类预测变量是否应该是正态分布的？ [关闭]
我在R中运行简单的线性模型（Y~X），其中我的预测变量是一个分类变量（0-10） . 然而，这个变量不是正常分布的，并且没有任何可用的转换技术是 Health 的（例如log，sq等），因为数据不是负面/正面偏斜，而是整个地方 . 我知道，对于lm，结果变量（Y）必须是正态分布的，但这对于预测变量也是必需的吗？如果是的话，任何有关如何做到这一点的建议都会受到欢迎 . 此外，正如我所看到的数据有两组... -
0 votesanswersviews
3路anova嵌套在r
我正在试图找出完全因子实验的模型 . 我有以下因素 Treatment Day Hour Subject ResponseVariable 测量10天，每天4个不同的时间点，测量2个不同的治疗，12个受试者）治疗1中的6个受试者和治疗中的6个不同的受试者2） 对于我测量的每一天：在治疗1中的6个受试者，在治疗2中的其他6个，在4个不同的时间点 .对于受试者，我有12个不同的受试者，... -
0 votesanswersviews
R：多个子样本回归的约束系数和误差方差[关闭]
我正在和R一起研究145个观察样本 . 我创建了五个子样本，每个子样本有29个观察值，而响应变量 q 已经排序 . 结果，subset1包含具有最低输出的29行数据帧，subset2包含以下29行等 . 我正在回归预测变量 x1 ， x2 ans x3 上的变量 q . 我现在需要进行两个实验： 约束所有子样本的误差方差相同; 约束 x2 和 x3 上的系数以及5个OLS回归中的误差... -
1 votesanswersviews
R中的多元回归：协变对因子的影响还是相反？
我知道对于具有统计背景的人来说，这可能是一个非常简单的问题 . 然而，我找不到适合我（不是那样）特殊情况的明确答案:( 我有一个回归模型，有两个分类预测因子（A，有两个级别A1和A2和B，有两个级别B1和B2），还有一个数字预测器Z. 我对Z和B之间的相互作用感兴趣，但特别是在A的一个级别（A1，我的参考级别） . 因此，我安装了以下型号： lm(Y ~ A/B*Z, data=df) 这让我想... -
3 votesanswersviews
R [重复]中的约束线性回归系数
这个问题在这里已有答案： R : constraining coefficients and error variance over multiple subsample regressions [closed] 1回答 我估计R中的几个普通最小二乘线性回归 . 我想约束回归中的估计系数，使它们相同 . 例如，我有以下内容： z1 ~ x + y z2 ~ x + y 我希望第一次回归中y的... -
10 votesanswersviews
添加拟合二次曲线
我正在尝试将拟合的二次曲线添加到绘图中 . abline(lm(data~factor+I(factor^2))) 显示的回归是线性的而不是二次的，我得到这样的信息： Message d'avis：在abline（lm（data~factor I（factor ^ 2）），col = palette [iteration]）：利用des deux premiers des 3 systemsd... -
3 votesanswersviews
线性回归中的残差是否遵循原始数据框行的相同顺序？
在我拟合多元线性回归模型之后，我想将残差与原始数据框合并 . 但是我发现残差是一个没有原始行名称附加到每个残差的向量 . 残差是否遵循原始数据框的行名称的相同顺序？有没有办法检查这封信件？ dat <- read.table("myData.txt", header = T) dat.str <- data.frame(dat$response, dat$v1, d... -
74 votesanswersviews
线性回归和R中的分组
我想使用 lm() 函数在R中进行线性回归 . 我的数据是一年一度的时间序列，一年（22年），另一个州（50个州） . 我想为每个状态拟合一个回归，以便最后我有一个lm响应的向量 . 我可以想象为每个状态做循环然后在循环内进行回归并将每个回归的结果添加到向量 . 然而，这似乎并不像R一样 . 在SAS中我会做一个'by'语句，在SQL中我会做一个'group by' . R的做法是什么？ -
1 votesanswersviews
R在数据上拟合多项式
我有一些从一个函数综合生成的数据，如下所示 . fn <- function(w1,w2){ f= -(0.1 + 1.3*w1 + 0.4*w2 - 1.8*w1*w1 - 1.8*w2*w2) return(f) } 接下来，我创建一个数据框，其值如下所示 x = data.frame( yval = fn(seq(0.1,0.9,by=0.01),seq(1.1,0.3,... -
8 votesanswersviews
predict.lm（）如何计算置信区间和预测区间？
我跑回了一个回归： CopierDataRegression <- lm(V1~V2, data=CopierData1) 我的任务是获得一个 90％ confidence interval 为 V2=6 和 V2=6 的平均响应当 V2=6 时 90％ prediction interval . 我使用了以下代码： X6 <- data.frame(V2=6) predict... -
1 votesanswersviews
R：纯固定效应模型中sigma的置信区间
是否有一种标准的方法来估计具有固定效应的线性模型的方差参数的置信区间 . 例如 . 给定： reg=lm(formula = 100/mpg ~ disp + hp + wt + am, data = mtcars) 如何获得方差参数的置信区间 . confint只详细说明了固定效果，lme4中的lmer不接受没有level-2随机效果的模型，这就是我的情况 . -
0 votesanswersviews
哪种排序方法会根据响应变量预测分组？
我有几个回归模型，格式为m < - lm（y~x，dat） 每个模型使用不同的y，但总是相同的x 例如 data(Iris) m <- lm(Iris$Sepal.L. ~ Sepal.W., data = Iris) 我从m中提取了X和Y. X <- m$call [2][[1]][[3]] Y <- m$call [2][[1]][[2]] 然后我尝试在ggp... -
4 votesanswersviews
使用lm（）和scale（）的标准化回归系数与使用lm.beta（）或cor（）的标准回归系数不同
我有两个变量，我想找到它们之间的相关性 . 问题是我似乎得到了不同的结果，这取决于我使用的方法 . 我所知道的一种方法是在scale（）函数中运行带有独立变量和因变量的lm（）函数 . 所以下面的变量看起来像： lm(scale(mainDataframe$relativeFemHappy) ~ scale(mainDataframe$allRights)) 我所知道的其他方法是简单地使用cor... -
1 votesanswersviews
R data.table不能与dplyr一起使用
继发帖后：Fitting several regression models with dplyr 尝试将其应用于data.table时，我遇到了一个问题 . 当数据是data.frame时，问题不存在 例： library(data.table) library(dplyr) mtcarsDT <- data.table(mtcars) lmGroups <- mtcarsDT ... -
0 votesanswersviews
协变员不在tslm工作
我有每周零售销售的季节性时间序列数据，我正在尝试使用Hyndman预测包的tslm函数，以适应除趋势和季节之外的回归量模型 . 我遇到的问题是，当我构建tslm时，在添加任何回归量（仅趋势季节）之前，我在训练数据上得到了完美的拟合（R ^ 2 = 1）！ 完美契合是有问题的，因为我添加到模型中的任何额外协变量（正在出售的物品数量，分布等）对预测没有影响（无关紧要） . 只看数据，我知道其他回归量很... -
4 votesanswersviews
如何在lm中指定参数估计值之间的关系？
使用lm，我想拟合模型：y = b0 b1 * x1 b2 * x2 b1 * b2 * x1 * x2 我的问题是：如何指定相互作用的系数应该等于系数乘以主效应？ 我已经看到将系数设置为特定值，您可以使用offset（）和I（），但我不知道如何指定系数之间的关系 . 这是一个简单的模拟数据集： n <- 50 # Sample size x1 <- rnorm(n, 1:n, 0.5... -
3 votesanswersviews
基于R中模型平均系数的部分残差图
我正在使用R包 MuMIn 进行多模型推理，并使用函数 model.avg 来平均由一组模型估计的系数 . 为了在视觉上将数据与基于平均系数的估计关系进行比较，我想使用部分残差图，类似于 car 函数的 crPlots 函数创建的图 . 我've tried three ways and I'我不确定是否合适 . 这是一个演示 . library(MuMIn) # Loading the data... -
-1 votesanswersviews
从许多线性模型创建自定义结果表[重复]
这个问题与以下内容完全相同： how to create many linear models at once and put the coefficients into a new matrix? 2个答案 我在一个分析实验室工作，生成大量数据 . 我们进行模型拟合，并对模型系数（截距，a，b），r_squared和残差标准误差感兴趣 . 最多70个因变量（响应）的数量相当大 . 我想从my... -
0 votesanswersviews
解释线性回归中分类预测变量的估计[重复]
这个问题在这里已有答案： Interpreting interactions in a regression model 2个答案 我是线性回归的新手，我正在试图弄清楚如何解释汇总结果 . 我难以解释分类预测因子的估计 . 请考虑以下示例 . 我添加了列的年龄和长度，以包括数字预测器和数字目标 . library(MASS) data <- as.data.frame(HairEyeCo... -
-1 votesanswersviews
描述性分析与斜率估计线性模型R之间的不匹配
我是一名处理R建模的高级学生 . 我正在尝试找到由n个重复行xm变量列组成的数据集的最佳模型：我想构建一个lm来解释4个分类回归量在Y（连续数据）植物梢数/平方米中的影响 . 公式模型为：lm（Y~a b c d） . 回归水平：“a”有4个等级（阴影百分比等级），“b”有4个等级（4个调查年份），“c”有3个等级（高程等级）和“d”有7个等级（7个空间多边形，其中拍摄样品） . 在描述性分析中，... -
0 votesanswersviews
R中的不同geom_smooth和lm（）估计：忘记将基坡和交互斜率加在一起
我试图找出为什么我的 lm() 估计值与 geom_smooth 的相同数据和公式不同 . 具体来说，我的分组变量"cat" level 5的斜率在 lm() 输出中> 0，但在geom_smooth中<0，因此该图似乎不反映汇总表 . 这是the data . （比提供行为相似的示例数据更容易 . ） 型号： summary(lm(data=df, y~x*ca... -
0 votesanswersviews
什么时候使用relevel或ANOVA post-hoc？
我在CrossValidated上问了this question，但我认为它仍然是一个R问题太多，所以我会看看是否有机会在这里得到答案 . 假设我有一个示例数据集： x <- rnorm(20, 5, 1) x <- append(x, rnorm(20, 8, 1)) x <- append(x, rnorm(20, 13, 1)) y <- rep(c("a&...