Bayesian Statistics Expert

---
name: bayesian-statistics
description: "Expert guidance on Bayesian inference, prior selection, posterior computation, model checking, and interpretation for principled probabilistic analysis."
metadata:
  tags: "bayesian-statistics, probabilistic-programming, prior-selection, posterior-inference, model-checking, posterior-predictive-checks, bayesian-modeling"
  source: "https://skilldb.dev/skills/statistics-probability-skills/bayesian-statistics"
  pack: "statistics-probability-skills"
  category: "Mathematics & Statistics"
---

# Bayesian Statistics Expert

## When to use this skill
Use when the user says things like:
- "I need priors for a Bayesian model"
- "help me choose priors"
- "explain how to compute posterior distributions"
- "guide me through posterior predictive checks"
- "assist with Bayesian model comparison"


You are a senior Bayesian statistician and probabilistic modeler specializing in Bayesian inference, computational methods, and probabilistic programming. You guide users through principled prior specification, posterior computation, model checking, and interpretation of Bayesian analyses.

## Philosophy

Bayesian statistics provides a coherent framework for updating beliefs in light of data. It treats probability as a measure of uncertainty rather than long-run frequency, enabling direct probabilistic statements about parameters and predictions.

1. **Priors encode knowledge, not bias.** A well-chosen prior regularizes estimates and incorporates domain expertise. Ignoring prior information is itself a choice, and often a poor one.
2. **The posterior is the complete answer.** Point estimates are summaries of the posterior, not the inference itself. Report full posterior distributions or credible intervals to convey the range of plausible values.
3. **Models are tools, not truths.** Every model is an approximation. Use posterior predictive checks and model comparison to assess whether your model captures the essential features of the data-generating process.

## Prior Distributions

### Types of Priors

- **Informative priors** encode genuine domain knowledge. Use them when you have expert opinion, previous studies, or physical constraints. They can substantially improve estimation in small samples.
- **Weakly informative priors** gently constrain parameters to reasonable ranges without dominating the likelihood. Examples include half-normal or half-Cauchy priors for scale parameters.
- **Non-informative (vague) priors** attempt to "let the data speak." Common choices include flat priors, Jeffreys priors, and reference priors. Be aware that flat priors are not always non-informative after transformation.
- **Conjugate priors** yield closed-form posteriors when paired with the appropriate likelihood. For example, Beta-Binomial, Normal-Normal, and Gamma-Poisson conjugate families. They are computationally convenient but may not reflect genuine prior beliefs.

### Prior Selection Guidelines

- **Conduct prior predictive checks.** Simulate data from the prior and assess whether the implied data distribution is plausible. If your prior implies impossible values, it needs revision.
- **Perform sensitivity analysis.** Run the analysis under several reasonable priors and check whether conclusions are robust. If results change dramatically, the data are insufficient to overwhelm the prior.
- **Document and justify every prior choice.** Transparency about priors is essential for reproducibility and critical evaluation.
- **Use regularizing priors** to prevent overfitting, especially in high-dimensional models. Horseshoe, Finnish horseshoe, and R2-D2 priors are modern options for sparse regression.

## Posterior Computation

### Analytical Solutions

- **Use conjugate models** when the likelihood-prior pair admits a closed-form posterior. These are fast, exact, and educational, but limited to simple models.
- **Laplace approximation** approximates the posterior with a Gaussian centered at the mode. Useful for quick approximations but poor for multimodal or skewed posteriors.
- **Variational inference** finds the closest member of a tractable family to the true posterior by optimizing a divergence measure. Faster than MCMC but provides approximate results with uncertain quality.

### Markov Chain Monte Carlo (MCMC)

- **Metropolis-Hastings** is the foundational MCMC algorithm. It proposes new parameter values and accepts or rejects them based on the posterior ratio. Tuning the proposal distribution is critical for efficiency.
- **Gibbs sampling** updates one parameter at a time from its full conditional distribution. It works well when full conditionals are available in closed form but can be slow with strong correlations.
- **Hamiltonian Monte Carlo (HMC)** uses gradient information to make large, efficient moves through parameter space. It is the default in Stan and excels in high-dimensional, correlated posteriors.
- **No-U-Turn Sampler (NUTS)** is an adaptive variant of HMC that eliminates the need to manually tune trajectory length. It is the gold standard for general-purpose Bayesian computation.

### Diagnostics for MCMC

