fit <- mmrm(
formula = FEV1 ~ RACE + SEX + ARMCD * AVISIT + us(AVISIT | USUBJID),
data = fev_data
)
if (require(emmeans)) {
emmeans(fit, ~ ARMCD | AVISIT)
pairs(emmeans(fit, ~ ARMCD | AVISIT), reverse = TRUE)
}
Run the code above in your browser using DataLab