胡纯严，胡良平，2*

（1. 军事科学院研究生院，北京 100850；2. 世界中医药学会联合会临床科研统计学专业委员会，北京 100029*通信作者：胡良平，E-mail：lphu927＠163.com）

在分析“人-时间资料”时，若基于“人-年数”，在所考察的处理因素分别处于“暴露”与“非暴露”水平下，求出两组受试对象各自的“发病密度”之后，再求出两“发病密度”之比，就获得了“率比”［1］，它类似于分析通常四格表资料所采用的效应指标“相对危险度”或“优势比”［2-5］；但在文献［6］中，直接将“率比”称为“相对危险度”。本文将介绍与“率比分析”有关的内容以及基于SAS 软件实现计算的方法。

1 分层人-时间资料中第i 层率比的点估计和置信区间估计

1.1 分层人-时间资料的列表格式

分层人-时间资料的列表格式见表1。

表1 分层人-时间资料的列表格式

1.2 第i层率比的点估计及置信区间估计

第i层“率比”是该层中“暴露水平”下的“发病密度（α1i/t1i）”与“非暴露水平”下的“发病密度（α2i/t2i）”之比值，相当于通常定性资料中的“相对危险度”或“优势比”［2-5］。计算方法见式（1）：

第i层率比RRi的100（1-α）%置信区间的计算见式（2）：

2 不同层间率比的齐性检验

为了将各层资料合并起来估计“合并率比”，资料需满足的前提条件是各层率比相等，这就需要进行不同层间率比的齐性检验［1，6］。设分层因素共有K个水平，齐性检验的检验假设如下：

H0：RR1=RR2=…=RRk；H1：至少两个层的RRi不同；α=0.05。

检验统计量见式（5）：

式（5）中，RRi为第i层率比［见式（1）］；Wi为第i层权重，RR为合并率比，其计算分别见式（6）、式（7）：

在式（6）中，Var［ln（RRi）］为ln（RRi）的方差，计算方法见式（8）：

在式（7）中，指数C的计算方法见式（9）：

3 合并率比的点估计与置信区间估计

当各层率比满足齐性要求时，可以对合并率比进行点估计和置信区间估计［1，6］。……