How to generate random number from inverse gamma distribution in an extreme conditions

46 次查看（过去 30 天）
chendong zhang
chendong zhang 2019-10-23
回答： Matteo Iacopini 2020-12-27
hello everyone
I want to get the sample from the inverse gamma distribution as below
x~IG(54,0.004)
but when I try to use the reciprocal relationship between gamma distribution and inverse gamma distribution, like
first, generate a random number from t~G(54,0.004), then set x=1./t, and the result is:
3.66281673846745 4.15049653026671 5.59965910607058
the matlab code is:
tmp=gamrnd(54,0.004,1,3);
res=1./tmp;
the result seems to have some problems, I also test the sample in R, the result is below
7.432166e-05 6.758658e-05 6.506619e-05
the R code is:
rinvgamma(3,54,0.004)
Is this the internal algorithm problem with Matlab in handling such extreme condition?

请先登录，再回答此问题。

回答（2 个）

Shep Bryan
Shep Bryan 2020-1-24
Accoding to wikipedia, if X ~ InvGamma(a,b) then 1/X ~ Gamma(a,b), therefore to sample from the inverse gamma you simply take the inverse of a random variable sampled from the gamma distribution.
Matteo Iacopini
Matteo Iacopini 2020-12-27
I don't you whether you have already solved the issue, however you can sample from an Inverse Gamma by drawing values from a Gamma and then takinng the inverse. Note that if X ~ Ga(a, b), where a= SHAPE and b= RATE (i.e., inverse of SCALE), then Y = 1/X ~ IG(a, b).
The point here is that MATLAB implements the SHAPE-SCALE parametrization of the Gamma distribution, not the SHAPE-RATE (which is the one you are looking for, see the relationship above).
Therefore, all you need to do to sample from IG(a,b) is
x = gamrnd(a, 1./b);
y = 1./x;
the first line generates from a Gamma with SHAPE a and RATE b, the second inverts to obtain Inverse Gamma random variates.

类别

在 Help Center 和 File Exchange 中查找有关 Gamma Distribution 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by