在职研究生医学统计学作业一峰王.docx

在职研究生医学统计学作业一峰王.docx

《在职研究生医学统计学作业一峰王.docx》由会员分享，可在线阅读，更多相关《在职研究生医学统计学作业一峰王.docx（25页珍藏版）》请在冰豆网上搜索。

在职研究生医学统计学作业一峰王

2013在职研究生医学统计学作业-一峰-王

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