BUGS est un acronyme pour l'inférence bayésienne utilisant l'échantillonnage de Gibbs; BUGS est également un progiciel pour cela.
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...
Je vais expliquer mon problème avec un exemple. Supposons que vous souhaitiez prédire le revenu d'un individu en fonction de certains attributs: {âge, sexe, pays, région, ville}. Vous avez un ensemble de données de formation comme ça train <- data.frame(CountryID=c(1,1,1,1, 2,2,2,2, 3,3,3,3),...
Bonjour gourous statistiques et assistants de programmation R, Je m'intéresse à la modélisation des captures d'animaux en fonction des conditions environnementales et du jour de l'année. Dans le cadre d'une autre étude, j'ai dénombré des captures sur environ 160 jours sur trois ans. Chaque jour,...
Cette question est un suivi technique de cette question . J'ai du mal à comprendre et à reproduire le modèle présenté dans Raftery (1988): Inférence pour le paramètre binomial NNN : une approche bayésienne hiérarchique dans WinBUGS / OpenBUGS / JAGS. Il ne s'agit pas seulement de code, donc il...
Je me demandais quel progiciel statistique de logiciels recommandez-vous pour effectuer l'inférence bayésienne. Par exemple, je sais que vous pouvez exécuter openBUGS ou winBUGS de manière autonome ou vous pouvez également les appeler à partir de R. Mais R a également plusieurs de ses propres...
Fermé. Cette question est hors sujet . Il n'accepte pas actuellement les réponses. Voulez-vous améliorer cette question? Mettez à jour la question afin qu'elle soit sur le sujet pour la validation croisée. Fermé il y a 12 mois . Je suis un cours sur les statistiques bayésiennes utilisant BUGS et...
J'ai un ensemble de données de série chronologique auquel j'essaie d'adapter un modèle de Markov caché (HMM) afin d'estimer le nombre d'états latents dans les données. Mon pseudo-code pour ce faire est le suivant: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ......
Il existe maintenant plusieurs approches différentes pour effectuer une méta-analyse de réseau ou une comparaison de traitements mixtes. Les plus utilisés et les plus accessibles sont probablement les suivants: dans un cadre bayésien : approche d'interaction conception par traitement dans WinBUGS...
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...
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...
Je suis un nouvel utilisateur de WinBUGS et j'ai une question pour votre aide. Après avoir exécuté le code suivant, j'ai obtenu les paramètres de beta0through beta4(statistiques, densité), mais je ne sais pas comment obtenir la prédiction de la dernière valeur de h, que j'ai définie NApour...
Que se passe-t-il lorsque vous n'avez aucune idée de la distribution des paramètres? Quelle approche devrions-nous utiliser? La plupart du temps, nous visons à sous-estimer si une certaine variable a une influence sur la présence / absence d'une certaine espèce, et si la variable est acceptée ou...
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...
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...
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...
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
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 à...
Je souhaite installer une ANOVA bayésienne à deux facteurs dans des BUGS ou utiliser un package R. Malheureusement, j'ai du mal à trouver des ressources sur ce sujet. Aucune suggestion? Même un article décrivant l'approche serait