Hey there,

I try to plot a multinomial regression. However, something does no work the way I wished...

*The formula is: mod.fak_diff_1996.7 <- multinom(factor(fac.num)~ factor(revenu), data=voxit_destill_1996)*

fac.num is a nominal variable with 3 values, revenu is a ordinal scale which is a measure for income. As I would like to visualise the probability of being part of fac.num(1) having an income value 1, 2, 3...

I tried with this plot, but there is something wrong/missing...

*pred.probs <- predict(mod.fak_diff_1996.7, type = "probs")*

*plot.probs.1996 <- ggplot(pred.probs, aes(x = revenu, color = level)) + eom_point(aes(y=mean)) + geom_errorbar(aes(ymin = lower, ymax = upper)) + theme_minimal() + ylab("Vorausgesagte Wahrscheinlichkeit") + xlab("Einkommen") + scale_y_continuous(labels = percent)*

Would be very greatful if someone could help me

Thanks,

Johanna