I want to run a mediation analysis with two mediators using brms. It is a multilevel model too. The dependent variable is coded as 1 or 0. There were 60 questions, and participants were randomly assigned 4 questions. I can’t use the mediation() function.