58066_zh-TW

liliconcom

JMP中的Dirichlet分佈（適用於Bayes之前）

嗨，

 

我正在嘗試一些貝葉斯引導，並希望在JMP Pro中進行。我想要使用的先前分佈之一是Dirichlet。我正在考慮使用列數據類型Expression來保存權重的向量並使用壓縮到標籤地址索引列中的唯一值，以便我知道向量的哪一行/列包含要應用的正確權重。

 

我可以使用R連接器（下面）生成Dirichlet向量，但是想在公式列中進行，這樣我就可以利用Pro中的Bootstrap功能。有什麼建議？

 

Names Default To Here( 1 );
R Init();

//input list, outputlist, Rcode
rc = R Execute(
{nSamples,pollResult},
{x},
 "
 library (DirichletReg)
 
 x <- rdirichlet(nSamples, c(728,584,138))"
);
Wait( 3 );
//get dataframe from R
resultMatrix = R Get(x);

//close connection
R Term();

//export results
dt = As Table( resultMatrix );

Re：JMP中的Dirichlet分佈（適用於Bayes之前）

結合此信息與展開矢量列加載項給了我我需要的東西。謝謝。 

Re：JMP中的Dirichlet分佈（適用於Bayes之前）

嗨斯蒂芬，

 

您可以將下面的代碼粘貼到列公式中，它應該從Dirichlet進行採樣。我建議你在使用它之前說服自己這是準確的，因為我沒有測試它。我是根據它寫的 https://en.wikipedia.org/wiki/Dirichlet_distribution#Random_number_generation

 

alpha = [1 2 3 4 5];
gams = Random Gamma(alpha)`;
gams/sum(gams); // values sampled from dirichlet with concentration parameters alpha 

 

米洛

Re：JMP中的Dirichlet分佈（適用於Bayes之前）

Dirichlet分佈不能直接在JMP中使用，但在這種情況下您不需要調用R. 您可以創建一個列公式，在該公式中使用JMP中可用的函數。例如，見“ 什麼是Dirichlet分佈的直觀解釋？ “ 欲獲得更多信息。