I am working with ecological count data in order to analyze differences/any contrast in species composition between warm and cold year communities. The abundances of species were recorded from multiple sites, with two trips per site. For example in September of 1999, species were collected/recorded twice from Site 1. I have columns separated into years, site, trip (replicates), and group (warm or cold). I need to test if there is any significance in counts due to site effect or replicates within sites (the trips). Basic ANOVA was what was going to be used. This is essentially how I am attempting to setup the model:
aov(Count ~ Site * Trip + Years, data)
Like most count data it is not normal and has lot of zeros so its over-dispersed and negative binomial will be used as well (using the mvabund
package). Some of my other parts of the analysis requires data to be fourth root transformed due to dominant species. But I am unsure if data should be transformed for ANOVA and GLM as well? When an ANOVA was performed with transformed data, there appeared to be significant $p$ values for the site and the interaction categories, but was neither was significant when applied on just raw counts. I'm just unsure if transformations are appropriate for this analysis.