Improving the reliability of district-level estimates ofteenage pregnancy in Malawi using small area estimation

技術要約：小領域推定法を用いたマラウイにおける地区レベルのティーンエイジ妊娠率推計の信頼性向上

問題提起
マラウイは、意思決定の分権化へと移行しており、地区レベルの計画策定や資源配分において、信頼性の高いローカルな健康指標を必要としている。人口保健調査（DHS）などの家庭訪問調査を通じて、全国および地域レベルのティーンエイジ妊娠率の推計値は得られているが、地区レベルの直接的な推計値はしばしば信頼性に欠ける。この信頼性の低さは、詳細な行政単位に対して設計された調査規模では不十分な、サンプルサイズの小ささとサンプリングの変動性の高さに起因している。その結果、直接的な推計値は大きな変動係数（CV）を示し、エビデンスに基づく政策立案や標的を絞った介入の妨げとなっている。本研究は、調査に基づく指標との整合性を失うことなく、これらの推計値を安定させるための統計的枠組みの必要性に対処するものである。

手法
本研究では、マラウイの32の行政地区におけるティーンエイジ妊娠率の推計精度の向上を図るため、エリアレベルのフェイ＝ヘリオット（FH）モデルを用いたモデルベースの小領域推定（SAE）手法を採用している。

• データソース： 主要データは、2024年マラウイ人口保健調査（MDHS）から取得し、15歳から19歳の女性に焦点を当てている。ティーンエイジ妊娠は、この年齢層の女性のうち、現在妊娠中であるか、または妊娠経験（流産を含む）がある割合と定義されている。補助変数は、2018年マラウイ人口住宅センサスおよび統合家計調査（IHS5）から導出された。
• 変数選択： 18個の地区レベルの補助変数の初期プールを評価した。10分割交差検証を用いた最小絶対収縮選択演算子（LASSO）回帰を使用し、最終的な連結モデルのための4つの主要な予測変数を選定した。
  1. 農村部に居住する女性の割合。
  2. 18歳未満で結婚した女性の割合。
  3. 女性世帯主の割合。
  4. 対数変換された人口密度。
• 統計モデル：
  • 直接推計： 設計ベースの直接推計値（$\hat{p}^{DIR}_i$）およびそのサンプリング分散（$\psi_i$）を、調査ウェイトとテイラー展開を用いて算出した。
  • 変換： FHモデルの線形性と正規性の仮定を満たすため、直接推計値にロジット変換を適用し、境界値（0または1）に対しては連続性補正（0.5）を行った。
  • フェイ＝ヘリオット（FH）モデル： 以下の2段階モデルを構築した。
    • サンプリングモデル： $\hat{p}_i = p_i + e_i$, ここで $e_i \sim N(0, \psi_i)$。
    • 連結モデル： $p_i = \mathbf{x}^\top_i \boldsymbol{\beta} + u_i$, ここで $u_i \sim N(0, \sigma^2_u)$ はエリア固有のランダム効果を表す。
  • 推定： 分散成分を推定するために、制限最大尤度法（REML）を用いてモデルを適合させた。得られた経験的最良線形不偏予測量（EBLUP）は、直接推計値と回帰ベースの合成推計値の加重合成として機能する。
  • 逆変換： 元の確率スケールに戻すために、バイアス補正を伴う逆変換を適用した。
  • 精度評価： 平均二乗誤差（MSE）および変動係数（CV）を、200回の反復によるパラメトリック・ブートストラップ法を用いて推定した。

主な結果

• モデルの性能： 連結モデルは、エリア間の変動の大部分を説明しており、調整済み$R^2$は0.666、FH特有の$R^2$は0.878であった。ランダム効果の推定分散は、ロジットスケールで0.00785であった。
• 安定性と精度： モデルベースのEBLUP推計値は、直接的な調査推計値と比較して変動が減少した。
  • 直接推計値： 平均 = 0.3070、標準偏差 = 0.0752、範囲 = 0.130–0.421。
  • EBLUP推計値： 平均 = 0.3066、標準偏差 = 0.0628、範囲 = 0.151–0.385。
  • EBLUPの四分位範囲はより狭くなっており、これはより高い安定性を示している。
• 誤差の減少： EBLUP推計値のMSEは、全32地区において直接推計値よりも低かった。変動係数（CV）もすべての地区で一貫してモデルベースの推計値の方が低く、効率性の向上が確認された。
• 診断： 残差分析（Q-Qプロット、シャピロ＝ウィルク検定）により、正規性と独立性の仮定が満たされていることが確認された。直接推計値とEBLUP推計値の間には強い正の相関（$r=0.81$）が観察され、モデルベースのアプローチがノイズを低減しつつ、調査データとの整合性を維持していることが示された。
• 具体的な調整： モデルは極端な直接推計値を緩和した。例えば、チカワ地区の直接推計値は42.1%から33.5%に低下し、一方でゾンバ地区の推計値は26.3%から32.0%に上昇した。これは、補助変数および近隣エリアからの「情報の共有（borrowing of strength）」を反映している。

意義と主張
本論文は、エリアレベルのフェイ＝ヘリオット・モデルの適用が、地区レベルの調査における小規模サンプルの限界を効果的に解決していると主張している。調査データと補助的なセンサス情報を統合することで、本研究はより安定し信頼性の高いティーンエイジ妊娠率の推計を実現している。

著者らは、これらの改善された推計値が、マラウイの分権化された保健システムのためのより強固なエビデンス基盤を提供すると断言している。具体的には、本研究の結果は、「国家青少年フレンドリー・ヘルスサービス戦略」や「児童婚終結に関する国家戦略」といった国家政策の計画、モニタリング、および評価を支援するものである。本研究は、SAE手法が、直接的な小規模下限推計には不十分なサンプルサイズを持つ低所得国において特に価値があり、より一貫した標的設定とリプロダクティブ・ヘルス介入の配分を可能にすると結論付けている。また、モデルは精度を向上させる一方で、その正確性は補助変数の質に依存することに注意を促しており、今後の研究ではこの手法を多変量時空間モデルへと拡張できる可能性があるとしている。