重复测量的方差20160428研究.ppt

printe 该选项给出两个球性检验结果，一个是针对repeated语句中指定的转换形式的变量（如contrast），另一个是针对正交转换的变量。 针对正交转换的变量的检验结果是验证是否满足球性假定的关键，因为球性假定是正交转换后的球性检验。如果repeated语句指定polynomial，由于polynomial为正交转换，这时两个球性检验的结果是一致的。当repeated语句指定其它转换形式时，两个球性检验结果不一致，以针对正交转换的变量的结果作为判断是否满足球性的标准 * * 若资料不满足球对称性，应对自由度进行效正，用Greenhouse-Geisser (G-G)法或Huynh-Feldt ( H-F )法的“球对称”系数Epsilon乘以受试对象内各变异的自由度，得效正“自由度”（调小），再查F界值表获得P值。 proc glm 选项; class 分组变量; model 分析变量=分组变量; lsmeans 分组变量; means 分组变量; repeated 重复测量变量名 测量水平 转换项 /选项。 三、实例及SAS实现 repeated 重复测量变量名 测量水平 转换项 /选项; 【repeated】后必须指定重复测量变量名，测量水平可选，转换项也是可选的。如： repeated time 4; 表示指定重复测量变量名为time，测量水平（即重复次数）为4次，中间均以空格隔开。其中time由读者自行命名，也可为其它名称，如repeated t 4，repeated measure 4等均可。 repeated语句后的转换项主要用于各次测量之间的比较，常用的有两种转换方式：contrast和polynomial。 contrast表示指定一个参照水平，其它均与参照水平比较。 如repeated time 4 contrast(1); 表示共有4次重复测量，以第1次测量为参照，第2、3、4次测量分别与第1次测量水平进行比较。 polynomial表示正交多项式对比，即正交转换后判断曲线的n阶水平是否有统计学意义。如：repeated time 4 polynomial; 表示4次重复测量，分别对一次、二次和三次曲线分析是否有统计学意义，也就是判断测量值随时间变化大致呈何种趋势。 printe 该选项可输出一些相关矩阵（不是我们所关心的），更重要的是给出球性检验结果。 summary 输出重复测量各水平间的比较结果。如果没有该选项，即使前面指定了contrast或polynomial，也不会显示比较结果 例 1：试分析ALT变化趋势。 data alt; input id t1 t2 t3 t4; cards; 1 160 105 147 135 2 415 371 258 182 3 327 94 36 51 4 174 113 63 50 5 201 26 55 20 6 289 20 17 21 7 85 44 56 62 8 176 165 136 83 9 76 215 34 81 10 75 94 51 59 ; run; proc glm; model t1 t2 t3 t4=/nouni; /* model语句指定4个因变量，即t1-t4， nouni表示不输出单变量分析结果 */ repeated time 4 contrast(1) /printe summary; /*repeated语句指定有4个时间点; contrast(1)表示后3个时间点分别与第1时 间点比较; printe输出球性检验; summary输出后3个时间点分别与第1时间 点比较的结果*/ run; 第一部分是重复测量的基本信息及部分矩阵信息。 给出了误差SSCP矩阵的偏相关系数阵以及误差SSCP（交叉乘积矩阵平方和）矩阵等。 第二部分是