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. Apply bonferrini

data(iris)
iris <- within(iris, {Species <- factor(Species)})

model <- lm(Sepal.Length ~ Species * Sepal.Width, data = iris)

TukeyHSD(model, "Species")

TukeyHSD(model, "Sepal.Width")

TukeyHSD(model, "Species:Sepal.Width")