我想估计几乎是逻辑的加权模型,但需要估计饱和度(S)。目前我正在使用通用的非线性优化方法,如下所示:
nls.model <- y ~ S/(1 + exp(-1 * (b1*x1+b2*x2+b3*x3)))
...
nlsGR.out <- nls(nls.model, data=x, start=nls.init.values, trace=F, weight=weight, control=nls.ctrl)
其中S是我感兴趣的饱和度.y是范围内的双倍[0.0,1.0]
我想知道我可以使用更高效,更健壮或恰当的套餐例程吗?