Bayesian spectrum analysis for parameter estimation is a rigorous statistical (non-Fourier-based) method. Herein the Bayesian quadrature NMR model is introduced and applied to analysis of 31P NMR time domain data from in vivo rat brain. Immunity to both the brain spectrum "baseline hump" and the phase twist is demonstrated.