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

公開日: 更新日:

【2022年10月1週】 【A000】生物統計学 【A072】統計的推定

この記事をシェアする
  • B!
サムネイル画像

本稿は、ジョン・ラチン(2020)『医薬データのための統計解析』の「問題2.2」の自作解答例です。検定の逆変換にもとづく母比率の信頼区間に関する問題です。

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

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

問題2.2

二項分布の正規近似により、標本成功確率は漸近的に、 pdN[π,π(1π)n] 標本成功確率を標準化すると、 n(pπ)π(1π)=ZN(0,1) このとき、母数 π100(1α)% 信頼区間は、 Z0.5αn(pπ)π(1π)Z0.5α このとき、この不等式を π について解いたときの値が上下の信頼限界となる。
略記のために Z0.5α=z とおいて両辺を2乗すると、 z2=n(pπ)2π(1π)π22pπ+p2(ππ2)nz2=0(z2n+1)π2(z2n+2p)π+p2=0 2次方程式の解の公式 ax2+bx+c=0x=b±b24ac2a より、 π=(z2n+2p)±(z2n+2p)24(z2n+1)p22(z2n+1)=(z2n+2p)±4p2+4pz2n+z4n24p24p2z2n2(z2n+1)=(z2n+2p)±4pz2n+z4n24p2z2n2(z2n+1)=(z2n+2p)±z2n(4p+z2n4p2)2(z2n+1)=(z2n+2p)±z2n{z2n+4p(1p)}2(z2n+1)=(z22n+p)±z24n{z2n+4p(1p)}z2n+1 したがって、100(1α)% 信頼区間は、 (Z0.5α22n+p)Z0.5α24n[Z0.5α2n+4p(1p)]Z0.5α2n+1π(Z0.5α22n+p)+Z0.5α24n[Z0.5α2n+4p(1p)]Z0.5α2n+1

通常、信頼区間を求める場合には、標準誤差 π(1π)n の母数 π をその一致推定量である p で置き換えるが、この方法の場合は、π のまま計算し、π に関する2次方程式の解として、上下の信頼限界を直接求める点に特徴がある。

サンプルサイズとイベント数が問題文で与えられた条件のとき、 p=246=0.04 π に対する信頼限界は、 θ^L=(1.962246+0.04)1.962446[1.96246+40.04(10.04)]1.96246+1=0.012005θ^U=(1.962246+0.04)+1.962446[1.96246+40.04(10.04)]1.96246+1=0.14533

参考文献

  • ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.18-19, p.80
  • Wilson, E.B.. Probable Inference, the Law of Succession, and Statistical Inference. Journal of the American Statistical Association. 1927;22(158):209-212, doi: https://doi.org/10.2307/2276774

関連記事

自己紹介

自分の写真

yama

大学時代に読書の面白さに気づいて以来、読書や勉強を通じて、興味をもったことや新しいことを学ぶことが生きる原動力。そんな人間が、その時々に学んだことを備忘録兼人生の軌跡として記録しているブログです。

このブログを検索

ブログ アーカイブ

QooQ