Binary Emax 모델의 편향 감소 방법 평가

초록

본 논문은 소규모 임상시험에서 이진 Emax 모델의 최대우도 추정이 보이는 편향과 발산 문제를 해결하기 위해 Cox‑Snell 보정, Firth 점수 수정, Jeffreys 사전 기반 최대 페널라이즈드 우도 추정(MPLE) 세 가지 방법을 비교한다. 시뮬레이션과 실제 TURANDOT 임상 데이터 분석 결과, MPLE가 가장 안정적이며 분산이 가장 작아 실무 적용에 유리함을 확인하였다.

상세 분석

이진 Emax 모델은 용량‑반응 관계를 로짓 형태로 표현하는 비선형 회귀이며, 파라미터 (E₀, E_max, ED₅₀) 추정에 최대우도법(ML)이 흔히 사용된다. 그러나 샘플 크기가 작거나 데이터가 완전·준완전 분리(separation) 현상을 보이면 ML 추정치는 무한대로 발산하거나 편향이 크게 나타난다. 논문은 이러한 문제를 해결하기 위해 세 가지 편향 감소 기법을 수학적으로 전개한다. 첫 번째인 Cox‑Snell 보정은 ML 추정값에 1/n 차수의 편향 항 B(θ̂) 를 빼는 사후 보정 방식으로, Fisher 정보 행렬의 역행렬과 3차 도함수의 결합을 필요로 한다. 하지만 정보 행렬이 병렬하거나 분리 현상이 있으면 역행렬이 불안정해 실제 적용이 제한된다. 두 번째인 Firth 방법은 점수 함수에 ½ tr