为了从R中的反 Gamma 分布采样,请按照以下正确方法进行操作:

#I want to sample an inverse-gamma(a,b)

a = 4
b = 9

x = 1/rgamma(1,a,b)

最佳答案

尽管@Dason和@Stephane已经评论过您的方法是有效的,但是R中有几个软件包可以执行此操作(在r inverse gamma上进行了谷歌搜索:

  • rigamma {pscl}
  • rinvgamma {MCMCpack}

  • 另请参见Wikipedia页面上的gamma distributioninverse gamma distribution,以了解两种分布的概率密度函数:

    Gamma 分布与:

    反 Gamma 。

    关于r - 从R中的反 Gamma 分布采样,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/18395368/

    10-12 21:51