In my master thesis, I'm investigating how the rating of the participants face recognition skill over a period of time of training. There are 15 sessions in total and each session there are 15 trials. Each trial the participant rate pictures of faces. There are correct and incorrect rating. The correct ratings will be summed together for each session.
I have carried out the following mixed model:
mod = lmer(sumright ~ session * IQc + (session|subject), data=df, na.action=na.omit)
- the dependent variable: sumright (continous scaled) ranges from 1,2,3... to 15. -> the participant can get max 15 correct answers per session
- independent variables: session (continous scaled) ranges from 1 to 15 -> there are 15 sessions in total; and
IQc (continous scaled) is centred -> IQc = IQ - mean(IQ)
There are 30 subjects, and 445 observations (rows of the dataframe).
We can clearly see that there are gaps between the residuals. I’m not sure whether there is a voilation of the homogenity assumption. I also checked whether the residuals are normal and they are. also the residual of the random effects (Intercept and session) are normaly distributed. The only thing is the residual plot, which I’m not sure whether it voilates the homogenity assumption. What do you think of the residual pattern above?
I also tried with a linear regression model. but the gaps between the residuals remain. So, very similar to the mixed model one.
As suggested in the comment, I run a CLMM and the problem described above is solved. The residualplot no longer has any gaps. So, the dependent variable has to be coded as a ordered factor.