I want to check the normality assumption of many different response variables, so I wrote function to create a histogram of each variable, faceted by each predictor. It work
