How to code up a bespoke probability density in Stan

How to code up a bespoke probability density in Stan

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

Mots-clés

Stanlog-probabilitédistribution personnaliséeHamiltonian Monte Carloinférence bayésienne

Résumé

Cette vidéo de Ben Lambert explique comment utiliser Stan pour échantillonner à partir d’une distribution de probabilité qui n’est pas incluse dans la bibliothèque standard de Stan. L’auteur commence par rappeler que Stan fonctionne sur la log-probabilité du modèle, et que toute distribution peut être utilisée si l’on peut exprimer sa log-densité. Il illustre cela avec un exemple simple : une distribution de la forme P(x|a) = (1-a) * x^(-a), avec a entre 0 et 1 et x entre 0 et 1. Il montre comment coder cette distribution dans Stan en utilisant la fonction ’target’ pour incrémenter la log-probabilité, et comment définir une fonction personnalisée dans le bloc ‘functions’. Il explique également comment utiliser le suffixe ‘_lpdf’ pour traiter la distribution comme une fonction de densité standard. Enfin, il mentionne que les constantes peuvent être omises de la log-probabilité sans affecter l’inférence. La vidéo inclut une démonstration pratique avec R et Stan, montrant comment générer des données et inférer le paramètre a.

162 mots

Évaluation critique

Valeur des informations & solidité de l’argumentation

La vidéo apporte une valeur pédagogique certaine en expliquant un aspect avancé de Stan : la création de distributions personnalisées. L’argumentation est solide, basée sur des principes mathématiques clairs et illustrée par un exemple concret. L’auteur justifie chaque étape et montre comment la méthode s’intègre dans le fonctionnement de Stan. La démonstration pratique renforce la crédibilité de l’explication.

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

La rigueur scientifique est bonne : les explications sont correctes et conformes à la documentation de Stan. L’auteur mentionne son livre et son site pour approfondir, mais ne cite pas de sources externes spécifiques dans la vidéo. Le titre est en adéquation parfaite avec le contenu. Aucun commentaire n’est fourni pour analyser les tendances du public.

129 mots

Adéquation titre / contenu

Le titre correspond parfaitement au contenu : il explique comment coder une densité de probabilité personnalisée dans Stan.

Qualité & fiabilité

8/10

Explication claire et pédagogique, avec exemples concrets et démonstration pratique. La méthode est correcte et conforme aux principes de Stan. Quelques approximations (équivalence des formulations) mais sans impact majeur.

Moments clés

Sources citées

Sources concordantes

  • Documentation de Stan sur les fonctions utilisateur — Confirme la méthode de définition de fonctions personnalisées dans Stan.

Apport & nouveautés

Cette vidéo apporte une explication claire et pratique de la création de distributions personnalisées dans Stan, une fonctionnalité peu documentée. Elle montre comment utiliser la fonction target et le bloc functions pour définir une log-densité, et comment l’utiliser dans un modèle. L’originalité réside dans la démonstration pas à pas avec un exemple concret.

Pour aller plus loin :

  • Documentation de Stan sur les fonctions utilisateur — Documentation officielle sur la définition de fonctions personnalisées dans Stan.
  • Hamiltonian Monte Carlo — Article Wikipédia sur l’algorithme HMC, base de Stan.
  • No-U-Turn Sampler — Article de Hoffman et Gelman présentant l’algorithme NUTS utilisé par Stan.

102 mots

Profil radar

Le profil radar montre des scores élevés en qualité et fiabilité, avec un niveau technique modéré, indiquant une vidéo pédagogique fiable mais nécessitant des bases en statistiques bayésiennes.

Fiabilité 8/10