重复测量设计资料的方差分析.pptVIP

例：1976年某单位报告了果胶驱铅的疗效观察。30名铅中毒的工人脱离现场后住院治疗，治疗前测得尿铅的均数为0.116mg/L，血铅均数为1.81mg/L。服用20天后再测，尿铅均数降为0.087mg/L，血铅均数降为0.73mg/L。说明果胶有较好的驱铅作用。 重复测量设计资料的方差分析 温医公卫学院黄陈平 重复测量设计资料的方差分析 重复测量资料的数据特征 重复测量数据的统计分析 第一节 重复测量资料的数据特征 前后测量设计 r①② =0.963, r①③ =-0.602, r②③ =-0.794 配对设计 r②④ =0.280, P=0.434; r③④ =-0.307, P=0.434; r②③ =0.828, P=0.003 前后测量设计与配对设计的区别 前后测量设计： 不同期观察试验结果，前后处理非随机分配。推论处理是否有效是有条件的，即假定测量时间对观察结果没有影响。 前后两次观察结果通常与差值不独立（存在相关）。 配对设计： 同对两实验单位可以随机分配处理，两个实验单位可同期观察试验结果 同对两实验单位观察结果分别与差值相互独立。 重复测量设计 当前后测量设计的重复测量次数≥3时，称重复测量设计或重复测量数据。 随机区组设计资料的方差分析 重复测量设计与随机区组设计的区别 重复测量设计区组内实验单元安排不随机、观察结果彼此不独立。 重复测量数据欲用随机区组方差分析比较处理组间差异，前提条件是重复测量误差的协方差矩阵要满足“球对称”（sphericity) 。 重复测量数据“球对称”检验Mauchlys Test of Sphericity 不满足“球对称”，可采用对F界值进行校正： “球对称”系数（0.528）× F界值的自由度（v1=3, v2=21）= F(1.58, 11.09) H0：服从球对称 H1：不服从~ α=0.10 第二节 重复测量数据统计分析 例12-3 将手术要求基本相同的15名患者随机分3组，在手术过程中分别采用A，B，C三种麻醉诱导方法，在T0 (诱导前) 、T1 、 T2 、 T3 、 T4 测量患者的收缩压，数据记录见表12-17。试进行方差分析。 SS组间(个体间)=SS干预+SS组间误差 SS组内(个体内)=SS时间+SS交互作用+SS组内误差 SS组间(个体间)=(5772+5912+…+6392)/5-C SS干预= (29922+31122+32052)/25-C SS时间=(18422+17762+…+19332)/15-C H0：服从球对称 H1：不服从~ α=0.10 接受H0，服从球对称 A，B，C三种麻醉诱导方法在不同时相的收缩压有差别 注意事项 1、各组例数不等时，需加权校正，统计软件处理已包含此项。 2、无平行对照的单组重复测量数据分析注意与随机区组方差分析的区别与联系。 （1）符合“球对称”，与随机区组方差分析等价。 （2）不符“球对称”，可选择“球对称”系数对F界值的自由度进行校正，或直接采用多变量方差分析。 THE END * * “球对称”指协方差矩阵的球形性质，即矩阵的主对角线元素（方差）相等、非主对角线元素（协方差）为零。 如果做重复测量资料分析时，发现各时间点观察值的协方差阵具有球性，那说明该资料各观察时间点的数据是独立的，所以这时候数据严格来说就算不上重复测量数据，是独立观察数据。如果不具有球性，说明各时间点观察值不独立，因此这个时候就需要对检验统计量的F的分子分母做相应调整。要不然更容易范I类错误。 *