Un type spécial d'algorithme de Markov Chain Monte Carlo (MCMC) utilisé pour simuler à partir de distributions de probabilités complexes. Il est validé par la théorie des chaînes de Markov et offre un large éventail d'implémentations possibles.
Je suis tombé sur cet article où il est dit que dans l'échantillonnage de Gibbs, chaque échantillon est accepté. Je suis un peu confus. Comment se fait-il que chaque échantillon qu'il a accepté converge vers une distribution stationnaire. En général, nous acceptons l'algorithme Metropolis comme min (1, p (x *) …
J'ai essayé de simuler à partir d'une densité bivariée utilisant les algorithmes de Metropolis dans R et je n'ai pas eu de chance. La densité peut être exprimée comme , où est la distribution de Singh-Maddalap ( x , y)p(X,y)p(x,y)p ( y| x)p(x)p(y|X)p(X)p(y|x)p(x)p ( x )p(X)p(x) p ( x ) …
Dans l'algorithme Metropolis – Hastings d'échantillonnage d'une distribution cible, supposons: πjeπi\pi_{i} soit la densité cible à l'état ,jeii πjπj\pi_j la densité cible à l'état proposé ,jjj hje jhijh_{ij} la densité de proposition pour la transition vers l'état étant donné l'état actuel ,jjjjeii uneje jaija_{ij} soit la probabilité d'acceptation de l'état …
J'ai essayé d'avoir une idée des différents problèmes dans les environnements fréquentistes où MCMC est utilisé. Je sais que MCMC (ou Monte Carlo) est utilisé pour ajuster les GLMM et peut-être dans les algorithmes EM Monte Carlo. Y a-t-il des problèmes plus fréquentistes lorsque MCMC est utilisé?
Dans la chaîne de Metropolis-Hastings Markov Monte Carlo, la distribution de la proposition peut être n'importe quoi, y compris le gaussien (selon Wikipedia). Q: Quelle est la motivation pour utiliser autre chose que le gaussien? Le gaussien fonctionne, c'est facile à évaluer, c'est rapide et tout le monde le comprend. …
Dans une implémentation MCMC de modèles hiérarchiques, avec des effets aléatoires normaux et un Wishart antérieur pour leur matrice de covariance, l'échantillonnage de Gibbs est généralement utilisé. Cependant, si nous changeons la distribution des effets aléatoires (par exemple, vers Student's-t ou un autre), la conjugaison est perdue. Dans ce cas, …
Il existe plusieurs versions des algorithmes adaptatifs de Metropolis Hastings. L'une est implémentée dans la fonction Metro_Hastingsde Rpackage MHadaptive, voir ici . La référence qui y figure, Spiegelhalter et al. (2002), ne contient malheureusement pas de description d'aucun algorithme adaptatif, à ma connaissance. Cependant, l' Metro_Hastingsalgorithme fonctionne très bien dans …
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.