我在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减半因素之外,还有一些问题。有没有人知道可能导致此错误的原因?