We want to determine if there are statistical differences in the structure of the microbial communities of a lake. For this, samples were taken, monthly, for 10 years, from two depths and considering two fractions (particle associated and free-living). In this way, we have a Species-Sites table where in the columns we have the OTUs and in the rows, for each date, the depth-fraction combinations. Normaly we perform the statistical analysis using the function adonis of package Vegan (Permanova). However, we would like to know if it is correct to assume that each one of the dates as repetitions? Is there any other test e.g. mixed model or LASSO that can determine significance in the differences between depths or fractions considering the time series? Any suggestion on how to perform this analysis will be appreciated.