Generate random number from inverse gamma distribution
5 views (last 30 days)
Show older comments
Hi, I want to get the sample from the inverse gamma distribution with two parameter a and b and I am wondering how I can do that in matlab? Please help me ! Thanks
0 Comments
Answers (1)
Tom Lane
on 27 May 2012
If I understand the inverse gamma correctly, you can generate from x=gamrnd(a,1./b) and then take y=1./x. This would require the Statistics Toolbox.
3 Comments
Zhiyong Yang
on 27 Nov 2018
Hi, Jonathan
Have you solved this problem?I have compared the pdf curves obtained by ksdensity method of samples draw by the mentioned way using gamrnd function and by analytical solution using the equation of inverse-gamma pdf. It is shown as the alpha is smaller than 2, the mentioned way cannot sample well. Please see the attached figure where symbols represent pdf from ksdensity method and lines represent pdf from analytical solutoin. ![untitled.bmp](https://www.mathworks.com/matlabcentral/answers/uploaded_files/197080/untitled.bmp)
![untitled.bmp](https://www.mathworks.com/matlabcentral/answers/uploaded_files/197080/untitled.bmp)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!