R软件极大似然估计maxlik计算的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

代码;
ll<-function(theta){
  
  a <- theta[1]
  b <- theta[2]

  ll<-N*log(a)+sum(log(exp(-b*x2)-exp(-b*x1)))-a*(1-exp(-b*337))
  return(-ll)

}
x1 <- c(9,21,32,36,43,45,50,58,63,70,71,77,78,87,91,92,95,98,104,105,116,149,156,247,249,250,337)
x2 <- c(0,9,21,32,36,43,45,50,58,63,70,71,77,78,87,91,92,95,98,104,105,116,149,156,247,249,250)
y <- 0:26
z <- 1:27
N <- length(x1);
res<-nlminb(c(1,5),ll,lower=c(0,0),upper=c(100,100))
print(res)
遇到的问题:
$par
[1] 1 5

$objective
[1] Inf

$convergence
[1] 0

$iterations
[1] 1

$evaluations
function gradient
       2        2

$message
[1] "X-convergence (3)"

Warning message:
In nlminb(c(1, 5), ll, lower = c(0, 0), upper = c(100, 100)) :
  NA/NaN function evaluation
> 麻烦帮忙解决一些

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：MAXLIK 极大似然估计 似然估计 极大似然 max 软件

初始值选的不好，使得ll等于无穷大

qoiqpwqr 发表于 2015-10-31 06:34
初始值选的不好，使得ll等于无穷大

如何选取合适的初始值呢？