Statistics for A multi-model approach to X-ray pulsars: Connecting spectral and timing models to pin down the intrinsic emission characteristics of magnetized, accreting neutron stars