- 연구방법: 2006-2017년 12개년도의 전국 250개 시군구 자료를 이용하였다. 연도별 각 시군구 쯔쯔가무시증 발생건수는 질병관리본부 감염병 감시자료를 활용하였으며 산림파괴수준은 위성자료를 전처리한 자료(Landsat)를 이용하였다. 포아송 베이지안 회귀모형으로 산림파괴수준과 쯔쯔가무시증의 연관성을 평가하였고, 과대산포(over-dispersion)와 영과잉(excessive zero)을 고려하여 negative binomial, zero-inflated Poisson, and zero-inflated negative binomial 모형들도 비교하였다. 시공간자기상관성을 고려하기 위해 위계모형을 활용하였고, deviance information criterion 값으로 최종 모형을 선택하였다. 모형에서 인구사회학적 요인, 경제학적 요인, 지리적 요인, 기상학적 요인을 보정하였다.

- 결과: 최종모형에서 산림파괴수준의 사분위범위(interquartile range) 증가에 따른 쯔쯔가무시증 발생의 상대위험도는 1.20였으며, 95% 신용구간(credible interval)은 1.15-1.24로 나타났다.

- 결론: 산림파괴는 쯔쯔가무시증 발생의 위험도 증가와 연관성이 있는 것으로 관찰되었다.