model1 <- lme(categorical_FOD_FODs~Condition_FODs*Gender_FODs,random=~1|ID/Ch,na.action=na.omit, method = "ML", data=long1.RSPL,control = lmeControl(opt = "optim"))
We used linear mixed effect model to analyse the interaction effect of condition and gender on the categorical FODs data.