统计分析思路

112
112 121 120 121 129 115 114 119
119
119 127 118 119 126 111 116 118
124
126 133 131 129 135 123 123 135
118
120 126 137 133 142 131 133 129
C
C Cቤተ መጻሕፍቲ ባይዱC
12
13 14 15
B法
A法
+ -
合计 172 12
184
+ 合计
8 68
76
180 80
260
例13 某抗癌新药的毒理研究中，将78只大鼠 按性别、窝别、体重、年龄等因素配成39对， 每对大鼠经随机分配分别接受甲剂量和乙剂 量注射，实验结果见下表。试分析该新药两 种不同剂量的毒性有无差异。
甲剂量 乙 剂 量 死亡(+) 生存(-) 合计
例5 8名受试对象在相同实验条件下分别接 受4种不同频率声音的刺激，他们的反应率 （%）资料见表。问4种频率声音刺激的反 应率是否有差别？
8名受试对象对4种不同频率声音刺激的反应率（%）
受试号 1 2 3 4 频率A 反应率 8.4 11.6 9.4 9.8 秩 1 1 2 2 频率B 反应率 9.6 12.7 9.1 8.7 秩 2 4 1 1 频率C 反应率 9.8 11.8 10.4 9.9 秩 3 2 4 3 频率D 反应率 11.7 12.0 9.8 12.0 秩 4 3 3 4
例8 将手术要求基本相同的15名患者随机分3 组，在手术过程中分别采用A，B，C三种 麻醉诱导方法，在T1（诱导前）、T2、T3、 T3、T4， 五个时相测量患者的收缩压，数 据记录见表。试进行方差分析。 分析方法 重复测量设计：一个处理因素有3个水 平，一个重复因素有5个水平，可分析交互 作用。（相应设计的协方差分析更合理）
n! C X !(n X )!
X n
例11 某医师研究用兰芩口服液与银黄口服液 治疗慢性咽炎疗效有无差别，将病情相似 的80名患者随机分成两组，分别用两中药 物治疗结果如下。
两种药物疗效比较
药物 兰芩口 服液 银黄口 服液 合 计 有效 无效存 合计
有效率 （%） 91.11 68.57
A，B两药联合运用的镇痛时间（min）
A药剂量
1.0 mg
2.5 mg
5.0 mg
5μg 105 80 65 75 115 80 85 120 125
B药剂量 15μg 115 105
80 125 130 90 65 120 100
30μg 75 95 85 135 120 150 180 190 160
试验组和对照组空腹血糖下降值(mmol/L)
试验组 X1 (n1=20) 对照组 X2 (n2=20)
-0.70 -5.60 2.00 2.80 0.70 3.50 4.00 5.80 7.10 -0.50
2.50
-1.60
1.70
3.00
0.40
4.50
4.60
2.50
6.00
-1.40
3.70
III
IV
99
90
104
91
203
181
48.8
50.0
合计
355
406
761
46.6
分析方法：多个样本率（构成比）的 比较， 2检验
例15 三种疗法治疗慢性喉炎的结果
例 治疗方法 综合治疗 电子治疗仪 数
两种闪烁液测定血浆中3H-cGMP的交叉试验
受试者 1 2
阶 段
Ⅰ A（760） B（860） Ⅱ B（770） A（855）
3
4 5 6 7 8
A（568）
A（780） B（960） B（940） A（635） B（440）
B（602）
B（800） A（958） A（952） B（650） A（450）
不同麻醉诱导时相患者的收缩压（mmHg）
方法
A A 1 2
序号
麻 醉 诱 导 时 相 T0 120 118 T1 108 109 T2 112 115 T3 120 126 T4 117 123
A
A A B B B B B C
3
4 5 6 7 8 9 10 11
119
121 127 121 122 128 117 118 131
B（74）
D（73）
E（76）
C（77）
D（73）
D（64）
E（75）
C（64）
D（77）
F（72）
A（82）
E（76）
F（61）
B（70）
C（82）
A（81）
B（61）
例6 下表是A、B两种闪烁液测定血浆中3HcGMP的交叉试验结果。第I阶段1、3、4、 7、9号用A测定，2、5、6、8、10号用B测 定；第II阶段1、3、4、7、9号用B测定，2、 5、6、8、10号用A测定。试对交叉试验结 果进行方差分析。 分析方法 二阶段交叉设计：处理因素，受试者， 试验阶段，不分析交互作用。
该设计为样本均数与总体均数的比较 分析方法： 如果服从正态分布，用单样本t检验 （one sample t test) 如果不服从正态分布，非参数检验方法 （ Wilcoxon signed rank test Wilcoxon one sample test）
例2 为比较两种方法对乳酸饮料中脂肪含量 测定结果是否不同，随机抽取了10份乳酸 饮料制品，分别用脂肪酸水解法和哥特里 －罗紫法测定其结果如下表第(1)~(3)栏。问 两法测定结果是否不同？
9
10
A（528）
B（800）
B（530）
A（803）
例7 观察A,B两种镇痛药物联合运用在产 妇分娩时的镇痛效果。A药取3个剂量： 1mg，2.5mg，3mg；B药也取3个剂量：5， 15，30；共9个处理组。将27名产妇随机等 分9组，记录分娩时的镇痛时间，见表。试 分析A,B两药联合运用的镇痛效果。 分析方法 析因设计：两个处理因素，各有3个水平， 可分析交互作用。
129
123 123 125
128
123 121 124
121
120 116 118
148
143 145 142
132
136 126 130
结果解释： 有交互效应 患者的收缩压在不同的诱 导方法下不同诱导时相变化的趋势不同。 注意：也可只有一个处理组
计数资料分析思路
例9 根据大量筛查得知某沿海省份成人 的α-地中海贫血基因携带率为7.6%。 某研究在本省某山区随机抽样调查成 人125例，检出α-地中海贫血基因携带 者12例，问该山区与本省一般成人的 α-地中海贫血基因携带率有无差异？
思路：单反应变量：脂肪含量 资料类型：计量资料 单因素：脂肪含量的测定方法 两份样本 配对设计 分析方法： 差值如果服从正态分布，用配对t检验 （paired t test) 如果不服从正态分布，非参数检验方法 （ Wilcoxon signed rank test Wilcoxon one sample test）
81.25
41(36.56） 24(28.44)
65
4（8.44） 11(6.56)
15
45 35
80
分析方法：两独立样本率比较，2检验， Fisher确切概率法，两样本率u检验
例12 有260份血清样品，每份样品一分为 二，用两种不同的免疫学检测方法检验 类风湿因子，结果如下。问两种检测结 果的阳性率是否相同。
分析方法：样本率与总体率的比较：单样本u 检验（条件）
例10 根据以往经验，新生儿染色体异常率 一般为1％，某医院观察了当地400名新生 儿，只有1例染色体异常，问该地新生儿染 色体异常是否低于一般。 分析方法：样本率与总体率的比较，发生 数X较小，且总体率π接近于0或1的资料， 用基于二项分布的确切概率法。 1 P( X k ) p ( X ) p( X ) CnX X (1 )n X X 0
例3 为研究国产四类新药阿卡波糖胶囊的降 血糖效果，某医院用40名Ⅱ型糖尿病病人 进行同期随机对照试验。试验者将这些病 人随机等分到试验组(用阿卡波糖胶囊)和对 照组(用拜糖平胶囊)，分别测得试验开始前 和8周后的空腹血糖，算得空腹血糖下降值 见表，能否认为该国产四类新药阿卡波糖 胶囊与拜糖平胶囊对空腹血糖的降糖效果 不同？
6.50
5.00
5.20
0.80
0.20
0.60
3.40
6.60
-1.10
6.00
3.80
2.00
1.60
2.00
2.20
1.20
3.10
1.70
-2.00
思路：单反应变量：血糖下降值 资料类型：计量资料 单因素：不同药物 两份样本 两独立样本设计
分析方法： 如果两样本都服从正态分布，且方差齐 性，用两样本t检验 （two sample t test)； 如果方差不齐，但服从正态分布，用t’检 验，或两样本秩和检验（Wilcoxon two sample test/Mann Whitney test)； 皆不满足，两样本秩和检验（Wilcoxon two sample test/Mann Whitney test)
例4 某医生为了研究一种降血脂新药的临床 疗效，按统一纳入标准选择120名高血脂患 者，采用完全随机设计方法将患者等分为4 组（安慰剂组和不同剂量组），进行双盲 试验。6周后测得低密度脂蛋白(mmol/L) ， 问4组患者的低密度脂蛋白含量有无差别?
分析方法： 若方差齐性，且各样本均服从正态分布， 用单向方差分析（one-way-ANOVA) 若方差不齐或某样本不服从正态分布， 用Kruskal-Wallis秩和检验 分析结果 若有统计学意义，进一步做多 重比较
统计分析思路 及案例分析
常用统计方法选择思路
反应变量是单变量，双变量还是多变量 资料类型 单因素还是多因素 单样本、两样本还是多样本 是否配对、配伍设计、重复测量 是否满足检验方法的前提条件
单反应变量计量资料的分析
例1 某医生测量了36名从事铅作业男性工人 的血红蛋白含量（原始数据略），算得其均 数为130.83g/L，标准差为25.74g/L。问从事 铅作业工人的血红蛋白是否不同于正常成年 男性平均值140g/L？ 思路：单反应变量：血红蛋白含量 资料类型：计量资料 单因素：是否从事铅作业工人 一份样本
家兔编号 （行区组） 1 1 2 3 4 5 6 C（87） B（73）
注射部位编号（列区组）
2 B（75） A（81） 3 E（81） D（87） 4 D（75） C（85） 5 A（84） F（64） 6 F（66） E（79）
F（73）
A（77）
E（73）
F（68）
B（74）
C（69）
A（78）
两种方法对乳酸饮料中脂肪含量的测定结果(%)
编号 (1) 1 2 3 4 5 6 7 8 9 10 哥特里－罗紫法 (2) 0.840 0.591 0.674 0.632 0.687 0.978 0.750 0.730 1.200 0.870 脂肪酸水解法 (3) 0.580 0.509 0.500 0.316 0.337 0.517 0.454 0.512 0.997 0.506 差值d (4)=(2)(3) 0.260 0.082 0.174 0.316 0.350 0.461 0.296 0.218 0.203 0.364
死亡(+)
生存(-) 合计
6（ a）
3（ c） 9
12（b）
18（d） 30
18
21 39
分析方法：配对设计率的比较： McNemar
检验
例14 黑色素瘤患者随访术后三年生存情况
治疗组 I II 三年生存 三年死亡 例数 例数 77 108 89 103 合计 185 192 三年生存 率（%） 41.6 46.4
Ri
5
6 7 8
8.3
8.6 8.9 7.8
2
1 1 1
8.0
9.8 9.0 8.2
1
3 2 2
8.6
9.6 10.6 8.5
3.5
2 3 3
8.6
10.6 11.4 10.8
3.5
4 4 4
思路： 单反应变量：反应率 资料类型：计量资料 两因素：不同频率声音， 不同受试者 四份样本 配伍组设计
分析方法： 若方差齐性，且各样本均服从正态分布， 用双向方差分析（two-way-ANOVA) 若不满足，用Friedman秩和检验 若有统计学意义，进一步做多重比较
例5 某研究者为了比较甲、乙、丙、丁、戊、 己 6种药物给家兔注射后产生的皮肤疱疹大 小（mm2），采用拉丁方设计，选用6只家 兔、并在每只家兔的6个不同部位进行注射， 试验结果如下。 分析方法 拉丁方设计：一个处理因素，两个区组因 素，用three-way-ANOVA，不能分析交互 作用
拉丁方设计与试验结果（皮肤疱疹大小，mm2）