我正在尝试估算R包nlsLM()
中minpack.lm
的复杂模型。它是这样的:
nlc <- nls.lm.control(maxiter = 10000, maxfev = 10000)
nlsLM(formula = log(demand+0.001) ~ Diff.model.log(m,pin,q,lambda1,lambda2,goal,length,time),
data = projects_long,
start = c(lambda1 = lambda1, lambda2 = lambda2, m = m, pin = pin, q = q),
upper = c(5,5,Inf,1,1), lower = c(-Inf,0,0,0,0),
algorithm = "port", control = nlc)
你可以看到我在我试图估计的五个参数中的四个上放置了下限0
,它们似乎确实具有约束力:我所获得的所有参数集都包括零。 ..
但事实上,在我的应用程序中,参数为零是不可解释的,而且我认为它们是严格正面的。我想知道如何在代码中指定绑定不可达?
非常感谢您的指导!