r - 95% CI 用于参数估计

Question

如何获得参数最大似然估计的95%置信区间gev.fit？ismev在 R 中的包。

的输出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

Answer 1

如果您可以提出一个连贯的论点，即 $mle 值与 $se 中的标准误差呈正态分布，那么您可以获得 95% 的 CI：

 LowerBounds <- gev.fit$mle - 1.96*gev.fit$se
 UpperBounds <- gev.fit$mle + 1.96*gev.fit$se