我有一个问题来计算AIC.实际上,我估计了我的3个模型的参数:"mod_linear",它是线性模型,"mod_exp"和"mod_logis"是两个非线性模型.
我使用了函数AIC():
AIC(mod_linear,mod_exp,mod_logis)
df AIC
mod_linear 4 3.015378
mod_exp 5 -11.010469
mod_logis 5 54.015746
但我试图用公式AIC = 2k + nlog(RSS/n)来计算AIC,其中K是参数的数量,n是样本的数量,RSS是剩余的平方和.
k=4
n=21
#Calcul of nls for the linear model:
mod_linear=nls(data$P~P_linear(P0,K0,a),data=data,
start=c(P0=4.2,K0=4.5,a=0.)
2*k+n*log(sum(residuals(mod_linear)^2)/n)
-56.58004
正如您所看到的,结果并不相同,而其他两个模型也是如此.有人可以帮帮我吗？
问候