# \donttest{
## fit a model
if(requireNamespace("cmdstanr")){
m <- brmcoda(complr = complr(data = mcompd, sbp = sbp,
parts = c("TST", "WAKE", "MVPA", "LPA", "SB"),
idvar = "ID", total = 1440),
formula = Stress ~ bz1_1 + bz2_1 + bz3_1 + bz4_1 +
wz1_1 + wz2_1 + wz3_1 + wz4_1 + (1 | ID),
chain = 1, iter = 500,
backend = "cmdstanr")
## extract population and group-level coefficients separately
fixef(m)
ranef(m)
## extract combined coefficients
coef(m)
}# }
Run the code above in your browser using DataLab