问题标签 [gamma-distribution]

计算 Gamma Inverse。我正在使用 phpExcel 库：

https://github.com/PHPOffice/PHPExcel/blob/develop/Classes/PHPExcel/Calculation/Statistical.php

我使用这个库并将其与 Excel 的结果进行比较。对于 PHP 和 Excel，我都用相同的值进行了测试，但得到了不同的结果：

PHP: GAMMAINV(0.999, 21.054, 1/0.266) 结果: 395.75187969925

在 Excel 中：=GAMMA.INV(0.999 , 21.054, 1/0.266) 结果：143.076

使用 BETA.INV 和 BETAINV 进行的测试运行良好。所以我不明白这里出了什么问题。有人可以解释和/或指出一个正确的图书馆吗？

尺度参数为 1 的 gamma 分布的对数似然可以写为：

(α−1)s−nlogΓ(α) 其中 alpha 是形状参数，s=∑logXi是充分的统计量。

随机抽取 n = 30 的样本，形状参数为 alpha = 4.5。使用newton_search和make_derivative，找到 alpha 的最大似然估计。使用 alpha 的矩估计，即 x 的平均值作为初始猜测。R中的对数似然函数是：

我创建了 make_derivative 函数，如下所示：

我还创建了一个newton_search包含该功能的make_derivative功能；但是，我需要使用newton_search对数似然函数的二阶导数，并且我不确定如何修复以下代码才能做到这一点：

提示：您必须应用对数似然的newton_search一阶和二阶导数（使用 数字导出）。make_derivative你的答案应该接近 4.5。

当我跑步时newton_search2(gllik(), 0.0001, mean(x), conv = 0.001)，我得到的答案应该是双倍的。

我正在尝试绘制一个 Inverse 的 ksdensity，gamma(0.001, 0.001)但该图只有一个点。我使用的命令是

的第一项f是实数，其他都是NaN。的第一项xi是实数，其他都是Inf。

你能帮我解决这个问题吗？

非常感谢。

我有两个向量：

我正在尝试使用非线性最小二乘函数来拟合数据，但我不断收到错误消息：

nlsModel(formula, mf, start, wts) 中的错误：初始参数估计处的奇异梯度矩阵

我的代码是：

我输入的参数是我知道的接近预期值的参数。有谁知道我在这里做错了什么？

我已经尝试了参数的各种初始值a，b和k，但我总是得到某种错误。

我对某个向量有疑问。我想知道它是否是伽马分布的，以及（如果是的话）参数（形状、速率）是什么。我的向量有 400 个条目，但让我们举个例子

我读了一些关于 fitdistr() 的东西。但我不太明白它的实际作用！我用我的真实（长）向量尝试了以下操作：

输出是什么意思？这些是我的拟合参数吗？我对它们进行了测试，但 KS-Test 给了我一个否定的结果：

那么您能否告诉我如何确定我的向量是否是伽马分布的以及参数是什么？

我想在观察数据的直方图上绘制从一组观察结果中得出的伽马密度函数。我能够为伽马拟合生成直方图和参数估计。这是针对来自主数据集的多个数据子集完成的。如何在此循环中创建的每个直方图上绘制伽马密度函数？

我目前有：

我在想我应该使用lines()，但是，不确定如何指定这个？

我有一个适合我的数据的伽马分布，使用libary(fitdistrplus). 我需要确定一种方法来定义可以“合理”预期的 x 值范围，类似于使用具有正态分布的标准偏差。

例如，与平均值相差两个标准差内的 x 值可以被认为是正态分布的预期值的合理范围。关于如何根据伽马分布的形状和速率参数定义类似范围的预期值有什么建议吗？

...也许类似于识别 x 的两个值，这两个值之间包含 95% 的数据？

具有形状参数 k 和尺度参数 theta 的 gamma 分布由=

在R中，如果我想找到 Gamma(10,0.5) 的伽马分布概率为 0.05 的分位数

我用了

但这不是我想要的价值。我使用时得到的期望值，

qgamma(0.05,10,0.5)那么和有什么区别qgamma(0.05,shape=10,scale=0.5)。
为什么我会得到两个完全不同的结果？

假设我们在 MATLAB 中有这个伽马分布：

我希望这部分分布具有更高的密度（x 轴范围）。如何在 MATLAB 中提取它？histfit我使用函数拟合这个分布。

附言。我的代码：

我正在尝试使用以下代码从逆伽马分布中采样，获取修改我在网上找到的一些代码。我不是 C++ 编码方面的专家，所以我需要你的一些解释和帮助。

我得到的错误是：