MMRM （三）

 

前边的博文，小胖简单介绍了MMRM的基本思想，而MMRM的SAS程序则主要通过PROC MIXED来实现，下边我们看一下一个最基本的MMRM的SAS程序：

 

proc mixed;

  class treat time subj;

  model y=treat time treat*time/solution ddfm=kr;

  repeated time/ subject=subj type=un;

  lsmeans treat*time / diff cl;

run;

 

解释：

 

（1）class语句：纳入分析所用到的变量，在这个例子里，我们选择了最简单的即治疗组（treat），访视时间点（time）以及受试者（subj）；

 

（2）model语句：

-            y为反应变量也就是测量值，然后我们在model中纳入了treat，time以及两者的交互项treat*time；

-            solution选项为显示出固定效应的参数估计

-            ddfm选项为指定计算分母自由度的方法，因为我们下边的repeated语句中指定TYPE=UN，这时我们应规定ddfm=kr

 

（3）repeated语句：

-            time为repeated effect

-            subject选项指定你在混合模型里的subjects，这里为subject=subj

-            type选项指定covariance structure，在MMRM分析中我们指定为Unstructured，即type=un

 

（4）lsmeans语句：

-            lsmeans  treat*time表示检验不同时间点治疗组间Least Squares Means的不同

-            diff选项表示显示出Least Squares Means的差值

-            cl选项表示计算各Least Squares Means的可信区间，缺省值为95％可信区间