boxplot() will now interpret glms or bigglms as fitted (propensity) scores, and return plots of overlap on linear predictors.