在职研究生医学统计学作业一峰王.docx
《在职研究生医学统计学作业一峰王.docx》由会员分享，可在线阅读，更多相关《在职研究生医学统计学作业一峰王.docx（25页珍藏版）》请在冰豆网上搜索。
在职研究生医学统计学作业一峰王
2013级南昌大学医学院在职研究生（南昌半脱产）《医学统计学》作业
姓名：
学号：
得分：
题目：
以某治疗胃炎疼痛症状的临床试验数据库为例，评价A和B两种治疗胃炎疼痛症状药物的疗效是否有差别。
该临床试验采用完全随机设计，数据库包括人口学特征、体检及基线指标（共13项）和疗效指标（2项），其中，人口学特征、体检及疗前基线指标包括：
年龄、性别、呼吸、心率、收缩压、舒张压、体温、身高、体重、病程、过敏史、治疗史、治疗前胃部疼痛评分；疗效指标包括：
治疗前后胃部疼痛评分差值、综合疗效（痊愈、显效、有效、无效）。
要求：
1）拟一份统计分析计划（字数不限），内容包括分析哪些指标来评价A和B两种药物疗效的差别、采用何种统计指标或方法分析这些指标。
答：
在此次实验中实验者随机抽取了100例实验对象并采用了完全随机设计方法来评价A和B两种治疗胃炎疼痛症状药物的疗效是否有差别。
因此首先应该判定随机分组的两组病人各项人口学特征、体检及基线指标的可比性分析；分析判定各项指标是否具有可比性；分析各项指标的类别并选用不同的统计学检验方法描述并评价各项指标及疗效。
统计分析计划：
通过分析各项指标的类别及可比性后，对治疗前后胃部疼痛评分差值、综合疗效分别进行t检验或单因素方差分析及秩和检验分析A、B两组药物治疗效果有无差别，选取当P＜0.05时，认为两组药物在疗效上有统计学差异，当P＞0.05时，认为两组药物在疗效上无显著性差异。
2）对每一个须分析指标进行统计描述和推断，列出主要的描述性指标及统计推断结果（检验统计量及P值）。
答：
（一）、各处理因素的统计描述：
1、根据题目要求创建两个纵列分别命名为“group1”、“d”；
2、各项指标分别绘制成P-P图分析数据分布形态；
经统计SPSS软件绘制出的P-P图，我们可以得出以下结论：
符合近似正态分布的计数指标有：
age、heart、breathe、dp、temp、weight、height、tt_1、d；另有sp、bc不符合正态分布。
2、利用统计SPSS软件绘制图标描述各项指标：
各项计数/量指标的均数、方差、标准差、四分位数间距、中位数表格；
Statistics
group
age
heart
breathe
sp
dp
temp
weight
height
bc
tt_1
d
A
N
Valid
49
49
49
49
49
49
49
49
49
49
49
Missing
0
0
0
0
0
0
0
0
0
0
0
Mean
45.45
76.51
19.02
117.65
76.20
36.563
59.79
164.71
27.43
16.86
11.6939
Median
47.00
78.00
19.00
120.00
76.00
36.500
60.00
165.00
14.00
15.00
11.0000
Std.Deviation
12.075
8.778
1.601
8.059
8.818
.2138
9.858
6.252
47.185
5.408
5.60953
Percentiles
25
36.00
70.50
18.00
110.00
70.00
36.400
52.50
160.00
7.00
12.00
8.0000
50
47.00
78.00
19.00
120.00
76.00
36.500
60.00
165.00
14.00
15.00
11.0000
75
54.00
80.00
20.00
120.00
80.00
36.800
65.50
169.00
24.00
20.50
16.0000
B
N
Valid
50
50
50
50
50
50
50
50
50
50
50
Missing
0
0
0
0
0
0
0
0
0
0
0
Mean
46.48
72.86
19.24
116.10
75.00
36.532
58.06
162.52
26.30
18.50
13.1600
Median
45.50
73.50
20.00
120.00
75.00
36.500
59.00
161.50
13.00
19.00
13.0000
Std.Deviation
12.430
7.217
1.393
9.552
6.350
.1834
8.082
6.377
38.733
6.329
5.65093
Percentiles
25
36.00
69.50
18.00
110.00
70.00
36.400
50.00
158.00
8.50
12.75
8.7500
50
45.50
73.50
20.00
120.00
75.00
36.500
59.00
161.50
13.00
19.00
13.0000
75
57.25
79.00
20.00
120.00
80.00
36.700
62.25
168.00
24.00
24.25
17.2500
标准差描述，不符合近似正态分布指标用中位数及四分位数间距描述；
age
heart
breathe
dp
temp
weight
height
tt_1
d
A均数
45.45
76.51
19.02
117.65
76.20
36.563
59.79
164.71
27.43
B均数
46.48
72.86
19.24
116.10
75.00
36.532
58.06
162.52
26.30
A标准差
12.075
8.778
1.601
8.059
8.818
.2138
9.858
6.252
47.185
B标准差
12.430
7.217
1.393
9.552
6.350
.1834
8.082
6.377
38.733
A组中位数
B组中位数
A组四分位数间距
B四分位数间距
Sp
120.00
120.00
120-110=10
120-110=10
bc
14.00
13.00
24.00-8.25=15.75
24.00-9.00=15.00
分类变量指标的描述：
CaseProcessingSummary
Cases
Valid
Missing
Total
N
Percent
N
Percent
N
Percent
group_1*gender
99
100.0%
0
.0%
99
100.0%
group_1*gms
99
100.0%
0
.0%
99
100.0%
group_1*zls
99
100.0%
0
.0%
99
100.0%
group_1*zylx
99
100.0%
0
.0%
99
100.0%
group_1*genderCrosstabulation
gender
Total
1
2
group_1
1.00
Count
22
27
49
%withingroup_1
44.9%
55.1%
100.0%
2.00
Count
16
34
50
%withingroup_1
32.0%
68.0%
100.0%
Total
Count
38
61
99
%withingroup_1
38.4%
61.6%
100.0%
group_1*gmsCrosstabulation
gms
Total
1
2
group_1
1.00
Count
49
0
49
%withingroup_1
100.0%
.0%
100.0%
2.00
Count
44
6
50
%withingroup_1
88.0%
12.0%
100.0%
Total
Count
93
6
99
%withingroup_1
93.9%
6.1%
100.0%
group_1*zlsCrosstabulation
zls
Total
1
2
group_1
1.00
Count
40
9
49
%withingroup_1
81.6%
18.4%
100.0%
2.00
Count
40
10
50
%withingroup_1
80.0%
20.0%
100.0%
Total
Count
80
19
99
%withingroup_1
80.8%
19.2%
100.0%
group_1*zylxCrosstabulation
zylx
Total
1
2
3
4
group_1
1.00
Count
2
23
21
3
49
%withingroup_1
4.1%
46.9%
42.9%
6.1%
100.0%
2.00
Count
5
21
23
1
50
%withingroup_1
10.0%
42.0%
46.0%
2.0%
100.0%
Total
Count
7
44
44
4
99
%withingroup_1
7.1%
44.4%
44.4%
4.0%
100.0%
（二）、实验统计的推断分析：
1、利用t检验或单因素方差齐性检验分析各计数指标是否方差齐；
SPSS统计软件对计数指标进行t检验可得下表：
Levene'sTestforEqualityofVariances
t-testforEqualityofMeans
F
Sig.
t
df
Sig.（2-tailed）
MeanDifference
Std.ErrorDifference
95%ConfidenceIntervaloftheDifference
Lower
Upper
age
Equalvariancesassumed
.268
.606
-.419
97
.676
-1.031
2.464
-5.920
3.858
Equalvariancesnotassumed
-.419
96.993
.676
-1.031
2.463
-5.919
3.857
gender
Equalvariancesassumed
5.549
.021
-1.318
97
.191
-.129
.098
-.323
.065
Equalvariancesnotassumed
-1.317
96.310
.191
-.129
.098
-.323
.065
heart
Equalvariancesassumed
.384
.537
2.262
97
.026
3.650
1.614
.448
6.853
Equalvariancesnotassumed
2.258
92.777
.026
3.650
1.617
.439
6.861
breathe
Equalvariancesassumed
.005
.946
-.729
97
.468
-.220
.301
-.818
.379
Equalvariancesnotassumed
-.727
94.635
.469
-.220
.302
-.819
.380
dp
Equalvariancesassumed
.102
.750
.781
97
.437
1.204
1.542
-1.856
4.265
Equalvariancesnotassumed
.778
87.143
.438
1.204
1.547
-1.871
4.279
temp
Equalvariancesassumed
3.095
.082
.781
97
.436
.0313
.0400
-.0481
.1107
Equalvariancesnotassumed
.780
94.209
.437
.0313
.0401
-.0483
.1108
weight
Equalvariancesassumed
2.341
.129
.956
97
.342
1.730
1.810
-1.863
5.322
Equalvariancesnotassumed
.954
92.672
.343
1.730
1.814
-1.872
5.332
height
Equalvariancesassumed
.215
.644
1.729
97
.087
2.194
1.269
-.325
4.714
Equalvariancesnotassumed
1.729
97.000
.087
2.194
1.269
-.325
4.713
tt_1
Equalvariancesassumed
2.804
.097
-1.387
97
.169
-1.643
1.184
-3.993
.708
Equalvariancesnotassumed
-1.389
95.247
.168
-1.643
1.182
-3.990
.704
d
Equalvariancesassumed
.031
.860
-1.295
97
.198
-1.46612
1.13183
-3.71249
.78024
Equalvariancesnotassumed
-1.295
96.983
.198
-1.46612
1.13174
-3.71232
.78008
上表可知heart指数的P值为0.025＜α=0.05，故选用协方差对其进行校正分析。
heart指数协方差如下：
TestsofBetween-SubjectsEffects
DependentVariable:
d
Source
TypeIIISumofSquares
df
MeanSquare
F
Sig.
CorrectedModel
98.274（a）
2
49.137
1.541
.219
Intercept
49.031
1
49.031
1.538
.218
heart
37.525
1
37.525
1.177
.281
group1
81.496
1
80.496
2.528
.115
Error
3060.474
96
31.880
Total
18615.000
99
CorrectedTotal
3158.747
98
aRSquared=.031（AdjustedRSquared=.011）
从上表可知：
经协方差校正后heart指数的F值为1.541，P值为0.219＞α=0.05。
综上所述经t检验后可知：
age、heart、breathe、dp、temp、weight、height、tt_1、d的P值均＞α=0.05，按α=0.05水准，可认为A、B两组资料方差齐；A、B两组药物治疗前后胃部疼痛评分差值d无统计学差别。
2、无序分类变量需选用卡方检验分析：
利用SPSS统计软件对无序分类变量gender、gms、zls进行卡方检验得下表有：
CaseProcessingSummary
Cases
Valid
Missing
Total
N
Percent
N
Percent
N
Percent
group_1*gender
99
100.0%
0
.0%
99
100.0%
group_1*gms
99
100.0%
0
.0%
99
100.0%
group_1*zls
99
100.0%
0
.0%
99
100.0%
group_1*zylx
99
100.0%
0
.0%
99
100.0%
group_1*gender
Crosstab
gender
Total
1
2
group_1
1.00
Count
22
27
49
ExpectedCount
18.8
30.2
49.0
2.00
Count
16
34
50
ExpectedCount
19.2
30.8
50.0
Total
Count
38
61
99
ExpectedCount
38.0
61.0
99.0
Chi-SquareTests
Value
df
Asymp.Sig.（2-sided）
ExactSig.（2-sided）
ExactSig.（1-sided）
PointProbability
PearsonChi-Square
1.741（b）
1
.187
.218
.133
ContinuityCorrection（a）
1.238
1
.266
LikelihoodRatio
1.746
1
.186
.218
.133
Fisher'sExactTest
.218
.133
Linear-by-LinearAssociation
1.723（c）
1
.189
.218
.133
.070
NofValidCases
99
aComputedonlyfora2x2table
b0cells（.0%）haveexpectedcountlessthan5.Theminimumexpectedcountis18.81.
cThestandardizedstatisticis1.313.
group_1*gms
Crosstab
gms
Total
1
2
group_1
1.00
Count
49
0
49
ExpectedCount
46.0
3.0
49.0
2.00
Count
44
6
50
ExpectedCount
47.0
3.0
50.0
Total
Count
93
6
99
ExpectedCount
93.0
6.0
99.0
Chi-SquareTests
Value
df
Asymp.Sig.（2-sided）
ExactSig.（2-sided）
ExactSig.（1-sided）
PointProbability
PearsonChi-Square
6.259（b）
1
.012
.027
.014
ContinuityCorrection（a）
4.329
1
.037
LikelihoodRatio
8.577
1
.003
.027
.014
Fisher'sExactTest
.027
.014
Linear-by-LinearAssociation
6.196（c）
1
.013
.027
.014
.014
NofValidCases
99
aComputedonlyfora2x2table
b2cells（50.0%）haveexpectedcountlessthan5.Theminimumexpectedcountis2.97.
cThestandardizedstatisticis2.489.
group_1*zlsCrosstab
zls
Total
1
2
group_1
1.00
Count
40
9
49
ExpectedCount
39.6
9.4
49.0
2.00
Count
40
10
50
ExpectedCount
40.4
9.6
50.0
Total
Count
80
19
99
ExpectedCount
80.0
19.0
99.0
Chi-SquareTests
Value
df
Asymp.Sig.（2-sided）
ExactSig.（2-sided）
ExactSig.（1-sided）
PointProbability
PearsonChi-Square
.043（b）
1
.837
1.000
.520
ContinuityCorrection（a）
.000
1
1.000
LikelihoodRatio
.043
1
.837
1.000
.520
Fisher'sExactTest
1.000
.520
Linear-by-LinearAssociation
.042（c）
1
.837
1.000
.520
.197
NofValidCases
99
aComputedonlyfora2x2table
b0cells（.0%）haveexpectedcountlessthan5.Theminimumexpectedcountis9.40.
cThe
升级会员 