adapted from plot_clusgap from phyloseq.
plot_clusgap(clusgap, title = "Gap Statistic results")
plot gap stats