Hmc Checker Patched «Linux»

# 5. BFMI if hasattr(inference_data, "sample_stats") and hasattr(inference_data.sample_stats, "bfmi"): bfmi = inference_data.sample_stats.bfmi.values.mean() if bfmi < bfmi_threshold: results["failures"].append(f"BFMI = {bfmi:.2f} < {bfmi_threshold}") results["passed"] = False

# For demo, create dummy data import pymc as pm with pm.Model(): x = pm.Normal("x") trace = pm.sample(1000, chains=2, return_inferencedata=True) hmc checker

return results if name == " main ": # Load your inference data # idata = az.from_netcdf("my_model.nc") # 5. BFMI if hasattr(inference_data