From: Fitting parametric random effects models in very large data sets with application to VHA national data