关于AIC值计算的求助！！！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我有以下数据
weight              day
38                    17
78                     42
44                     20
（数据略）

要拟合模型

Richards

W=A/(1+Bexp(-Kt))M

W:weight       t:day
想看他的拟合优度，AIC,BIC值的计算！谢谢各位大侠！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：AIC Richard Weight Eight weigh 求助 AIC

你先求出似然值，然后根据公式就很容易得到所求的aic和bic

您好！

x <- c(0.02, 0.02, 0.06, 0.06, 0.11, 0.11, 0.22, 0.22, 0.56, 0.56,1.10, 1.10)

y <- c(76, 47, 97, 107, 123, 139, 159, 152, 191, 201, 207, 200)
fn <- function(p) sum((y - (p[1] * x)/(p[2] + x))^2)

plot(x, y)

xfit <- seq(.02, 1.1, .05)

yfit <- 200 * xfit/(0.1 + xfit)

lines(spline(xfit, yfit))
out <- nlm(fn, p = c(200, 0.1), hessian = TRUE)
sqrt(diag(2*out$minimum/(length(y) - 2) * solve(out$hessian)))

plot(x, y)

xfit <- seq(.02, 1.1, .05)

yfit <- 212.68384222 * xfit/(0.06412146 + xfit)

lines(spline(xfit, yfit))
以上是我要拟合的模型及其代码。
我不清楚AIC值要什么语句才能出来？
望指点

先拟合，再直接调用stats包的AIC（）函数
比如：
fit<-lm(formula,data)
AIC(fit)