본문으로 건너뛰기
#frontier-candidate#vaccination#public-health

영양보충을 예방접종과 묶자 홍역백신 접종률이 올랐다

나이지리아의 지역 단위 무작위 시험을 읽는다. 새로운 백신의 효과가 아니라 영양보충을 묶은 서비스의 접종률 변화를 평가했다.

증거 단계: 지역 단위 무작위 현장시험. 접종률을 평가했으며, 새로운 백신이나 홍역 감염·사망 감소를 입증한 시험은 아니다.

#30초 결론

나이지리아의 20개 지역에서 정기 예방접종과 영양보충을 함께 제공하는 전략을 시험했다. 종료 조사에서 카드로 확인한 첫 홍역백신 접종률은 대조군 41.2%, NutriVax군 55.3%였다. 주분석의 조정 오즈비는 2.08(95% 신뢰구간 1.30–3.35)이었다. 오즈는 ‘접종한 쪽:접종하지 않은 쪽’의 비이므로, 이 값은 접종 확률이 2.08배라는 뜻이 아니다. 또 다른 분석에서 보고한 20.1%포인트 증가는 시작 수준과 지역 구조를 반영한 추정값이며, 종료 비율끼리 뺀 14.1%포인트와 구분한다. [원문 Primary outcome]

핵심 질문은 보충제가 면역을 얼마나 높였느냐가 아니라, 같은 방문에서 필요한 서비스를 함께 제공하면 예방접종 이용이 늘어나는가다. 실제 감염·사망 감소와 다른 지역에서의 효과는 이 시험만으로 확정되지 않는다.

#같은 방문에서 두 가지 필요를 해결하면

NutriVax는 백신 제품명이 아니라 영양보충을 정기 예방접종 서비스와 연결한 운영 전략이다. 연구는 북부 나이지리아 20개 지역을 두 군으로 나눴다. 보충제 제공 대상은 생후 6–23개월이지만 주요 결과를 조사한 아동은 12–23개월이었다. 서비스를 받을 수 있는 연령과 효과 측정 대상의 연령은 같지 않다.

Methods의 실제 신규 제공 대상은 6–22개월이며, 23개월에 도달하거나 12개월간 보충을 마치면 종료했다. 급성 영양실조 아동에게 같은 예방용 보충제를 그대로 제공한 것이 아니라 치료 서비스로 연결했다. 따라서 이 연구를 모든 영유아에게 동일하게 적용할 처방처럼 읽으면 안 된다.

또 중요한 운영 조건이 있다. 보호자가 예정된 접종을 거절해도, 연령 등 자격을 갖춘 아동은 영양보충을 받을 수 있었다. “백신을 맞아야만 음식을 받았다”는 조건부 보상 시험이 아니다. 두 군 모두에 예방접종·영유아 영양 관련 안내가 제공됐고, NutriVax 군에 보충제 이용 안내가 추가됐다.

지역 단위로 배정한 두 서비스 전략과 아동의 접종카드 확인 과정을 연결한 설명도
서비스를 배정한 단위와 결과를 관찰한 단위를 나누어 읽는 설명도. 원문 Methods를 바탕으로 독립 제작했으며 연구 현장 사진이 아니다.

#1,604명을 각각 무작위 배정한 시험은 아니다

무작위 배정한 것은 20개 지역이며, 종료 시점 조사에서 관찰한 아동이 1,604명이다. 같은 지역의 아동은 같은 의료시설·거리·서비스 여건의 영향을 공유할 수 있다. 따라서 모든 아동을 서로 완전히 독립적인 실험 반복으로 취급하면 불확실성을 잘못 계산할 수 있다.

또 시작과 종료 시점은 각각 별도의 횡단면 가구 조사다. 같은 1,604명을 처음부터 끝까지 추적해 접종 여부가 바뀌는 것을 센 분석이 아니다. 이 구분은 그래프의 두 시점을 개인의 성장 곡선처럼 연결하지 않기 위해 필요하다.

