# NOT RUN {
fit <- brm(time | cens(censored) ~ age * sex + disease + (1|patient),
data = kidney, family = lognormal())
summary(fit)
# summarize via 'emmeans'
library(emmeans)
rg <- ref_grid(fit)
em <- emmeans(rg, "disease")
summary(em, point.est = mean)
epred <- emmeans(fit, "disease", dpar = "mean")
summary(epred, point.est = mean)
# }
Run the code above in your browser using DataLab