在以下程序中,我无法找到正确的方法来获得答案。我的意思是,找不到合适的参数值:
f.sigma <- function(sigma, x)
{
n <- length(x)
sum1 <- sum(x^2)
sum2 <- sum(log(1-exp(-1*(x/sigma)^2)))
sum3 <- sum(x^2/(exp((x/sigma)^2)-1))
ans <- -1*n + sum1/sigma^2 + (n/sum2+1)/sigma^2*sum3
return(ans)
}
burrx.mle1 <- function(x, lower, upper)
{
n <- length(x)
sigma.root <- uniroot(f.sigma, interval=c(lower,upper), x=x)
sigmahat <- **sigma.root$root**
thetahat <- -1*n/sum(log(1-exp(-1*(x/sigmahat)^2)))
ans <- list(theta.hat=thetahat, sigma.hat=sigmahat)
return(ans)
}
burrx.mle1(c(4,1,5),1,5)
另一件事是什么意思我的sigma.root $ root?基本上什么是$ root的用途? 在burrx.mle1(c(4,1,5),1,5)这是什么意思我的c(4,1,5)?