こんにちは、
私はいくつかのベイジアンブートストラップを試していますし、JMP Proでそれをやりたいと思います。私が使いたいと思う以前の配布の1つが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 );この情報をベクトル列アドインを展開私に必要なものだけを与えてくれました。ありがとう。
こんにちはスティーブン、
以下のコードを列式に貼り付けることができます。そして、それは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
ミロ
ディリクレ分布はJMPでは直接利用できませんが、この場合はRを呼び出す必要はありません。JMPで利用可能な機能を使用する列式を作成できます。例えば、「 ディリクレ分布の直感的な説明は何ですか? " 詳細については。