我是使用 R 的新手,我不熟悉编写 R 循环。我需要你帮助我完成以下 R 编程。R 代码的目标是为所有具有不同 df 的模型获取 AIC。
我写的 R 代码是:
library(splines)
library(dlnm)
For (i in 1:30)
argvar1 <-list(type="bs", df=2[i],cen=50)
arglag1<- list(type="ns",df=3)
cb1 <-crossbasis(AFH6w,lag=24,argvar=argvar1,arglag=arglag1)
argvar2 <-list(type="ns", df=11, cen=-2)
arglag2 <- list(type="ns",df=3)
cb2 <-crossbasis(OutTw,lag=24,argvar=argvar2,arglag=arglag2)
model1 <-lm(NH3cH6w~ cb1 +cb2+DenH6w+as.factor(Month))
AIC[i]=AIC(model1)
}
我认为最后一行有问题,但经过疲惫的挣扎后我无法弄清楚。
我将非常感谢您的帮助!
舒乐