我正在使用以下循环以适应多个模型:
lambda = seq(from=0,to=10.2,by=0.2)
for(i in lambda){
assign(paste("fitb",i, sep=""),
lmmlasso(y = mydata.B$Bloodpress,
x = x.matrix, z = z, grp = grp1,
lambda = i, pdMat = "pdIdent"))
}
所以我最终在名单中列出了很多模型:fitb0,fitb0.2等等
我需要从每个列表中提取AIC值并将它们放在一个向量中,而不是使用代码:
AIC=c(fitb0$aic, fitb0.2$aic ... fitb10$aic)
我该怎么办?
答案 0 :(得分:0)
也许在循环中抓住aic?
lambda = seq(from=0,to=10.2,by=0.2)
i=1 ; aic=NA
for(i in lambda){
wholelassoresults=lmmlasso(y = mydata.B$Bloodpress,
x = x.matrix, z = z, grp = grp1,
lambda = i, pdMat = "pdIdent"))
assign(paste("fitb",i, sep=""),wholelassoresults)
aic[i]=wholelassoresults$aic
}