Optim()函数错误
Optim()函数错误
提问于 2014-02-15 04:49:05
这是我的代码。
beta1 = function(a,b,t) { beta(a+(1/t),b) }
beta2 = function(a,b,t) { beta(a+(2/t),b) }
eb11 = function(a,b,t) { beta2(a,b,t)/beta(a,b) }
eb12 = function(a,b,t) { (beta1(a,b,t)-beta2(a,b,t))/beta(a,b) }
eb22 = function(a,b,t) { 1 + (beta2(a,b,t)-2*beta1(a,b,t))/beta(a,b) }
eb11r11 = function(a,b,t) { beta2(a,b,t)*beta(a,b)/beta1(a,b,t)^2 }
eb12r12 = function(a,b,t) { (beta1(a,b,t)-beta2(a,b,t))*beta(a,b)/beta1(a,b,t)/(beta(a,b)-beta1(a,b,t)) }
eb22r22 = function(a,b,t) { (beta(a,b)^2 + (beta2(a,b,t)-2*beta1(a,b,t))*beta(a,b))/(beta(a,b)-beta1(a,b,t))^2 }
gbetloglik = function(a,b,t) {
loglik = n1*log(eb11r11(a,b,t)) + n2*log(eb12r12(a,b,t)) + n3*log(eb22r22(a,b,t))
return(-loglik)
}
abt = optim(c(0.5,0.5,1),gbetloglik,lower=c(0.001,0.001,0.001),method="L-BFGS-B")$par我想要做的是找到a，b和t来最大化'gbetloglik‘函数。但是我得到了这个错误。
Error in 2/t : 't' is missing似乎缺少函数'beta2‘的第三个参数。当我在gbetloglik函数中直接输入三个数字时，它工作得很好。这个问题只出现在optim()函数中。有谁知道吗？
回答 1
Stack Overflow用户
发布于 2014-02-15 05:12:51
看起来您曲解了optim函数的第一个参数。第一个参数只是为正在优化的1个参数提供初始值。在您的示例中，这为gbetloglik的一个参数提供了3个初始猜测。此调用将起作用：
abt = optim(0.5,gbetloglik,lower=c(0.001,0.001,0.001),method="L-BFGS-B", b=0.5, t= 0.5)$par但不会在所有三个参数上都进行优化，它只会在给定b和t的情况下优化a。要优化所有参数，您需要从here安装一个外部包。下面是来自nlmrt的一个示例
ydat = c(6.308, 6.94, 9.638, 12.866, 17.069, 23.192, 31.443, 37.558, 51.156, 64.948, 77.995, 91.972)
tdat = seq_along(ydat)
start1 = c(b1=1, b2=1, b3=1)
eunsc = y ~ b1/(1+b2*exp(-b3*tt))
anlxb1g =try(nlxb(eunsc, start=start1, trace=FALSE, data=data.frame(y=ydat, tt=tdat)))
print(anlxb1g)
anlxb1g$coefficients 页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/21789261
复制相关文章
相似问题
腾讯云开发者
