Questions marquées «jags»

Je suis sur le point d'essayer un environnement de type BUGS pour estimer les modèles bayésiens. Y at-il des avantages importants à considérer dans le choix entre OpenBugs ou JAGS? L'un est-il susceptible de remplacer l'autre dans un avenir prévisible? Je vais utiliser le sampler choisi avec Gibbs...

J'ai trouvé quelques distributions pour lesquelles BUGS et R ont des paramétrisations différentes: Normal, log-Normal et Weibull. Pour chacun d'eux, je suppose que le deuxième paramètre utilisé par R doit être transformé inversement (1 / paramètre) avant d'être utilisé dans BUGS (ou JAGS dans mon...

Dans son article largement cité Prior distributions for variance parameters in hierarchical models (916 citation à ce jour sur Google Scholar) Gelman propose que de bonnes distributions a priori non informatives pour la variance dans un modèle bayésien hiérarchique soient la distribution uniforme...

Je pensais que je pourrais jouer avec une sélection de variables bayésiennes, à la suite d'un bel article de blog et des articles liés. J'ai écrit un programme dans rjags (où je suis plutôt une recrue) et récupéré des données de prix pour Exxon Mobil, ainsi que certaines choses qui sont peu...

Je viens de commencer à apprendre à utiliser Stan et rstan. À moins que je ne sois toujours confus sur le fonctionnement de JAGS / BUGS, je pensais que vous deviez toujours définir une distribution préalable d'une sorte pour chaque paramètre du modèle à partir duquel tirer. Il semble que vous...

J'essaie d'adapter un modèle hiérarchique à l'aide de jags et du package rjags. Ma variable de résultat est y, qui est une séquence d'essais bernoulli. J'ai 38 sujets humains qui se produisent sous deux catégories: P et M. D'après mon analyse, chaque locuteur a une probabilité de succès dans la...

Il y a plusieurs articles mathématiques qui décrivent le lasso bayésien, mais je veux un code JAGS correct et testé que je peux utiliser. Quelqu'un pourrait-il publier un exemple de code BUGS / JAGS qui implémente une régression logistique régularisée? Tout schéma (L1, L2, Elasticnet) serait...

J'essaie de mettre en place un modèle de poisson zéro gonflé dans R et JAGS. Je suis nouveau chez JAGS et j'ai besoin de conseils sur la façon de le faire. J'ai essayé avec ce qui suit où y [i] est la variable observée model { for (i in 1:I) { y.null[i] <- 0 y.pois[i] ~ dpois(mu[i]) pro[i] <-...

J'ai utilisé rjags pour exécuter MCMC sur un modèle, spécifié dans le langage JAGS. Existe-t-il un bon moyen d'extraire ce modèle et d'effectuer des prédictions avec lui (en utilisant les distributions postérieures de mes paramètres)? Je peux re-spécifier le modèle en R et brancher les modes de mes...

Il semble que les conditions complètes soient souvent assez difficiles à dériver, mais des programmes comme JAGS et BUGS les dérivent automatiquement. Quelqu'un peut-il expliquer comment il génère algorithmiquement des conditions complètes pour toute spécification de modèle...

Dans Ron peut "pondérer préalablement" une glmrégression via le paramètre des poids . Par exemple: glm.D93 <- glm(counts ~ outcome + treatment, family = poisson(), weights=w) Comment cela peut-il être accompli dans un modèle JAGSou BUGS? J'ai trouvé un article sur ce sujet, mais aucun ne fournit...

J'essaie d'adapter un modèle à temps discret dans R, mais je ne sais pas comment le faire. J'ai lu que vous pouvez organiser la variable dépendante dans différentes lignes, une pour chaque observation de temps, et utiliser la glmfonction avec un lien logit ou cloglog. En ce sens, j'ai trois...

Je construis un modèle bayésien hiérarchique assez complexe pour une méta-analyse utilisant R et JAGS. Pour simplifier un peu, les deux niveaux clés du modèle ont où est la ème observation de la (dans ce cas, les rendements des cultures GM vs non GM) dans l'étude , est l'effet pour l'étude , les s...

Je viens de (re) lire pourquoi nous n'avons généralement pas à nous soucier des comparaisons multiples . En particulier, la section «Résultats multiples et autres défis» mentionne l'utilisation d'un modèle hiérarchique pour les situations où il existe plusieurs mesures connexes de la même personne...

J'essaie de construire un modèle où la réponse est une proportion (c'est en fait la part des votes qu'un parti obtient dans les circonscriptions). Sa distribution n'est pas normale, j'ai donc décidé de le modéliser avec une distribution bêta. J'ai également plusieurs prédicteurs. Cependant, je ne...

Gelman & Hill (2006) disent: Dans Bugs, les résultats manquants dans une régression peuvent être facilement gérés en incluant simplement le vecteur de données, les NA et tout. Les bogues modélisent explicitement la variable de résultat, et il est donc trivial d'utiliser ce modèle pour, en e ff...

J'ai une question sur la façon de régler un problème de censure dans JAGS. J'observe un mélange bivarié normal où les valeurs X ont une erreur de mesure. Je voudrais modéliser les véritables «moyens» sous-jacents des valeurs censurées observées. ⌈ xt r u e+ ϵ ⌉ = xo b s e r v e d ϵ ∼ N( 0 , s d=...

Fermé. Cette question est hors sujet . Il n'accepte pas actuellement les réponses. Vous souhaitez améliorer cette question? Mettez à jour la question afin qu'elle soit sur le sujet pour la validation croisée. Fermé il y a 2 ans . R semble être capable de produire de jolis tracés récapitulatifs à...

Dans R, c (3,1,0) * c (2,0,1) == c (6,0,0). Ce n'est pas un produit scalaire et ce n'est pas un produit croisé. Premièrement, quel est le nom de ce produit, et deuxièmement, cela fonctionne-t-il dans WinBUGS, OpenBUGS et / ou

Problème Je voudrais faire une inférence sur un système analogue à mourir avec un nombre inconnu de côtés. Le dé est lancé plusieurs fois, après quoi je voudrais déduire une distribution de probabilité sur un paramètre correspondant au nombre de côtés du dé, θ. Intuition Si après 40 rouleaux vous...