用optim进行优化
用optim进行优化
提问于 2015-05-17 15:57:48
我不知道下面的代码出了什么问题，在经历了很长时间的痛苦之后，可能一整天.我真的希望有人能提前找到problem...Thanks，....I不知道为什么R给了我相同的初始值
lambda<-1
likelihood<- function(lambda) {
T6<- exp(lambda*2.82)/(exp(lambda*2.82)+exp(lambda*0.28))
T5<- exp(lambda*2.69)/(exp(lambda*2.69)+exp(lambda*((T6*0.38)+((1-T6)*2.92))))
T4<- exp(lambda*2.52)/(exp(lambda*2.52)+exp(lambda*((T5*0.51)+((1-T5)*2.82))))
T3<- exp(lambda*2.30)/(exp(lambda*2.30)+exp(lambda*((T4*0.68)+((1-T4)*2.69))))
T2<- exp(lambda*2.00)/(exp(lambda*2.00)+exp(lambda*((T3*0.90)+((1-T3)*2.52))))
T1<- exp(lambda*1.60)/(exp(lambda*1.60)+exp(lambda*((T2*1.20)+((1-T2)*2.30))))
l6<-(T6^0)*((1-T6)^281)
l5<-(T5^0)*((1-T5)^281)
l4<-(T4^2)*((1-T4)^(281-2))
l3<-(T3^19)*((1-T3)^(281-19))
l2<-(T2^93) *((1-T2)^(281-93))
l1<-(T1^166)*((1-T1)^(281-166))
-(l1*l2*l3*l4*l5*l6)
}
optim(lambda,likelihood)$par回答 1
Stack Overflow用户
回答已采纳
发布于 2015-05-17 16:16:14
代码没有任何问题(就编程错误而言)。如果您认为它不应该返回上限，那么您可能在代码中函数的定义上犯了错误。
首先，当您只有一个参数来估计您的方法时，应该是Brent，如下所示：
> optim(lambda,likelihood, method='Brent', lower=1.5, upper=200)$par
[1] 200它总是将上限显示为正确的上限的原因是，您的函数只定义在-100 - 300之间(大约)，并且它的值总是为零，如下所示(不管lambda的值如何)。当这种情况发生时，optim总是将上限显示为正确的值。您可以在下面看到您的函数：
plot(likelihood, from=-100, to=1000, xlab='lambda')
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/30288905
复制相关文章
相似问题
腾讯云开发者
