Produces a scatterplot matrix with histograms in diagonal, 2D kernel
density estimates and contours in the lower half and bivariate scatterplots
with lowess smooth curves and Pearson correlation values
in the upper half as a side effect.
Returns NULL invisibly.