如何获得R的95%置信区间？

我想在函数中为我的参数找到MLE的95%CI，但我不知道如何找到

给定的函数是一个具有

f（x）=Cx^（-mu）

我使用R中的bbmle包计算了mu的MLE

互联网上的一些人说使用profile-likelion可以做到这一点，但我不确定如何在R中使用profile-likelike，或者使用其他方法也可以得到同样的结果

非常感谢

更新：

load("fakedata500.Rda")
> library(stats4)
> library(bbmle)
> x<-fakedata500
> pl <- function(u){-length(x)*log(u-1)-length(x)*(u-1)*log(min(x))+u*sum(log(x))}

mle1<-mle2(pl, start=list(u=2), data=list(x))
> summary(mle1)
Maximum likelihood estimation

Call:
mle2(minuslogl = pl, start = list(u = 2), data = list(x))

Coefficients:
Estimate Std. Error z value     Pr(z)    
u  2.00510    0.04495  44.608 < 2.2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

-2 log L: 1300.166 

load（“fakedata500.Rda”）
>图书馆（stats4）
>图书馆（bbmle）
>x pl（从注释转换而来。）

如果您使用的是
bbmle
软件包中的
mle2
，您只需说出
confint（mle1）
即可获得95%的配置文件置信区间。请参阅
？confint.mle2
，或尝试
vignette（“mle2”，package=“bbmle”）
并搜索“confint”以获取更多信息。
请提供包括数据和代码在内的信息。@Thomas我已使用一些R代码更新了我的问题。希望这是有意义的。如果你正在使用
mle2
你应该只需要说
confint（mle1）
就可以得到95%的资料置信区间。PS尝试
vignette（“mle2”，package=“bbmle”）
并搜索“confint”…@BenBolker成功了！谢谢！！！嘿，本！你知道最近有没有关于如何找到MLL函数的置信区间的更新？这些方法似乎已经过时，我已经搜索了几个月，试图找到为MLL生成CI的方法。如果你今天有时间，如果你能看看我的问题，我将不胜感激