My biological question is the following: is the variance in behaviour similar in all regions studied, and if not, how does it change?
The variable 'behaviour' is a poisson-distributed count data, the experiment was reproduced in 4 different regions and the glmm I ran on the data gave a significant effect of region.
My first approach was to perform a single Fligner-killen test on the whole data set (fligner.test() in R) and then a series of post-hoc pairwise comparisons among regions with the same test. I used Bonferroni-adjusted signifificance threshold (alpha = 0.05/6 for six pairwise tests = 0.00833) but that looks really low to me, and probably too conservative.
I was also looking a bit over the coefficient of variation, but did not find any sound statistics to do with these.
Is there a better way to answer my question?