시작 조사는 2024년 5월에 1,597명, 종료 조사는 2025년 5월에 1,604명을 포함했다. 지역을 배정한 뒤 그 지역의 같은 연령대 집단을 두 차례 조사한 구조다. 따라서 “한 아이가 보충제를 먹은 뒤 접종했다”는 개인별 시간 순서를 이 두 비율에서 복원할 수 없다. [원문 Characteristics of clusters and survey participants]

그림을 읽고 확인하기: 위 설명도에서 무작위 배정의 화살표는 지역을 향해야 할까, 아동 한 명 한 명을 향해야 할까? 지역을 향해야 한다. 아동은 배정된 지역에서 결과를 조사한 단위다. 조사 아동 수를 늘리는 것과 무작위로 배정한 지역 수를 늘리는 것은 같은 설계 변화가 아니다.

#“두 배”는 무엇의 두 배인가

논문이 보고한 OR 2.08은 접종 오즈의 비교다. 접종 확률 자체가 2.08배라는 뜻은 아니다. 확률 p에 대응하는 오즈는 p/(1−p)라는 통계적 정의를 사용한다. 이 정의는 이 글의 설명용 선수 개념이지 연구가 새로 제안한 방정식이 아니다. 지역 내 비슷함과 시작 시점의 접종률을 고려한 분석값을 단순한 두 집단의 비율로 바꿔 읽지 않는다.

고등학생부터: 확률과 오즈를 숫자로 구분하기
접종한 쪽과 접종하지 않은 쪽의 비 odds=p1−p\mathrm{odds}=\frac{p}{1-p} 식이 길면 좌우로 이동해 읽으세요. 키보드는 수식에 초점을 둔 뒤 방향키를 사용합니다. 기호·연산·계산 과정 펼치기

p는 접종할 확률이고 1−p는 접종하지 않을 확률이다. 둘은 단위가 없는 0–1 사이 값이다. 오즈는 둘의 비다. 설명을 위한 가상값 p=0.5에서는 오즈가 0.5/(1−0.5)=1이다. 접종:미접종 오즈가 2:1인 또 다른 가상 상황에서는 확률이 2/(2+1)=2/3≈66.7%다. 이 숫자는 개념을 설명하기 위한 가상 예시이며 이번 시험의 결과가 아니다. 확률 50%가 두 배인 100%가 된다는 계산과 오즈가 두 배가 되는 계산은 다르다는 점만 보여준다.

오즈비 OR은 두 오즈의 비다. 이번 논문의 2.08은 단순히 종료 비율 두 개를 이 식에 대입한 값이 아니라 혼합효과 모형의 지역 조건부 추정이다. 가상 예시의 50%를 출발점으로 이 시험의 접종률을 예측해서는 안 된다. [선수 통계 정의; 원문 Fig. 2 및 Statistical analysis]

마찬가지로 20.1%포인트는 논문의 차이의 차이 분석 결과다. 최종 시점의 두 접종률을 빼거나 단순 전후 변화만 계산한 값과 섞어서는 안 된다. 각각 어떤 비교를 했는지 먼저 확인해야 한다.

원문 비율을 표로 나누면 차이가 보인다.

카드로 확인한 첫 홍역 접종시작 조사종료 조사단순 변화
대조군29.9%41.2%+11.3%포인트
NutriVax군25.0%55.3%+30.3%포인트

종료 시점만 빼면 14.1%포인트이고, 두 변화량을 빼면 19.0%포인트다. 어느 쪽도 논문이 모형으로 추정한 20.1%포인트와 같은 계산은 아니다. 논문은 지역별 차이와 시작 수준 등을 반영한 혼합효과 로지스틱 모형을 사용했다. 단순 비율의 산술 점검과 모형 추정값을 분리해 제시해야 한다.

