Updated EDA code with function for producing a multivariate Q-Q plot of Mahalanobis distance (based on the corresponding function in package 'semutils', but simplified and ported to use 'ggplot2'). TODO: add a call within EDA code to include building this plot into MVN assessment.