The relsa_reference functions calculates the default RELSA reference model - to keeps things simple. The reference set is the surgery data (transmitter-implanted, female B6 mice)

relsa_reference(
  vars = c("bwc", "burON", "hr", "hrv", "temp", "act"),
  normthese = c("hr", "burON", "hrv", "temp", "act"),
  turnthese = c("hr", "temp")
)

Value

list with baselines (bsl) and cluster levels (levels)