我正在尝试使用MuMIn软件包手动计算模型的AICc（校正的AIC）。MuMIn软件包中的AICc（模型）是否与下面的aicc相同？
model <- lm(Y~., df) #df consists of input vectors and Y.
n <- nrow(model$model)
k <- length(model$coefficients)+1
aicc <- -2*stats::logLik(model)[1]+2*k+2*k*(k+1)/(n-k-1)字符串
我尝试做这个比较，因为我的代码可以在R 3.6.3下运行，而MuMIn包不再支持R 3.6.3。很难让我的代码兼容MuMIn（R>4.2.0）。
1条答案
按热度按时间cpjpxq1n1#
似乎是这样：
字符串
更一般地说，您可以使用
n <- nobs(model); k <- length(coef(model))+1;大多数模型类型都有这些访问器方法。