医学统计学第四版各章例题SAS与STATA实现第三章.docx
《医学统计学第四版各章例题SAS与STATA实现第三章.docx》由会员分享,可在线阅读,更多相关《医学统计学第四版各章例题SAS与STATA实现第三章.docx(18页珍藏版)》请在冰豆网上搜索。
医学统计学第四版各章例题SAS与STATA实现第三章
医学统计学(第四版)各章例题SAS/STATA实现
(第三章)
例3J
若果巾1999年18岁男生身咼服从均数为167.7cm,标准岸为53cm的正态分布。
从该
正态分布NU67•7,5-32)cm总体中随机抽样100次即共抽取样本“二100个,每次样本含量川)
二10人得到每个样本均数Xj及标准差S/如图3-1和下表3-1所示。
图3-11999年某市18岁男生身aN(16775.32)抽样示意图
表3-1N(167•7•5.32)总体中100个随机样本的乂j.S丿和95%CI(«j=10)
样木号
乂j
Sj
9S%CZ
样本号
乂J
Sj
95%CZ
1
167.41
2.74
165.45
169.37
51
16&47
3・91
165,67
17L27
165.56
6.57
160.86
170.26
52
165.95
3・76
163.26
168.64
3
16S.20
5.36
164.37
172.03
53
168.87
5.77
164.74
173,00
4
166.67
4.81
163.24
170.11
*54
169.53
2.07
16S.05
171.00
5
164.89
5.41
161.02
168.76
55
166.10
5.58
162.11
170,10
6
166.36
4.50
16344
169.58
56
167.20
4.56
163.94
17047
7
166.16
4.04
163.27
169.05
57
170.50
7.66
165.02
175.98
S
169.11
5.71
165.02
173.19
58
166.44
4.93
162.91
169.97
9
167.17
8.26
16L27
173.08
59
16&68
4.52
16545
17L91
10
166.13
5.24
1623S
169.87
60
16&40
6.95
16343
173.37
11
167.71
6.42
163.12
172.31
61
171.21
630
166.70
175.72
12
16&68
5.93
164.44
172.92
62
170.33
4.34
167.23
173,44
13
166.83
3.69
16449
169.47
63
169.03
7.38
163.75
17431
14
169.62
4.81
166.18
173.06
64
16763
4.58
164.36
170,90
15
166.95
3.64
16435
169.56
65
16&66
3・33
166.27
171.04
16
170.29
4.91
166.78
173.80
66
16&84
2.78
166.85
170,83
17
169.20
5.72
165.11
173.30
67
169.31
5.31
165.51
173.11
1S
167.65
2.79
165.65
169.65
68
168.46
4.81
16302
171.90
19
166.51
5.39
162.65
170.36
69
168.60
5.4S
164,68
172.52
•20
163.28
3.19
16L00
165.57
70
168.47
5.05
164.86
172,09
21
166.29
4.95
162.75
169.84
71
165.6S
5.19
161.97
169.40
167.65
5.27
163S8
171.42
72
165.68
8.22
159.80
171.56
23
167.64
4.61
16435
170.94
73
168.03
4.89
164.53
171.53
24
172.61
7.74
167.07
178.15
74
169.37
5.00
16579
172.94
25
166.65
4.12
163.70
169.59
75
169.16
8.36
163,18
175,14
26
165.19
4.41
162.04
168.34
*76
171.27
4.99
167.71
174,84
27
168.80
7.68
16331
174.30
77
16&36
4.50
165,14
171.58
28
167如
2.58
166.14
169.83
78
168.50
3.55
165,96
171,04
29
168.41
3.43
165.95
170.86
79
168.08
5.33
164.27
171.90
30
167.75
7.53
162.36
173.13
80
165.51
4.71
162.14
168.88
♦31
164.25
4.30
161.17
167.33
S1
167.59
3.73
164.93
17026
32
166.42
5.19
16271
170.13
*82
171.12
4.40
167.98
174,27
33
166.90
4.41
163.74
170.05
83
165.92
5.11
162.26
169.58
34
166.77
4.34
163&
169&
84
16786
4.44
164.69
171,04
35
165.77
5.34
161.95
169.59
85
167.43
6.15
163.03
171.83
36
16442
6.63
15938
168.86
86
16790
6.13
163.51
172.28
37
169.83
4.20
166.82
172.84
87
167.59
633
163.06
172.12
38
165.16
4.01
162.29
168.02
88
167.74
4・60
16445
17L03
39
166.59
6.20
1623
171.03
89
167.40
8.27
161.49
173.32
40
165.65
3.56
163.10
168.20
90
167.1S
6.00
162.89
171.48
41
165.72
4.17
162.74
168.71
91
16643
3.87
163.66
169,21
42
166.22
7.44
1603
171.54
92
166.62
4.08
163.70
169.54
43
167.71
6.12
163.33
172.09
93
166.30
4.84
162.83
169.76
44
16725
5.24
163.50
170.99
94
169.70
5.26
165.94
17345
45
165.69
5.91
161.46
169.92
95
169.17
632
164.65
173.69
46
169.06
5.65
165.03
173.10
96
167.89
6.07
163.54
172.23
47
16&76
6.14
16436
173.15
97
167.48
6.03
163.16
171.79
4S
16&64
4.54
16539
171.89
98
169.93
4.80
166.50
173.37
49
167.72
3.82
164.99
170.45
99
16940
5.57
16342
17339
50
170.39
4.15
16742
173.35
100
165.69
5.09
162,06
16933
*:
表该样本资料算得的可信区间未包含已知总体均数167.7cm
例3•5
某医生测量了36洛从事铅作业男性工人的血红蛋白含量,算得加均数为130.83或L,标准
差为25-74g/L.问从事铅作业工人的血红蛋白是否不同于正常成年男性平均值140gzL?
⑴建立检验假设,确定检验水准
Ho:
严3=140g/L,即从事铅作业的男性工人平均血红蛋白含量与正常成年男性平均值
相等
H1:
/岸MF140弓L即从事铅作业的男性工人平均血红蛋白含量与正常成年男性平均值不
0=0.05
⑵讣算检验统il嗤
本例«=36>乂=130-83g/L,425-74g/L,“o=140g/L。
按公式(3-15)
仟36•1=35
J30.83一斗」38
25.74/后
⑶确定p值,做岀推断结论
以仟35、f=—2.138=2.138査附表2的r界值表,因"0,05/235<2.138V’0•02/2
-35♦
故双尾概率0.02
.接受Hi,有统讣学意义。
结合本题可认为从事铅作业的男性工人平均血红蛋白含量低于正常成年男性。
(1)SAS实现
SAS代码如下:
dataex3_5;
n=36;
s_m二130.83;
std=25・74;
p_m二140;
df=n-l:
t二abs(s_m-p_m)/(std-sqrt(a));p=(1-probt(t,df))*2;
procprint:
artp:
结果如下:
ObsP
■sumhb
(2)STATA实现
Variable
Obs
MeanSid,Dev,
Mm
Max
hb|
36
130.833325.74102
88
178
・tiesthb=140
Oiie~sanq>let
■test
\toable
1
Obs
Mean
Std・Err
Std.Dev
[95%Conf.Intenal]
hb|
36
130.8333
4.29017
25.74102
122.1238139.5428
mean=meaii(hb)-2・1367
Ho:
mean=140
例3•6
degreesoffreedom=35
Ha:
mean<140Ha:
mean!
=140Ha:
mean>140
Pr(T<1)二0.0198Pr(iT|>|tD=0.0397Pr(T>t)二0.9802
为比较两种方法对乳酸饮料屮脂肪含量测;4•结果是否不同,随机抽取了10份乳酸饮料
制品,分別用脂肪酸水解法和哥特里一罗紫法测定貝结果如表3•5第(1X3)栏。
问两法测崔结果是否不同?
表3-5两种方法对乳酸饮料屮脂肪含量的测崔结果(%)
編号
哥特里一罗紫法
脂肪酸水解法
差值d
<1)
(2)
(3)
(4X2)73)
1
0.840
0580
0.260
2
0.591
0509
0.0S2
3
0.674
0500
0.174
4
0.632
0316
0.316
L
0
0.687
0337
0.350
6
0978
0517
0.461
7
0.750
0,454
0.296
8
0730
0512
0.21S
9
L200
0997
0.203
10
0.870
0506
0.364
2.724
⑴建立检验假设,确定检验水准
Hz%即两种方法的测定结果相同
沟耳工0•即两种方法的测企结果不同
0=0.05
⑵汁算检验统订•量
本例《二10・加]2・724,"二0•8483.7二ZJ/N二2-724/10二0•2724。
0.8483-22空L
叫一二0.log?
10-1
按公式(3-⑹
v二iOT二9
尸。
27247.925,
0J087/710
(3)确企P值,做出推断结论
査附表2的r界值表得P<0.001o按0=0.05水准,拒绝尿.接受局,有统计学意义。
可认为两种方法对脂肪含量的测定结果不同,哥特里一罗紫法测总结果较高。
(1)SAS实现
SAS代码如下:
严用MEANS作配对资料两个样本均数比较的t检验*/
dataex3_6:
inputxlx2
d二xl•x2:
cards;
0.840
0.580
0.591
0.509
0.674
0.500
0.632
0.316
0.687
0.337
0.978
0.517
0.750
0.454
0.730
0.512
1.200
0.997
0.870
0.506
procmeanstprt:
vard
nm;
V
严用UNIVARIATE过程作配对资料两样本均数比较的t检验*/procunivariatedata二ex3_6;
vard:
nm;
MEANSPROCEDURE
分析变量:
d
t值Pr>|t!
7.93<0001
SAS系统
UNIVARIATEPROCEDURE变&矩
037900
四分位极差0.14700
位置检验:
MuO二0检验
Studentt符号符统株量…
t7.925976Pr>|t|
27.5
分位数(是义5)
分位数100%最大
估il值
值99%
0.4610
0.4610
95%
0.4610
90%
0.4125
75%Q3
03500
50%中位数
0.2780
25%QI
0.2030
10%
0.1280
5%
0.0820
1%
0.0820
0%最小值
0.0820
-p值
<0001
Pr>=|M0.0020
Pr>=|S|0.0020
(2)STATA实现
・ttestxl=x2
Ha:
mean<
Pairedttest
0
Ha:
mean!
=
0
Ha:
mean>0
'z/anable!
Obs
Mean
Std・Err
StdDe\•
(95%ConfInieml]
Xll
10
.7952
・0583003
・1843618
・6633155・9270845
x2
10
.5228
・0588125
・1859814
•389757•655843
diff
10
.2724
・034368
・1086812
1946542•3501458
niean(diff)=
mean(xl
■x2)
t=7,9260
Ho:
meaa(di3)
二0
degrees
offreedom=9
Ha:
<0
Ha
:
meaii(dif1£)
■
0
Ha:
mean(difi)>0
Pr(Tp论
>t|)=0,0000
PrCT>1)=0.0000
T
・gend=xl"x2
・ttestd二0
Oiie-sanq>lettest
例3•7
为研究国产四类新药阿卡波糖胶囊的降血糖效果,某医院用40名n型糖尿病病人进行同期随机对照试验。
研究考将这些病人随机等分到试验组(用阿卡波糖胶囊)和对照组(用拜唐苹胶囊),分别测得试验开始前和8周时的空腹血糖,算得空腹血糖卜•降值见表3•6,能否认
为该国产四类新药阿卡波糖胶囊与拜膚苹胶囊对空腹血糖的降糖效果不同?
表3-6阿卡波糖胶囊组和拜唐苹胶囊组空腹血糖卜•降值(mm&L)
阿卡波糖胶囊组必
-070
•5.60
2.00
2.80
0.70
3.50
4.00
5•80
7.10
•050
(hi=20)
2•50
•1.60
1.70
3・00
0.40
4.50
4.60
2•50
6.00
-1.40
拜唐苹胶囊组上
3・70
6.50
5.00
5.20
0.80
0.20
0.60
340
6.60
-1.10
(”2•20)
6.00
3.80
2.00
L60
2.00
2.20
1.20
3.10
1.70
-2.00
⑴建立检验假设,确定检验水准
Hz“LM,即阿卡波糖胶囊组与拜唐苹胶囊组空腹血糖下降值的总体均数相等
Hz“iHM,即阿卡波糖胶囊组与拜膺苹胶囊组空腹血糖下降值的总体均数不等
0=0.05
⑵讣算检验统il嗤
今算得阿卡波糖胶囊组空腹血糖下降值均数X|二2•0650(minoLT),标准差
Si=3-0601(uimol/L);拜唐苹胶囊组空腹血糖卞降值均数2.6250(nmiolL),标准差
S,=24205(inmol/L)a按公式(3•18)
乂一¥()650—2.6250—冋JSf+S;J3.060F+2.4205,
1(q-l)S:
+(n2-l)S;(1I1z、/、
v"i+n2-2二2(«-——1)二2(20-1)二38
20
⑶确定P值,做出推断结论
以1二38、t二
:
-0-642二0-642,査附表2的/界值表得P>0.50.按a=(M)5水准,
不拒绝%无统计学意义。
还不能认为阿卡波糖胶囊与拜腐苹胶囊对空腹血糖的降糖效果不同。
(1)SAS实现
SAS代码如下:
dataex3_7:
inputXif_n_〈21thenc=l;
elsec=2:
cards;
-070-5,602,002,800.703.504・005.807.10-0・50
2.50-L60L703.000404.504602.506.00-1.40
3・706.505.005200.800.200.603.406.60-1.10
6.003.802.001.602.002.201.203.101.70-2・00
procnesi;
var
X;
class
c;
run:
结果如下:
TheTTESTProcedure
「able:
X
Mean
StdErr
Minimum
2010650
3.0601
0.6S43H5.600
202,6250
2.4205
05412
■2.000
Method
■05600
2.7589
0.8724
Mean95%CLMean
2,06500632S3.4972
2,6250L492237578
Di£F(l-2)Pooled-05600H2.32621.20622-
7589
Di£f(l-2)SatterthTT-0560023293L2093
aite
Metiiod\*aiicesDFPr>|tl
Maximum
7,1000
6.6000
StdDev95%CLStdDev
3,06012.32724.4696
24205LS40S3,5353
225473.5556
Pooled
Equal
38-0,640.5248
Unequal36.086064
0.5250
EqualityofX*nances
Metiiod
MumDF
DenDF叭yuePr>F
FoldedF
19
19
1.600.3153
⑵STATA实现
・nestxl=x2jinp
Ha:
diff!
二0Ha:
diff>0
Pr(TT>ND=0.5248Pr(T>t)=0.7376
l\vo-samplettestwichequalvariances
Variable
Obs
Mean
Std・Err・
Std.Dev
[95%ConfImervalJ
xl|
20
2.065
・6842697
3,060147
•6328071
3.49刀93
x2
20
2.625
.54124
2420499
1.492172
3.757828
combmed
40
2.345
•4329231
2.738046
1.46933
3.22067
&ff
••56
•8724482
-2326179
1206179
06419
diff=mean(xl)-inean(x2)
・stackxlx2,
into(x)・rename
・stackgroup・ttestx・
by(group)
l\vo~samplettestuirhequal
variances