需要帮助R编码

时间:2013-09-30 20:28:36

标签: r

我是使用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)
}

我认为最后一行存在问题,但在经历了疲惫的挣扎之后,我无法弄明白。

我真的很感谢你的帮助!

舒乐

1 个答案:

答案 0 :(得分:1)

您缺少大括号括号,for应为小写:

for (i in 1:30) {
        argvar1 <-list(type="bs", df=2[i],cen=50)
        arglag1<- list(type="ns",df=3)
        ...
}

除此之外很难说。请在发布前阅读How to make a great R reproducible example?