本稿は、ジョン・ラチン(2020)『医薬データのための統計解析』の「問題8.7」の一部に対する自作解答例です。ポアソン回帰分析に必要なサンプルサイズと検出力に関する問題です。
なお、閲覧にあたっては、以下の点にご注意ください。
- スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。
- 曝露(発症)状況を表す右下の添え字は、「0」である場合($n_0,\pi_0$ など)や「2」である場合($n_2,\pi_2$ など)がありますが、どちらも「非曝露群(コントロール群)」を表しています。
- 漸近的な性質を用いる際は、①中心極限定理が成り立つ、②漸近分散を推定する際に、母数をその一致推定量で置き換えることができるということが成り立つと仮定しています。
- 著作権の関係上、問題文は、掲載しておりません。上述の参考書をお持ちの方は、お手元にご用意してご覧ください。
- この解答例は、筆者が自作したものであり、公式なものではありません。あくまでも参考としてご覧いただければ幸いです。
問題8.7.0:効果量・有意水準・検出力の関係式(ポアソンモデル)
2重同次ポアソンモデルの仮定より、
\begin{align}
E \left(d_{ij}\right)=V \left(d_{ij}\right)=\lambda_it_{ij}
\end{align}
期待値の性質より、
\begin{align}
E \left(d_{i\bullet }\right)&=E \left(\sum_{j=1}^{n_i}d_{ij}\right)\\
&=E \left(\sum_{l=1}^{n_i}{\lambda_it_{ij}}\right)\\
&=\lambda_iE \left(\sum_{l=1}^{n_i}t_{ij}\right)\\
&=\lambda_i \cdot t_{i\bullet }
\end{align}
ポアソン分布の正規近似により、漸近的に、
\begin{align}
d_{i\bullet } \sim \mathrm{N} \left(\lambda_i \cdot t_{i\bullet },\lambda_i \cdot t_{i\bullet }\right)
\end{align}
パラメータ $\lambda_i$ の最尤推定量 ${\hat{\lambda}}_i=\frac{d_{i\bullet }}{t_{i\bullet }}$ について、期待値と分散の性質より
\begin{align}
E \left({\hat{\lambda}}_i\right)&=E \left(\frac{d_{i\bullet }}{t_{i\bullet }}\right)\\
&=\frac{\lambda_i \cdot t_{i\bullet }}{t_{i\bullet }}\\
&=\lambda_i\\
V \left({\hat{\lambda}}_i\right)&=V \left(\frac{d_{i\bullet }}{t_{i\bullet }}\right)\\
&=\frac{\lambda_i \cdot t_{i\bullet }}{t_{i\bullet }^2}\\
&=\frac{\lambda_i}{t_{i\bullet }}
\end{align}
したがって、漸近的に、
\begin{align}
{\hat{\lambda}}_1 \sim \mathrm{N} \left(\lambda_1,\frac{\lambda_1}{t_{1\bullet }}\right) \quad {\hat{\lambda}}_2 \sim \mathrm{N} \left(\lambda_2,\frac{\lambda_2}{t_{2\bullet }}\right)
\end{align}
最尤推定量の差を $\hat{\delta}={\hat{\lambda}}_1-{\hat{\lambda}}_2$ とすると、正規分布の再生性より、
\begin{align}
{\hat{\delta}}_1 \sim \mathrm{N} \left(\lambda_1-\lambda_2,\frac{\lambda_1}{t_{1\bullet }}+\frac{\lambda_2}{t_{2\bullet }}\right)
\end{align}
帰無仮説 $H_0:\lambda_1=\lambda_2=\lambda$ のもとでは、
\begin{align}
{\hat{\delta}}_0 \sim \mathrm{N} \left[0,\lambda \left(\frac{t_{1\bullet }+t_{2\bullet }}{t_{1\bullet }t_{2\bullet }}\right)\right]
\end{align}
臨床的有意差・有意水準・検出力の関係式より、
\begin{gather}
\left|\delta_1-\delta_0\right|=Z_\alpha \cdot \sigma_0-Z_{1-\beta} \cdot \sigma_1\\
\left|\lambda_1-\lambda_2\right|=Z_\alpha \cdot \sqrt{\lambda \left(\frac{t_{1\bullet }+t_{2\bullet }}{t_{1\bullet }t_{2\bullet }}\right)}-Z_{1-\beta} \cdot \sqrt{\frac{\lambda_1}{t_{1\bullet }}+\frac{\lambda_2}{t_{2\bullet }}}
\end{gather}
したがって、2つのグループ間の率の差に対する $Z$ 検定の検出力は、曝露の合計患者年数の関数である。
また、$E \left(d_{i\bullet }\right)=\lambda_i \cdot t_{i\bullet }$ より、
\begin{align}
\sigma_0^2&=\frac{\lambda}{t_{1\bullet }}+\frac{\lambda}{t_{2\bullet }}\\
&=\frac{\lambda^2}{\lambda t_{1\bullet }}+\frac{\lambda^2}{\lambda t_{2\bullet }}\\
&=\frac{\lambda^2}{E \left(d_{1\bullet }\right)}+\frac{\lambda^2}{E \left(d_{2\bullet }\right)}\\
&=\lambda^2 \left\{\frac{E \left(d_{1\bullet }\right)+E \left(d_{2\bullet }\right)}{E \left(d_{1\bullet }\right)E \left(d_{2\bullet }\right)}\right\}
\end{align}
\begin{align}
\sigma_1^2&=\frac{\lambda_1}{t_{1\bullet }}+\frac{\lambda_2}{t_{2\bullet }}\\
&=\frac{\lambda_1^2}{\lambda_1t_{1\bullet }}+\frac{\lambda_2^2}{\lambda_2t_{2\bullet }}\\
&=\frac{\lambda_1^2}{E \left(d_{1\bullet }\right)}+\frac{\lambda_2^2}{E \left(d_{2\bullet }\right)}
\end{align}
したがって、
\begin{align}
\left|\lambda_1-\lambda_2\right|=Z_\alpha \cdot \sqrt{\lambda^2 \left\{\frac{E \left(d_{1\bullet }\right)+E \left(d_{2\bullet }\right)}{E \left(d_{1\bullet }\right)E \left(d_{2\bullet }\right)}\right\}}-Z_{1-\beta} \cdot \sqrt{\frac{\lambda_1^2}{E \left(d_{1\bullet }\right)}+\frac{\lambda_2^2}{E \left(d_{2\bullet }\right)}}
\end{align}
すなわち、2つのグループ間の率の差に対する $Z$ 検定の検出力は、各グループの期待事象数の関数でもある。
$\blacksquare$
参考文献
- ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.455-456
- ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.446-448
0 件のコメント:
コメントを投稿