Please also provide the following information in addition to your question:

- Operating System: Windows 10
- brms Version: 2.4.3

Hi All,

I’m trying to understand what exactly brms is doing with the following model:

bf_1 <- bf(y|trials(n)~z+x)+binomial()

bf_2 <- bf(x~z)+gaussian()

brms.out <- brm(bf_1+bf_2+set_rescor(FALSE),data=dat,

chains = 2, iter = 1000,

warmup = 500, thin = 1,

control=list(adapt_delta=0.9999,max_treedepth=30))

I’d really like to be able to write down this model in mathematical notation. But I can’t figure out from the output how this model is any different than the following:

brms.out2 <- brm(bf_1,data=dat,

chains = 2, iter = 1000,

warmup = 500, thin = 1,

control=list(adapt_delta=0.9999,max_treedepth=30))

brms.out3 <- brm(bf_2,data=dat,

chains = 2, iter = 1000,

warmup = 500, thin = 1,

control=list(adapt_delta=0.9999,max_treedepth=30))

Any guidance?