求关于repeat measurement的MANOVA和univariate的问题

fb_flora

看书的时候遇到的问题
看书的时候，知道GLM中的REPEATED语句默认会用两种方法（MANOVA和univariate）处理repeated measurement的数据，因此，我试着写了个univariate的程序，
求教：为什么这个GLM过程产生的univariate approach部分的结果是不一样的？
O(∩_∩)O谢谢！



数据及程序如下，吃药的数据 WD*指的是该*月的收集的因变量的值，一共是month0  month1 month2 month3 month4各收集一次，构成了repeated measurement
PAT是subject的编号

DATA WDVIS;
INPUT TREATMNT $ PAT WD0 WD1 WD2 WD3 WD4;
CENTER = INT(PAT/100);
DATALINES;
ACT 101 190 212 213 195 248
ACT 102 98 137 185 215 225
ACT 104 155 145 196 189 176
ACT 105 245 228 280 274 260
ACT 107 182 205 218 194 193
ACT 110 140 138 187 195 205
ACT 111 196 185 185 227 180
ACT 113 162 176 192 230 215
ACT 116 195 232 199 185 200
ACT 117 167 187 228 192 215
ACT 118 123 165 145 185 215
ACT 120 105 144 119 168 165
ACT 202 161 177 162 185 192
ACT 204 255 242 330 284 319
ACT 205 144 195 180 184 213
ACT 208 180 218 224 165 200
ACT 209 126 145 173 175 140
ACT 211 175 155 154 164 154
ACT 214 227 218 245 235 257
ACT 215 175 197 195 182 193
PBO 103 187 177 200 190 206
PBO 106 205 230 172 196 232
PBO 108 165 142 195 185 170
PBO 109 256 232 252 326 292
PBO 112 197 182 160 210 185
PBO 114 134 115 150 165 170
PBO 115 196 166 166 188 205
PBO 119 167 144 176 155 158
PBO 121 98 102 89 128 130
PBO 201 167 175 122 162 125
PBO 203 123 136 147 130 135
PBO 206 95 102 154 105 112
PBO 207 181 177 140 212 230
PBO 210 237 232 245 193 245
PBO 212 144 172 163 158 188
PBO 213 182 202 254 185 173
PBO 216 165 140 153 180 155
PBO 217 196 195 204 188 178
;
PROC GLM DATA = WDVIS;
CLASS TREATMNT CENTER;
MODEL WD0 WD1 WD2 WD3 WD4 = TREATMNT CENTER TREATMNT*CENTER / SS3;
REPEATED MONTH  / PRINTE SUMMARY;
QUIT;
RUN;

data wdvis1;
set wdvis;
do month=0 to 4;
        if month=0 then respose=wd0;
        else if month=1 then response=wd1;
        else if month=2 then response=wd2;
        else if month=3 then response=wd3;
        else if month=4 then response=wd4;
        output;
end;
run;
PROC GLM DATA = WDVIS1;
CLASS TREATMNT CENTER PAT MONTH;
MODEL response = TREATMNT CENTER TREATMNT*CENTER pat(TREATMNT*CENTER) month month*TREATMNT month*CENTER month*TREATMNT*CENTER/ SS3;
random pat(TREATMNT*CENTER);
test H=TREATMNT CENTER TREATMNT*CENTER E=pat(TREATMNT*CENTER) ;
QUIT;
RUN;

mono

你的第一个程序和第二个程序是不一样的东西。

你的第一段程序本身就有 manova 和univariate test.

使用manova 还是用univariate test 要根据 sphercity test 的结果来定。
Orthogonal Components 的p=0.0108 说明要用 Manova test 更加合适，其实采用univariate test 也是可以的。

Sphericity Tests
VariablesDFMauchly's CriterionChi-SquarePr > ChiSq
Transformed Variates90.163250958.754135<.0001
Orthogonal Components90.515955321.4512470.0108



你第二段程序加上了random effect 和第一段程序有了本质的区别。

fb_flora

mono 发表于 2014-8-29 21:26
你的第一个程序和第二个程序是不一样的东西。

你的第一段程序本身就有 manova 和univariate test.

谢谢，我试了一下（我把第二个的random语句去掉了），结果还是不一样。其实我是想看一下怎么自己写univariate 的模型才能和sas中直接给出的是一样的，换言之，我希望能弄明白sas这个repeated语句是计算了个什么样的univariate，跪求（ORZ）大神指引