Split findVG to computation and plots

[hypercompetent]

Currently, findVG both computes statistics for variable genes and can optionally output plots to a file as a side effect. Instead, we should split this up to 3 functions: 1 computes the statistics 1 generates the plots 1 saves the plots (or leave this up to the user to use ggsave() or cowplot).