第十八章病例对照研究实习指导(定)

?

m

i

x

i

2

?

(

< br>?

m

i

x

i

)

2

]

i

?

0

i

?

0

k

k

二、应用< /p>

成组设计病例

-

对照研究资料的分析步骤：

(1)

整理频数表。

(2)

计算优势比。

(3)

对优势比做假设检验。

(4)

估计优势比的置信区间。

(5)

如果危险因素是多水平等级变量变，可以进一步做趋势检验。

(6)

如果存在混杂因素，可以用分 层分析或后面介绍的多因素

logistic

回归模

第三节

匹配设计资料的分析

18-3

一、计算

优势比的估计：

OR

?

优势比的假设检验：

?

2

?

b

c

?

b

?

c

?

1

?

< br>

?

2

?

b

?

c

?

优势比的

95%

置信 区间估计用

Miettinen

法。

二、应用

配对设计病例对照研究资料的分析步骤：

(1)

整理频数表。

(2)

估计优势比。

(3)

对优势比做假设检验。

(4)

估计优势比的置信区间。

1: M

或

N: M

< br>配伍设计资料、多个水平危险因素的分析和混杂因素的校正用

后面介绍的条件

logistic

回归模型更为方便，这里不要求掌握。

[

案例讨论参考答案

]

案例

18-1

1.

该研 究设计是在现况调查的基础上作病例

-

对照研究资料分析，

但

该研究没有严格按病例

-

对照研究进行设计，特别是对照的选择不能代表产生病

例的人群，存在严重的偏 倚

,

例如，病例组与对照组的年龄、性别、肥胖和体力

劳动可能有明显差别。

2.

由于在设计上的缺陷，病例组与对照组在一些重要的混杂因素方面存在不均

衡 ，因此在分析中必须校正这些混杂因素的影响，如性别、年龄、职业或体力

劳动史等。若 混杂因素只是一、两个，可以用分层分析方法校正；若混杂因素

较多，可采用第

19

章的

logistic

回归。

[

电脑实验及结果解释

]

实验

18-1

Mantel- Haenszel

分层分析

18-4

程序

18-1

Mantel- Haenszel

分层分析及说明

行号

01

02

03

04

05

20

21

22

23

24

程

序

DATA

mh;

INPUT fre age drink id;

CARDS;

5

1

1

0

5

1

2

0

… …

;

PROC

FREQ

;

TABLE age*id*drink/CHISQ CMH

；

WEIGHT fre;

RUN

;

说

明

建立< /p>

SAS

数据集

mh

；

定义并输入变量；

调用频数表

FREQ

过程；

定义表

,

?

2

检验和分层分 析；

指定权重变量；

运行程序；

运行结果：

Output

窗口：

Table 1 of id by drink Table 2 of id by drink

Controlling for age=1 Controlling for age=2

id drink id drink

Frequency| Frequency|

Row Pct | 1| 2| Total Row Pct | 1| 2| Total

---------+--------+--------+ ---------+--------+--------+

0 | 5 | 5 | 10 0 | 25 | 21 | 46

| 50.00 | 50.00 | | 54.35 | 45.65 |

---------+--------+--------+ ---------+--------+--------+

1 | 35 | 270 | 305 1 | 29 | 138 | 167

| 11.48 | 88.52 | | 17.37 | 82.63 |

---------+--------+--------+ ---------+--------+--------+

Total 40 275 315 Total 54 159 213

Table 3 of id by drink Table 4 of id by drink

Controlling for age=3 Controlling for age=4

id drink id drink

Frequency| Frequency|

Row Pct | 1| 2| Total Row Pct | 1| 2| Total

---------+--------+--------+ ---------+--------+--------+

0 | 42 | 34 | 76 0 | 24 | 44 | 68

| 55.26 | 44.74 | | 35.29 | 64.71 |

---------+--------+--------+ ---------+--------+--------+

1 | 27 | 139 | 166 1 | 18 | 119 | 137

| 16.27 | 83.73 | | 13.14 | 86.86 |

---------+--------+--------+

---------+--------+--------+

Total 69 173 242 Total 42 163 205

以上是四个年龄层的频数表。

Summary Statistics for id by drink

Controlling for age

Cochran-Mantel-Haenszel Statistics (Based on Table Scores)

Statistic

Alternative Hypothesis

DF

Value

Prob

--- -------------------------------------------------- ---------------------------

1

Nonzero Correlation

1

86.5130

<.0001

2

Row Mean Scores Differ

1

86.5130

<.0001

3

General Association

1

86.5130

<.0001

以上是

M-H

分层

?

检验结果。< /p>

Estimates of the Common Relative Risk (Row1/Row2)

Type of Study

Method

Value

95% Confidence Limits

----------------------------------------- --------------------------------------------

------------

Case-Control

Mantel-Haenszel

5.2873

3.6505

7.6580

(Odds Ratio)

Logit

5.3480

3.6949

7.7406

18-5

2

Cohort

Mantel-Haenszel

3.1750

2.4829

4.0599

(Col1 Risk)

Logit

3.2533

2.5590

4.1359

Cohort

Mantel-Haenszel

0.6093

0.5298

0.7008

(Col2 Risk)

Logit

0.6388

0.5587

0.7305

以上是优势比的估计 ，选用

OR

值估计。

Breslow-Day Test for

Homogeneity of the Odds Ratios

--------------------------------------

Chi-Square

1.8676

DF

3

Pr > ChiSq

0.6003

以上是各层相对危险性的齐性检验，差异没有统 计学意义，说明各年龄组吸烟与疾病关联程度相同。

Total Sample Size = 975

实验

18-2

过度匹配

程序

18-2

过度匹配电脑实验的

SAS

程序及说明

程

序

行号

01

DATA

case1;

02

DO i=

1

TO

100

;

03

no=i; cough=

1

; smoke1=

1

; smoke2=

1

;

04

IF UNIFORM(

0