二変量のガンマ分布を持つ乱数を生成するR

Question

二変量のガンマ分布を持つ乱数を生成する方法。密度は：

F _{（X、 Y）}（ X、 Y）= α ^{のp + qを} X ^P-1（YX）^Q-1 E^{からα Y}/[Γ（ P）Γ（ Q）]、＆＃のx1d540。 _{0＆＃x2264; x＆＃x2264; 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) 

は、n個のパラメータpの二変量分布からの値、Q、アルファ