[논문 리뷰] Two-stage approaches to the analysis of occupancy data I: The homogeneous case

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

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에 대해 각각 독립적인 이항 우도로 분해되는 부분우도를 유도함으로써 명시적 해석 가능한 추정량을 가능하게 함.
• 이단계 접근법을 구현: 먼저 조건부우도를 통해 p를 추정하고, 이후 변환된 매개변수 η를 통해 ψ를 추정함으로써 동시 최대화를 피함.
• 시뮬레이션과 실제 frog 검출 데이터 적용을 통해 전체 우도와의 효율성 및 편향을 비교함.
• R 패키지 VGAM을 사용해 조건부우도 접근법을 구현하고 전체 우도 추정과 결과를 비교함.

실험 결과

연구 질문

• RQ1 동질적 점유 모델에서 직교 매개변수화가 최대우도 추정을 단순화하는 자연스러운 이단계 추론 절차를 가능하게 할 수 있는가?
• RQ2 부분우도 접근법은 전체 우도 대비 효율성과 계산 가능성 측면에서 어떻게 비교되는가?
• RQ3 작은 검출 확률과 적은 수의 조사 기회가 추정량의 편향과 표준오차 정확도에 어떤 영향을 미치는가?
• RQ4 많은 공변량을 포함한 모델에서 이단계 방법이 계산 복잡성을 감소시키면서도 높은 효율성을 유지할 수 있는가?

주요 결과

• η = ψθ와 p의 직교 매개변수화는 점차적 독립적인 최대우도 추정량을 도출하여 추론을 단순화하고 계산 복잡성을 감소시킴.
• 조건부우도에서 유도된 검출 확률 추정량 p는 표준 이항 추정량과 동일하여 간편한 구현이 가능함.
• 부분우도 접근법은 ψ와 p에 대해 명시적 해석 가능한 추정량을 도출하며, ψ의 추정량은 전체 우도 대비 90% 이상의 상대적 효율성을 달성함.
• 시뮬레이션 결과 두 방법 모두 유의미한 편향이 없었고, 추정된 표준오차는 신뢰할 수 있었지만, 작은 p와 τ 값에서는 약간의 과대추정이 발생함.
• 27개의 서식지와 τ=4인 실제 데이터 적용에서 부분우도는 p̃=0.889와 ψ̃=0.556을 도출하였고, 전체 우도는 p̂=0.780과 ψ̂=0.557을 도출함으로써 검출 확률 추정치에 뚜렷한 차이가 있음.
• p 추정치의 격차는 부분우도가 전체 우도에서 사용된 첫 번째 방문 정보를 무시하기 때문에, 조사 기회 간 검출 확률의 이질성이 존재할 가능성을 시사함.