model <- multinom(categorical_FOD_FODs ~ Condition_FODs + Language_used_FODs + Condition_FODs*Language_used_FODs, indvar_FODs) summary(model) write code for the post-hoc analysis, treating each level of the variables equally
posthoc_FODs <- glht(model, linfct = mcp(Condition_FODs = "Tukey", Language_used_FODs = "Tukey", Condition_FODs*Language_used_FODs = "Tukey")) summary(posthoc_FODs)