医学统计学ch王一任.pptx

《医学统计学ch王一任.pptx》由会员分享，可在线阅读，更多相关《医学统计学ch王一任.pptx（83页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 讲述内容讲述内容第一节第一节 方差分析的基本思想及其应用条件方差分析的基本思想及其应用条件第二节第二节 完全随机设计资料的方差分析完全随机设计资料的方差分析第三节第三节 随机区组设计资料的方差分析随机区组设计资料的方差分析 第四节第四节 拉丁方设计资料的方差分析拉丁方设计资料的方差分析 第五节第五节 两阶段交叉设计资料的方差分析两阶段交叉设计资料的方差分析第六节第六节 多个样本均数间的多重比较多个样本均数间的多重比较第七节第七节 多样本方差比较的多样本方差比较的BartlettBartlett检验检验 和和LeveneLevene检验检验第1页/共83页第一节第一节 方差分析的基方差分析的基

2、本思想及其应用条件本思想及其应用条件第2页/共83页目的：目的：推断多个总体均数是否有差别。推断多个总体均数是否有差别。也可用于两个也可用于两个方法：方法：方差分析，即多个样本均数比较方差分析，即多个样本均数比较的的F检验。检验。基本思想：基本思想：根据试验设计的类型，将全根据试验设计的类型，将全部测量值总的离均差平方和及其自由度部测量值总的离均差平方和及其自由度分解为两个或多个部分，除随机误差作分解为两个或多个部分，除随机误差作用外，每个部分的变异可由某个因素的用外，每个部分的变异可由某个因素的作用作用(或某几个因素的交互作用或某几个因素的交互作用)加以解加以解释释。第3页/共83页应用条件

3、：应用条件：总体总体正态且方差相等正态且方差相等 样本样本独立、随机独立、随机设计类型：设计类型：完全随机设计资料的方差分析完全随机设计资料的方差分析随机区组设计资料的方差分析随机区组设计资料的方差分析拉丁方设计资料的方差分析拉丁方设计资料的方差分析两阶段交叉设计资料的方差分析两阶段交叉设计资料的方差分析第4页/共83页完全随机设计资料的方差分析的基本思想完全随机设计资料的方差分析的基本思想 合计 N S 第i个处理组第j个观察结果第5页/共83页变异（变异（离均差平方和离均差平方和）分解）分解：离离均均差差平平方方和和(sum of squares of deviations from me

4、an，SS)：第6页/共83页1.1.总变异总变异:所有测量值之间总的变所有测量值之间总的变异程度，异程度，计算公式为计算公式为校正系数校正系数：第7页/共83页 2组组间间变变异异：各各组组均均数数与与总总均均数数的的离离均均差差平平方方和和，表表示示处处理理因因素素的的作作用用和和随随机机误误差差的的影影响响，计算公式为：计算公式为：如果如果有有g个总体均数个总体均数有有差别差别 如果如果无无g个总体均数个总体均数无无差别差别第8页/共83页 3 3组内变异：组内变异：各组内各测量值各组内各测量值X Xijij与其所与其所在组的均数的差值的平方和，表示随机误差在组的均数的差值的平方和，表示

5、随机误差的影响。计算公式为的影响。计算公式为:第9页/共83页三种变异的关系三种变异的关系：第10页/共83页 均方差，均方均方差，均方(mean square，MS)。第11页/共83页检验统计量：检验统计量：如果如果 ，则，则 都为随都为随机误差机误差 的估计，的估计，F F值应接近于值应接近于1 1。如果如果 不全相等，不全相等，F F值将明显大于值将明显大于1 1。用用F F界值（单侧界值）确定界值（单侧界值）确定P P值。值。第12页/共83页第二节第二节完全随机设计资料的方差分析完全随机设计资料的方差分析第13页/共83页 将将全全部部试试验验对对象象分分配配到到g个个处处理理组组

6、（水水平平组组），各各组组分分别别接接受受不不同同的的处处理理，试试验验结结束束后后比比较较各各组组均均数数之之间间的的差差别别有有无无统统计计学学意意义义，推推论处理因素的效应。论处理因素的效应。一、完全随机设计一、完全随机设计第14页/共83页 例例4-14-1 某某医医生生为为了了研研究究一一种种降降血血脂脂新新药药的的临临床床疗疗效效，按按统统一一纳纳入入标标准准选选择择120名名患患者者，采采用用完完全全随随机机设设计计方方法法将将患患者者等等分分为为4组组进进行行双盲试验。问如何进行分组？双盲试验。问如何进行分组？第15页/共83页（1 1）完全随机分组方法：）完全随机分组方法：1

