数据文件：
下面是我自己编的程序，用极大似然估计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))