[論文レビュー] Two-stage approaches to the analysis of occupancy data I: The homogeneous case

Natalie Karavarsamis, Richard Huggins|arXiv (Cornell University)|Mar 30, 2018

Spatial and Panel Data Analysis参考文献 10被引用数 4

ひとこと要約

本稿では、直交パラメータ化と部分尤度を用いた二段階推定手法を提案し、同質的占有モデルの最尤推定を簡素化する。検出確率と占有確率の間で直交性を確保するためのパラメータ変換により、条件付き尤度と部分尤度を別々に用いた推定が可能となり、計算負荷が最小限で、かつ完全尤度法に比べてわずかな効率の損失しか生じない、効率的で解析的に扱いやすい推定量が得られる。

ABSTRACT

Occupancy models are used in statistical ecology to estimate species dispersion. The two components of an occupancy model are the detection and occupancy probabilities, with the main interest being in the occupancy probabilities. We show that for the homogeneous occupancy model there is an orthogonal transformation of the parameters that gives a natural two-stage inference procedure based on a conditional likelihood. We then extend this to a partial likelihood that gives explicit estimators of the model parameters. By allowing the separate modelling of the detection and occupancy probabilities, the extension of the two-stage approach to more general models has the potential to simplify the computational routines used there.

研究の動機と目的

• 検出確率と占有確率の同時モデル化を伴う、複雑な完全尤度推定を簡素化すること。
• 多くの共変数を含むモデルでは完全尤度が、可能なモデル数の指数的増加により計算不能になるため、その計算負荷を軽減すること。
• 検出確率と占有確率の別々のモデル化を可能にする二段階推定手順を開発し、計算の安定性と効率性を向上させること。
• 直交パラメータ化と部分尤度の活用により、占有確率と検出確率の明示的で解析的に取り扱いやすい推定量を導出すること。

提案手法

• 検出確率と占有確率の間の直交性を確保するため、パラメータの直交変換を適用し、η = ψθ（占有済みの場所で検出される確率）を定義し、p（訪問回数ごとの検出確率）を保持する。これにより、漸近的に独立な最尤推定量が得られる。
• 占有済みの場所における検出回数に基づく条件付き尤度を用い、検出確率の二項尤度が得られる。
• ηとpの二つの独立な二項尤度に分解される部分尤度を導出し、明示的な解析的推定量が得られる。
• 二段階アプローチを実装する：まず条件付き尤度によりpを推定し、次に変換されたパラメータηを用いてψを推定する。これにより、同時最大化を回避できる。
• シミュレーションと実データ（カエルの検出データ）への応用を通じて、完全尤度との比較により、効率性とバイアスを検証する。
• RパッケージVGAMを用いて条件付き尤度アプローチを実装し、完全尤度推定との結果を比較する。

実験結果

リサーチクエスチョン

• RQ1 同質的占有モデルにおける直交パラメータ化は、最尤推定を簡素化する自然な二段階推定手順をもたらすか？
• RQ2 部分尤度アプローチは完全尤度と比較して、効率性と計算可能性の点でどのように異なるか？
• RQ3 小さな検出確率と少ない調査回数が、推定量のバイアスと標準誤差の正確性に与える影響は何か？
• RQ4 多くの共変数を含むモデルにおいて、二段階法は計算複雑性を軽減しつつも高い効率性を維持できるか？

主な発見

• η = ψθとpのパラメータの直交変換により、漸近的に独立な最尤推定量が得られ、推定が簡素化され、計算複雑性が低減する。
• 条件付き尤度から得られる検出確率推定量pは、標準的な二項推定量と等価であり、実装が容易である。
• 部分尤度アプローチにより、ψとpの両方の明示的解析的推定量が得られ、ψの推定量は完全尤度と比較して90％以上の相対効率を達成する。
• シミュレーションでは両手法に顕著なバイアスは認められず、推定された標準誤差も信頼性があるが、pとτが小さい場合にはわずかに過大評価が生じる。
• 27の調査地点（τ=4）における実データ応用では、部分尤度によりp̃=0.889、ψ̃=0.556が得られた一方、完全尤度ではp̂=0.780、ψ̂=0.557が得られ、検出確率の推定値に顕著な差が生じた。
• pの推定値の差異は、部分尤度が完全尤度で用いられる初回訪問の情報を無視するため、調査回数に応じた検出確率の非定常性（異質性）の可能性を示唆している。