Hi all,

I have some multivariate data (species counts) and want to investigate the effect of some treatment over time.

I am using the mvabund package in R, which fits a separate GLM to each species.
Specifically, my data are invertebrate taxa collected from pitfall traps under two treatments. These traps have been repeatedly measured each summer over ~15 years. Thus my data are repeated measures.

I'm trying to figure out how to restrict permutations to account for this lack of independence between years using the permute package.

I can't theoretically figure out how to achieve this.

Thank you in advance for any advice.