I have a set of 9 different factor levels from my independent variable to be compared against each other. Here are the results of the different assumption tests in R. I'm just going through my methodology in the hope that someone can guide me to the right conclusion.
- The Shapiro-Wilk test is significant and shows that the my dependent variable is not normally distributed. I read somewhere that the assumption of normality doesn't always have to be met if the size of the dataset is substantial as it will almost always deviate from normality when it's the case.
- Mauchly's test for sphericity is significant, therefore the assumption of sphericity has been violated. My Greenhouse-Geiser correction confirmed it, but it's p-value is significant, which indicates that the levels of my independent variable significantly affects the dependent variable.
Where do I go from here?
The example I'm basing my methodology on in a textbook goes on to use these methods on the same dataset every time (though that dataset is different than mine)
- Repeated measures ANOVA
- A multilevel approach
- Robust test
Based on what I know, I think I have to go for a robust method given the assumptions above, but I have no idea if this is the right call.
lmer(). With regards to normality, what you are worried about in MLM is the normality of the residuals. You can easily look at these after running yourlmermodel. See https://www.ssc.wisc.edu/sscc/pubs/MM/MM_DiagInfer.html – Erik Ruzek Dec 29 '20 at 20:58I stumbled upon this article by Wilcox in which he writes that when normality and homoscedasticity are not met, classic inferential methods based on means (e.g., the ANOVA, F-test) are not the way to go, and you should switch to a robust method.
– bolleke Dec 30 '20 at 08:48