How to code up a model with discrete parameters in Stan

How to code up a model with discrete parameters in Stan

🎙 Ben Lambert 👥 148K 📅 14 novembre 2018 ⏱ 21 min 👁 8K 📄 tutoriel 🧭 2026-08-17
Disponible en : Français (actuel) English

Mots-clés

StanBayesianMCMCMarginalisationLog-sum-exp

Résumé

Cette vidéo de Ben Lambert, issue d’un cours sur la statistique bayésienne, explique comment utiliser Stan pour inférer des modèles contenant des paramètres discrets. Stan ne permet pas de déclarer directement des paramètres discrets, mais il est possible de les marginaliser dans la log-vraisemblance. L’auteur illustre la méthode avec un exemple concret : une série d’expériences de lancers de pièce où le nombre de lancers N (discret, uniforme entre 5 et 8) et la probabilité de pile θ (continue, uniforme entre 0 et 1) sont inconnus. Il détaille la dérivation mathématique : écrire la distribution jointe, marginaliser sur N, puis utiliser la fonction log_sum_exp pour une stabilité numérique. Il montre ensuite l’implémentation dans Stan, avec les blocs data, transformed data, parameters, transformed parameters, model et generated quantities. Dans le bloc transformed parameters, il calcule la log-probabilité non normalisée pour chaque valeur possible de N. Dans le bloc model, il ajoute à la cible le log_sum_exp de ces valeurs. Enfin, dans generated quantities, il calcule les probabilités postérieures normalisées de chaque valeur de N en utilisant l’astuce exponentielle et log_sum_exp. Il mentionne le bénéfice de la Rao-Blackwellisation. Les résultats montrent une forte probabilité postérieure pour N=5. La vidéo est claire et bien structurée, adaptée à un public ayant des bases en statistique bayésienne et en programmation Stan.

216 mots

Évaluation critique

Valeur des informations & solidité de l’argumentation

La valeur des informations est élevée : la vidéo fournit une méthode pratique et éprouvée pour traiter les paramètres discrets dans Stan, un problème courant en modélisation bayésienne. L’argumentation est solide : l’auteur explique chaque étape mathématique, justifie la nécessité de marginaliser, et montre comment implémenter la solution. Il utilise un exemple concret et reproductible, ce qui renforce la crédibilité. La démonstration est progressive et bien illustrée par le code.

Rigueur scientifique, qualité des sources, adéquation du titre

La rigueur scientifique est bonne : l’auteur est un enseignant reconnu en statistique bayésienne, et la méthode présentée est conforme aux pratiques standards. Les sources citées dans la description (site personnel et playlist) sont pertinentes pour approfondir. Le titre est exact et reflète le contenu. Aucune source externe n’est citée dans la vidéo elle-même, mais la méthode est bien établie.

145 mots

Adéquation titre / contenu

Le titre correspond parfaitement au contenu : la vidéo explique comment coder un modèle avec paramètres discrets dans Stan.

Qualité & fiabilité

8/10

Explication claire et pédagogique de la marginalisation des paramètres discrets dans Stan, avec code et exemples concrets. La méthode est standard et bien présentée, mais sans vérification indépendante des résultats.

Moments clés

Sources citées

Sources concordantes

  • Stan User's Guide — Documentation officielle de Stan, qui décrit la marginalisation des paramètres discrets comme une technique recommandée.

Apport & nouveautés

Cette vidéo apporte une explication claire et pratique de la marginalisation des paramètres discrets dans Stan, une technique essentielle pour les modèles bayésiens. Elle comble un manque en montrant pas à pas comment implémenter cette approche, avec un exemple concret et du code. L’accent mis sur la stabilité numérique via log_sum_exp est particulièrement utile.

Pour aller plus loin :

  • Stan User’s Guide — Documentation officielle de Stan, chapitre sur les paramètres discrets et la marginalisation.
  • Rao-Blackwell theorem — Concept mentionné dans la vidéo, amélioration de l’efficacité de l’estimation.
  • Log-sum-exp — Fonction mathématique utilisée pour la stabilité numérique.

97 mots

Profil radar

Le profil radar montre des scores élevés et équilibrés dans toutes les dimensions, indiquant une vidéo de qualité homogène : bonne quantité d'information, contenu fiable, niveau technique adapté et fiabilité globale solide.

Fiabilité 8/10