R solo alternative a BUGS [chiuso]


13

Sto seguendo un corso sulle statistiche bayesiane usando BUGS e R. Ora, conosco già BUGS, è fantastico ma non mi piace molto usare un programma separato piuttosto che solo R.

Ho letto che ci sono molti nuovi pacchetti bayesiani in R. Esiste un elenco o un riferimento su quali pacchetti ci sono per le statistiche bayesiane e cosa fanno? E esiste un'alternativa al pacchetto R per la flessibilità di BUGS?

Risposte:


8

Puoi dare un'occhiata al pacchetto MCMCglmm fornito con vignette molto belle. C'è anche una bayesglm()funzione per il montaggio di modelli lineari generalizzati bayesiani nel pacchetto braccio , di Andrew Gelman. Ho anche sentito parlare di una versione blmer / bglmerfunzioni future per la modellazione gerarchica nello stesso pacchetto.


2
Nel braccio del pacchetto ci sono i bug di fucntion, che ti permettono di chiamare bug da R. Questo è quello che uso nella mia ricerca. Nel blog di Gelman c'è un esempio di calling bug di R.
Manoel Galdino,


6

Alcune persone che conosco hanno usato JAGS . La sintassi di JAGS è simile a BUGS.


(+1 ma penso che l'OP stia cercando una soluzione R pura.) Funziona benissimo con il pacchetto rjags , ma dobbiamo ancora specificare il nostro modello nella sintassi BUGS in un file esterno.
chl

5

Secondo la vista delle attività bayesiane. Aggiungerei un voto per MCMCpack , un pacchetto maturo che offre una varietà di modelli. Per la maggior parte è anche abbastanza ben documentato.


0

Le prestazioni sono il motivo principale per cui le persone usano WinBUGS / OpenBUGS / JAGS rispetto a pacchetti come MCMglmm. È molto difficile non pratico scrivere un efficiente campionatore Gibbs in nativo R. Esistono pacchetti che consentono di eseguire modelli BUGS da uno script R, in particolare RBUGS e BUGSParallel .


1
MCMCglammè un cattivo esempio perché "[a] ll simulazione viene eseguita in C / C ++ usando la libreria CSparse per sistemi lineari sparsi" (vedi abstract ).
Bernd Weiss,

1
-1; vedi @Bernd. La maggior parte dei pacchetti maturi usa il codice compilato. Il motivo principale è che BUGS et al sono più flessibili in quanto possono adattarsi a più modelli. Sebbene ciò possa portare a un calcolo più efficiente poiché un pacchetto R - anche con codice compilato - deve essere più generale, potrebbe non esserlo.
JMS

1
MCMCpack utilizza il codice C / C ++ compilato, ottimizzato per l'attività a portata di mano, quindi in realtà è più veloce rispetto a fare qualcosa in un pacchetto generalizzato come JAGS (per un'attività specifica).
Wayne,
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.