7、.编编号号：120名名高高血血脂脂患患者者从从1开开始始到到120，见表，见表4-2第第1行；行；2.取取随随机机数数字字：从从附附表表15中中的的任任一一行行任任一一列列开开始始，如如第第5行行第第7列列开开始始，依依次次读读取取三三位位数数作作为为一一个个随随机机数数录录于于编编号下，见表号下，见表4-2第第2行；行；第16页/共83页3.3.编编序序号号：按按数数字字从从小小到到大大 (数数据据相相同同则则按按先先后后顺顺序序）编序号，见表编序号，见表4-24-2第第3 3行。行。4.4.事事先先规规定定：序序号号1-301-30为为甲甲组组，序序号号31-6031-60为为乙乙组组，序

8、序号号61-9061-90为丙组，序号为丙组，序号91-12091-120为丁组，见表为丁组，见表4-24-2第四行。第四行。第17页/共83页（2 2）统计分析方法选择：）统计分析方法选择：1.对于正态分布且方差齐同的资料，常采用单因素方差分析的单向分类的方差分析或成组资料的t 检验（g=2）；2.对于非正态分布或方差不齐的资料，可进行数据变换或采用Wilcoxon秩和检验。第18页/共83页二、变异分解二、变异分解 第19页/共83页 例例4-2某某医医生生为为了了研研究究一一种种降降血血脂脂新新药药的的临临床床疗疗效效，按按统统一一纳纳入入标标准准选选择择120名名高高血血脂脂患患者者，

9、采采用用完完全全随随机机设设计计方方法法将将患患者者等等分分为为4组组（具具体体分分组组方方法法见见例例4-1），进进行行双双盲盲试试验验。6周周后后测测得得低低密密度度脂脂蛋蛋白白作作为为试试验验结结果果，见见表表4-3。问。问4个处理组患者的低密度脂蛋白含量总体均数有无差别个处理组患者的低密度脂蛋白含量总体均数有无差别?第20页/共83页表表4-3 44-3 4个处理组低密度脂蛋白测量值个处理组低密度脂蛋白测量值(mmol/L)(mmol/L)第21页/共83页三、分析步骤H0：即4个试验组总体均数相等 H1：4个试验组总体均数不全相等 2.计算检验统计量:1.建立检验假设，确定检验水准:

10、第22页/共83页第23页/共83页表4-5 完全随机设计方差分析表列方差分析表第24页/共83页3.确定P值，作出推断结论：按 水准，拒绝H0，接受H1，认为4个试验组ldl-c总体均数不相等，即不同剂量药物对血脂中ldl-c降低影响有差别。第25页/共83页注意：注意：方方差差分分析析的的结结果果拒拒绝绝H0，接接受受H1，不不能能说说明明各各组组总总体体均均数数间间两两两两都都有有差差别别。如如果果要要分分析析哪哪些些两两组组间间有有差差别别，可可进进行行多多个个均均数数间间的的多多重重比比较较（见见本本章章第第六六节节）。当当g=2时时，完完全全随随机机设设计计方方差差分分析析与与成成

11、组组设设计计资资料料的的t检检验验等等价价，有有。第26页/共83页第三节第三节随机区组设计资料的方差分析随机区组设计资料的方差分析第27页/共83页一、随机区组设计配伍组设计配伍组设计(randomizedblockdesign)先先按按影影响响试试验验结结果果的的非非处处理理因因素素（如如性性别别、体体重重、年年龄龄、职职业业、病病情情、病病程程等等）将将受受试试对对象象配配成成区区组组，再再分分别别将将各各区区组组内内的的受受试试对对象象随随机机分分配配到到各各处处理理或对照组。或对照组。（1 1）随机分组方法：第28页/共83页（2 2）随机区组设计的特点 随随机机分分配配的的次次数数

12、要要重重复复多多次次，每每次次随随机机分分配配都都对对同同一一个个区区组组内内的的受受试试对象进行，且各个处理组受试对象数量相同，对象进行，且各个处理组受试对象数量相同，区组内均衡区组内均衡。在在进进行行统统计计分分析析时时，将将区区组组变变异异离离均均差差平平方方和和从从完完全全随随机机设设计计的的组组内内离离均均差差平平和和中分离出来，从而中分离出来，从而减小组内离均差平方和减小组内离均差平方和（误差平方和），提高了统计检验效率。（误差平方和），提高了统计检验效率。第29页/共83页例4-3如何按随机区组设计，分配如何按随机区组设计，分配5 5个个区组的区组的1515只小白鼠接受甲、乙、丙

13、三种抗癌药只小白鼠接受甲、乙、丙三种抗癌药物？物？分分组组方方法法：先先将将小小白白鼠鼠按按体体重重编编号号，体体重重相相近近的的3只只小小白白鼠鼠配配成成一一个个区区组组，见见表表4-6。在在随随机机数数字字表表中中任任选选一一行行一一列列开开始始的的2位位数数作作为为1个个随随机机数数，如如从从第第8行行第第3列列开开始始纪纪录录，见见表表4-6；在在每每个个区区组组内内将将随随机机数数按按大大小小排排序序；各各区区组组中中内内序序号号为为1的的接接受受甲甲药药、序序号号为为2的的接接受受乙乙药药、序序号号为为3的接受的接受丙药丙药，分配结果见表，分配结果见表4-6。第30页/共83页第3

14、1页/共83页（3 3）统计方法选择：1.正正态态分分布布且且方方差差齐齐同同的的资资料料，应应采采用用单单因因素素方方差差分分析析的的双双向向分分类类方方差差分分析析或或配配对对t检验（检验（g=2）；）；2.当当不不满满足足方方差差分分析析和和t检检验验条条件件时时，可可对对数数据据变变换换或或采采用用随随机机区区组组设设计计资资料料的的FriedmanM检验检验。第32页/共83页 表4-7 随机区组设计的试验结果 第33页/共83页二、变异分解二、变异分解(1)总变异：总变异：反映所有观察值之间的变异反映所有观察值之间的变异,记为记为SS总总。(2)处处理理间间变变异异：由由处处理理因

15、因素素的的不不同同水水平平作作用用和和随随机机误差产生的变异，记为误差产生的变异，记为SS处理处理。(3)区区组组间间变变异异：由由不不同同区区组组作作用用和和随随机机误误差差产产生生的的变异，记为变异，记为SS区组区组.(4)误误差差变变异异：完完全全由由随随机机误误差差产产生生的的变变异异，记记为为SS误误差差。对总离均差平方和及其自由度的分解，有对总离均差平方和及其自由度的分解，有:第34页/共83页 表4-8 随机区组设计资料的方差分析表 第35页/共83页三、分析步骤 例例4-4某某研研究究者者采采用用随随机机区区组组设设计计进进行行实实验验，比比较较三三种种抗抗癌癌药药物物对对小小

16、白白鼠鼠肉肉瘤瘤抑抑瘤瘤效效果果，先先将将15只只染染有有肉肉瘤瘤小小白白鼠鼠按按体体重重大大小小配配成成5个个区区组组，每每个个区区组组内内3只只小小白白鼠鼠随随机机接接受受三三种种抗抗癌癌药药物物（具具体体分分配配方方法法见见例例4-3），以以肉肉瘤瘤的的重重量量为指标，试验结果见表为指标，试验结果见表4-9。问三种不同的药物的抑瘤效果有无差别？。问三种不同的药物的抑瘤效果有无差别？第36页/共83页 表表4-9 不同药物作用后小白鼠肉瘤重量不同药物作用后小白鼠肉瘤重量（g g）第37页/共83页H0：，即即三三种种不不同同药药物物作用后作用后小白鼠肉瘤重量的小白鼠肉瘤重量的总体均数相等总

17、体均数相等 H1：三种不同药物作用后小白鼠肉瘤重：三种不同药物作用后小白鼠肉瘤重量的量的总体均数不全相等总体均数不全相等第38页/共83页第39页/共83页第40页/共83页 据据 1=2、2=8查查附附表表3的的F界界值值表表，得得在在=0.05的的水水准准上上，拒拒绝绝H0，接接受受H1，认认为为三三种种不不同同药药物物作作用用后后小小白白鼠鼠肉肉瘤瘤重重量量的的总总体体均均数数不不全全相相等等，即即不不同同药药物物的的抑抑瘤瘤效效果果有差别。同理可对区组间的差别进行检验。有差别。同理可对区组间的差别进行检验。第41页/共83页注意：注意：方方差差分分析析的的结结果果拒拒绝绝H0，接接受受

18、H1，不不能能说说明明各各组组总总体体均均数数间间两两两两都都有有差差别别。如如果果要要分分析析哪哪些些两两组组间间有有差差别别，可可进进行行多多个个均均数数间间的的多多重重比比较较（见见本本章章第第六六节节）。当当g=2时时，随随机机区区组组设设计计方方差差分分析析与与配配对对设设计计资资料料的的t检检验验等等价价，有有。第42页/共83页 随随机机区区组组设设计计确确定定区区组组应应是是对对试试验验结结果果有有影影响响的的非非处处理理因因素素。区区组组内内各各试试验验对对象象应应均均衡衡，区区组组之之间间试试验验对对象象具具有有较较大大的的差差异异为为好好，这这样样利利用用区区组组控控制制

19、非非处处理理因因素素的的影影响响，并并在在方方差差分分析析时时将将区区组组间的变异从组内变异中分解出来。间的变异从组内变异中分解出来。因因此此，当当区区组组间间差差别别有有统统计计学学意意义义时时，这这种种设设计计的的误误差差比比完完全全随随机机设设计计小小，试试验效率得以提高。验效率得以提高。第43页/共83页第四节拉丁方设计资料的方差分析（不讲）第44页/共83页第五节两阶段交叉设计资料的方差分析（不讲）第45页/共83页第六节第六节多个样本均数间的多重比较多个样本均数间的多重比较 （multiple comparisonmultiple comparison）第46页/共83页多重比较不

20、能用两样本均数比多重比较不能用两样本均数比较的较的 t t 检验检验！若用两样本均数比较的若用两样本均数比较的t 检验检验进行多重比较，将会加大犯进行多重比较，将会加大犯类错误类错误（把本无差别的两个总体均数判为有（把本无差别的两个总体均数判为有差别）的概率。差别）的概率。第47页/共83页 例例 如如，有有 4个个 样样 本本 均均 数数，两两 两两 组组 合合 数数 为为 ，若若用用 t 检检验验做做6次次比比较较，且且每每次次比比较较的的检检验验水水准准定定为为=0.05，则则每每次次比比较较不不犯犯类类错错误误的的概概率率为为（10.05），6次次均均不不犯犯类类错错误误的的概概率率为

21、为 ，这这时时，总总的的检检验验水水准准变变为为 ，远远比比0.05大大。因因此此，样样本本均均数数间间的的多多重重比比较不能用两样本均数比较的较不能用两样本均数比较的 t 检验。检验。第48页/共83页适用条件适用条件：当当方方差差分分析析的的结结果果为为拒拒绝绝H0，接接受受H1时时，只只说说明明g个个总总体体均均数数不不全全相相等等。若若想想进进一一步步了了解解哪哪些些两两个个总总体体均均数数不不等等，需需进进行行多多个个样样本本均均数间的两两比较或称多重比较。数间的两两比较或称多重比较。第49页/共83页一、LSD-t检验（leastsignificantdifference）适用范围

22、：一对或几对在专业上有特殊意义的样本均数间的比较。第50页/共83页检验统计量检验统计量t的计算公式为的计算公式为式中 第51页/共83页注意：第52页/共83页第53页/共83页 例例4-7 对例对例4-2资料，问高血脂资料，问高血脂患者的降血脂新药患者的降血脂新药2.4g组、组、4.8g组、组、7.2g组与安慰剂组的低密度组与安慰剂组的低密度脂蛋白含量总体均数有无差别？脂蛋白含量总体均数有无差别？第54页/共83页，即即降降血血脂脂新新药药2.4g组组与安慰剂与安慰剂组的低密度脂蛋白含量总体均数相等组的低密度脂蛋白含量总体均数相等，即即降降血血脂脂新新药药2.4g组组与安慰剂与安慰剂组的低

23、密度脂蛋白含量总体均数不等组的低密度脂蛋白含量总体均数不等=0.05降血脂新药降血脂新药2.4g2.4g组与安慰剂组的比较：组与安慰剂组的比较：第55页/共83页第56页/共83页新药新药4.8g组组VS安慰剂组安慰剂组:LSD-t为为-4.297.2g组组VS安慰剂组安慰剂组:LSD-t为为-8.59。同理：按水准，降血脂新药4.8g组、7.2g组与安慰剂组间差别有统计学意义。第57页/共83页二、Dunnett-t 检验 适用条件：g-1个实验组与一个对照组均数差别的多重比较，检验统计量Dunnett-t有专门的界值表(附表5)。第58页/共83页式中 计算公式为：Dunnett-第59页

24、/共83页例例4-8对例对例4-2资料，问高血脂资料，问高血脂患者的三个不同剂量降血脂新药组与安慰患者的三个不同剂量降血脂新药组与安慰剂组的低密度脂蛋白含量总体均数是否有剂组的低密度脂蛋白含量总体均数是否有差别？差别？H0：i=0，即各实验组与安慰剂组的低密度即各实验组与安慰剂组的低密度脂蛋白含脂蛋白含量总体均数相等量总体均数相等H1：i 0，即各实验组与安慰剂组的低密度即各实验组与安慰剂组的低密度脂蛋白含量总体均数不等脂蛋白含量总体均数不等=0.05第60页/共83页Dunnett-Dunnett-Dunnett-第61页/共83页第62页/共83页三、SNK-q检验（Student-New

25、man-Keuls）适用于多个样本均数两两之间的全面比较。第63页/共83页检验统计量q的计算公式为第64页/共83页例4-9对例4-4资料，问三种不同药物的抑瘤效果两两之间是否有差别？H0：A=B，即任两对比较组的总体均数相等H1：AB，即任两对比较组的总体均数不相等=0.05第65页/共83页将三个样本均数由小到大排列，并编组次：第66页/共83页 列列出出对对比比组组，并并计计算算两两对对比比组组的的均均数数之差，写出两对比组包含的组数之差，写出两对比组包含的组数a a。已已知知=8=8和和a a，查查附附表表4 4的的q q界界值值，得得出出相应的相应的q q界值。界值。以以实实际际的

26、的q q值值和和相相应应的的q q界界值值作作比比较较，确定对应的确定对应的P P值值 。第67页/共83页表4-15 多个均数两两比较值 第68页/共83页结论：结论：可认为可认为A A药和药和B B药、药、C C药的抑瘤药的抑瘤 效果有差别，还不能认为效果有差别，还不能认为B B药和药和C C药药的的抑瘤效果有差别。抑瘤效果有差别。第69页/共83页第七节第七节多样本方差比较的Bartlett检验和Levene检验第70页/共83页 在在进进行行方方差差分分析析时时要要求求所所对对比比的的各各组组即即各各样样本本的的总总体体方方差差必必须须是是相相等等的的，这这一一般般需需要要在在作作方方

27、差差分分析析之之前前，先先对对资资料料的的的的方方差差齐齐性性进进行行检检验验，特特别别是是在在样样本本方差相差悬殊时，应注意这个问题。方差相差悬殊时，应注意这个问题。对对两两样样本本方方差差进进行行齐齐性性检检验验的的方方法法前前已已介介绍绍。本本节节介介绍绍多多样样本本（也也适适用用于于两两样样本本）方方差差齐齐性性检检验验的的BartlettBartlett检检验验法法和和LeveneLevene检验法检验法。第71页/共83页一、一、Bartlett检验检验第72页/共83页第73页/共83页第74页/共83页第75页/共83页表4-17 例4-2的方差齐性检验表 第76页/共83页 第77页/共83页二、二、Levene检验检验资料要求：可不具有正态性。资料要求：可不具有正态性。检验统计量：检验统计量：W计算公式：计算公式：第78页/共83页式中式中第79页/共83页检验步骤：检验步骤：第80页/共83页第四章练习题一、最佳选择题一、最佳选择题全做全做二、计算分析题二、计算分析题 第第1 1、3 3、7 7题题第81页/共83页第82页/共83页谢谢您的观看！第83页/共83页