我很抱歉在这里问这个....但我确实需要你的帮助。
我正在创建一个表单。代码是:
laglist <- list(list(fun="ns", degree=2,int=TRUE),
list(fun="ns", degree=3,int=TRUE),
list(fun="bs",degree=2,int=TRUE),
list(fun="bs",degree=3, int=TRUE))
laglist <- rep(laglist,2)
varlist <- list(list(fun="ns", degree=2, knots=c(52.19, 82.51,237.02),int=TRUE),
list(fun="ns", degree=3, knots=c(52.19, 82.51,237.02),int=TRUE),
list(fun="ns",degree=4,knots=c(52.19, 82.51,237.02),int=TRUE))
varlist <- varlist[c(rep(1,4),3,2,3,3)]
for(i in 1:4) {
tab2a[i,1:2] <- c(argvar=varlist[[i]],arglag=laglist[[i]])
cb1 <- crossbasis(AFH6,lag=c(2,24), argvar=varlist[[i]],arglag=laglist[[i]])
tab2a[i,3] <- ncol(cb1)
model<- lm(NH3cH6~cb1+ns(DenH6,15)+ns(HenageH6,11)+as.factor(LightH6)+as.factor(Time))
tab2a[i,4] <- formatC(AIC(model),digits=1,format="f")
tab2a[i,5] <- formatC(BIC(model),digits=1,format="f") }
但总有人说:“tab2a [i,1:2]&lt; - c(argvar = varlist [[i]],arglag = laglist [[i]]”中的错误“要替换的项目数是不是替换长度的倍数。
我不知道如何修复它。希望你能帮我!非常感谢!!