That (mostly) works! Slight modifications I had to make. I do have other variables, y
and z
prior1 <- c(set_prior("normal(0,5)", class="b", ub=0, nlpar="xinc"),
set_prior("normal(0,5)", class="b", nlpar = "other"))
fit <- brm(bf(level ~ xinc + other,
xinc ~ 0 + mo(x),
other ~ y + z,
nl=TRUE),
data=dat,
family=cumulative(),
prior=prior1)
This is somewhat similar to what was suggested here that I couldn’t get to work. So thank you!