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