Implementasi Dirichlet cdf?

Saya perlu menghitung CDF Dirichlet , tetapi saya hanya dapat menemukan implementasi dari PDF .

Apakah Anda tahu ada perpustakaan (terutama di R) yang mengimplementasikannya?

r matlab pdf cdf dirichlet-distribution Ricky Robinson
sumber

Tidak secara langsung mengetahui adanya. Tetapi mungkin ada sesuatu yang bisa dilakukan. Apa yang perlu Anda lakukan dengannya?

Glen_b -Reinstate Monica

Saya perlu mengambil pelengkap CDF dan menganggapnya sebagai nilai p saya.

Ricky Robinson

Hmm. Jadi ya, jika Anda memerlukan , Anda perlu melakukan cdf. Ide simulasi Zen jelas merupakan cara untuk melakukannya (dan semakin tinggi jumlah dimensi, semakin baik mulai terlihat), tetapi jika Anda melakukannya, gunakan salah satu paket dengan implementasi bawaan dari . Jika hanya 3-variate atau mungkin 4-variate (komponen terakhir, tentu saja, redundan) mungkin ada baiknya mencoba numerik quadrature.

1 - P (X_{1} \leq x_{1}, X_{2} \leq x_{2}, . . ., X_{k} \leq x_{k})

$1-P(X_1\leq x_1,X_2\leq x_2, ...,X_k\leq x_k)$ rdirichlet

Glen_b -Reinstate Monica

Jawaban:

Ingat itu, jika independen , untuk , maka $Y_i$ $\mathrm{Gamma}(a_i,b)$ $i=1,\dots,k$

(X_{1}, \dots, X_{k}) = (\frac{Y_{1}}{\sum_{j = 1}^{k} Y_{j}}, \dots, \frac{Y_{k}}{\sum_{j = 1}^{k} Y_{j}}) \sim D i r i c h l e t (a_{1}, \dots, a_{k}) .

$(X_1,\dots,X_k) = \left(\frac{Y_1}{\sum_{j=1}^k Y_j}, \dots, \frac{Y_k}{\sum_{j=1}^k Y_j} \right) \sim \mathrm{Dirichlet}(a_1,\dots,a_k) \, .$

Buktinya dapat ditemukan pada halaman 594 dari Luc Devroye ini buku .

Oleh karena itu, satu kemungkinan adalah untuk menghitung perkiraan Monte Carlo untuk dimulai dengan gammas. Di , coba ini:

F_{X_{1}, \dots, X_{k}} (t_{1}, \dots, t_{k}) = P {X_{1} \leq t_{1}, \dots, X_{k} \leq t_{k}},

$F_{X_1,\dots,X_k}(t_1,\dots,t_k)=P\left\{X_1\leq t_1,\dots, X_k\leq t_k\right\} \, ,$ R

pdirichlet <- function(a, t) {
    N <- 10000
    rdirichlet <- function(a) { y <- rgamma(length(a), a, 1); y / sum(y) }
    x <- replicate(N, rdirichlet(a), simplify = FALSE)
    sum(sapply(x, function(x) prod(x <= t))) / N
}

Saya tidak memeriksa kode. Gunakan dengan hati-hati. Jika Anda menemukan bug, beri tahu kami.

Zen
sumber

Ada fungsi rdirichlet vektor di R sudah - pada kenyataannya beberapa dari mereka (di gtools, MCMCpackdan dirmultmisalnya).

Glen_b -Reinstate Monica

@ Zen saya punya a <- c(6, 20,2)cara untuk mendapatkan cdf Drichelt? Apakah matriks t 2 x 2?

skor 324

Saya menggunakan kode di atas, tetapi itu melemparkan kesalahan.

skor 324

@ score324 Vektor Dirichlet biasanya memiliki setiap elemen dalam dan jumlahnya adalah , jadi berada di luar dukungan

[0, 1]

$[0,1]$

1

$1$

(6, 20, 2)

$(6, 20,2)$

Henry

Ada perpustakaan? Mathematica memilikinya. Berikut kode untuk contoh plot CDF Dirichlet dari dokumentasi:

Plot3D[CDF[DirichletDistribution[{1, 3, 2}], {x, y}], {x, 0, 1}, {y, 0, 1}]

Mike Z.
sumber

Bagaimana cara mendapatkan ekspresi untuk CDF [DirichletDistribution [{1, 3, 2}]]?

AIB