第七章两独立样本定量资料的统计分析的Stata实现Word文档格式.docx

第七章两独立样本定量资料的统计分析的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，可以认为甲镇原发性肝癌的死亡率比乙镇高。