On the magnitude, sign and ranking of recanting-twin path-specific effects

기술 요약: 회차하는 쌍둥이 경로 특이적 효과(Recanting-Twin Path-Specific Effects)의 크기, 부호 및 순위에 관하여

문제 제기
인과 매개 분석은 노출이 결과에 미치는 총 효과를 특정 인과 경로로 분해하는 것을 목표로 한다. 그러나 중간 혼란 변수(노출에 의해 영향을 받으면서 매개체-결과 관계를 혼란시키는 변수)가 존재하는 경우, 표준적인 자연 경로 특이적 효과(natural path-specific effects, PSEs)는 식별 불가능해진다. "회차하는 목격자(recanting witness)" 문제로 알려진 이 문제는, 특정 자연 PSE를 정의하기 위해 어떤 반사실적 변수가 서로 다른 노출 조건 하에서 동시에 특정 값을 가져야 한다는 모순이 발생하며, 이로 인해 관측 데이터로부터의 식별이 불가능해지기 때문에 발생한다.

이를 해결하기 위해, Díaz(2023)는 회차하는 목격자를 그 분포로부터의 무작위 추출값으로 대체함으로써 관측 데이터로부터 식별 가능한 클래스인 "회차하는 쌍둥이(recanting-twin)" PSE를 제안하였다. 이러한 효과들은 단조성(monotonicity) 하에서는 자연 PSE와 방향성이 일치하고 경로 특이적 날카로운 귀무가설(sharp null criteria)을 충족하지만, 단조성이 위배되는 경우(즉, 비선형, 비단조적 설정) 자연 PSE와 비교하여 크기, 부호 및 상대적 순위 측면에서 얼마나 잘 근사하는지는 여전히 불분명하다. 구체적으로, 회차하는 쌍둥이 프레임워크가 식별 불가능한 자연적 대응물과 비교하여 효과의 방향이나 경로 중요도의 순위를 신뢰성 있게 보존하는지는 알려져 있지 않다.

방법론
저자들은 자연 PSE와 비교하여 회차하는 쌍둥이 PSE의 거동을 규명하기 위해 분석적 유도와 시뮬레이션 연구를 결합하여 사용하였다.

1. 분석적 유도:

   • 선형 설정: 저자들은 선형 결과 모델 하에서 자연 PSE와 회차하는 쌍둥이 PSE에 대한 폐쇄형 표현식(closed-form expressions)을 먼저 유도하였다. 이들은 결과 모델에서 매개체와 중간 혼란 변수 사이에 상호작용이 없거나, 데이터 생성 과정에 관한 특정 강력한 가정(예: 일정한 개별 수준 효과)이 있는 경우 두 추정치가 정확히 일치함을 입증하였다.
   • 비모수적 경계: 이진 결과를 갖는 일반적인 설정에 대해, 저자들은 회차하는 쌍둥이 PSE와 자연 PSE 사이의 절대적 차이에 대한 식별 가능한 비모수적 상한 및 하한을 유도하였다. 이러한 경계는 총 간접 효과가 두 프레임워크 모두에서 식별 가능하다는 사실에 기초한다. 하한은 삼각 부등식을 통해 유도되었으며, 상한은 이진 결과의 제약 조건과 잠재적 결과 간의 구조적 관계를 활용하였다.

2. 시뮬레이션 연구:

   • 이진 또는 연속형 매개체($Z, M$)와 결과($Y$)를 포함하는 네 가지 시나리오에 걸쳐 종합적인 시뮬레이션 연구가 수행되었다.
   • 연구는 노출이 중간 혼란 변수에 미치는 직접 효과, 혼란 변수가 매개체에 미치는 효과, 혼란 변수가 결과에 미치는 효과, 그리고 결정적으로 중간 혼란 변수의 반사실적 값들($Z(0)$과 $Z(1)$) 사이의 상관관계($\rho$)를 변화시키며 320개의 수치적 설정을 평가하였다.
   • 저자들은 두 추정치 유형 간의 부호(효과가 양수인지 음수인지) 및 순위(네 가지 인과 경로의 상대적 크기 순서)에서의 불일치 빈도를 평가하였다.

