我正在尝试在R中运行一系列非线性回归。我试图拟合的模型描述了病毒群体在三个不同温度下的时间进度。我正在尝试将函数nlsList()与groupedData对象一起使用:
surv <- groupedData(LOG10N ~ t | Temp, data.frame(cbind(rbind(surv4, surv22,surv56),
Temp=rep(c(4, 22, 56),
each=8))))
其中LOG10N是总体的log10,t是时间,温度是温度。使用nlstools包我试图适应mafart模型,所以我使用mafart模型公式来创建自我启动:
mafart.self <- selfStart(~ LOG10N0-(t/delta)^p, surv,
parameters=c('p', 'delta', 'LOG10N0'))
getInitial(LOG10N~mafart.self, data=surv56)
当我尝试获取初始值时,我收到以下错误:
Error in object[[3L]][[1L]] : object of type 'symbol' is not subsettable
如果我跳过getInitial()步骤并使用如下的nlsList()函数,我会得到同样的错误:
> curve.maf <- nlsList(LOG10N~mafart.self, surv)
Error in model[[3]][[1]] : object of type 'symbol' is not subsettable
我尝试使用nls2()分别运行回归,但结果完全令人失望,因为它们返回的值与初学者相同。