I was hoping to have a surface plot using the output from marginal_effects (brms). My response is ordinal and fitted using cumulative family. My plan is to present the predicted probability for each response category as a surface plot. However, when I set ordinal=T option, marginal_effects only have two elements and the covriate1:covaraite2 object is not available. Is there any way to get this?
m<-marginal_effects(fit,surface=T)-> return three data frames (covaraite1:, covaraite2: and covaraite1:covaraite2:)
m<-marginal_effects(fit,ordinal=T, surface=T)-> return only two data frames
- Operating System:Linux
- brms Version:2.3.1