표를 읽고 확인하기: NutriVax군이 30.3%포인트 올랐으니 그 전부가 전략의 효과일까? 그렇게 단정할 수 없다. 같은 기간 대조군도 11.3%포인트 올랐다. 단순 산술로는 두 변화의 차이인 19.0%포인트를 구할 수 있지만, 이것도 지역 구조를 고려한 논문의 조정 추정과 같지는 않다. ‘이전보다 증가’와 ‘대조군보다 더 증가’를 나누어 읽는 연습이다.

대학생부터: 두 번 빼는 산술과 조정 모형의 차이
표의 관측 비율로만 계산한 변화의 차이 Δraw=(55.3−25.0)−(41.2−29.9)=19.0  pp\Delta_{\mathrm{raw}}=(55.3-25.0)-(41.2-29.9)=19.0\;\mathrm{pp} 식이 길면 좌우로 이동해 읽으세요. 키보드는 수식에 초점을 둔 뒤 방향키를 사용합니다. 기호·연산·계산 과정 펼치기

첫 괄호는 NutriVax군의 전후 변화, 두 번째는 대조군의 전후 변화다. 모든 숫자는 위 표의 백분율이며, 차이의 단위 pp는 퍼센트포인트다. 이 식은 관측 비율의 산술을 보여 주기 위한 재표현이지, 논문 회귀모형을 대신하는 새 분석식이 아니다.

두 시점의 조사 대상이 서로 달라 같은 아동의 접종 전후 변화를 뜻하지 않는다. 지역 구조와 시작 수준 등을 반영한 모형은 다른 추정 과정을 사용한다. 원문의 조정 추정 20.1 pp에는 95% 신뢰구간 13.7–26.5 pp도 함께 제시된다. 19.0과 20.1의 차이를 반올림만으로 설명하거나, 같은 수치가 나와야 한다고 요구하면 안 된다. [원문 Primary outcome 및 Table 3]

주분석의 OR 2.08과 추가 차이의 차이 분석의 OR 2.39도 서로 다른 분석값이다. 둘 중 큰 값을 골라 “효과”라고 쓰면 안 된다. 접종률 차이, 오즈비, 불확실성 구간과 분석 종류를 함께 보아야 한다.

대학원부터: 같은 모형에서 오즈 척도와 확률 척도를 읽기

주분석은 종료 시점의 아동 접종 여부를 결과로 두고, 서비스 배정군과 지역별 시작 접종률을 고정효과로, 지역을 임의 절편으로 반영했다. 임의 절편은 지역마다 출발 수준이 다를 수 있음을 모형에 담는 장치다. 지역 안의 아동을 모두 독립적인 반복으로 보는 계산과 다르다. 보고된 지역 조건부 OR은 2.08(95% 신뢰구간 1.30–3.35, p=0.004)이었다. [원문 Statistical analysis 및 Fig. 2]

추가 DiD 분석에는 군·시점·군×시점 항이 들어간다. 군×시점 항은 두 군의 전후 변화 양상이 다른지를 표현한다. 이 항의 OR은 2.39(95% 신뢰구간 1.76–3.25)이고, 같은 모형에서 예측한 확률의 차이로 표현한 값이 20.1%포인트다. 보충표 4는 확률 척도의 신뢰구간을 5,000회 모수적 부트스트랩으로 구했다고 설명한다. 이는 개별 아동을 5,000번 새로 관찰했다는 뜻이 아니라, 적합 모형을 바탕으로 불확실성을 계산하는 반복 절차다. [보충표 4]

공개 심사 답변에는 시작 접종률 보정을 뺀 비교도 나온다. 저자들이 보고한 DiD는 보정 모형에서 20.1%포인트, 보정을 뺀 모형에서 20.2%포인트였다. 이는 이 비교에서 결론이 해당 보정 하나에 크게 좌우되지 않았다는 근거이지, 모든 통계 가정이나 편향 문제가 해소됐다는 뜻은 아니다. 저자는 시작 수준 보정과 지역 임의 절편이 각각 시작 시점 차이와 지역 내 상관을 다룬다고 설명했다. 이 설명과 재분석은 저자의 답변이며 여기서 원자료로 독립 재현한 결과는 아니다. [공개 Peer Review File, PDF 21쪽, 통계 심사 의견 13에 대한 답변]

