Expand description
ρ-posterior certificate / escalation DATA types (contract-down #1521).
These are the plain-data carriers that a fit result STORES
(UnifiedFitResult::rho_posterior_{certificate,escalation}) and that the
gam-solve REML evaluator returns. The COMPUTATION that produces them — the
PSIS certificate, the Tier-1 Gauss-Hermite quadrature, and the Tier-2 NUTS
escalation (which pulls the gam-inference hmc_io sampler) — stays UP in the
monolith inference::rho_posterior, which re-exports these types so its
construction sites name them unchanged. Contract-downed here (the neutral
criterion-contract crate) so gam-solve can store/return them without a
back-edge into gam-inference.
Structs§
- RhoMixture
Node - One node of the criterion-closure Tier-1 mixture (#938): a
ρlocation, its normalized posterior mass, and the exact profiled criterion value there. - RhoPosterior
Certificate - The Tier-0
ρ-uncertainty certificate for a fit. - RhoPosterior
Mixture - Tier-1 deliverable (#938):
π(ρ|y)as a discrete mixture of conditional Gaussians, with the posterior moment summary ofρitself. - RhoPosterior
Samples - Tier-2 deliverable (#938):
π(ρ|y)draws from NUTS with the exact profiled gradient, whitened by the exact outer Hessian atρ̂.
Enums§
- RhoCertificate
- Reliability tier read off the Pareto tail-shape
k̂of theρ-importance weights. - RhoPosterior
Escalation - The auto-selected escalation outcome when the Tier-0 certificate reads
RhoCertificate::Escalate(#938): Tier 1 (deterministic quadrature) forK ≤ 4, Tier 2 (NUTS overρ) forK ≤ 16, and an HONEST report that escalation is unavailable beyond that — never a silently-degraded answer.
Constants§
- ESCALATE_
K_ HAT k̂above which the self-normalized importance estimator’s central limit theorem no longer applies and the proposal must be replaced rather than reweighted (Vehtari, Simpson, Gelman, Yao & Gabry, Pareto smoothed importance sampling, JMLR 25(72), 2024, §3): fork > 0.7the practical pre-asymptotic convergence rate of PSIS collapses and the estimate is declared unreliable. How finely a given draw count can resolve this boundary is NOT a property of the cutoff:gam_solve::psis::shape_resolutionreports the standard error of the fitted shape, and a verdict against this cutoff only means something when the truth sits several such standard errors away from it.- PLUG_
IN_ CERTIFIED_ K_ HAT k̂below which the proposal’s importance weights have finite variance (k < 1/2⇒ the generalized-Pareto tail has a second moment), so the plug-in answer plus the first-order correction needs no reweighting.
Traits§
- RhoPosterior
Escalator - The gam-inference-tier producer of the Tier-0
ρ-certificate and the auto-selected Tier-1/Tier-2 escalation (trait-inversion #1521).
Functions§
- rho_
posterior_ escalator - The registered
ρ-posterior certificate/escalation producer, orNonewhen the sampler tier is not linked / not yet initialized (gam-solve then declines the certificate and escalation — a safe no-op leaving plug-in intervals). - set_
rho_ posterior_ escalator - Register the monolith’s
hmc_io-backedρ-posterior certificate/escalation producer. Called once at process init by the gam-inference tier. First writer wins; a later call is ignored (returnsErrwith the boxed value) so a re-init can never swap a live producer mid-run.