I am using the mvabund package to run the model-based analysis of multivariate abundance data.

My species abundance is expressed as Number.of.Individuals / 0.0625 m2. I have run the model to test the significance of factors -> anova (model1) and also to see which species are significantly contributing to the changes observed -> anova(model1, p.uni='adjusted').

As I wanted to present the data per m2 instead of per 0.0625 I have transformed all of my species abundances to densities (to obtain everything in m2). Then, I have re-run the model. And while the overall effect were the same, the p.uni test was very different ( much less species had a p