Fatal编程技术网
  • r/
  • R can'；找不到输入结果的正确方法？

R can'；找不到输入结果的正确方法？

r

R can'；找不到输入结果的正确方法？,r,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)/sigm

在下面的程序中，我找不到正确的方法来得到答案。我的意思是，找不到正确的参数值：

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)

f.sigma请阅读以下内容：关于第二个问题：查看
？uniroot
，其中
值是：
至少包含四个组件的列表…
。因此
uniroot
返回一个列表，您可以使用
uniroot（…）$root
访问
root
元素。。。。对于您的编辑：它与此相同：
burrx.mle1（x=c（4,1,5），lower=1，upper=5）
。但是c（4,1,5）是什么意思呢？我想您需要开始更多地阅读关于R的内容。
c（1,2,3）
是一个向量。