umx 4.5: 트윈·경로 기반 SEM 확장과 새로운 분석 도구

초록

umx 4.5는 OpenMx와 연동된 R 패키지로, 교차시차 패널 모델, MR‑DoC, 정의 변수 직접 지정, Ωnyx 경로 가져오기, 검열 데이터 처리 등 5가지 주요 기능을 추가한다. 이를 통해 장기·인과 트윈 설계와 유전역학 분석을 보다 직관적이고 재현 가능하게 만든다.

상세 분석

umx 4.5는 기존 umx(1.8.0)에서 제공하던 간결한 문법과 자동 라벨링을 유지하면서, 최신 OpenMx 기능을 포괄한다. 가장 눈에 띄는 확장은 네 가지 고수준 함수이다. 첫째, umxCLPM은 전통적인 CLPM과 RI‑CLPM을 모두 구현하며, 무작위 절편을 도입해 개인 간 차이를 분리한다. 또한 도구변수(IV) 옵션을 제공해 교차시차 경로에 외생 변수를 삽입함으로써 인과 추론의 자유도를 늘린다. 둘째, umxMRDoC은 다중 트윈 설계에 폴리제닉 스코어(PRS)를 도구변수로 활용하는 MR‑DoC와 MR‑DoC2 모델을 구현한다. 이는 유전적 배경 혼란을 제어하면서 양방향 인과 효과와 직접적인 플리오트로피 경로를 동시에 추정할 수 있게 한다. 셋째, umxRAM()에 정의 변수를 직접 지정할 수 있게 함으로써, 행별 고유값을 이용한 평균·공분산 조정이 간편해졌다. 기존에는 정의 변수를 algebra 형태로 수동 코딩해야 했지만, 이제 umxPath(defn=‘defVar’)만으로 구현 가능하다. 넷째, Ωnyx와의 연동이 강화돼 Ωnyx에서 추출한 OpenMx RAM 코드와 경로를 umxTwinMaker()가 자동으로 ACE 구조로 변환한다. 이는 그래픽 인터페이스를 선호하는 연구자가 복잡한 트윈 모델을 손쉽게 구축하도록 돕는다.

추가적으로, 검열 데이터 분석을 위한 ICU(Integrated Censored‑Uncensored) 메서드가 xmu_make_bin_cont_pair_data() 함수로 제공된다. 이는 검출 한계 이하의 값을 이진 변수, 그 이상을 연속 변수로 변환해 하나의 모델에 동시에 포함시키며, 기대 공분산 행렬을 자동으로 트리밍한다. 정의 변수와 결합된 경우 행별 결측이 발생할 때는 xmu_update_covar()가 99999라는 자리표시자를 삽입해 데이터 손실을 최소화한다.

데이터 전처리 측면에서는 umx_residualize()가 넓은 형식(wide)과 긴 형식(long) 트윈 데이터 모두에 적용 가능하도록 설계돼, 공변량을 사전 제거하거나 정의 변수 대신 사용할 때 발생할 수 있는 외생성 가정 위반을 명시적으로 경고한다. 또한 성별 제한 모델링이 다섯 개의 트윈 그룹(예: MZ 남, MZ 여, DZ 남, DZ 여, 이성쌍)까지 확장돼 정량적·정성적 차이를 동시에 추정한다.

전반적으로 umx 4.5는 복잡한 유전·환경 모델링을 R 사용자에게 친숙한 인터페이스와 자동 보고서(umxSummary) 기능으로 제공한다. 모델 적합도 지표와 해석 가이드를 LaTeX·HTML로 출력해 논문 작성 과정을 단축하고, 파워 분석(umxPower)과 ACE, CP, simplex 등 기존 트윈 모델에 대한 공통화된 공변량 처리 기능을 강화한다. 이러한 통합은 재현 가능성, 투명성, 그리고 분석 파이프라인의 효율성을 크게 향상시킨다.