我正在努力分析数据，可能需要一些帮助。.我正在研究am割对大猩猩社会行为发展的影响。我观察到11位男性，其中4位被cast割，另外7位完好无损。我已经观察了他们四年，因此对于每个男性，我们观察到四个数据点：2015, 2016, 2017 and 2018。我想比较cast割和完整的雄性多年来的发展，但也要考虑到年龄和动物园的可能影响（因为雄性生活在不同的动物园中）。
到目前为止，我认为我将使用以下代码对ANOVA进行重复测量：
lme_aggression = lme(Aggression ~ Group + Age + Zoo, data=aggressionlong, random = ~1|Focal)
anova(lme_aggression)
summary(lme_aggression)
我真的不知道那么多统计数据，所以我想知道这是否正确？
如果是，那么下一个问题是，对于许多行为而言，残差的正态性假设无法满足，并且使用Friedman检验作为ANOVA的替代方法在R中不起作用...
感谢您的帮助！！ 谢谢。