用R中的雙變量伽馬分佈生成隨機數

Question

如何用雙變量伽馬分佈生成隨機數。密度爲：

˚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上找到任何包，而沒有在文獻中做任何包。

Answer 1

這是簡單的：

1. 生成X〜伽瑪（P，阿爾法）（阿爾法是在製劑中的速率參數）

2. 生成W〜伽瑪（Q，阿爾法），獨立X

   的

3. 計算Y = X + W

4. （X，Y）具有所需的二元分佈。

在R（假定P，Q，α-和n已定義）：

x <- rgamma(n,p,alpha) 
y <- x + rgamma(n,q,alpha) 

從與參數P，Q，α-