fixef.brmcoda_pivot | R Documentation |
Estimate pivot balance coordinates
## S3 method for class 'brmcoda_pivot'
fixef(object, summary = TRUE, ...)
object |
An object of class |
summary |
Should summary statistics be returned instead of the raw values? Default is |
... |
currently ignored. |
A list of brmcoda
for each pivot balance coordinate.
if(requireNamespace("cmdstanr")){
cilr <- complr(data = mcompd, sbp = sbp,
parts = c("TST", "WAKE", "MVPA", "LPA", "SB"), idvar = "ID",
total = 1440)
# inspects ILRs before passing to brmcoda
names(cilr$between_logratio)
names(cilr$within_logratio)
names(cilr$logratio)
# model with compositional predictor at between and within-person levels
m1 <- brmcoda(complr = cilr,
formula = Stress ~ bilr1 + bilr2 + bilr3 + bilr4 +
wilr1 + wilr2 + wilr3 + wilr4 + (1 | ID),
chain = 1, iter = 500,
backend = "cmdstanr")
m_pb <- fixef.brmcoda_pivot(m1)
m_pb_raw <- fixef.brmcoda_pivot(m1, summary = FALSE)
brms::posterior_summary(m_pb_raw)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.