用R中的双变量伽马分布生成随机数
问题描述：
如何用双变量伽马分布生成随机数。密度为：用R中的双变量伽马分布生成随机数
˚F（ X，ÿ）（ X，ÿ）= α P + Q X P-1（YX）Q-1ë- αÿ/[Γ（ p）Γ（ q）]，＆＃x1d540; 0＆＃x2264; X＆＃x2264; y
其中y> x> 0，α> 0，p> 0且q> 0。
我没有在R上找到任何包，而没有在文献中做任何包。
答
这是简单的：
生成X〜伽玛（P，阿尔法）（阿尔法是在制剂中的速率参数）
-
生成W〜伽玛（Q，阿尔法），独立X
的 计算Y = X + W
（X，Y）具有所需的二元分布。
在R（假定P，Q，α-和n已定义）：
x <- rgamma(n,p,alpha)
y <- x + rgamma(n,q,alpha)
从与参数P，Q，α-
二元麦凯分布二元分布产生n个值在VGAM库中实现，检查bivgamma.mackay –
似乎这个包只使用最大似然估计来估计双变量伽马分布的参数。 – fsbmat