我有几行 R code 不 太熟悉 是 hierachical modeling 算 varaince， covariance 的。fixed model。 恳请高手， 前辈转换成SAS code 用 PROC MIX。 不胜感激。

library(nlme)
hlm1 <- lme(fixed = fatigue ~ week + trt + trt.week,
            random = ~ 1 | id,
            correlation = corAR1(0, ~ week | id),
            data = Schutte,
            method = "REML")
summary(hlm1)