[論文レビュー] Parameter estimation of beta-geometric model with application to human fecundability data

Brijesh Singh, P. S. Pudir+1 more|arXiv (Cornell University)|May 25, 2014

Statistical Methods and Bayesian Inference参考文献 3被引用数 3

ひとこと要約

本稿では、第3回国際家族計画調査（NFHS-III）のデータを用いて、人間の不妊能耐性を推定するためのベータ幾何分布モデルを提案する。古典的（モーメント法および最尤法）およびベイズ推定法を適用し、モデルがカップル間での不妊能耐性の不均一性を効果的に捉えていることを示す。シミュレーション結果から、実世界の出生前健康データに対する推定手法の性能が安定していることが明らかになった。

ABSTRACT

The present study deals with the estimation of the mean value of fecundability by fitting a theoretical distribution from the observed month of first conception of the married women who did not use any contraceptive method before their first conception. It is assumed that fecundability is fixed for a given couple, but across couples it varies according to a specified distribution. Under the classical approach, methods of moment and maximum likelihood are used while for Bayesian approach we use the above two estimates as prior for fecundability parameter. A real data analysis from the third National Family Health Survey (NFHS-III) is analyzed as an application of model. Finally, a simulation study is performed to access the performance of the several of methods used in this paper

研究の動機と目的

• 統計的分布を用いてカップル間での不妊能耐性のばらつきをモデル化すること。
• 避妊を行わない女性の初回受胎月の観察データを用いて、平均不妊能耐性を推定すること。
• 導出された事前分布を用いたベイズ的手法と、古典的推定法（モーメント法および最尤法）を比較すること。
• シミュレーションおよびNFHS-IIIの実データ解析を通じて、モデルの性能を評価すること。
• 集団レベルの出生前健康研究における不妊結果の評価に向けた堅牢な統計的フレームワークを提供すること。

提案手法

• 不妊能耐性を、カップルごとに異なるベータ分布に従う潜在変数としてモデル化する。
• 受胎までの月数をベータ幾何分布でモデル化し、成功確率がベータ分布に従う幾何分布待ち時間の仮定を採用する。
• 古典的手法として、モーメント法および最尤推定法（MLE）を用いてモデルのパラメータを推定する。
• ベイズ推論においては、古典的推定値を不妊能耐性パラメータの事前分布として用いる。
• 推定手法の効率性と正確性を比較するためのシミュレーションスタディを実施する。
• NFHS-IIIの実データを解析し、モデルが人間の不妊データへの適用可能性を検証する。

実験結果

リサーチクエスチョン

• RQ1 ベータ幾何分布モデルは、避妊を行わないカップルの集団における初回受胎までの期間の分布をどれほど適切に捉えているか？
• RQ2 モーメント法と最尤法のうち、どちらの推定手法が不妊能耐性のパラメータ推定においてより正確か？
• RQ3 古典的推定値を事前分布として用いるベイズ的手法は、古典的手法に比べてパラメータ推定をどのように改善するか？
• RQ4 第3回国際家族計画調査（NFHS-III）の実世界データに適用した際、モデルの実効的性能はいかがなものか？
• RQ5 さまざまなシミュレーション条件下で、バイアス、分散、平均二乗誤差の観点から、推定手法はどのように比較されるか？

主な発見

• ベータ幾何分布モデルは、カップル間での不妊能耐性の未観測の不均一性を効果的に説明できる。
• シミュレーションスタディにおいて、最尤推定法がモーメント法よりもより精度の高いパラメータ推定を達成した。
• 古典的推定値を事前分布として用いるベイズ的手法は、推定の正確性を向上させるとともに分散を低減した。
• モデルはNFHS-IIIデータに良好に適合し、推定された平均不妊能耐性は既知の出生率パターンと整合的であった。
• シミュレーション結果から、特に最尤法およびベイズ的手法が、さまざまなサンプルサイズ下でも安定した性能を示した。
• 本研究は、不妊能耐性に分布的仮定を組み込むことで、集団調査における不妊推定の信頼性が向上することを示した。