"JAGS è solo un altro Gibbs Sampler. È un programma per l'analisi di modelli gerarchici bayesiani che utilizza la simulazione Markov Chain Monte Carlo (MCMC) non del tutto diversa da BUGS". (Http://mcmc-jags.sourceforge.net/)
Sto per provare un ambiente in stile BUGS per stimare i modelli bayesiani. Ci sono dei vantaggi importanti da considerare nella scelta tra OpenBugs o JAGS? È probabile che uno sostituisca l'altro nel prossimo futuro? Userò il Gibbs Sampler scelto con R. Non ho ancora un'applicazione specifica, ma piuttosto sto …
Ho trovato alcune distribuzioni per le quali BUGS e R hanno parametrizzazioni diverse: normale, log-normale e Weibull. Per ognuno di questi, ho capito che il secondo parametro utilizzato da R deve essere trasformato inverso (1 / parametro) prima di essere utilizzato in BUGS (o JAGS nel mio caso). Qualcuno sa …
Nel suo ampiamente citato documento Distribuzioni precedenti per parametri di varianza in modelli gerarchici (916 citazione finora su Google Scholar) Gelman propone che buone distribuzioni precedenti non informative per la varianza in un modello bayesiano gerarchico siano la distribuzione uniforme e la distribuzione della mezza t. Se capisco bene le …
Ho pensato di poter giocare con un po 'di selezione delle variabili bayesiane, seguendo un bel post sul blog e i relativi documenti. Ho scritto un programma in rjags (dove sono piuttosto un principiante ) e ho recuperato i dati sui prezzi per Exxon Mobil, insieme ad alcune cose che …
Ho appena iniziato a imparare a usare Stan e rstan. A meno che non sia sempre stato confuso su come funzionava JAGS / BUGS, ho pensato che dovevi sempre definire una distribuzione precedente di qualche tipo per ogni parametro nel modello da cui attingere. Sembra che non sia necessario farlo …
Esistono diversi documenti matematici che descrivono il lazo bayesiano, ma voglio un codice JAGS testato e corretto che posso usare. Qualcuno potrebbe pubblicare un codice BUGS / JAGS di esempio che implementa la regressione logistica regolarizzata? Qualsiasi schema (L1, L2, Elasticnet) sarebbe fantastico, ma è preferito Lasso. Mi chiedo anche …
Sto cercando di adattare un modello gerarchico usando jags e il pacchetto rjags. La mia variabile di risultato è y, che è una sequenza di prove di bernoulli. Ho 38 soggetti umani che si esibiscono in due categorie: P e M. Sulla base della mia analisi, ogni oratore ha una …
Ho usato rjags per eseguire MCMC su un modello, specificato nel linguaggio JAGS. C'è un buon modo per estrarre quel modello ed eseguire previsioni con esso (usando le distribuzioni posteriori dei miei parametri)? Posso specificare nuovamente il modello in R e collegare le modalità dei miei parametri posteriori; Mi chiedo …
Sto cercando di impostare un modello poisson gonfiato a zero in R e JAGS. Sono nuovo di JAGS e ho bisogno di alcune indicazioni su come farlo. Ho provato con quanto segue in cui y [i] è la variabile osservata model { for (i in 1:I) { y.null[i] <- 0 …
Sembra che i condizionali completi siano spesso piuttosto difficili da ricavare, ma programmi come JAGS e BUGS li derivano automaticamente. Qualcuno può spiegare come generano algoritmicamente i condizionali completi per qualsiasi specifica di modello arbitraria?
Sto costruendo un modello bayesiano gerarchico piuttosto complesso per una meta-analisi usando R e JAGS. Semplificando un po ', i due livelli chiave del modello hanno dove è la esima osservazione del endpoint (in questo caso, rese colturali GM / non GM) nello studio , è l'effetto per lo studio …
In R"peso precedente" possiamo glmregredire tramite il parametro pesi . Per esempio: glm.D93 <- glm(counts ~ outcome + treatment, family = poisson(), weights=w) Come può essere realizzato in un JAGSo BUGSmodello? Ho trovato alcuni articoli che ne discutono, ma nessuno di loro fornisce un esempio. Sono interessato principalmente agli esempi …
Ho appena (ri) letto il perché di Gelman (di solito) non dobbiamo preoccuparci di confronti multipli . In particolare, la sezione "Risultati multipli e altre sfide" menziona l'uso di un modello gerarchico per situazioni in cui vi sono più misure correlate della stessa persona / unità in tempi / condizioni …
Sto cercando di costruire un modello in cui la risposta è proporzionale (in realtà è la percentuale di voti che un partito ottiene nei collegi elettorali). La sua distribuzione non è normale, quindi ho deciso di modellarlo con una distribuzione beta. Ho anche diversi predittori. Tuttavia, non so come scriverlo …
Ho una domanda su come adattare un problema di censura in JAGS. Osservo una miscela bivariata normale in cui i valori X presentano errori di misurazione. Vorrei modellare i veri "mezzi" sottostanti dei valori censurati osservati. ⌈ xt r u e+ ϵ ⌉ = xo b s e r v …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.