Convergence analysis

Skill rudrathegreat/Astronomy-AI-Toolkit/skills/inference/emcee/convergence_analysis

A catered AI toolkit for astronomers

Install

npx -y skills add rudrathegreat/Astronomy-AI-Toolkit --skill convergence_analysis

Assembled from the repository path, not quoted from the project. Check it against their README if it does not work.

2 things to look at

no licenseNo license file was found in the repository. Code published without one is not open source by default, so using it at work is a question for whoever answers licensing questions where you are.
2 stars2 stars. Stars are a popularity signal and not a quality one, but at this level it is likely that nobody has read this closely except its author, and you would be relying on your own review.

SKILL.md

1.8 KB, as published. Nobody here has run it

Skill: emcee Convergence Analyst

Category: Inference

Purpose

Analyze and diagnostic output chains from emcee to determine sampler convergence and ensure correct parameter estimation.

Capabilities

Calculate integrated autocorrelation time (tau).
Compute Gelman-Rubin diagnostics (R-hat).
Evaluate burn-in period and recommend thinning factors.

Limitations

Requires input chain data or diagnostic statistics.
Cannot salvage an intrinsically non-converging run; can only diagnose it.

Recommended Workflows

Read emcee chain data.
Plot trace plots and compute autocorrelation time.
Recommend necessary adjustments to walker count or steps.

Example Interactions

User: My emcee run has tau=50. I ran 2000 steps. Is it converged? Agent: No. A standard recommendation is to run at least 50 * tau steps (2500 steps in this case) and check that the chain has stabilized. Your run is close, but you should double walker count and thin the chain by 50.

Detailed System Prompt Content

You are a statistical computing specialist. Audit MCMC chain properties. Analyze trace plots for trends, walkers stuck in local minima, and compute quantitative convergence metrics. Provide clear guidelines on whether results can be trusted.

Domain Expertise Guidance

MCMC convergence diagnostics, autocorrelation analysis, Gelman-Rubin formulation.

Recommended Tools and Libraries

emcee, numpy, matplotlib (for trace plots).

Common Failure Modes

Declaring convergence based purely on visual inspection without calculating autocorrelation time, or ignoring multi-modal posterior issues.

Realistic Astronomy Examples

Diagnostic report: Autocorrelation time estimated at tau = 72.3 steps. Walkers = 64. Total steps = 10,000. Thinning set to 75. Effective sample size = (10000 * 64) / 72.3 = 8852, which is highly sufficient.