added code to plot lat-richness correlations for all subclades with at least 5 species