我在 nlme 中运行一个自启动的逻辑模型,并收到了经典的错误消息:
Error in nlme.formula(model = MASS.CUL ~ SSlogis(GD10C, Asym, xmid, scal), :
step halving factor reduced below minimum in PNLS step
我首先输入代码:
shortG <- groupedData(MASS.CUL ~ GD10C | Expu2, data = short)
nlsL <- nlsList(MASS.CUL ~ SSlogis(GD10C, Asym, xmid, scal) , data = shortG)
shortG.nlme <- nlme(nlsL, random=pdDiag(Asym + xmid + scal ~ 1))
一切正常。然后我去用代码标准化方差:
shortG.nlmeb <- update(shortG.nlme, weights = varPower())
这是我收到错误消息的地方。现在,通常我可以通过简单地增加 pnlsTol 来在圆孔中安装一个方形钉。但是,在将其首先增加到 0.01,然后增加到 0.1,然后增加到 1,然后增加到 10,然后增加到 10,000,只是为了好玩,我不断收到错误消息。这让我相信除了 pnls 减半因素之外还有一些问题。有没有人知道可能导致此错误的原因?