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
all.equal(model$coefficients, coef(model)) tidy_fod_mult <- tidy(model, conf.int = TRUE) tidy_fod_mult