本稿では、平均発生率を指標とするコホート研究における平均発生率比の信頼区間を行っています。
なお、閲覧にあたっては、以下の点にご注意ください。
- スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。
- 曝露(発症)状況を表す右下の添え字は、「0」である場合($n_0,\pi_0$ など)や「2」である場合($n_2,\pi_2$ など)がありますが、どちらも「非曝露群(コントロール群)」を表しています。
- 漸近的な性質を用いる際は、①中心極限定理が成り立つ、②漸近分散を推定する際に、母数をその一致推定量で置き換えることができるということが成り立つと仮定しています。
【定理】平均発生率比の信頼区間
【定理】
平均発生率比の信頼区間
Confidence Interval for Incidence Rate Ratio
平均発生率を指標とするマッチングなしのコホート研究における平均発生率比と対数平均発生率比を \begin{gather} \delta=\mathrm{IRR}=\frac{\lambda_1}{\lambda_0} \quad \hat{\delta}=\mathrm{\widehat{IRR}}=\frac{{\hat{\lambda}}_1}{{\hat{\lambda}}_0}\\ \theta=\log{\mathrm{IRR}}=\log{\frac{\lambda_1}{\lambda_0}} \quad \hat{\theta}=\log{\mathrm{\widehat{IRR}}}=\log{\frac{{\hat{\lambda}}_1}{{\hat{\lambda}}_0}} \end{gather} とするとき、 母集団の平均発生率比の $100 \left(1-\alpha\right)\%$ 信頼区間は、漸近的に \begin{gather} 100 \left(1-\alpha\right)\%\ \mathrm{C.I.}= \left[\delta_L,\delta_U\right]\\ \delta_L=\mathrm{exp} \left(\theta_L\right) \quad \delta_U=\mathrm{exp} \left(\theta_U\right) \end{gather} で与えられる。 ただし、 \begin{gather} \theta_L=\hat{\theta}-\hat{\sigma} \cdot Z_{0.5\alpha} \quad \theta_U=\hat{\theta}+\hat{\sigma} \cdot Z_{0.5\alpha}\\ {\hat{\sigma}}^2=\frac{1}{d_1}+\frac{1}{d_0} \quad \hat{\sigma}=\mathrm{\widehat{S.E.}} \left(\hat{\theta}\right) \end{gather}
導出
標本対数平均発生率比の漸近分布は、 \begin{align} \hat{\theta}\xrightarrow[]{d}\mathrm{N} \left[\log{\frac{\lambda_1}{\lambda_0}},\frac{1}{\lambda_1T_1}+\frac{1}{\lambda_0T_0}\right] \end{align} 標本対数平均発生率比の分散の一致推定量は、 \begin{align} \hat{V} \left(\hat{\theta}\right)=\frac{1}{d_1}+\frac{1}{d_0} \end{align} 標本対数平均発生率比の漸近的な標準誤差 $\mathrm{\widehat{S.E.}} \left(\hat{\theta}\right)=\hat{\sigma}$ は、 \begin{gather} \hat{\sigma}=\sqrt{\frac{1}{d_1}+\frac{1}{d_0}} \end{gather} 母集団の対数平均発生率比の $100 \left(1-\alpha\right)\%$ 信頼区間は、 \begin{gather} \hat{\theta}-\hat{\sigma} \cdot Z_{0.5\alpha} \le \log{\delta} \le \hat{\theta}+\hat{\sigma} \cdot Z_{0.5\alpha} \end{gather} ここで、スペースの節約のために以下のようにおいて、 \begin{gather} \theta_L=\hat{\theta}-\hat{\sigma} \cdot Z_{0.5\alpha} \quad \theta_U=\hat{\theta}+\hat{\sigma} \cdot Z_{0.5\alpha} \end{gather} 逆変換を行うと、母集団の平均発生率比の $100 \left(1-\alpha\right)\%$ 信頼区間は、 \begin{gather} \mathrm{exp} \left(\theta_L\right) \le \delta \le \mathrm{exp} \left(\theta_U\right) \end{gather} $\blacksquare$
参考文献
- ケネス・ロスマン 著, 矢野 栄二, 橋本 英樹, 大脇 和浩 監訳. ロスマンの疫学. 篠原出版新社, 2013, p.236-238
- 丹後 俊郎, 松井 茂之 編集. 医学統計学ハンドブック. 朝倉書店, 2018, p.509
0 件のコメント:
コメントを投稿