r - 在 lm() 调用中更新数据

Question

是否有与 lm 调用对象的数据部分的更新等效的方法？

例如，假设我有以下模型:

dd = data.frame(y=rnorm(100),x1=rnorm(100))
Model_all <- lm(formula = y ~ x1, data = dd)

有没有一种方法可以对 lm 对象进行操作以产生等效的效果:

Model_1t50 <- lm(formula = y ~ x1, data = dd[1:50,])

我正在尝试从样本预测测试中构建一些 psudo，拥有一个 lm 对象并简单地滚动数据会非常方便。

Answer 1

我相当确定 update 实际上会做您想要的!

example(lm)
dat1 <- data.frame(group,weight)
lm1 <- lm(weight ~ group, data=dat1)
dat2 <- data.frame(group,weight=2*weight)
lm2 <- update(lm1,data=dat2)
coef(lm1)
##(Intercept)    groupTrt 
##      5.032      -0.371 
coef(lm2)
## (Intercept)    groupTrt 
##     10.064      -0.742 

如果您希望从中获得效率提升，您会感到失望——R 只是替换新参数并重新评估调用(参见 update.default 的代码) .但它确实使代码更清晰......