如何加入保存类的glm对象("glm“"lm")

Question

我试图估计一些logistic回归模型的Nagelkerke R^2。因为我有很多模型，所以我尝试用for自动化这个过程。但是，当我将glm模型加入到列表中时，来自package fmsb的NagelkerkeR2()函数给了我一个错误。

##Calculate Nagelkerke's R Squared
x<- list(m01,m02,m03,m04,m05,m06,m07,m08,m09,m10)
Nag_R2 <- NA

for (i in 1:10){
   NagelkerkeR2(x[1])[2] <- Nag_R2 
}

之所以出现此错误，是因为该函数只处理"glm“"lm”对象。那么，是否有一种方法将glm对象作为单个"glm“"lm”对象加入？

Answer 1

我想你是在找这样的东西：

Nag_R2 <- numeric(length(x))  ## pre-allocate vector
for (i in 1:10) {
   Nag_R2[i] <- NagelkerkeR2(x[[i]])[2]
}

这是假设NagelkerkeR2返回的第二个组件是要保存的数值(我还没有检查)。

请注意，如果您希望返回一个对象(而不是length-1 list )，则需要使用[[]]访问列表的元素：

​

​