4  MMRM模型sas和R结果对比

ep1_adeff_mmrm_womac <- read.xlsx("Data/2023053_ep1_adeff_mmrm_womac.xlsx")%>%
    mutate(AVISITN = factor(AVISITN),
         KLRGRN = factor(KLRGRN),
         TRT01P = fct_relevel(TRT01P, "安慰剂组"))

fit_womac_mmrm <- mmrm(
  formula = CHG ~ TRT01P*AVISITN + BASE+ KLRGRN + us(AVISITN | USUBJID),
  data = ep1_adeff_mmrm_womac,
  control =  mmrm_control(method = "Kenward-Roger" , vcov ="Kenward-Roger") 
)
emmeans(fit_womac_mmrm,~TRT01P*AVISITN, at = list(AVISITN="11"))

   TRT01P AVISITN emmean   SE  df lower.CL upper.CL
 安慰剂组 11       -40.1 3.57 202    -47.1    -33.0
 试验组   11       -37.1 3.58 207    -44.2    -30.1

Results are averaged over the levels of: KLRGRN 
Confidence level used: 0.95 

ep1_adeff_mmrm_womac_pair_lsmean <- pairs(emmeans(fit_womac_mmrm, ~ TRT01P*AVISITN, at = list(AVISITN="11")), reverse = T)
ep1_adeff_mmrm_womac_pair_lsmean%>%as_tibble()

# A tibble: 1 × 6
  contrast                              estimate    SE    df t.ratio p.value
  <chr>                                    <dbl> <dbl> <dbl>   <dbl>   <dbl>
1 试验组 AVISITN11 - 安慰剂组 AVISITN11     2.96  4.87  180.   0.607   0.545

include_graphics("Data/image/Screenshot 2025-12-31 145316.png")