Top. MIRI consistency plot. Posterior predictive distribution $p(f(\theta )+ϵ|x_{\text{MIRI}})$ of noisy simulations spectra for the 99.7%, 95% and 68.7% quartiles (hues of blue) obtained from the retrieval on MIRI data alone extended to near-infrared wavelengths, and overlaid on top of the WFC3+GNIRS+MIRI observation x_obs (black line). Bottom. Residuals of the posterior predictive samples, normalized by the inflated standard deviation of the noise distribution for each spectral channel and a horizontal line at 0 for reference (in black).