running-bayesian-workflow

Diagnose Bayesian regression and hierarchical-model workflows with PyMC, NumPyro, Stan, brms, and rstanarm.

2|Updated May 23, 2026
One-click install
npx skills add https://github.com/rocklambros/rcs --skill running-bayesian-workflow
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: running-bayesian-workflow
Source: https://github.com/rocklambros/rcs/tree/main/skills/ml-datasci/running-bayesian-workflow
Command: npx skills add https://github.com/rocklambros/rcs --skill running-bayesian-workflow

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill includes references (resource) components.

What problem does it solve?

It helps you fit Bayesian regression and hierarchical models without shipping unreliable results, by forcing a disciplined workflow for priors, sampling, diagnostics, and uncertainty reporting.

Core Features & Use Cases

  • Prior discipline: Choose weakly informative priors that are scaled to the data instead of vague defaults that cause funnels and divergences.
  • End-to-end validation: Run prior-predictive checks, NUTS sampling, posterior-predictive checks, and model comparison with LOO or WAIC.
  • Diagnostic gatekeeping: Stop interpretation until r_hat, ESS bulk, ESS tail, divergences, and E-BFMI all pass threshold.
  • Use case: A data scientist fitting a hierarchical logistic regression in PyMC can use this Skill to fix divergences, reparameterize safely, and report credible intervals only after the fit is clean.

Quick Start

Use the running-bayesian-workflow skill to diagnose my Bayesian model, tighten the priors, verify the sampling diagnostics, and report credible intervals only after the posterior passes all checks.

Frequently Asked Questions about running-bayesian-workflow

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I fix divergences when fitting hierarchical models in PyMC?

Resolve divergences by applying non-centered parameterizations to hierarchical models and scaling weakly informative priors to your data, which eliminates the funnel dynamics that destabilize NUTS sampling.

What Bayesian model diagnostics are required before interpreting posterior results?

Required Bayesian diagnostics include checking r_hat, ESS bulk, ESS tail, divergences, and E-BFMI, which must all pass their threshold limits to form a hard diagnostic gate before any credible interval interpretation.

How do I compare Bayesian models using LOO or WAIC?

Compare Bayesian models using LOO or WAIC by running posterior-predictive validation after NUTS sampling, which provides the standardized metrics needed to evaluate relative predictive accuracy across different model fits.

Does this Bayesian workflow support Stan and NumPyro fits?

Yes, the Bayesian workflow supports Stan, NumPyro, brms, and rstanarm fits, applying prior-predictive checks, NUTS sampling diagnostics, and credible interval reporting uniformly across these different frameworks.

Why do vague default priors cause issues in Bayesian regression?

Vague default priors cause issues in Bayesian regression because they induce funnel geometries and divergences during NUTS sampling, whereas weakly informative priors scaled to your data stabilize the posterior exploration.

What is the best way to validate posterior-predictive checks for hierarchical models?

The best way to validate posterior-predictive checks is to enforce a hard diagnostic gate on ESS bulk, ESS tail, and divergences first, ensuring the posterior distribution is clean before comparing simulated data against observed outcomes.