统计分析思路统计分析思路 及案例分析及案例分析常用统计方法选择思路常用统计方法选择思路 反应变量是单变量，双变量还是多变量反应变量是单变量，双变量还是多变量 资料类型资料类型 单因素还是多因素单因素还是多因素 单样本、两样本还是多样本单样本、两样本还是多样本 是否配对、配伍设计、重复测量是否配对、配伍设计、重复测量 是否满足检验方法的前提条件是否满足检验方法的前提条件单反应变量计量资料的分析单反应变量计量资料的分析 例例1 1 某医生测量了某医生测量了3636名从事铅作业男性工人名从事铅作业男性工人 的血红蛋白含量（原始数据略），算得其均的血红蛋白含量（原始数据略），算得其均 数为数为130.83g/L130.83g/L，标准差为，标准差为25.74g/L25.74g/L。问从事。问从事 铅作业工人的血红蛋白是否不同于正常成年铅作业工人的血红蛋白是否不同于正常成年 男性平均值男性平均值140g/L140g/L？ 思路：思路：单反应变量：血红蛋白含量单反应变量：血红蛋白含量资料类型：计量资料资料类型：计量资料单因素：是否从事铅作业工人单因素：是否从事铅作业工人一份样本一份样本 该设计为样本均数与总体均数的比较该设计为样本均数与总体均数的比较 分析方法：分析方法： 如果服从正态分布，用单样本如果服从正态分布，用单样本t t检验检验 （one sample t test)one sample t test) 如果不服从正态分布，非参数检验方法如果不服从正态分布，非参数检验方法 （ Wilcoxon signed rank testWilcoxon signed rank testWilcoxon one sample test Wilcoxon one sample test）例例2 2 为比较两种方法对乳酸饮料中脂肪含量为比较两种方法对乳酸饮料中脂肪含量 测定结果是否不同，随机抽取了测定结果是否不同，随机抽取了1010份乳酸份乳酸 饮料制品，分别用脂肪酸水解法和哥特里饮料制品，分别用脂肪酸水解法和哥特里 罗紫法测定其结果如下表第罗紫法测定其结果如下表第(1)(3)(1)(3)栏。问栏。问 两法测定结果是否不同？两法测定结果是否不同？两种方法对乳酸饮料中脂肪含量的测定结果(%)编号编号 (1)(1)哥特里罗紫法哥特里罗紫法 (2)(2)脂肪酸水解法脂肪酸水解法 (3)(3)差值差值d d (4)=(2)(4)=(2) (3)(3)1 10.8400.8400.5800.5800.2600.2602 20.5910.5910.5090.5090.0820.0823 30.6740.6740.5000.5000.1740.1744 40.6320.6320.3160.3160.3160.3165 50.6870.6870.3370.3370.3500.3506 60.9780.9780.5170.5170.4610.4617 70.7500.7500.4540.4540.2960.2968 80.7300.7300.5120.5120.2180.2189 91.2001.2000.9970.9970.2030.20310100.8700.8700.5060.5060.3640.364思路思路：单反应变量：脂肪含量：单反应变量：脂肪含量资料类型：计量资料资料类型：计量资料单因素：脂肪含量的测定方法单因素：脂肪含量的测定方法 两份样本两份样本配对设计配对设计 分析方法：分析方法：差值差值如果服从正态分布，用配对如果服从正态分布，用配对t t检验检验 （paired t test)paired t test)如果不服从正态分布，非参数检验方法如果不服从正态分布，非参数检验方法 （ Wilcoxon signed rank testWilcoxon signed rank testWilcoxon one sample test Wilcoxon one sample test）例例3 3 为研究国产四类新药阿卡波糖胶囊的降为研究国产四类新药阿卡波糖胶囊的降 血糖效果，某医院用血糖效果，某医院用4040名名型糖尿病病人型糖尿病病人 进行同期随机对照试验。试验者将这些病进行同期随机对照试验。试验者将这些病 人随机等分到试验组人随机等分到试验组( (用阿卡波糖胶囊用阿卡波糖胶囊) )和对和对 照组照组( (用拜糖平胶囊用拜糖平胶囊) )，分别测得试验开始前，分别测得试验开始前 和和8 8周后的空腹血糖，算得空腹血糖下降值周后的空腹血糖，算得空腹血糖下降值 见表，能否认为该国产四类新药阿卡波糖见表，能否认为该国产四类新药阿卡波糖 胶囊与拜糖平胶囊对空腹血糖的降糖效果胶囊与拜糖平胶囊对空腹血糖的降糖效果 不同？不同？ 试验组试验组试验组试验组X X1 1 ( (n n1 1=20)=20)-0.70-0.70-5.60-5.602.002.002.802.800.700.703.503.504.004.005.805.807.107.10-0.50-0.502.502.50-1.60-1.601.701.703.003.000.400.404.504.504.604.602.502.506.006.00-1.40-1.40对对对对照照组组组组X X2 2 ( (n n2 2=20)=20)3.703.706.506.505.005.005.205.200.800.800.200.200.600.603.403.406.606.60-1.10-1.106.006.003.803.802.002.001.601.602.002.002.202.201.201.203.103.101.701.70-2.00-2.00试验组和对照组空腹血糖下降值(mmol/L) 思路：思路：单反应变量：血糖下降值单反应变量：血糖下降值资料类型：计量资料资料类型：计量资料单因素：不同药物单因素：不同药物 两份样本两份样本两独立样本设计两独立样本设计分析方法：分析方法：如果两样本都服从正态分布，且方差齐如果两样本都服从正态分布，且方差齐 性，用两样本性，用两样本t t检验检验 （two sample t test)two sample t test)；如果方差不齐，但服从正态分布，用如果方差不齐，但服从正态分布，用t t检检 验，或两样本秩和检验（验，或两样本秩和检验（Wilcoxon two Wilcoxon two sample test/Mann Whitney test)sample test/Mann Whitney test)；皆不满足，两样本秩和检验（皆不满足，两样本秩和检验（Wilcoxon Wilcoxon two sample test/Mann Whitney test)two sample test/Mann Whitney test)例例4 4 某医生为了研究一种降血脂新药的临床某医生为了研究一种降血脂新药的临床 疗效，按统一纳入标准选择疗效，按统一纳入标准选择120120名高血脂患名高血脂患 者，采用完全随机设计方法将患者等分为者，采用完全随机设计方法将患者等分为4 4 组（安慰剂组和不同剂量组），进行双盲组（安慰剂组和不同剂量组），进行双盲 试验。试验。6 6周后测得低密度脂蛋白周后测得低密度脂蛋白(mmol/L) (mmol/L) ， 问问4 4组患者的低密度脂蛋白含量有无差别组患者的低密度脂蛋白含量有无差别? ? 分析方法：分析方法：若方差齐性，且各样本均服从正态分布若方差齐性，且各样本均服从正态分布 ，用单向方差分析（，用单向方差分析（one-way-ANOVA)one-way-ANOVA)若方差不齐或某样本不服从正态分布，若方差不齐或某样本不服从正态分布， 用用Kruskal-WallisKruskal-Wallis秩和检验秩和检验分析结果分析结果 若有统计学意义，进一步做多若有统计学意义，进一步做多 重比较重比较例例5 85 8名受试对象在相同实验条件下分别接名受试对象在相同实验条件下分别接 受受4 4种不同频率声音的刺激，他们的反应率种不同频率声音的刺激，他们的反应率 （%）资料见表。问）资料见表。问4 4种频率声音刺激的反种频率声音刺激的反 应率是否有差别？应率是否有差别？8 8名受试对象对名受试对象对4 4种不同频率声音刺激的反应率（种不同频率声音刺激的反应率（%）受受试试试试 号号频频频频率率A A频频频频率率B B频频频频率率C C频频频频率率D D反反应应应应率率秩秩反反应应应应率率秩秩反反应应应应率率秩秩反反应应应应率率秩秩1 18.4 8.41 19.6 9.6 2 2 9.8 9.83 311.7 11.7 4 42 211.611.61 112.7 12.7 4 4 11.811.82 212.0 12.0 3 33 39.4 9.42 29.1 9.1 1 1 10.410.44 49.8 9.8 3 34 49.8 9.82 28.7 8.7 1 1 9.9 9.93 312.0 12.0 4 45 58.3 8.32 28.0 8.0 1 1 8.6 8.6 3.5 3.5 8.6 8.6 3.5 3.56 68.6 8.61 19.8 9.8 3 3 9.6 9.62 210.6 10.6 4 47 78.9 8.91 19.0 9.0 2 2 10.610.63 311.4 11.4 4 48 87.8 7.81 18.2 8.2 2 2 8.5 8.53 310.8 10.8 4 4思路：思路：单反应变量：反应率单反应变量：反应率资料类型：计量资料资料类型：计量资料两因素：不同频率声音，两因素：不同频率声音， 不同受试者不同受试者 四份样本四份样本配伍组设计配伍组设计 分析方法：分析方法：若方差齐性，且各样本均服从正态分布若方差齐性，且各样本均服从正态分布 ，用双向方差分析（，用双向方差分析（two-way-ANOVA)two-way-ANOVA)若不满足，用若不满足，用FriedmanFriedman秩和检验秩和检验若有统计学意义，进一步做多重比较若有统计学意义，进一步做多重比较例例5 5 某研究者为了比较甲、乙、丙、丁、戊某研究者为了比较甲、乙、丙、丁、戊 、己、己 6 6种药物给家兔注射后产生的皮肤疱疹种药物给家兔注射后产生的皮肤疱疹 大小（大小（mmmm2 2），采用拉丁方设计，选用），采用拉丁方设计，选用6 6只只 家兔、并在每只家兔的家兔、并在每只家兔的6 6个不同部位进行注个不同部位进行注 射，试验结果如下。射，试验结果如下。 分析方法分析方法拉丁方设计：一个处理因素，两个区组因拉丁方设计：一个处理因素，两个区组因 素，用素，用three-way-ANOVAthree-way-ANOVA，不能分析交互，不能分析交互 作用作用家兔家兔编编编编号号 （行区（行区组组组组）注射部位注射部位编编编编号（列区号（列区组组组组）1 12 23 34 45 56 61 1C C（8787）B B（7575）E E（8181）D D（7575）A A（8484）F F（6666）2 2B B（7373）A A（8181）D D（8787）C C（8585）F F（6464）E E（7979）3 3F F（7373）E E（7373）B B（7474）A A（7878）D D（7373）C C（7777）4 4A A（7777）F F（6868）C C（6969）B B