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") )
list with baselines (bsl) and cluster levels (levels)