我希望您能够将每个lm对象的系数和标准误差提取出来,并将它们组合成一个数据框,缺失的预测变量用NA填充。
任何建议将不胜感激!
set.seed(12345)
x<-matrix(rnorm(1000),nrow=100,ncol=10)
colnames(x)<-paste("x",1:10,sep="")
df<-data.frame(y=rnorm(100),x)
m1<-vector('list', 10)
for ( i in 2:11){
eqn <- as.formula(paste("y ~", paste(colnames(df)[2:i], collapse='+')))
m1[[i-1]] <- lm(eqn, df)
}
任何建议将不胜感激!
m1
的方式引起的。(你第一个模型拟合的结果被放入了m1[2]
,导致m1[1]
的值为NULL
。)我修改了你的问题以纠正这个问题(通过在倒数第二行中分配到m1[i-1]
),现在我的所有代码都应该可以正常工作了。 - Josh O'Brien