> chisq.test(table(indvar_FODs$Gender_FODs, indvar_FODs$categorical_FOD_FODs)) Pearson's Chi-squared test data: table(indvar_FODs$Gender_FODs, indvar_FODs$categorical_FOD_FODs) X-squared = 1.1716, df = 4, p-value = 0.8828 > chisq.test(table(indvar_FODs$Condition_FODs, indvar_FODs$categorical_FOD_FODs)) Pearson's Chi-squared test data: table(indvar_FODs$Condition_FODs, indvar_FODs$categorical_FOD_FODs) X-squared = 105.87, df = 6, p-value < 2.2e-16 analyse the interaction effect between the Condition_FODs and Gender_FODs on the categorical_FOD_FODs

model <- lm(categorical_FOD_FODs ~ Condition_FODs*Gender_FODs, data=indvar_FODs)
anova(model)