brmcoda_pivot: Estimate pivot balance coordinates

In florale/multilevelcoda: Estimate Bayesian Multilevel Models for Compositional Data

Estimate pivot balance coordinates

Description

Estimate pivot balance coordinates

Usage

brmcoda_pivot(object, ...)

Arguments

object	An object of class brmcoda.
...	Further arguments passed to brm.

Value

A list of brmcoda for each pivot balance coordinate.

Examples

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_1 <- brmcoda_pivot(m1)
  }

florale/multilevelcoda documentation built on Oct. 14, 2024, 4:51 a.m.