2024년 통계분석계획서는 DiD를 민감도·강건성 분석으로 계획했고, 출판 논문은 추가 분석으로 제시한다. 주분석과 계획된 보조 분석이 같은 질문·척도라고 섞어 쓰지 않는다. 원문은 모형 항을 문장으로 설명하므로, 여기에 확인하지 않은 회귀계수·분포 가정을 붙여 ‘논문에 있는 회귀식’으로 만들지 않았다. [통계분석계획서 v1.0, 19–20쪽]

주분석의 불확실성을 계산하는 절차도 별도로 구분해야 한다. 통계분석계획서 v1.0은 Kenward–Roger 보정을 계획했지만, 출판 논문 Methods는 주분석의 신뢰구간과 p값을 프로파일 우도로 설명한다. 같은 Methods에서 우도비 검정을 명시한 곳은 추가 DiD 분석의 지역 조건부 OR이다. 따라서 두 분석의 검정 절차를 같은 것으로 합치지 않는다. OR 2.08을 읽을 때 ‘모든 분석 절차가 사전계획 그대로였다’고 단정하지 않는 이유다. 이번 열람 범위에서는 변경 이유나 정식 계획 변경 기록까지 확인하지 못했다. 절차가 다르다는 관측만으로 결과의 오류를 확정하는 것도 아니다.

보충표 14의 각주는 주분석 OR에도 프로파일 우도 신뢰구간과 우도비 검정을 명시한다. 또한 부평가변수들을 여러 번 검정할 때 우연한 양성 결과가 늘어나는 문제를 줄이기 위해 Benjamini–Hochberg 절차로 거짓발견률을 5% 수준에서 관리했다고 설명한다. 이는 각각의 결과가 틀릴 확률이 5%라는 뜻이 아니다. 주평가변수는 이 부평가변수 묶음과 별도로 검정했다. 따라서 여러 접종 지표의 개선을 모두 독립된 확증 결과처럼 더해서 읽지 않는다. 이 보충표 설명도 사전계획과 달라진 이유까지 알려주는 것은 아니다. [보충표 14, PDF 19쪽·인쇄 18쪽]

#기록이 없으면 어떻게 세었나

주요 분석은 예방접종 카드로 확인한 접종을 사용했다. 카드가 없거나 카드에 첫 홍역 접종 기록이 없는 아동은 미접종으로 분류했다. 따라서 서비스가 실제 접종과 카드 보관·기록에 미친 영향을 구분해서 해석할 필요가 있다. 이는 면역 수준을 직접 검사한 결과도 아니다.

카드 보유 인원은 조사 시점을 나누어 읽어야 한다. 시작 조사에서는 대조군 796명 중 384명(약 48%), NutriVax군 801명 중 318명(약 40%)이었다. 종료 조사에서는 보충표 7의 두 행정구역 인원을 합하면 대조군 474/801명(약 59%), NutriVax군 552/803명(약 69%)이다. 이는 서로 다른 횡단면 표본의 집계이지 같은 아동이 카드를 새로 보관하게 된 수가 아니다. 카드 미보유는 누락값을 대치한 것이 아니라 주평가변수의 정의에 따라 실패로 세었다. 따라서 이 결과는 우선 ‘카드로 확인한 접종’에 관한 것이며, 기록의 확보 정도가 두 군에서 같았다고 가정해서는 안 된다. [본문 Characteristics; 보충표 7, PDF 12쪽·인쇄 11쪽]

