Is there a way in R to fit a mediation model with a binary mediator and a multicategorical (3) outcome

Question

So I\'ve been trying to fit the model mentioned in the title description. So far, I\'ve tried the package mediate, but it does not work. I know the model for m requires a lo