Is there a way to use ggpredict and get the standard errors (or confidence intervals) for predicted probabilities of a multinomial logistic regression model using multinom
and ggpredict
? Thank you.
var1<-sample(c('A', 'B', 'C'), size=1000, replace=T)
var2<-rnorm(n=1000)
var3<-rnorm(n=1000)
df<-data.frame(var1, var2, var3)
library(nnet)
mod1<-multinom(var1~var2+var3, data=df)
library(ggeffects)
preds<-ggpredict(mod1, terms="var2")
data.frame(preds)