공간 인과 텐서 완성: PFAS 오염과 다중 노출·다중 결과 분석

본 논문은 퍼- 및 폴리플루오로알킬 물질(PFAS)이라는 복합 화학물질군이 여러 건강 결과에 미치는 인과 효과를 정확히 추정하기 위한 새로운 통계적 프레임워크를 제시한다. PFAS는 물, 기름, 열에 대한 저항성 때문에 1950년대 이후 광범위하게 사용되어 왔으며, ‘영원한 화학물질’이라는 별명을 가질 정도로 환경에 장기 잔류한다. 실제 환경에서는 PFOA와 PFOS와 같이 여러 PFAS가 동시에 존재하며, 각각의 화학적 특성과 노출 경로가 다르기 때문에 개별 효과를 분리하기가 어렵다. 또한, 노출과 건강 결과 모두 지리적 패턴을 보이며, 측정되지 않은 공간적 교란 요인(예: 미측정 산업 활동, 지역적 환경 요인)이 존재할 가능성이 높다. 이러한 두 가지 난제—다중 노출·다중 결과와 공간 교란—를 동시에 해결하는 것이 본 연구의 핵심 목표이다. 기존 연구는 (1) 혼합 노출 분석을 위해 가중량합 회귀, 베이지안 커널 머신, 양자화 기반 g-컴퓨테이션 등을 사용했지만, 대부분 단일 결과에만 적용되며 공간 교란을 고려하지 못한다. (2) 공간 회귀 방법은 하나의 노출과 하나의 결과에만 적용 가능해 복합 오염 물질 연구에 부적합하다. (3) 최근 인과 텐서 완성 접근법은 잠재 결과를 저차원 구조에 맞춰 보간하지만, 노출 간 상호작용을 무시하거나 공간 정보를 활용하지 않는다. 이에 저자들은 저차원 Tucker 텐서를 기반으로 한 “공간 인과 텐서 완성” 모델을 제안한다. 관측된 데이터는 N × L × O 형태의 3차 텐서 Y로 구성한다(N: 지역 단위, L: 노출 조합(2^K), O: 결과 수). 각 지역 i는 K 차원의 이진 노출 벡터 A_i를 가지고, L은 가능한 노출 조합 수이다. 잠재 결과 텐서는 Y* = G ×₁ U₁ ×₂ U₂ ×₃ U₃ 로 표현되며, G는 핵심 텐서, U₁은 공간 요인, U₂는 노출 요인, U₃는 결과 요인을 나타낸다. 공간 요인 U₁은 측정된 공변량 Z와 미측정 공간 요인 S의 선형 결합(Z η_Z + S η_S)으로 모델링된다. 여기서 S는 그래프 라플라시안 고유벡터(스펙트럴 베이스)로 근사되며, 이는 인접 지역 간 교란이 유사하다는 평활성 가정을 수학적으로 구현한다. 추정 절차는 세 단계로 진행된다. 첫째, 초기 Tucker 분해를 통해 G와 U₁, U₂, U₃를 얻는다. 둘째, 라플라시안 고유벡터를 이용해 S와 η_S를 추정한다. 셋째, 투사 경사 하강법(S‑PGD)을 사용해 전체 파라미터를 동시에 최적화한다. 이 알고리즘은 각 단계에서 저차원성 제약과 정규화를 투사함으로써 수렴성을 보장한다. 이론적 기여는 두 가지 주요 결과에 집중한다. (1) 추정된 잠재 결과 텐서 \(\hat Y\)의 Frobenius norm 오류가 O_p(√((r₁r₂r₃)/N)) 로 수렴한다는 강력한 오류 한계를 증명한다. 여기서 r₁, r₂, r₃는 각각 공간, 노출, 결과 차원의 저차원 순위이다. (2) 평균 처리 효과(ATE) 추정에 대해 이중 강건(doubly robust) 분산 추정기를 구축해, 모델이 부분적으로 잘못 지정되더라도 일관된 표준 오차를 제공한다. 이러한 결과는 기존 텐서 완성 기반 인과 추정이 제공하지 못했던 불확실성 정량화를 가능하게 한다. 시뮬레이션에서는 12가지 시나리오(공간 교란 강도, 노출‑결과 상호작용 복잡도, 누락된 잠재 결과 비율 등)를 설정해 제안 방법을 검증하였다. 모든 경우에서 제안 모델은 평균 제곱 오차와 편향 면에서 기존 방법(단일 노출 회귀, 공간 라플라시안 보정 없는 텐서 완성)보다 우수하였다. 실제 데이터 적용에서는 EPA의 PFAS 모니터링 데이터와 CDC의 질병 유병률 데이터를 결합했다. 5,495개의 공공 급수 시스템(PWS)을 대상으로, PFOA와 PFOS 두 물질을 이진 노출 변수로 정의하고, 13개의 만성 질환(관절염, 고혈압, 암, 천식 등)을 결과로 설정하였다. 인구통계학적 변수(Z)를 보정하고, 라플라시안 기반 공간 교란 보정을 적용하였다. 결과적으로, 기존 연구에서 보고된 다수의 PFAS–질환 연관성이 공간 교란 보정 후 크게 약화되었으며, 오직 PFOS 노출이 고혈압, 치아 손실, 천식, 비만과의 연관성만이 통계적으로 유의하게 남았다. 이는 공간 교란을 무시할 경우 위양성 결과가 과대 평가될 위험을 강조한다. 결론적으로, 이 논문은 (1) 다중 노출·다중 결과 인과 추정을 위한 저차원 텐서 프레임워크, (2) 미측정 공간 교란을 스펙트럴 방식으로 보정하는 새로운 알고리즘, (3) 강력한 이론적 보증과 실증적 검증을 제공함으로써, 환경 역학 및 공중보건 분야에서 복합 오염 물질의 효과를 보다 정확히 추정할 수 있는 도구를 제시한다.