Sto usando PyMC3 per eseguire modelli bayesiani sui miei dati.
Sono nuovo nella modellazione bayesiana ma secondo alcuni post di blog , Wikipedia e QA di questo sito, sembra essere un approccio valido per utilizzare il fattore Bayes e il criterio BIC per essere in grado di scegliere quale modello rappresenti meglio i miei dati (quello che genera i miei dati).
Per calcolare il fattore Bayes, ho bisogno della probabilità relativa per i modelli che voglio confrontare. Forse è un po 'confuso per me, ma penso che ci sia due modi per ottenere la probabilità (correggimi se sbaglio):
il modo algebrico quando i modelli sono semplici: vedi la pagina del fattore Bayes dell'esempio di Wikipedia
il modo numerico: questo è ciò che fa PyMC3 con gli algoritmi MCMC
Come posso accedere alla probabilità e quindi confrontare i miei modelli in PyMC3? Ho trovato un model.logp
metodo che secondo il documento è la "funzione di densità della probabilità di registro". Posso usarlo per ottenere la probabilità?
Domanda bonus: quando vengono confrontati due modelli, viene calcolato il rapporto tra entrambe le probabilità. Cosa succede se si desidera confrontare diversi modelli?
Un esempio concreto di PyMC3 sarebbe molto utile!