如何将SAS混合模型代码转换为R？

Question

我试图将在SAS中使用的类似混合模型应用于R中的其他分析。我需要了解如何用R编写初始代码。 这将是一个具有重复测量和随机因素的混合模型。

这是我在SAS中使用的：

proc sort data=comp1; by age order subj light clock;
  PROC MIXED DATA = comp1;by age;
  CLASS  subj order light clock;

  MODEL mel = light clock light*clock/DDFM = KR;
  RANDOM subj order/ subject=subj ;
  Repeated clock /TYPE = ar(1)  ;


  LSMEANS light * clock/PDIFF CL;
ods output  diffs=diffcm1 ;
run;

您知道如何用R编写吗？ 谢谢