희귀 수생조류 종의 존재와 풍부성을 희소 데이터에서 추정하는 새로운 통계 모델

초록

본 논문은 관측이 결측되고 영(0)값이 과다한 희귀 수생조류 데이터에 대해, 저차원 잠재 가우시안 구조와 영-인플레이션을 결합한 Poisson‑Log‑Normal 모델(ZI‑PLN)을 제안한다. 변분 EM 알고리즘으로 파라미터를 추정하고, 결측값의 예측 구간, 연도별 추세 및 변곡점 검출 절차를 제공한다. 시뮬레이션과 국제 수생조류 조사(IWC) 데이터를 통해 모델의 정확성과 실용성을 검증한다.

상세 분석

이 연구는 희귀 조류와 같이 관측이 드물고 영값이 다수인 다변량 카운트 데이터를 다루기 위해 기존 Poisson‑Log‑Normal(PLN) 모델에 영‑인플레이션(Bernoulli)과 저차원 잠재 구조를 동시에 도입한 ZI‑PLN‑PCA 모델을 설계하였다. 먼저, 각 사이트‑연도(i,j)에서 종의 존재 여부 Uᵢⱼ를 로짓 선형 모델(logit πᵢⱼ = xᵢⱼᵀγ)로 정의하고, 존재할 경우에만 Poisson(λᵢⱼ) 분포를 따르게 하여 영‑인플레이션을 구현한다. λᵢⱼ는 고정 효과 xᵢⱼᵀβ와 잠재 가우시안 Zᵢⱼ의 지수합으로 표현되며, Zᵢ = C Wᵢ (Wᵢ∼N(0,I_q)) 형태의 저차원 구조를 통해 연도 간 상관관계를 효율적으로 추정한다. 여기서 C는 p × q 차원의 로딩 행렬이며, Σ = C Cᵀ 로 정의되는 공분산 행렬은 저계수(rank ≤ q) 제약을 받아 과적합을 방지한다.

모델 식별성은 두 가지 주요 가정에 기반한다. 첫째, 모든 연도 쌍 (j,k)에 대해 최소 하나의 사이트가 두 연도 모두 관측된 집합 Q가 전체 연도 쌍을 포괄한다(