如何获得gev.fit
中参数最大似然估计的95%置信区间? R中的ismev
包。
gev.fit
的输出:
$conv
[1] 0
$nllh
[1] 194.6467
$mle
[1] 93.33096665 17.53896061 -0.02969443
$se
[1] 2.92704621 2.06566895 0.09614705
答案 0 :(得分:1)
如果您可以做出一个连贯的论证,即$ mle值通常以$ se中的标准错误分布,那么您可以获得95%CI:
LowerBounds <- gev.fit$mle - 1.96*gev.fit$se
UpperBounds <- gev.fit$mle + 1.96*gev.fit$se