> table(indvar_FODs$categorical_FOD_FODs) ascending descending identity 1683 1432 206 > table(indvar_FODs$Condition_FODs) A B C D 880 880 880 880 > table(indvar_FODs$Language_used_FODs) English German Hungarian Italian Turkish 880 160 160 160 2160 categorical_FOD_FODs is the dependent variable. please apply multinomial logistic regression. treat each level of the Condition_FODs variable (4 levels) and Language_used_FODs variable (5 levels) equally. include the interaction effect of variables in the analysis.

# multinomial logistic regression with interaction effect
model_mlr_interaction <- mlr(categorical_FOD_FODs ~ Condition_FODs * Language_used_FODs, data=indvar_FODs, family=multinomial("identity"))  
summary(model_mlr_interaction)