PyMC Bayesian Modeling
Guides Claude through the complete PyMC 5 Bayesian workflow — priors, NUTS sampling, convergence diagnostics, LOO model comparison, and prediction.
Data & AnalyticsAdvanced★ 367⑂ 36AI score 9/10Last updated: Aug 29, 2026
What it does
Turns Claude into a disciplined Bayesian modeling assistant built around PyMC 5.
- 8-step workflow: standardize data → define model and priors → prior predictive check → NUTS sampling → convergence diagnostics → posterior predictive check → LOO/WAIC comparison → prediction
- Diagnostic thresholds: R-hat < 1.01, ESS > 400, zero divergences, Pareto-k < 0.7, with prescribed fixes for each failure
- Prior/distribution selection table: Normal, HalfNormal, StudentT, Beta, Dirichlet, NegativeBinomial and when each applies
- Ready recipes: hierarchical models with non-centered parameterization, logistic regression, Gaussian Processes
- Troubleshooting table: divergent transitions, low ESS, slow sampling,
logp = -inf, max tree depth
Who it's for
- Researchers and analysts who need full uncertainty quantification (HDIs), not just point estimates
- Anyone modeling grouped or nested data (patients within hospitals, students within schools)
- Biostatistics, clinical trial, psychometrics, and Bayesian A/B testing practitioners
- PyMC users who get stuck interpreting diagnostics or comparing candidate models
Examples
- "Fit a hierarchical model of recovery time across 20 hospitals" → produces non-centered parameterization code plus an
az.summarywalkthrough - "My sampler reported 150 divergences" → raises
target_acceptto 0.95, suggests reparameterization and tighter priors, then re-checks R-hat/ESS - "Which is better, the linear or quadratic model?" → refits with
log_likelihood=True, runsaz.compare(ic="loo")and validates with Pareto-k diagnostics
· · · Install guide · · ·
Try it now, no install
Paste this into Claude to use the skill without installing anything.
Read the instructions in this file and follow them to help me: https://raw.githubusercontent.com/jaechang-hits/SciAgent-Skills/HEAD/skills/biostatistics/pymc-bayesian-modeling/SKILL.md What I want: (describe your task here)
If Claude can't open the link, open it yourself and paste the contents instead.
↓ If it works for you, download the ZIP below and install it. Then it runs on its own — no pasting each time.
Install in the Claude app (no terminal)
- Download the ZIP with the button below.
- In Claude, open Settings → Capabilities and turn on 'Code execution and file creation'. (one time)
- Go to Customize → Skills → + → 'Upload a skill' and upload the ZIP.
Install in Claude Code
Let Claude do it — paste this into Claude Code
Install the skill I found on Claude Skill Mart. Copy the skills/biostatistics/pymc-bayesian-modeling folder from the GitHub repo jaechang-hits/SciAgent-Skills into my ~/.claude/skills/pymc-bayesian-modeling/. When it's done, tell me in one line what this skill can do.
Install with a command instead
git clone https://github.com/jaechang-hits/SciAgent-Skills.git && mkdir -p ~/.claude/skills && cp -r SciAgent-Skills/skills/biostatistics/pymc-bayesian-modeling ~/.claude/skills/⚠ This is a third-party skill. Check the source repository before installing.
- Open your terminal.
- Clone the repository:
git clone https://github.com/jaechang-hits/SciAgent-Skills.git - Create the skills directory:
mkdir -p ~/.claude/skills - Copy this skill:
cp -r SciAgent-Skills/skills/biostatistics/pymc-bayesian-modeling ~/.claude/skills/ - Install the Python dependencies:
pip install "pymc>=5" arviz numpy matplotlib - Restart Claude Code and run
/skillsto confirmpymc-bayesian-modelingis listed. - Try a prompt such as "Build a Bayesian logistic regression in PyMC for this dataset and check convergence."
View source on GitHub ↗License: NOASSERTION