将滞后变量添加到lm模型中？
我在一个时间序列上使用了lm，实际上效果非常好，并且速度超级超级快。 假设我的模型是：> formula <- y ~ x 我在训练集上训练它：> train <- data.frame( x = seq(1,3), y = c(2,1,4) ) > model...
如何从lm结果中获得RMSE？
我知道$sigma和均方根误差（RMSE）的概念之间有一点小差异。那么，我想知道在R中从lm函数中获取RMSE的最简单方法是什么？res<-lm(randomData$price ~randomData$carat+ randomData$cut...
当在lm（）函数中编写回归公式时，R中的“/”表示交互作用。
公式就像这样。我不是很理解符号“/”的用法。它似乎通常用于虚拟变量。但是我对它的用法不确定。lm(y~x/z,data = data.frame(x = rnorm(6), y = rnorm(6), z = rep(0:1,each=3)))
在特定范围内绘制拟合线R。
使用R语言，我想画出两个变量之间的线性关系图，但是我希望拟合直线仅出现在数据范围内。 例如，如果我的代码如下，我希望该直线仅存在于x和y值为1:10的范围内（使用默认参数，该直线超出数据点的范围）。x <- 1:10 y <- 1:10 plot(x,y) abline(lm(y...
在lm()中使用列编号而非列名
我希望可以通过列编号而不是名称来指定模型的列，例如lm(bp~height+age, data=mydata)。 我尝试过lm(mydata[[1]]~mydata[[2]]+mydata[[3]])但问题在于：拟合模型中的系数被命名为mydata[[2]]，mydata[[3]]等，而我想...
使用lm()对象列表进行预测
我有一些数据，我经常运行回归分析。每个“块”数据都会配合不同的回归方程。例如每个州可能有一个不同的函数来解释因变量。这似乎是一个典型的“拆分-应用-组合”问题，因此我正在使用 plyr 包。我可以轻松创建 lm() 对象的列表，这很有效。但我还不太确定如何稍后使用这些对象来预测另一个数据框中的...
plot.lm(): 提取诊断 Q-Q 图中标记的数字
对于下面的简单示例，您可以看到在随后的图中有某些点被识别出来。我如何提取在这些图中被标识的行号，特别是正态 Q-Q 图？set.seed(2016) maya <- data.frame(rnorm(100)) names(maya)[1] <- "a" maya$b <- ...
使用BIC准则运行逐步线性模型
能否将逐步线性模型设置为使用BIC准则而不是AIC？我一直在尝试，但它仍然使用AIC值计算每个步骤。null = lm(data[,1] ~ 1) full = lm(data[,1] ~ age + bmi + gender + group) step(null, scope = list(...
我该如何使用ggplot绘制lm()的残差图？
我希望能够得到一个关于 lm() 模型残差的漂亮图形。目前我使用 plot(model$residuals)，但我想要更好看的图形。如果我尝试使用 ggplot 绘制，就会收到以下错误信息： ggplot2 不知道如何处理数字类别的数据