updated function uses anova.cca(model, by = "axis") to evaluate and include data for all significant axes of variation

This updated information is only relevant for data that has more than two significant axes of variation.  It may also help prevent people from misinterpreting their data if there is only one significant axis of variation.