私は 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)
}
最後の行に問題があると思いますが、苦労した後はわかりません。
私は本当にあなたの助けに感謝します!
シューレ