问题标签 [gamma-distribution]

我想从 C++ 中的伽马分布生成样本。我正在超级计算机上编译我的代码，它没有最新版本的 g++（也无法更新）。我编写了代码，它使用 chrono 和 random 从 gamma 分布中生成样本。不幸的是，g++ 版本不允许我编译代码（我也不能使用标志 -std=c++11 或 -std=c++0x）。我也不能使用 boost/TR1，因为编译器没有它。

除了回到基础知识（从指数总和生成伽马等）之外，还有另一种方法可以从这个分布中采样吗？我正在考虑纯 c++ 而不是实现 MCMC。

如何从 Java 中的伽马分布生成数字？

这是我在这里的第一篇文章，我希望我会遵守社区的所有规则。

我正在尝试使用来自镶嵌包的函数 antiD 计算 R 中形状参数 2 和比例参数 3 的伽马分布方差。我使用的R代码如下

问题是我得到的数字没有意义，因为这些参数的 Gamma 分布的方差应该等于 2*3^2 = 18 （关于 Gamma 分布的 Wiki 页面）。此外，如果我将 10^4 作为variance() 的上限（默认下限为0），它将返回以下内容：

从 10^4 到 10^5 的积分将是：

有谁知道为什么variance(10^5)在这种情况下会产生荒谬的结果？我也将不胜感激对帖子风格的任何其他评论。

我在 SAS 中有以下代码，我需要为 R 中的相同数据集重新创建此模型。我无法找到合适的代码来重新创建模型，尤其是将数据的分布更改为伽马分布。

SAS 代码

任何帮助，将不胜感激

matlab 函数nbinfit返回负二项式的值r和p。matlab 中是否有等效的 MLE 函数，它返回负二项式的“生态”或“Polya”参数化的 mu（均值）和 theta（色散参数）的值？

（类似于fitdistrR 中的东西）

我正在尝试创建一个函数来计算二项式公式，但我不断收到错误“伽马 fn 的值超出范围”

我使用 scipy 尝试了Saullo Castro 的示例代码来获取 gamma 分布的参数。

给我熵结果，这是没有意义的。我在这里错在哪里？

在 scipy 文档中，它指出

在哪里

我试图估计特威迪（或复合泊松伽马）分布的有限混合。我搜索了我能想到的任何资源，但没有找到任何关于如何做到这一点的资源。

我目前正在尝试使用 R 中的flexmix包编写不同的 M-step 驱动程序，如第 12-14 页的flexmix小插图中所述。这是我的代码，它依赖于cplm包：

但是，这会导致以下错误：

dtweedie(y, xi = p, mu = mu, phi = phi) 中的错误：对不一致数组的二元运算

有没有人做过或见过 tweedie 分布的有限混合？您能指出我使用flexmix还是其他方式来完成此任务的正确方向吗？

我想使用 Matlab 中的“randg”函数将乘法伽马噪声添加到图像中并去除该噪声。我们必须记住，噪声应该具有均值 1 和级别 4。它应该遵循 Gamma 定律（使用 Gamma 概率分布函数）。添加噪声后的图像变为

f=u*v; 其中 f = 噪声图像，u = 原始图像，v = 噪声图像。

伽马定律为： gv(v)=L^L/(Γ(L)) v^(L-1) exp(-Lv) 1_(v≥0)

其中 L 是噪声水平，v 是噪声。

这是我尝试过的代码：

但我无法得到预期的结果。请给我一个方法。

我发布这个部分是为了得到一个问题的答案，部分是为了让其他人以后可能会找到它并让他们免于头痛。

我正在尝试按照Epix Analytics的公式在 NetLogo 中创建 beta-pert 分布

我在 Beta 分布 wiki 页面上使用 gamma distbutiotn 的转换创建了一个 Beta(Alpha1, alpha2) 分布：

我创建的代码是：

我的问题是关于维基百科页面从 Gamma 分布到 Beta 分布的转换。它有一个 theta 作为伽马分布的第二部分，但从未定义那是什么。我把它当作占位符，只要 X 和 Y 的 theta 相同，就不会影响结果。我玩过不同的 theta 值，分布没有明显变化。这是一个正确的假设吗？