无法找到输入结果的正确方法?

时间:2015-02-17 03:32:21

标签: r

在以下程序中,我无法找到正确的方法来获得答案。我的意思是,找不到合适的参数值:

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)?

0 个答案:

没有答案