Distribuzioni di Dirichlet in JMP (per Bayes precedente)Ciao,
Sto provando alcuni bootstrap bayesiani e mi piacerebbe farlo in JMP Pro. Una delle precedenti distibuzioni che vorrei usare è Dirichlet. Stavo pensando di usare il tipo di dati colonna Expression per contenere i vettori dei pesi e usare il Comprimi nell'etichetta Addin per indicizzare i valori univoci nella colonna, quindi so quale riga / colonna del vettore contiene il peso corretto da applicare.
Posso generare i vettori di Dirichlet usando il connettore R (sotto) ma vorrei farlo in una colonna formula in modo da poter sfruttare la funzione Bootstrap in Pro. Eventuali suggerimenti?
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: Distribuzioni di Dirichlet in JMP (per Bayes precedente)Combinando queste informazioni con Espandi Componente aggiuntivo colonna vettoriale mi ha dato proprio quello di cui avevo bisogno. Grazie.
Re: Distribuzioni di Dirichlet in JMP (per Bayes precedente)Ciao Stefano,
Puoi incollare il codice qui sotto in una formula di colonna, e dovrebbe campionare da un Dirichlet. Ti consiglierei di convincerti che è accurato prima di usarlo, dal momento che non l'ho testato. L'ho scritto basandomi su 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
Milo
Re: Distribuzioni di Dirichlet in JMP (per Bayes precedente)La distribuzione di Dirichlet non è disponibile direttamente in JMP, ma in questo caso non è necessario chiamare R. È possibile creare una formula di colonna in cui si utilizzano le funzioni disponibili in JMP. Vedi, ad esempio, " Qual è una spiegazione intuitiva della distribuzione di Dirichlet? " per maggiori informazioni.