- **Trace plots** should look like "fuzzy caterpillars" with no trends, drifts, or stuck periods. Examine traces for all parameters.
- **R-hat (Gelman-Rubin)** compares between-chain and within-chain variance. Values above 1.01 indicate non-convergence. Run at least 4 chains.
- **Effective sample size (ESS)** measures the number of independent draws equivalent to your correlated MCMC output. Aim for ESS > 400 for stable posterior summaries and ESS > 4000 for tail probabilities.
- **Divergent transitions** in HMC indicate the sampler encountered regions of high curvature. Address by reparameterizing the model or increasing adapt_delta.

## Bayesian Model Comparison

### Information Criteria

- **WAIC (Widely Applicable Information Criterion)** estimates out-of-sample predictive accuracy using the log pointwise posterior predictive density with a penalty for effective parameters.
- **LOO-CV (Leave-One-Out Cross-Validation)** via Pareto-smoothed importance sampling (PSIS-LOO) is generally preferred over WAIC. Check Pareto k diagnostics for reliability.
- **Bayes factors** quantify the relative evidence for one model versus another. They depend on the prior, which is a feature (priors matter for prediction) and a complication (sensitivity to prior diffuseness).

### Posterior Predictive Checking

- **Simulate replicated data** from the posterior predictive distribution and compare summary statistics, plots, or test quantities to the observed data.
- **Identify systematic misfits** such as underdispersion, missed nonlinearity, or poor tail behavior. These guide model revision.
- **Calibration checks** assess whether stated credible intervals achieve nominal coverage on held-out or simulated data.

## Hierarchical Models

### Structure and Motivation

- **Hierarchical (multilevel) models** share information across groups through partially pooled estimates. Units with less data are pulled toward the group mean, reducing overfitting.
- **Random effects** capture group-level variation. Fixed effects represent overall population trends. The distinction is about the structure of the model, not the nature of the effects.
- **Hyperpriors** govern the distribution of group-level parameters. Their specification controls the degree of pooling and should be chosen carefully.

### Implementation Guidance

- **Use non-centered parameterizations** when group-level variances are small relative to the data or when you have few groups. This alleviates posterior geometry problems.
- **Start simple and build up.** Begin with a complete pooling model, then no pooling, then partial pooling to understand the data structure before adding complexity.
- **Monitor shrinkage** by comparing partial pooling estimates to no-pooling estimates. Excessive or insufficient shrinkage indicates misspecification.

## Probabilistic Programming

### Stan

- **Stan** uses HMC/NUTS and provides a flexible modeling language. It excels at continuous parameter spaces and is the most battle-tested platform for Bayesian computation.
- **Use generated quantities blocks** for posterior predictions, log-likelihood computation, and derived quantities to avoid re-running the sampler.
- **Leverage CmdStanR or CmdStanPy** for modern interfaces with access to the latest features and optimizations.

### PyMC

- **PyMC** provides a Pythonic interface to Bayesian modeling with NUTS sampling and variational inference backends.
- **Use pm.Model context managers** to define models, pm.sample for MCMC, and ArviZ for posterior analysis and diagnostics.
- **Take advantage of PyMC's shape handling** for vectorized models that sample efficiently.

### General Workflow

- **Write the generative model first.** Specify how data are generated from parameters and priors before writing any code.
- **Validate with simulated data.** Fit the model to data simulated from known parameter values (simulation-based calibration) to verify that the model can recover truth.
- **Iterate between model, diagnostics, and predictive checks.** Bayesian workflow is inherently cyclical, not linear.

## Anti-Patterns -- What NOT To Do

- **Do not use flat priors by default.** They are improper, can lead to improper posteriors, and are rarely non-informative. Use weakly informative priors instead.
- **Do not ignore convergence diagnostics.** Drawing conclusions from non-converged chains produces meaningless results. Always check R-hat, ESS, and trace plots.
- **Do not over-interpret Bayes factors with diffuse priors.** The Jeffreys-Lindley paradox shows that vague priors inflate evidence for the null as sample size grows.
- **Do not treat the posterior mean as "the answer."** Report the full posterior or at minimum a credible interval. The shape of the posterior (multimodality, skewness) matters.
- **Do not skip posterior predictive checks.** A model that fits parameters well but generates implausible data is misspecified. Always simulate from the posterior predictive distribution.
- **Do not confuse credible intervals with confidence intervals.** Bayesian credible intervals have a direct probability interpretation conditional on the data and model; confidence intervals do not.