주요 기여 및 결과

• 동등 조건: 본 논문은 매개체-중간 혼란 변수 간의 상호작용이 없는 선형 모델에서는 자연 PSE와 회차하는 쌍둥이 PSE가 동일함을 확립하였다. 이러한 상호작용이 존재하는 경우, 두 값은 서로 달라지며, 그 차이는 상호작용의 강도, 처치가 혼란 변수에 미치는 효과, 그리고 반사실적 혼란 변수 값들 사이의 상관관계에 따라 달라진다.
• 식별 가능한 경계: 두 효과 유형 사이의 절대적 차이에 대한 비모수적 상한 및 하한의 유도는 응용 연구자들에게 실질적인 도구를 제공한다. 이러한 경계는 관찰되지 않은 반사실적 상관관계를 식별하지 못하더라도, 회차하는 쌍둥이 PSE와 자연 PSE 사이의 크기 차이가 상당한지 여부를 평가할 수 있게 해준다.
• 불일치에 관한 시뮬레이션 결과:
  • 순위: 비선형, 비단조적 설정에서 두 프레임워크 간의 경로 특이적 효과 순위는 빈번하게 불일치하였다. 불일치율은 네 가지 시나리오 전반에 걸쳐 15%에서 56% 사이로 나타났다.
  • 부호: 효과의 부호에 대한 불일치는 **5%에서 43%**의 설정에서 발생하였다.
  • 불일치의 동인: 부호와 순위 모두에서 발생하는 불일치 빈도는 중간 혼란 변수의 반사실적 값들($Z(0)$과 $Z(1)$) 사이의 관찰 불가능한 상관관계($\rho$)에 의해 강력하게 영향을 받는다. 구체적으로, $\rho$가 더 음수값이 될수록 순위 불일치가 더 빈번하게 발생한다.
  • 중간 혼란의 예측 가치: 중간 혼란의 존재를 식별하는 파라미터인 $\psi_{P2 \lor P3}$는 순위 또는 부호 불일치의 정도와 관련이 없는 것으로 나타났다. 즉, 중간 혼란을 탐지하는 것이 두 PSE 유형 간의 불일치 정도를 예측하지는 못한다.
  • 불일치의 크기: 특정 시나리오(예: 이진 $Z$, 연속형 $M$, 이진 $Y$)에서 부호 불일치가 빈번하게 발생했지만, 이러한 불일치 발생 시 절대적인 크기 차이는 대개 작았다(통상 0.05 미만).

의의 및 한계
본 논문은 회차하는 쌍둥이 PSE가 중간 혼란이 존재하는 상황에서 자연 PSE에 대한 실행 가능하고 식별 가능한 대안을 제공하지만, 비단조적 설정에서는 주의해서 해석되어야 한다고 결론짓는다. 연구 결과는 인과 경로의 방향이나 상대적 중요성을 결정하기 위해 회차하는 쌍둥이 PSE에만 의존하는 것은, 특히 반사실적 혼란 변수 값들 사이의 상관관계가 음수인 경우 잘못된 결론을 도출할 수 있음을 시사한다.

저자들은 이러한 불일치를 유발하는 파라미터들(구체적으로 반사실적 값들 사이의 상관관계)이 실제로는 일반적으로 식별 불가능하기 때문에, 연구자들이 회차하는 쌍둥이 근사가 충분한 시점이 언제인지 쉽게 판단할 수 없음을 강조한다. 본 논문은 두 프레임워크가 일치하는 결론을 내는 구체적인 조건을 식별하기 위한 후속 연구를 촉구한다.

저자들이 인정한 한계점은 다음과 같다:

• 시뮬레이션 파라미터가 특정 임상적 또는 정책적 맥락과 연결되지 않아, 실질적인 영향력을 평가하는 데 한계가 있다.
• 제안된 비모수적 경계는 매끄럽거나 경로별로 미분 가능하지 않으므로, 표준적인 비모수적 추정의 점근적 성질(asymptotic properties)을 확립하는 것을 어렵게 만들지만, 파라미터적 접근법을 위한 부트스트랩 절차는 여전히 유효하다.