第七章两独立样本定量资料的统计分析的Stata实现Word文档格式.docx
《第七章两独立样本定量资料的统计分析的Stata实现Word文档格式.docx》由会员分享，可在线阅读，更多相关《第七章两独立样本定量资料的统计分析的Stata实现Word文档格式.docx（11页珍藏版）》请在冰豆网上搜索。
7.86
例7-2为研究接触某重金属对人体血胰岛素水平有无影响，研究者从接触某重金属的职业工人中随机抽取14人，从非接触工人中随机抽取14人，测量每个工人的血胰岛素水平（ulu/ml）（表7-2），试分析上述两个人群的血胰岛素平均水平有无差异。
表7-2某重金属接触工人和非接触工人血胰岛素水平（ulu/ml）
接触重金属工人
7.84
8.11
8.45
9.38
10.22
12.28
11.57
11.38
8.23
8.92
3.02
4.54
7.11
6.04
非接触工人
17.94
7.42
9.12
9.45
16.48
24.12
8.57
6.74
9.23
16.29
5.73
例7-3随机抽取30名在校大学生，男性16人，女性14人，检测血红蛋白含量（g/ul），问男性和女性大学生的血红蛋白平均含量有无差别？
（数据略）
例7-4例7-1中要检验两种饮水小鼠的肝脏锌平均含量有无差别，需要先检验实验组和对照组所来自的总体方差是否相等，即对例7-1中资料进行方差齐性检验。
例7-1的Stata数据格式如下：
group
x
1
2
3
7.1
4
5
6
7
8
9.3
9
10
11
12
13
7.2
14
15
16
17
18
19
20
正态性检验的程序和结果如下：
swilkxifgroup=="
实验组"
对实验组资料进行正态性检验
Shapiro-WilkWtestfornormaldata
Variable|ObsWVzProb>
z
-------------+-------------------------------------------------
x|100.979840.311-1.8030.96432
P=0.96432>
，按正态性假定进行统计分析。
对照组"
对对照组资料进行正态性检验
Shapiro-WilkWtestfornormaldata
x|100.935021.0010.0020.49903
P=0.49903>
方差齐性检验程序和结果如下：
sdtestx,by（group）方差齐性检验
Varianceratiotest
------------------------------------------------------------------------------
Group|ObsMeanStd.Err.Std.Dev.[95%Conf.Interval]
---------+--------------------------------------------------------------------
对照组|106.968.2352625.74396536.4357997.500201
实验组|108.196.38405211.2144797.3272149.064786
combined|207.582.26054561.1651957.0366728.127328
Ho:
sd（对照组）=sd（实验组）
F（9,9）observed=F_obs=0.375
F（9,9）lowertail=F_L=F_obs=0.375
F（9,9）uppertail=F_U=1/F_obs=2.665
Ha:
sd
（1）<
sd
（2）Ha:
sd
（1）!
=sd
（2）Ha:
sd
（1）>
sd
（2）
P<
F_obs=0.0802P<
F_L+P>
F_U=0.1604P>
F_obs=0.9198
方差齐性检验的P=0.1604，可以按方差齐性的假定作进一步统计分析。
例7-5对例7-1资料，试比较两种饮水小鼠肝脏平均锌含量有无差别？
相应的Stata程序为
ttestx,by（group）两独立样本t检验
输出结果为:
Two-samplettestwithequalvariances
diff|-1.228.4503826-2.174219-.2817813
Degreesoffreedom:
18
mean（对照组）-mean（实验组）=diff=0
diff<
0Ha:
diff!
=0Ha:
diff>
0
t=-2.7266t=-2.7266t=-2.7266
t=0.0069P>
|t|=0.0138P>
t=0.9931
t检验的P=0.0138<
0.05，可以认为饮用污染的水使的小鼠肝脏中锌平均含量高于饮用未污染水的小鼠。
例7-6对例7-2中的资料进行方差齐性检验，得F＝4.13，P<
0.02（具体计算略），可认为两总体方差不等。
试比较接触重金属工人和非接触工人的血胰岛素平均水平有无差别？
例7-2的Stata数据为：
21
22
23
24
25
26
27
28
由于方差不齐，所以采用方差不齐的t检验。
程序和结果为：
ttestx,by（group）unequal两独立样本t检验，方差不齐
Two-samplettestwithunequalvariances
非接触工|1411.171.4296815.3493758.08136314.25864
接触重金|148.363571.70326992.6313956.8442499.882894
combined|289.766786.82708214.3765078.06975311.46382
diff|2.8064291.593291-.5290436.1419
Satterthwaite'
sdegreesoffreedom:
18.9433
mean（非接触工）-mean（接触重金）=diff=0
t=1.7614t=1.7614t=1.7614
t=0.9528P>
|t|=0.0943P>
t=0.0472
t检验的P=0.0943>
0.05，可以认为接触重金属工人和非接触工人的血胰岛素平均水平没有差别。
例7-7为研究孕妇患有妊娠合并症对葡萄糖耐受水平有无影响，对17名孕妇作葡萄糖耐受水平试验，其中有9名孕妇曾患有妊娠合并症，而另外8名没有患妊娠合并症（表7-4中第
（1）、（3）列）。
试比较患与未患妊娠合并症的孕妇葡萄糖平均耐受水平是否相同？
表7-4患与未患妊娠合并症的孕妇葡萄糖耐受水平试验比较
患妊娠合并症
未患妊娠合并症
葡萄糖耐受水平
（1）
秩次
（2）
（3）
（4）
110
120
4.5
119
140
133
162
127
184
141
132
117
128
135
177
143
181
n1=8
T1=47.5
n2=9
T2=105.5
例7-7的Stata数据如下：
Stata程序和结果如下：
ranksumx,by（g）
Two-sampleWilcoxonrank-sum（Mann-Whitney）test
g|obsranksumexpected
-------------+---------------------------------
1|847.572
2|9105.581
combined|17153153
unadjustedvariance108.00
adjustmentforties-0.13
----------
adjustedvariance107.87
Ho:
x（g==1）=x（g==2）
z=-2.359
Prob>
|z|=0.0183
t检验的P=0.0183<
0.05，可以认为患与未患妊娠合并症的孕妇葡萄糖平均耐受水平不同。
未患妊娠合并症的孕妇葡萄糖平均耐受水平比患的高。
例7-8例7-2资料经方差齐性检验认为两总体方差不等，为比较重金属接触工人和非接触工人的血胰岛素水平有无差别，也可以用Wilcoxon秩和检验。
例7-9为检测两种饮料被大肠杆菌污染的情况，分别从每种饮料中各随机抽取10份样品，每份样品取出1ml作细菌培养。
甲饮料中共得到大肠杆菌菌落450个，乙饮料中共得到大肠杆菌菌落230个，试问两种饮料的平均菌落数有无差别？
Stata数据如下：
w
450
230
分析程序和结果如下：
poissonxgroup
PoissonregressionNumberofobs=2
LRchi2
（1）=72.47
chi2=0.0000
Loglikelihood=-7.612088PseudoR2=0.8264
x|Coef.Std.Err.zP>
|z|[95%Conf.Interval]
-------------+----------------------------------------------------------------
group|-.6711683.0810558-8.280.000-.8300348-.5123018
_cons|6.780416.115050958.930.0006.554927.005912
Poisson回归分组的P<
0.001，可以认为甲饮料中的大肠杆菌比乙饮料多。
例7-10为了解2005年某地区甲乙两镇原发性肝癌的死亡情况，该地区慢性病控制中心在甲乙两镇进行了抽样调查。
甲镇调查了123443人，死于原发性肝癌的病例97人，乙镇调查了180274人，死于原发性肝癌的63人。
试问甲乙两镇原发性肝癌的死亡率有无差别？
Stata的数据为
97
123443
63
180274
分析程序和结果为：
poissonxgroup,exposure（w）
LRchi2
（1）=25.86
Loglikelihood=-6.1989818PseudoR2=0.6759
group|-.8102747.1618094-5.010.000-1.127415-.493134
_cons|5.174376.238977321.650.0004.705995.642763
W|（exposure）
0.001，可以认为甲镇原发性肝癌的死亡率比乙镇高。
升级会员 