SAS proc freq 风险比

Question

我正在尝试使用 PROC FREQ 计算风险比率和相关的置信限。这在一般情况下工作正常，但在两个实验组的所有测试单元都遇到事件的情况下，风险比没有计算出来，而应该是一个。

考虑以下事项：

Treated event: 2
Treated no event: 0
Not-treated event: 2
Not-Treated no event: 0

Risk ratio = 2/(2+0) / 2/(2+0)  
       = 2/2  
       = 1

PROC FREQ 然而返回注意：grp by 表的行或列总和为零。没有为此表计算统计数据。

我使用了 zeros 选项，但这并不能解决问题。

  DATA temp;
      grp = 1;
      occur = 0;
      count = 2;
      OUTPUT;
      occur = 1;
      count = 0;
      OUTPUT;
      occur = 0;
      grp = 2;
      count = 2;
      OUTPUT;
      occur = 1;
      count = 0;
      OUTPUT;
    RUN;
    PROC FREQ DATA = temp ORDER = internal NOPRINT  ;
      WEIGHT  count /zeros;
      TABLES grp*occur   / RELRISK        ;
      OUTPUT OUT = odds RELRISK       ;
    RUN;