95％CI用于参数估计
95% CI for parameters estimation

发表于 2013-04-16 04:05:57

活跃于 2013-04-16 06:13:16

查看 584 次

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

2回复

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

2回复

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

1回复

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

1回复

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

1回复

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

1回复

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

1回复

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

1回复

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

1回复

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

2回复

请问如何用R写跑;