보충제를 한 번이라도 받은 비율도 배정과 다르다. NutriVax군은 386/803명, 대조군도 33/801명이 수령했다고 보고됐다. 원문이 약 4%의 지역 간 혼입으로 설명한 부분이다. 무작위 배정에 따른 분석은 서비스를 함께 제공하도록 배정한 전략의 효과를 평가하며, 실제 수령자만의 효과를 분리한 분석은 아니다.

특히 보충제를 실제 받은 아동과 받지 않은 아동의 사후 비교는 무작위 배정 비교가 아니다. 보건소를 더 쉽게 방문할 수 있거나 원래 서비스 이용 의향이 높은 가정이 더 많이 수령했을 수도 있다. 원문도 이 분석을 사후·비무작위 분석으로 표시한다. “수령한 아이의 접종률이 더 높았다”는 사실만으로 보충제 수령 자체의 인과 효과를 분리할 수 없다.

조건을 바꿔 생각하기: 두 지역의 실제 접종률은 같지만 한쪽에서 접종카드를 더 잘 보관한다면, 카드로 확인한 접종률도 같을까? 반드시 그렇지는 않다. 이 연구의 주평가변수는 카드 확인을 요구하므로 기록 확보의 차이가 측정값에 영향을 줄 수 있다. 이는 가능한 측정 경로를 설명한 가정이지, 이번 결과가 전부 카드 보관으로 생겼다는 주장은 아니다.

#무작위 배정과 눈가림은 별개다

지역은 공개 추첨으로 배정됐다. 두 행정구역 안에서 각각 대조군 5개·NutriVax군 5개 지역을 뽑아 지역 구성을 맞췄지만, 배정 결과를 숨긴 시험은 아니다. 자료 수집 담당자는 시작 조사까지 배정을 몰랐고 서비스 제공팀과 독립적이었다. 그러나 시행 뒤에는 보호자와 현장 인력이 어떤 서비스를 제공하는지 알 수 있었다. 무작위 배정은 비교의 출발점을 만드는 절차이고, 눈가림은 기대나 행동이 관찰에 미칠 영향을 줄이는 절차이므로 하나를 했다고 다른 하나도 했다고 읽지 않는다. [원문 Randomization and masking]

#다음에 확인할 것

연구 기간에는 ALIMA와 현지 당국의 조율로 백신·보충제 품절이 발생하지 않았다고 보고했다. 공급이 자주 끊기는 환경에서도 같은 결과가 나올지는 별도로 검증해야 한다. 비용효율과 세부 실행 장벽·촉진 요인은 이번 결과 논문에서 보고하지 않고 후속 분석 대상으로 남겼다. 접종률 증가만으로 비용 대비 효과까지 확인됐다고 쓰지 않는 이유다. [원문 Discussion]

이 결과를 다른 지역에 적용하려면 접근성, 공급 체계와 비용을 살펴야 한다. 서비스 방문이 늘어도 보충제 공급이 끊기거나 대기 시간이 늘면 같은 효과를 기대하기 어렵다. 접종 기록 보관과 실제 접종을 구분한 결과, 다른 지역에서의 반복 시험, 장기 운영 비용도 다음 판단에 필요하다. 이 시험의 첫 홍역백신 카드 확인 접종률 55.3%만으로 지역 전체의 감염 위험이 충분히 낮아졌다고 판단할 수는 없다. [원문 Discussion]

원문과 열람 범위

Cazes et al., Nature Medicine, DOI 10.1038/s41591-026-04675-1. 초록, Results 일부, Methods의 Statistical analysis 및 권리 조항을 확인했다. 전체 보충자료와 도판 픽셀은 미검토다. 원본 도판은 삽입하지 않았다.

추가로 역할 검토에서 공식 보충자료의 표 4·6과 2024년 통계분석계획서 19–20쪽을 대조했다. 전체 분석 재현이나 모든 자료의 검수가 완료된 것은 아니다.

공개 심사 기록에서는 통계 심사 의견 13과 저자 답변(PDF 20–21쪽)을 추가 확인했다. 전체 52쪽의 심사 기록을 검토한 것은 아니다.

Connect