重复测量方差分析 精选文档.ppt

时间点间 误差 4. 计算 F 值 球对称性通常采用 Mauchlys test 检验来判断 其结果按 α ＝ 0.1 水准检验，不满足球对称性， 对系数 ε 进行校正，其结果如下： 结果显示：治疗前与治疗后不同时间转氨酶平均水平不同。 5. 单组重复测量方 差分析数据结构 6. 基本程序格式 Nouni 不打印输出单变量分析结果 Printe 产生 Mauchly 球性检验的统计量 χ 2 、 P 值 分析： 单组重复测量数据分析的缺陷只能分析 观察对象的观测值在不同时间点的差别。 专业认为： 不同时间点上的观测值变化可能是“处 理”的作用，也可能是患者病情的自然变化，与“处 理”无关，如果要分析“处理”效应，必须设立一个 平行对照组，通过组间差别的大小说明“处理”组效 应的大小。两独立样本重复测量设计是将 N 个受试对象 随机地等分两组，一组作为实验组，另一组作为对照 组。 为研究国产某药品与同类型进口药品对慢 性乙型肝炎患者谷丙转氨酶（ ALT ）水平的影 响，将 20 名慢性乙型肝炎患者随机分为两组， 一组服用国产药作为实验组，一组服用进口药 作为对照组。对每一患者在治疗前、治疗后 12 周、 24 周、 36 周重复 4 次测量 ALT 水平，实验结 果见下表。问两种药物对慢性乙型肝炎患者的 ALT 水平影响是否不同？ 5. 举例：两组重复测量数据方差分析 表 国产药 进口药 分析结果 同一患者重复测量值就是组间差异 患者之间的个体差异就是组内误差 两样本重复测量数据方差分析结果 程序如下 两组重复测 量数据结构 相关分析结果 球对称性检验结果 方差分析结果 结果显示汇总 总结论为： 不同药物对 ALT 水平未见不同，但两药物在不同时间 点 ALT 有差别，药物与时间的交互作用无意义。 问题？ 从以上分析可见：无论单组或两 组分析，我们仅算出总的有差 别，两两比较未知的？ 分析方法与一般方差分析不同。 Contrast(1) 产生不同水 平间同括弧内指定的参考 水平对比检验 Summary 定义每个对 比组生成方差分析表 后多个重复测量时间点与前第 1 时间点比 原程序如下： 结果分析 两组比较 数据分析 所关注的结果如下 Contrast(1) 结果比较结构 有时试验仅对先后两个时间点的比较感兴趣， 也希望比较水平 1 与 2 、 2 与 3 、 3 与 4 比 可用选择项 profile ，程序如下： proc glm data=SASfeng.p178; model time0 time12 time24 time36 =A /nouni; repeated time profile /printm summary ; run; 提示： 显示不同水平比较时的结构 方 差 分 析 单因素 ：完全随机设计 两因素：随机区组设计 多因素：？？？ 重复测量设计 单组 多组 单因素方差设计只涉及一个处理因素 该因素 至少 有两个水平 有两水平时：称为两样本均数比较 两水平以上：称多个样本均数比较的方差分析 知识点：多重比较时有 特定的方法，不能用两 样本均数比较，此时容 易加大Ⅰ类错误（把本 无差别的两个总体判为 有差别）的概率。 为什么？ Why? 举例： 有 4 个样本均数 ? ? 6 4 2 ? 如果用 t 检验每次比较选 α ＝ 0.05 ， 不犯Ⅰ错误的概率 1 － α 6 次不犯Ⅰ错误（ 1 － α ） 6 总的水准： 1 － ( 1 － α) 6 ＝ 1 － ( 1 － 0.05) 6 =0.26 比 0.05 大多了！！ 比较的次数越多犯Ⅰ错误的概率越大！！ 把无差别的结果判为有差别 完全随机设计 如比较 4 种饲料对小鼠体重增加量的影响，处理 因素是饲料，有 4 个水平（不同饲料）。 完全随机设计是将 n 个小鼠随机分为 4 组。 随机区组设计 是将 n 个小鼠按 出生体重相近 的原则， 4 个一组 相配（称为区组）后，再随机分不同的水平组， 称为两因素方差分析。 非处理因素 应用条件 1. 各样本是相互独立的随机样本，均服 从正态分布。 2. 相互比较的各样本的总体方差相等， 具有方差齐性。 重复测量设计 一、重复测量资料的数据特征 当对同一受试对象在不同时间重复测量次数 p≥3 时，称为重复测量设计或重复测量数据。