# Fit a mixed‑effects model fit <- xw_lme( outcome ~ time + treatment + (1 | patient_id), data = df, method = "REML" )