如何获得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

===============>>#1 票数:1 已采纳

如果您可以提出一个连贯的论点,即$ mle值在$ se中正态分布且带有标准错误,则可以通过以下方式获得95%的CI

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

  ask by Wee Wang Wang translate from so

未解决问题?本站智能推荐:

2回复

R仿真中的参数估计

我对R相当陌生,并且正在探索仿真以估计参数n(整数) 1)Z是从N(0,1)开始的n个绘制的向量 2)max(Z)> 4的概率等于0.25 R中估计参数n以满足这两个条件的最佳方法是什么? 尝试避免代码中的循环或穷举搜索时,我陷入了困境。 谢谢!
2回复

lm参数估计

错误变量长度不同 我对此错误感到困惑,而且我不知道该怎么办。 编辑:我正在这种情况下,其中n = 20,参数b0 = 0,且b = 1为true,且独立和误差均值均值为0和sd = 1。 这可能吗? 非常感谢!
1回复

R中非线性优化的自举参数估计:为什么与常规参数估计不同?

这是我的问题的简短版本。 代码如下。 我使用optim()计算了R中非线性von Bertalanffy生长方程的参数,现在我试图通过自举将95%置信区间添加到von B生长系数K中。 对于至少一年的数据,当我总结增长系数K的自举输出时,自举的均值和中值参数估计与估计的参数完全不同:
1回复

威布尔分布参数估计误差

我使用以下函数来估计三参数威布尔分布。 但是,当我运行上述功能时,出现以下错误。 有什么办法可以克服这个问题? 谢谢!
1回复

参数估计中的optim()函数错误

以下代码是找到广义beta分布的3个参数。 然后,它给出如下错误并带有许多警告。 警告消息是 n1,n2和n3由某些随机数生成方案生成,并且(n1,n2,n3)决定是否发生错误。 (对于某些(n1,n2,n3)不会发生此错误) 我不知道此错误的原因是什么,“ be
1回复

初始参数估计时的奇异梯度矩阵

我将Lotka-Volterra模型拟合到一些数据,而我的简单代码是: 但它返回: 奇异梯度” 如果我从公式中删除乘以a1的“ PL”,那么它将起作用。 我不能跳过它,请帮助
1回复

另一:“在初始参数估计单数梯度矩阵”

我正在尝试使用nls函数拟合Gompertz曲线,这总是会给我“奇异梯度”错误。 我知道有很多问题专门针对此错误,但是我仍然无法弄清我的问题。 我尝试了多种公式变体,起始值,nlsLM函数,...我感觉好像缺少了一些非常简单的东西,但我一直在搜索数小时。 这是我的数据和功能: 感谢
1回复

如何在R中对参数估计施加约束?

我正在使用mlogit软件包在R上拟合logit模型。 虚拟变量的所有参数(var1,var2,var3)均呈正态分布。 如何对一个变量(var1)施加约束以使其参数为零均值? 当前代码为var1返回非零均值。 谢谢! 我知道在STATA中对参数施加约束相对容易。 对于我的模型
1回复

初始参数估计时的nls奇异梯度矩阵

我试图将第一个零阶动力学曲线拟合到某些数据,当我绘制数据并将一些初始参数添加到模型时,我得到了很好的视觉拟合,但是当我运行并使用nls时,我得到了错误: 我一直在处理参数,但没有成功,因此我根据模型模拟了一些数据,但仍然遇到相同的错误。 任何建议可能有什么问题:
2回复

什么是R中的SAS参数估计

请问如何用R写 跑;