R语言 t分布参数的极大似然估计问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

数据文件： glass fibre.txt (326 Bytes)

2015-3-21 15:19:45 上传

下面是我自己编的程序，用极大似然估计t分布的参数：
问题是：运行之后，产生很多warnings, 比如：

1: In fnOrig(theta, ...) : 'gammafn'里的值在范围外
……
29: In log(v * pi * sigma^2) : 产生了NaNs……等等等等。
虽然能够得到参数的估计值，但与数据严重不符，而在被人的文献中，这些数据是比较符合的。

这是什么问题呀？哪里出错了？
求助


x <- scan("glass fibre.txt")

loglik2 <- function(par)
{
        v <- par[1]
        mu <- par[2]
        sigma <- par[3]
        t <- (x-mu)/sigma

        ll <- n*log( gamma((v+1)/2)) - n*log( gamma(v/2) ) - n*0.5*log(v*pi*sigma^2)  - 0.5*(v+1)* sum( log( 1+ ((x-mu)/sigma)^2/v ) )
       
      ll
}
res2 <- maxLik(loglik2, start=c(2,0,1))

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：极大似然估计 似然估计 极大似然 R语言 t分布 sigma 程序

请问楼主这个问题解决了吗，我也遇到了相同的问题，十分困惑