Choose your reading route
eduKateSG · Why Mathematics? · Energy Diagnostics and E-BFMI for HMC Reliability Models: Kinetic Energy, Marginal Energy, Exploration and Failure Warnings
Start with the question, audit the worked example, or jump to the assumptions and practice. Every chapter stays open for uninterrupted reading.
Understand the question
Define the quantity and follow the mechanism.Route 2 · Chapters 6–10
Reproduce the example
Audit the numbers one step at a time.Route 3 · Chapters 11–15
Check assumptions
See where design and calibration matter.Route 4 · Chapters 16–20
Explore applications
Connect the method to a concrete investigation.Route 5 · Chapters 21–25
Study and verify
Practise, inspect failure cases and read further.
Full chapter index · Practice and checking · Mathematics Learning Hub
The importance of mathematics becomes visible when Hamiltonian Monte Carlo must move between energy levels of a Bayesian reliability posterior. At each transition, momentum refreshment supplies kinetic energy. If that energy distribution cannot efficiently move through the marginal energy distribution of the target, the chain can explore poorly even when trace plots look busy. Energy diagnostics and the empirical Bayesian fraction of missing information, E-BFMI, help reveal this mismatch.
For a teaching calculation, suppose the mean squared change in total energy between successive draws is 1.8 and the sample variance of total energy is 6.0. Using the common empirical ratio, E-BFMI=1.8/6.0=0.30. Software conventions and finite-sample details should be checked; a value near a warning threshold is a prompt to investigate energy histograms, divergences, parameterisation and effective sample size rather than a universal pass–fail verdict.
Full chapter index
Chapters 1–5 · Understand the question
Chapters 6–10 · Reproduce the example
Chapters 11–15 · Check assumptions
Chapters 16–20 · Explore applications
Chapters 21–25 · Study and verify
Begin with total Hamiltonian energy. Total energy combines potential energy from the negative log density with kinetic energy from momentum. A reliability draw has position and auxiliary momentum. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Energy is algorithmic, not physical component energy. Keep meanings separate.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Save per-chain energy values, squared transitions, variance convention and software version. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use hamiltonian energy as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Did You Know? Energy diagnostics can flag weak exploration even when parameter trace plots look active and all chains have similar means. That is why a careful report includes the construction, not only the final scalar.
Now compare refreshed energy changes with marginal energy spread. The posterior density determines potential energy up to an additive constant. Low posterior density corresponds to high potential energy. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Improper targets make the geometry ill-defined. Verify the posterior.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Inspect energy histograms alongside the scalar empirical ratio. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use potential energy as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
A low ratio warns that exploration may be incomplete. The momentum distribution supplies kinetic energy at each refresh. A Gaussian momentum uses the mass matrix. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Poor metric choice shapes available energy changes. Inspect adaptation.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Keep E-BFMI distinct from model fit, ESS and divergence counts. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use kinetic energy as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Begin with total Hamiltonian energy. Integrating over position and momentum gives the energy levels the chain must explore. A hierarchical funnel has a broad energy distribution. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Local movement can miss remote levels. Compare energy histograms.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Investigate metric, step size, centring, priors and weak identification. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use marginal energy distribution as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Primary research: https://arxiv.org/abs/1604.00695. Read the authors’ definition and conditions before carrying the method into a censored reliability model.
Now compare refreshed energy changes with marginal energy spread. Successive total energies show how momentum refreshment moves the chain across levels. Small changes repeat for many iterations. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Busy parameter traces can still have weak energy movement. Plot energy differences.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Recompute the ratio and preserve every diagnostic column. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use energy transitions as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
A low ratio warns that exploration may be incomplete. Average squared differences between successive energy values. The teaching mean squared change is 1.8. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Autocorrelation and finite length affect the estimate. Use post-warmup draws.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Save per-chain energy values, squared transitions, variance convention and software version. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use empirical numerator as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Did You Know? Energy diagnostics can flag weak exploration even when parameter trace plots look active and all chains have similar means. That is why a careful report includes the construction, not only the final scalar.
Begin with total Hamiltonian energy. Estimate the sample variance of total energy. The teaching variance is 6.0. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Different variance conventions can slightly change finite-sample ratios. State software convention.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Inspect energy histograms alongside the scalar empirical ratio. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use empirical denominator as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Now compare refreshed energy changes with marginal energy spread. Divide 1.8 by 6.0. The illustrative empirical ratio is 0.30. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. A rounded value near a threshold is not decisive alone. Inspect the full panel.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Keep E-BFMI distinct from model fit, ESS and divergence counts. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use worked e-bfmi 0.30 as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
| Quantity | Value |
|---|---|
| Mean squared energy change | 1.8 |
| Sample variance of energy | 6.0 |
| Illustrative E-BFMI | 0.30 |
| Interpretation | Investigate with the full diagnostic panel |
A low ratio warns that exploration may be incomplete. Low E-BFMI suggests momentum refreshment explores the marginal energy distribution inefficiently. The chain remains in a narrow range of energies. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. The diagnostic does not identify one unique cause. Investigate geometry.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Investigate metric, step size, centring, priors and weak identification. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use interpretive intuition as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Begin with total Hamiltonian energy. Compare the distribution of transition energy changes with the marginal energy distribution. One histogram is much narrower. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Histogram bins can affect appearance. Show several views.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Recompute the ratio and preserve every diagnostic column. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use energy histograms as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Now compare refreshed energy changes with marginal energy spread. Some software guidance uses values below roughly 0.3 as a warning. The example lies near that boundary. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. A threshold is heuristic and implementation-specific. Report value and context.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Save per-chain energy values, squared transitions, variance convention and software version. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use warning thresholds as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Did You Know? Energy diagnostics can flag weak exploration even when parameter trace plots look active and all chains have similar means. That is why a careful report includes the construction, not only the final scalar.
A low ratio warns that exploration may be incomplete. Energy mismatch and divergent transitions can occur together in difficult geometry. A funnel shows both warnings. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. One diagnostic can be good while the other fails. Check both.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Inspect energy histograms alongside the scalar empirical ratio. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use divergences as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Second authoritative source: https://mc-stan.org/docs/reference-manual/analysis.html. Compare its reference distribution, diagnostics and finite-simulation cautions with this worked example.
Begin with total Hamiltonian energy. Long NUTS trajectories may accompany poor energy exploration. Many draws reach maximum depth. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Depth and E-BFMI describe different mechanisms. Do not merge them.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Keep E-BFMI distinct from model fit, ESS and divergence counts. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use tree depth as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Now compare refreshed energy changes with marginal energy spread. Low energy efficiency can reduce ESS even when R-hat is near one. All chains explore the same narrow energy band. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Agreement does not prove completeness. Use energy diagnostics too.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Investigate metric, step size, centring, priors and weak identification. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use r-hat and ess as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
A low ratio warns that exploration may be incomplete. Scale parameters can create energy distributions hard for a fixed metric to traverse. Batch variance approaches zero. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. More draws may reproduce the same limitation. Try justified non-centring.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Recompute the ratio and preserve every diagnostic column. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use hierarchical funnels as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Begin with total Hamiltonian energy. Heavy tails broaden energy behaviour and challenge Euclidean HMC. Frailty scale has rare huge values. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. A finite run may underestimate tail difficulty. Use rank and tail checks.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Save per-chain energy values, squared transitions, variance convention and software version. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use heavy-tailed posteriors as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Did You Know? Energy diagnostics can flag weak exploration even when parameter trace plots look active and all chains have similar means. That is why a careful report includes the construction, not only the final scalar.
Now compare refreshed energy changes with marginal energy spread. The metric controls kinetic energy and coordinate scaling. A dense matrix handles correlated parameters. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Adaptation from a poorly explored warmup can be misleading. Use adequate warmup.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Inspect energy histograms alongside the scalar empirical ratio. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use mass matrix as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
A low ratio warns that exploration may be incomplete. Step size controls integration accuracy but not the entire energy-refresh mismatch. Smaller steps remove divergences but E-BFMI remains low. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Tuning epsilon alone may not solve exploration. Reparameterise.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Keep E-BFMI distinct from model fit, ESS and divergence counts. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use step size as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Begin with total Hamiltonian energy. Weak information creates broad ridges and varying curvature. Few failures leave shape and scale uncertain. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Computational tuning cannot create observations. Strengthen model and design.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Investigate metric, step size, centring, priors and weak identification. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use censored reliability as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Now compare refreshed energy changes with marginal energy spread. Coordinates that regularise geometry can improve energy movement. Non-centred batch effects raise E-BFMI. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Changing coordinates must preserve the model. Compare posterior predictions.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Recompute the ratio and preserve every diagnostic column. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use reparameterisation as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
A low ratio warns that exploration may be incomplete. Energy diagnostics should be examined for each chain. One chain has E-BFMI 0.18 and others 0.45. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. A pooled value can hide failure. Report the minimum and distribution.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Save per-chain energy values, squared transitions, variance convention and software version. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use per-chain reporting as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Did You Know? Energy diagnostics can flag weak exploration even when parameter trace plots look active and all chains have similar means. That is why a careful report includes the construction, not only the final scalar.
Begin with total Hamiltonian energy. Test energy diagnostics across simulated weak and strong data regimes. Sparse groups trigger warnings as expected. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Simulation from the model does not validate field realism. Add predictive checks.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Inspect energy histograms alongside the scalar empirical ratio. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use simulation benchmarks as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Now compare refreshed energy changes with marginal energy spread. Save energy values, step sizes, metrics, flags and estimator definition. A reviewer recomputes 1.8/6.0. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Only reporting a warning message prevents audit. Archive diagnostic columns.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Keep E-BFMI distinct from model fit, ESS and divergence counts. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use reproducibility as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
A low ratio warns that exploration may be incomplete. E-BFMI connects variance, squared differences, ratios, distributions and mechanics. Verify 1.8/6=0.30. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. The lesson is whether refreshed energy can traverse the target. Compare two histograms.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Investigate metric, step size, centring, priors and weak identification. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use parent and student learning route as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Useful next reading: Bukit Timah Mathematics curriculum overview, cross-subject learning approach, previous Why Mathematics article, Mathematics Learning Hub, eduKate Mathematics Learning System, approach to learning mathematics. Choose the destination that matches the next question rather than treating links as a checklist.
Begin with total Hamiltonian energy. E-BFMI is not model fit, posterior probability or a universal pass score. It complements divergences, R-hat, ESS and MCSE. This step matters because Hamiltonian energy movement, marginal exploration and diagnostic interpretation require the mathematical object, prediction target and reliability question to stay aligned. Good energy diagnostics cannot prove scientific validity. Read Betancourt on cotangent disintegrations.
A reproducible audit should record variables, units, observation boundaries, censoring codes, model version and software settings. Recompute the ratio and preserve every diagnostic column. Rebuild the displayed number from saved inputs before interpreting it. A result that cannot be reproduced is not ready to guide reliability work.
Try a small learning experiment: divide mean squared energy change 1.8 by energy variance 6.0 to obtain an illustrative E-BFMI of 0.30. Use common misconceptions and next reading as the focus. Predict the direction of change before using software, alter one defensible input, and finish with two sentences: what the calculation supports, and what remains uncertain. This habit turns mathematics education into transferable problem-solving rather than formula recall.
Quick answers
What must be fixed before checking a model or algorithm? Define the dataset unit, censoring mechanism, time scale, parameterisation, simulation design and diagnostic rule.
Does a satisfactory diagnostic prove the model is true? No. It supports a particular check under stated assumptions and simulation conditions.
Which mathematics matters? Probability defines the reference distribution, simulation makes it observable, proportions summarise tail events, variance describes uncertainty and graphs expose structured mismatch.
How should a student practise? Preserve intermediate quantities, predict directions before calculating, change one input at a time and explain the result in ordinary language.
What follows an internal check? Sensitivity analysis, independent computation, held-out or later-data prediction and a separate decision analysis when cost or safety matters.
Extended practice laboratory
Practice case 1: Hamiltonian energy
Setup: Use this chapter’s illustrative case: A reliability draw has position and auxiliary momentum.
Task: Keep meanings separate.
Checking target: Energy is algorithmic, not physical component energy.
Practice case 2: Potential energy
Setup: Use this chapter’s illustrative case: Low posterior density corresponds to high potential energy.
Task: Verify the posterior.
Checking target: Improper targets make the geometry ill-defined.
Practice case 3: Kinetic energy
Setup: Use this chapter’s illustrative case: A Gaussian momentum uses the mass matrix.
Task: Inspect adaptation.
Checking target: Poor metric choice shapes available energy changes.
Practice case 4: Marginal energy distribution
Setup: Use this chapter’s illustrative case: A hierarchical funnel has a broad energy distribution.
Task: Compare energy histograms.
Checking target: Local movement can miss remote levels.
Practice case 5: Energy transitions
Setup: Use this chapter’s illustrative case: Small changes repeat for many iterations.
Task: Plot energy differences.
Checking target: Busy parameter traces can still have weak energy movement.
Practice case 6: Empirical numerator
Setup: Use this chapter’s illustrative case: The teaching mean squared change is 1.8.
Task: Use post-warmup draws.
Checking target: Autocorrelation and finite length affect the estimate.
Practice case 7: Empirical denominator
Setup: Use this chapter’s illustrative case: The teaching variance is 6.0.
Task: State software convention.
Checking target: Different variance conventions can slightly change finite-sample ratios.
Practice case 8: Worked E-BFMI 0.30
Setup: Use this chapter’s illustrative case: The illustrative empirical ratio is 0.30.
Task: Inspect the full panel.
Checking target: A rounded value near a threshold is not decisive alone.
Practice case 9: Interpretive intuition
Setup: Use this chapter’s illustrative case: The chain remains in a narrow range of energies.
Task: Investigate geometry.
Checking target: The diagnostic does not identify one unique cause.
Practice case 10: Energy histograms
Setup: Use this chapter’s illustrative case: One histogram is much narrower.
Task: Show several views.
Checking target: Histogram bins can affect appearance.