Question

我的任务是遇到以下问题：

在服用阿司匹林的1000名受试者中，80名患有心脏病，65名患有中风。在安慰剂总共2,000名受试者中，240名患有心脏病，165名患有中风。

我被问到阿司匹林治疗心脏病和中风是否有显着益处。对于两种结果中的每种结果，使用阿司匹林的RR是多少？

我的主要问题是设置数据线。这是我到目前为止，但我的输出窗口看起来不正确。

另一个问题是弄清楚如何解释不同的样本量以及某人可能患有心脏病和中风的事实。

DATA ODDS;
   INPUT OUTCOME $ EXPOSURE $ COUNT;
DATALINES;
HeartAttack 1-Yes 80
HeartAttack 2-No 240
Stroke 1-Yes 65
Stroke 2-No 165
;
PROC FREQ DATA=ODDS;
   TITLE "Odds Ratio Aspirin";
   TABLES EXPOSURE*OUTCOME / CHISQ CMH;
   WEIGHT COUNT;
RUN;

编辑1：

DATA ODDS;
   INPUT OUTCOME $ EXPOSURE $ COUNT;
DATALINES;
HeartAttack 1-Yes 80
HeartAttack 2-No 240
NoHeartAttack 1-Yes 920
NoHeartAttack 2-No 1760
Stroke 1-Yes 65
Stroke 2-No 165
NoStroke 1-Yes 935
NoStroke 2-No 1835
;
PROC FREQ DATA=ODDS;
   TITLE "Odds Ratio Aspirin";
   TABLES EXPOSURE*OUTCOME / CHISQ CMH;
   WEIGHT COUNT;
RUN;

Answer 1

安慰剂为您提供预期的分配。我们必须分别处理中风和心脏病发作，因为没有相互作用的数据。（如果没有相互作用，我们预计会有少数患者同时进行，但可能存在负面或掩蔽的相互作用，例如，如果心脏病发作是致命的，或者可能存在累积的相互作用，心脏病发作通常先发生中风或反之亦然）。我们无法回答任何问题。

一旦你达到了预期，那就是两个卡方测试，每个测试有两个箱子。没有一个卡方测试有4个箱子。

（如果你想从头开始编写卡方显着性测试代码，我会在我的书中插入一个插件，而不使用任何连接表）。

Answer 2

这就是您的数据和代码应该是什么样子。您可能需要在TABLES语句中翻转顺序，以便根据您的情况适当计算相对风险。我不打算检查是否是这种情况，因为如果需要你可以轻松改变。

DATA HeartAttack;
   INPUT OUTCOME $ EXPOSURE $ COUNT;
DATALINES;
HeartAttack 1-Yes 80
HeartAttack 2-No 240
NoHeartAttack 1-Yes 920
NoHeartAttack 2-No 1760
;
PROC FREQ DATA=HeartAttack;
   TITLE "Heart Attack Odds Ratio Aspirin";
   TABLES EXPOSURE*OUTCOME / CHISQ CMH Relrisk;
   WEIGHT COUNT;
RUN;

DATA Stroke;
   INPUT OUTCOME $ EXPOSURE $ COUNT;
DATALINES;
Stroke 1-Yes 65
Stroke 2-No 165
NoStroke 1-Yes 935
NoStroke 2-No 1835
;
PROC FREQ DATA=Stroke;
   TITLE "Stroke Odds Ratio Aspirin";
   TABLES EXPOSURE*OUTCOME / CHISQ CMH Relrisk;
   WEIGHT COUNT;
RUN;

如果你有原始数据我会建议使用它而不是首先创建数据。它留下了错误的空间，然后你也可以处理交互。

确定阿司匹林是否明显优于安慰剂的程序

2 个答案: