ジョン・ラチン（2020）『医薬データのための統計解析』 問題2.5 解答例

本稿は、ジョン・ラチン（2020）『医薬データのための統計解析』の「問題2.5」の自作解答例です。リスク差・対数相対リスク・対数オッズ比の漸近的性質に関する問題です。

なお、閲覧にあたっては、以下の点にご注意ください。

• スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。
• 曝露（発症）状況を表す右下の添え字は、「0」である場合（$n_0,{\pi }_0$ など）や「2」である場合（$n_2,{\pi }_2$ など）がありますが、どちらも「非曝露群（コントロール群）」を表しています。
• 漸近的な性質を用いる際は、①中心極限定理が成り立つ、②漸近分散を推定する際に、母数をその一致推定量で置き換えることができるということが成り立つと仮定しています。
• デルタ法を用いる際、剰余項（2次の項）が漸近的に無視できる（$0$に確率収束する）と仮定しています。
• 上述の参考書では、標準正規分布の上側 $100\alpha \mathrm{\%}$ 点を $Z_{1-\alpha }$ と表記していますが、本サイトでは、$Z_{\alpha }$ としています。そのため、参考書に載っている式の形式と異なる部分があります。
• 著作権の関係上、問題文は、掲載しておりません。上述の参考書をお持ちの方は、お手元にご用意してご覧ください。
• この解答例は、筆者が自作したものであり、公式なものではありません。あくまでも参考としてご覧いただければ幸いです。

 

問題2.5.1：曝露効果の指標の漸近分布（対立仮説）

二項分布の正規近似により、標本比率は漸近的に $$\begin{array}{c}{p}_1\sim N[{\pi }_1,\frac{{\pi }_1(1-{\pi }_1)}{n_1}]\\ p_2\sim N[{\pi }_2,\frac{{\pi }_2(1-{\pi }_2)}{n_2}]\end{array}$$

リスク差

標本リスク差を $\widehat{\mathrm{RD}}=p_1-p_2$ とすると、正規分布の再生性より、 $$\begin{array}{r}\widehat{\mathrm{RD}}\sim \mathrm{N}\{{\pi }_1-{\pi }_2,\frac{{\pi }_1(1-{\pi }_1)}{n_1}+\frac{{\pi }_2(1-{\pi }_2)}{n_2}\}\end{array}$$ 母比率を一致推定量である標本比率で置き換えると、漸近分散の一致推定量は、 $$\begin{array}{r}{\hat{\sigma }}_1^2=\frac{p_1(1-p_1)}{n_1}+\frac{p_2(1-p_2)}{n_2}\end{array}$$ $◼$

 

対数リスク比

まず、 $$\begin{array}{c}g\left({\pi }_i\right)=\log {\pi }_i\\ g\left(p_i\right)=\log p_i\end{array}$$ と変数変換する。 デルタ法を用いて、$g\left(p_i\right)$ を期待値 $E\left(p_i\right)={\pi }_i$ まわりでテイラー展開すると、$g\left(p_i\right)$ の1階微分は、 $$\begin{array}{r}{g}^{\prime }\left(p_i\right)=\frac{1}{p_i}\end{array}$$ よって、デルタ法における期待値と分散の公式より、 $$\begin{array}{rl}E\left\{g\left(p_i\right)\right\}& \cong E\left[g\left({\pi }_i\right)\right]\\ & =\log {\pi }_i\end{array}$$ $$\begin{array}{rl}V\left[g\left(p_i\right)\right]& \cong {\left\{g^{\prime }\left({\pi }_i\right)\right\}}^{2}V\left(p_i\right)\\ & =\frac{1}{{\pi }_i^2}\cdot \frac{{\pi }_i(1-{\pi }_i)}{n_i}\\ & =\frac{1-{\pi }_i}{n_i{\pi }_i}\end{array}$$ スラツキーの定理より、 $$\begin{array}{r}\log p_i\stackrel{d}{\to }\mathrm{N}[\log {\pi }_i,\frac{1-{\pi }_i}{n_i{\pi }_i}]\end{array}$$ このとき、標本対数リスク比は、 $$\begin{array}{r}\log \widehat{\mathrm{RR}}=\log \frac{p_1}{p_2}=\log p_1-\log p_2\end{array}$$ したがって、正規分布の再生性より、 $$\begin{array}{r}\log \widehat{\mathrm{RR}}\stackrel{d}{\to }\mathrm{N}[\log {\pi }_1-\log {\pi }_2,\frac{1-{\pi }_1}{n_1{\pi }_1}+\frac{1-{\pi }_2}{n_2{\pi }_2}]\end{array}$$ $◼$

 

対数オッズ比

まず、 $$\begin{array}{c}g\left({\pi }_i\right)=\log \frac{{\pi }_i}{1-{\pi }_i}\\ g\left(p_i\right)=\log \frac{p_i}{1-p_i}\end{array}$$ と変数変換する。 デルタ法を用いて、$g\left(p_i\right)$ を期待値 $E\left(p_i\right)={\pi }_i$ まわりでテイラー展開すると、$g\left(p_i\right)$ の1階微分は、 $$\begin{array}{rl}{g}^{\prime }\left(p_i\right)& =\frac{1-p_i}{p_i}\cdot \frac{d}{d{p}_i}\left(\frac{p_i}{1-p_i}\right)\\ & =\frac{1-p_i}{p_i}\cdot \frac{(1-p_i)-1}{{(1-p_i)}^2}\\ & =\frac{1}{p_i(1-p_i)}\end{array}$$ よって、デルタ法における期待値と分散の公式より、 $$\begin{array}{rl}E\left\{g\left(p_i\right)\right\}& \cong E\left[g\left({\pi }_i\right)\right]\\ & =\log \frac{{\pi }_i}{1-{\pi }_i}\end{array}$$ $$\begin{array}{rl}V\left[g\left(p_i\right)\right]& \cong {\left\{g^{\prime }\left({\pi }_i\right)\right\}}^{2}V\left(p_i\right)\\ & =\frac{1}{{\pi }_i^2{(1-{\pi }_i)}^2}\cdot \frac{{\pi }_i(1-{\pi }_i)}{n_i}\\ & =\frac{1}{n_i{\pi }_i(1-{\pi }_i)}\end{array}$$ スラツキーの定理より、 $$\begin{array}{r}\log p_i\stackrel{d}{\to }\mathrm{N}[\log \frac{{\pi }_i}{1-{\pi }_i},\frac{1}{n_i{\pi }_i(1-{\pi }_i)}]\end{array}$$ このとき、標本対数オッズ比は、 $$\begin{array}{r}\log \widehat{\mathrm{OR}}=\log \frac{p_1}{1-p_1}-\log \frac{p_2}{1-p_2}\end{array}$$ したがって、正規分布の再生性より、 $$\begin{array}{r}\log \widehat{\mathrm{OR}}\stackrel{d}{\to }\mathrm{N}[\log \frac{{\pi }_1}{1-{\pi }_1}-\log \frac{{\pi }_2}{1-{\pi }_2},\frac{1}{n_1{\pi }_1(1-{\pi }_1)}+\frac{1}{n_2{\pi }_2(1-{\pi }_2)}]\end{array}$$ $◼$

 

問題2.5.2：対数リスク比の漸近分散の一致推定量

母比率を一致推定量である標本比率で置き換えると、期待値と分散の一致推定量は、 $$\begin{array}{c}\hat{E}[\log \widehat{\mathrm{RR}}]=\log p_1-\log p_2\\ \hat{V}(\log \widehat{\mathrm{RR}})=\frac{1-p_1}{n_1{p}_1}+\frac{1-p_2}{n_2{p}_2}\end{array}$$ 分散についてさらに式を整理すると、$a=n_1{p}_1,b=n_2{p}_2$ より、 $$\begin{array}{rl}{\hat{\sigma }}_1^2& =\frac{1-p_1}{a}+\frac{1-p_2}{b}\\ & =\frac{1}{a}-\frac{p_1}{a}+\frac{1}{b}-\frac{p_2}{b}\\ & =\frac{1}{a}-\frac{1}{n_1}+\frac{1}{b}-\frac{1}{n_2}\end{array}$$ $◼$

 

問題2.5.3：対数オッズ比の漸近分散の一致推定量

母比率を一致推定量である標本比率で置き換えると、期待値と分散の一致推定量は、 $$\begin{array}{c}\hat{E}[\log \widehat{\mathrm{OR}}]=\log \frac{p_1}{1-p_1}-\log \frac{p_2}{1-p_2}\\ \hat{V}(\log \widehat{\mathrm{OR}})=\frac{1}{n_1{p}_1(1-p_1)}+\frac{1}{n_2{p}_2(1-p_2)}\end{array}$$ 分散についての式を部分分数分解すると、 $$\begin{array}{rl}{\hat{\sigma }}_1^2& =\frac{1}{n_1{p}_1}+\frac{1}{n_1(1-p_1)}+\frac{1}{n_2{p}_2}+\frac{1}{n_2(1-p_2)}\end{array}$$ ここで、各セルの値について、 $$\begin{array}{c}a=n_1{p}_1\\ c=n_1(1-p_1)\\ b=n_2{p}_2\\ d=n_2(1-p_2)\end{array}$$ したがって、 $$\begin{array}{r}{\hat{\sigma }}_1^2=\frac{1}{a}+\frac{1}{b}+\frac{1}{c}+\frac{1}{d}\end{array}$$ $◼$

 

問題2.5.4：リスク比・オッズ比の漸近分散

リスク比

二項分布の正規近似により、標本比率は漸近的に $$\begin{array}{c}{p}_1\sim N[{\pi }_1,\frac{{\pi }_1(1-{\pi }_1)}{n_1}]\\ p_2\sim N[{\pi }_2,\frac{{\pi }_2(1-{\pi }_2)}{n_2}]\end{array}$$ 標本比率ベクトルを $$\begin{array}{r}\mathit{p}=\left(\begin{array}{c}{p}_1\\ p_2\end{array}\right)\end{array}$$ 期待値ベクトルを $$\begin{array}{r}\mathit{\pi }=\left(\begin{array}{c}{\pi }_1\\ {\pi }_2\end{array}\right)\end{array}$$ 分散・共分散行列を $$\begin{array}{r}\mathbf{\Omega }=\left[\begin{array}{cc}\frac{{\pi }_1(1-{\pi }_1)}{n_1}& 0\\ 0& \frac{{\pi }_2(1-{\pi }_2)}{n_2}\end{array}\right]\end{array}$$ として、 $$\begin{array}{c}G\left(\mathit{\pi }\right)=\mathrm{RR}=\frac{{\pi }_1}{{\pi }_2}\\ G\left(\mathit{p}\right)=\widehat{\mathrm{RR}}=\frac{p_1}{p_2}\end{array}$$ と変数変換する。 多変量のデルタ法を用いて $G\left(\mathit{p}\right)$ を期待値 $E\left(\mathit{p}\right)=\mathit{\pi }$ まわりでテイラー展開すると、偏導関数ベクトルは、 $$\begin{array}{r}\mathit{H}\left(\mathit{\pi }\right)=\left(\begin{array}{c}\frac{G\left(\mathit{\pi }\right)}{\partial {\pi }_1}\\ \frac{G\left(\mathit{\pi }\right)}{\partial {\pi }_2}\end{array}\right)=\left[\begin{array}{c}\frac{1}{{\pi }_2}\\ -\frac{{\pi }_1}{{\pi }_2^2}\end{array}\right]\end{array}$$ 多変量のデルタ法の期待値と分散の公式より、 $$\begin{array}{r}E\left[G\left(\mathit{p}\right)\right]\cong G\left(\mathit{\pi }\right)=\frac{{\pi }_1}{{\pi }_2}\end{array}$$ $$\begin{array}{rl}V\left[G\left(\mathit{p}\right)\right]& =\left[\begin{array}{cc}\frac{1}{{\pi }_2}& -\frac{{\pi }_1}{{{\pi }_2}^2}\end{array}\right]\left[\begin{array}{cc}\frac{{\pi }_1(1-{\pi }_1)}{n_1}& 0\\ 0& \frac{{\pi }_2(1-{\pi }_2)}{n_2}\end{array}\right]\left[\begin{array}{c}\frac{1}{{\pi }_2}\\ -\frac{{\pi }_1}{{{\pi }_2}^2}\end{array}\right]\\ & =\left[\begin{array}{cc}\frac{1}{{\pi }_2}& -\frac{{\pi }_1}{{{\pi }_2}^2}\end{array}\right]\left[\begin{array}{c}\frac{{\pi }_1(1-{\pi }_1)}{n_1{\pi }_2}\\ -\frac{{\pi }_1(1-{\pi }_2)}{n_2{\pi }_2}\end{array}\right]\\ & =\frac{1}{{\pi }_2}\frac{{\pi }_1(1-{\pi }_1)}{n_1{\pi }_2}+\frac{{\pi }_1}{{{\pi }_2}^2}\frac{{\pi }_1(1-{\pi }_2)}{n_2{\pi }_2}\\ & =\frac{{\pi }_1(1-{\pi }_1)}{n_1{\pi }_2^2}+\frac{{\pi }_1^2(1-{\pi }_2)}{n_2{\pi }_2^3}\end{array}$$ いっぽう、 $$\begin{array}{rl}{\mathrm{RR}}^{2}V(\log \widehat{\mathrm{RR}})& =\frac{{\pi }_1^2}{{\pi }_2^2}(\frac{1-{\pi }_1}{n_1{\pi }_1}+\frac{1-{\pi }_2}{n_2{\pi }_2})\\ & =\frac{{\pi }_1(1-{\pi }_1)}{n_1{\pi }_2^2}+\frac{{\pi }_1^2(1-{\pi }_2)}{n_2{\pi }_2^3}\end{array}$$ したがって、漸近的に $$\begin{array}{r}V\left(\widehat{\mathrm{RR}}\right)\cong {\mathrm{RR}}^{2}V(\log \widehat{\mathrm{RR}})\end{array}$$ $◼$

 

オッズ比

二項分布の正規近似により、標本比率は漸近的に $$\begin{array}{c}{p}_1\sim N[{\pi }_1,\frac{{\pi }_1(1-{\pi }_1)}{n_1}]\\ p_2\sim N[{\pi }_2,\frac{{\pi }_2(1-{\pi }_2)}{n_2}]\end{array}$$ 標本比率ベクトルを $$\begin{array}{r}\mathit{p}=\left(\begin{array}{c}{p}_1\\ p_2\end{array}\right)\end{array}$$ 期待値ベクトルを $$\begin{array}{r}\mathit{\pi }=\left(\begin{array}{c}{\pi }_1\\ {\pi }_2\end{array}\right)\end{array}$$ 分散・共分散行列を $$\begin{array}{r}\mathbf{\Omega }=\left[\begin{array}{cc}\frac{{\pi }_1(1-{\pi }_1)}{n_1}& 0\\ 0& \frac{{\pi }_2(1-{\pi }_2)}{n_2}\end{array}\right]\end{array}$$ として、 $$\begin{array}{c}G\left(\mathit{\pi }\right)=\mathrm{OR}=\frac{{\pi }_1}{1-{\pi }_1}\cdot \frac{1-{\pi }_2}{{\pi }_2}\\ G\left(\mathit{p}\right)=\widehat{\mathrm{OR}}=\frac{p_1}{1-p_1}\cdot \frac{1-p_2}{p_2}\end{array}$$ と変数変換する。 多変量のデルタ法を用いて $G\left(\mathit{p}\right)$ を期待値 $E\left(\mathit{p}\right)=\mathit{\pi }$ まわりでテイラー展開すると、偏導関数ベクトルは、 $$\begin{array}{r}\mathit{H}\left(\mathit{\pi }\right)=\left(\begin{array}{c}\frac{G\left(\mathit{\pi }\right)}{\partial {\pi }_1}\\ \frac{G\left(\mathit{\pi }\right)}{\partial {\pi }_2}\end{array}\right)=\left[\begin{array}{c}\frac{1-{\pi }_2}{{\pi }_2{(1-{\pi }_1)}^2}\\ -\frac{{\pi }_1}{{\pi }_2^2(1-{\pi }_1)}\end{array}\right]\end{array}$$ 多変量のデルタ法の期待値と分散の公式より、 $$\begin{array}{r}E\left[G\left(\mathit{p}\right)\right]\cong G\left(\mathit{\pi }\right)=\frac{{\pi }_1}{1-{\pi }_1}\cdot \frac{1-{\pi }_2}{{\pi }_2}\end{array}$$ $$\begin{array}{rl}V\left[G\left(\mathit{p}\right)\right]& =\left[\begin{array}{cc}\frac{1-{\pi }_2}{{\pi }_2{(1-{\pi }_1)}^2}& -\frac{{\pi }_1}{{\pi }_2^2(1-{\pi }_1)}\end{array}\right]\left[\begin{array}{cc}\frac{{\pi }_1(1-{\pi }_1)}{n_1}& 0\\ 0& \frac{{\pi }_2(1-{\pi }_2)}{n_2}\end{array}\right]\left[\begin{array}{c}\frac{1-{\pi }_2}{{\pi }_2{(1-{\pi }_1)}^2}\\ -\frac{{\pi }_1}{{\pi }_2^2(1-{\pi }_1)}\end{array}\right]\\ & =\left[\begin{array}{cc}\frac{1-{\pi }_2}{{\pi }_2{(1-{\pi }_1)}^2}& -\frac{{\pi }_1}{{\pi }_2^2(1-{\pi }_1)}\end{array}\right]\left[\begin{array}{c}\frac{{\pi }_1(1-{\pi }_2)}{n_1{\pi }_2(1-{\pi }_1)}\\ -\frac{{\pi }_1(1-{\pi }_2)}{n_2{\pi }_2(1-{\pi }_1)}\end{array}\right]\\ & =\frac{1-{\pi }_2}{{\pi }_2{(1-{\pi }_1)}^2}\frac{{\pi }_1(1-{\pi }_2)}{n_1{\pi }_2(1-{\pi }_1)}+\frac{{\pi }_1}{{\pi }_2^2(1-{\pi }_1)}\frac{{\pi }_1(1-{\pi }_2)}{n_2{\pi }_2(1-{\pi }_1)}\\ & =\frac{{\pi }_1{(1-{\pi }_2)}^2}{n_1{\pi }_2^2{(1-{\pi }_1)}^3}+\frac{{\pi }_1^2(1-{\pi }_2)}{n_2{\pi }_2^3{(1-{\pi }_1)}^2}\end{array}$$ いっぽう、 $$\begin{array}{rl}{\mathrm{OR}}^{2}V(\log \widehat{\mathrm{OR}})& =\frac{{\pi }_1^2}{{(1-{\pi }_1)}^2}\cdot \frac{{(1-{\pi }_2)}^2}{{\pi }_2^2}[\frac{1}{n_1{\pi }_1(1-{\pi }_1)}+\frac{1}{n_2{\pi }_2(1-{\pi }_2)}]\\ & =\frac{{\pi }_1{(1-{\pi }_2)}^2}{n_1{\pi }_2^2{(1-{\pi }_1)}^3}+\frac{{\pi }_1^2(1-{\pi }_2)}{n_2{\pi }_2^3{(1-{\pi }_1)}^2}\end{array}$$ したがって、漸近的に $$\begin{array}{r}V\left(\widehat{\mathrm{OR}}\right)\cong {\mathrm{OR}}^{2}V(\log \widehat{\mathrm{OR}})\end{array}$$ $◼$

上述の参考書の問題文（p.81）では、 $$\begin{array}{}\text{(2.164)}& V\left(\widehat{\mathrm{OR}}\right)\cong {\mathrm{OR}}^{2}V(\log \widehat{\mathrm{RR}})\end{array}$$ となっていますが、 これは、おそらく誤植だと思われます。  

問題2.5.5：曝露効果の指標の漸近分布（帰無仮説）

リスク差

${\pi }_1={\pi }_2=\pi$ を代入すると、帰無仮説のもとで、 $$\begin{array}{r}\widehat{\mathrm{RD}}|H_0\stackrel{d}{\to }N[0,\pi (1-\pi )(\frac{1}{n_1}+\frac{1}{n_2})]\end{array}$$ 共通の母比率 $\pi$ の一致推定量は、 $$\begin{array}{r}\hat{\pi }=\frac{a+b}{n_1+n_2}=\frac{m_1}{N}\end{array}$$ これを用いて、母比率を一致推定量である標本比率で置き換えると、漸近分散の一致推定量は、 $$\begin{array}{rl}{\hat{\sigma }}_0^2& =\frac{m_1}{N}\cdot \frac{m_2}{N}\cdot \frac{n_1+n_2}{n_1{n}_2}\\ & =\frac{m_1}{N}\cdot \frac{m_2}{N}\cdot \frac{N}{n_1{n}_2}\\ & =\frac{m_1{m}_2}{N{n}_1{n}_2}\end{array}$$ $◼$

対数リスク比

帰無仮説 $H_0:{\pi }_1={\pi }_2(=\pi )$ のもとでは、 $$\begin{array}{r}\log \widehat{\mathrm{RR}}\stackrel{d}{\to }\mathrm{N}[0,\frac{1-\pi }{\pi }(\frac{1}{n_1}+\frac{1}{n_2})]\end{array}$$ 共通の母比率 $\pi$ の一致推定量は、 $$\begin{array}{r}\hat{\pi }=\frac{a+b}{n_1+n_2}=\frac{m_1}{N}\end{array}$$ これを用いて、母比率を一致推定量である標本比率で置き換えると、漸近分散の一致推定量は、 $$\begin{array}{rl}{\hat{\sigma }}_0^2& =\frac{1}{\hat{\pi }(1-\hat{\pi })}(\frac{1}{n_1}+\frac{1}{n_2})\\ & =\frac{N}{m_1}\cdot \frac{N}{m_2}\cdot \frac{N}{n_1{n}_2}\\ & =\frac{N^3}{m_1{m}_2{n}_1{n}_2}\end{array}$$ $◼$

対数オッズ比

帰無仮説 $H_0:{\pi }_1={\pi }_2(=\pi )$ のもとでは、 $$\begin{array}{r}\log \widehat{\mathrm{OR}}\stackrel{d}{\to }\mathrm{N}[0,\frac{1}{\pi (1-\pi )}(\frac{1}{n_1}+\frac{1}{n_2})]\end{array}$$ 共通の母比率 $\pi$ の一致推定量は、 $$\begin{array}{r}\hat{\pi }=\frac{a+b}{n_1+n_2}=\frac{m_1}{N}\end{array}$$ これを用いて、母比率を一致推定量である標本比率で置き換えると、漸近分散の一致推定量は、 $$\begin{array}{rl}{\hat{\sigma }}_0^2& =\frac{1-\hat{\pi }}{\hat{\pi }}\cdot \frac{N}{n_1{n}_2}\\ & =\frac{N}{m_1}\cdot \frac{m_2}{N}\cdot \frac{N}{n_1{n}_2}\\ & =\frac{m_2}{m_1}\cdot \frac{N}{n_1{n}_2}\end{array}$$ $◼$

 

問題2.5.6：リスク比の定義による影響（分散）

対数の性質より、 $$\begin{array}{c}\log {\widehat{\mathrm{RR}}}_{2:1}=\log p_2-\log p_1\\ \log {\widehat{\mathrm{RR}}}_{1:2}=\log p_1-\log p_2\end{array}$$ したがって、 $$\begin{array}{r}\log {\widehat{\mathrm{RR}}}_{2:1}=-\log {\widehat{\mathrm{RR}}}_{1:2}\end{array}$$ しかし、分散の性質より、 $$\begin{array}{c}V(\log {\widehat{\mathrm{RR}}}_{2:1})=V(\log p_2)+V(\log p_1)\\ V(\log {\widehat{\mathrm{RR}}}_{1:2})=V(\log p_1)+V(\log p_2)\end{array}$$ したがって、 $$\begin{array}{r}V(\log {\widehat{\mathrm{RR}}}_{2:1})=V(\log {\widehat{\mathrm{RR}}}_{1:2})\end{array}$$ $◼$

問題2.5.7：リスク比の定義による影響（信頼区間）

ここで、 $$\begin{array}{c}\hat{\theta }=\log {\widehat{\mathrm{RR}}}_{2:1}=-\log {\widehat{\mathrm{RR}}}_{1:2}\\ {\hat{\sigma }}^2=V(\log {\widehat{\mathrm{RR}}}_{2:1})=V(\log {\widehat{\mathrm{RR}}}_{1:2})\end{array}$$ とすると、 対数リスク比の信頼区間は、${\hat{\theta }}_L\le \log RR\le {\hat{\theta }}_U$ なので、
（i）$\log {\widehat{\mathrm{RR}}}_{2:1}$ の上下信頼限界 $$\begin{array}{rl}{\hat{\theta }}_{L21}& =\log {\widehat{\mathrm{RR}}}_{2:1}-Z_{0.5\alpha }\cdot \hat{\sigma }\\ & =\hat{\theta }-Z_{0.5\alpha }\cdot \hat{\sigma }\end{array}$$ $$\begin{array}{r}{\hat{\theta }}_{U21}=\log {\widehat{\mathrm{RR}}}_{2:1}+Z_{0.5\alpha }\cdot \hat{\sigma }\\ {\hat{\theta }}_{U21}=\hat{\theta }+Z_{0.5\alpha }\cdot \hat{\sigma }\end{array}$$ （ii）$\log {\widehat{\mathrm{RR}}}_{1:2}$ の上下信頼限界 $$\begin{array}{rl}{\hat{\theta }}_{L12}& =\log {\widehat{\mathrm{RR}}}_{1:2}-Z_{0.5\alpha }\cdot \hat{\sigma }\\ & =-\hat{\theta }-Z_{0.5\alpha }\cdot \hat{\sigma }\\ & =-(\hat{\theta }+Z_{0.5\alpha }\cdot \hat{\sigma })\\ & =-{\hat{\theta }}_{U21}\end{array}$$ $$\begin{array}{rl}{\hat{\theta }}_{U12}& =\log {\widehat{\mathrm{RR}}}_{1:2}+Z_{1-0.5\alpha }\cdot \hat{\sigma }\\ & =-\hat{\theta }+Z_{1-0.5\alpha }\cdot \hat{\sigma }\\ & =-(\hat{\theta }-Z_{0.5\alpha }\cdot \hat{\sigma })\\ & =-{\hat{\theta }}_{L21}\end{array}$$ また、$\exp \left({\hat{\theta }}_L\right)\le \mathrm{RR}\le \exp \left({\hat{\theta }}_U\right)$ より、 $$\begin{array}{c}\exp \left({\hat{\theta }}_{L21}\right)=\exp (-{\hat{\theta }}_{U12})=\frac{1}{\exp \left({\hat{\theta }}_{U12}\right)}\\ \exp \left({\hat{\theta }}_{U21}\right)=\exp (-{\hat{\theta }}_{L12})=\frac{1}{\exp \left({\hat{\theta }}_{L12}\right)}\end{array}$$ $◼$

 

参考文献

• ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.80-82
• ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.23-28
• Woolf, B.. On estimating the relation between blood group and disease. Annals of human genetics. 1955;19(4):251-253, doi: https://doi.org/10.1111/j.1469-1809.1955.tb01348.x

関連記事

• コホート研究
• ケース・コントロール研究