치료가 손상의 원인이 되는 순간, 그리고 방패가 문을 함께 닫는다는 것.
세균성 수막염은 보통 "균을 빨리 죽이고 스테로이드를 같이 준다"로 가르쳐집니다. 이 모델은 대신 곱 하나에서 출발합니다.
손상 flux = k_kill(t) × N(t) × Y_lysis × (1 − E_dex(t))그리고 하나의 사실만 덧붙입니다: k_kill 은 첫 항생제 투여에서 0 → E_max 로 점프하는데, N 은 정확히 그 순간 최대다.
따라서 이 곱은 치료의 첫 몇 시간에 최대가 됩니다. 항생제는 뇌를 구하는 동시에, 그 순간 가장 많은 염증 화물을 쏟아붓습니다.
시뮬레이션이 시간마다 출력하는 값입니다 (S02: 세프트리악손 단독, 스테로이드 없음):
t[h] CEF_CSF C/MIC kkill/h log10Nc 용해flux CW PLY Mg TNF
0.5 0.82 27 1.37 6.97 12.78 3.9 1.90 0.69 530
1.0 2.57 86 1.40 6.79 8.74 9.3 4.14 0.82 627
2.0 4.51 150 1.40 6.45 3.93 14.9 5.57 0.93 729
3.0 5.24 175 1.40 6.10 1.77 17.0 5.17 0.95 728
4.0 5.36 179 1.40 5.75 0.80 17.4 4.25 0.95 683
8.0 3.99 133 1.40 4.39 0.03 15.0 1.36 0.93 515
24.0 4.33 144 1.40 1.55 0.00 5.9 0.01 0.86 408
세포벽 부하는 투여 전 1.3 CWU/mL 에서 3–4 시간 뒤 17.4 로 13배 오릅니다. 용소(pneumolysin)는 2 시간에 정점을 지납니다. TNF 는 530 → 최고 737 pg/mL. 균은 사라지는데 염증은 그 직후에 최대가 됩니다 — 문헌에 기록된 "항생제 투여 후 CSF 사이토카인 상승" 이 곱의 산술적 결과입니다.
같은 항생제, 시작 시점만 다르게:
| 항생제 지연 | AUC_lysis | 최고 락테이트 | 멸균 시간 | 청력 dB | 사망 % |
|---|---|---|---|---|---|
| 0 h | 22.3 | 5.2 | 30.4 h | 0.1 | 11.3 |
| 3 h | 116.7 | 5.2 | 40.2 h | 0.2 | 16.6 |
| 6 h | 313.6 | 7.3 | 47.7 h | 8.9 | 20.5 |
| 12 h | 610.4 | 8.4 | 58.1 h | 19.2 | 25.2 |
12 시간 늦으면 쏟아지는 화물의 총량이 27배가 됩니다. N 이 μ ≈ 0.85 /h 로 자라 있기 때문이고, 그 N 을 결국 다 죽여야 하기 때문입니다. "빨리 주라"와 "천천히 죽여라"는 충돌하지 않습니다 — 최적해는 방패를 올린 상태로 최대한 빨리 입니다.
혈액-CSF 장벽 투과성 Pb 는 염증이 열고 덱사메타손이 닫습니다. 열리면
친수성 항생제가 들어오고, 동시에 알부민·호중구·물이 들어옵니다. 그래서
스테로이드는 한 팔로 돕고 다른 팔로 방해합니다. 어느 쪽이 이기는지는
그 항생제의 C/MIC 여유가 결정합니다.
| 시나리오 | VAN CSF AUC | Pb 최고 | CSF 멸균 | 사망 % |
|---|---|---|---|---|
| 감수성균 · CEF+VAN+DEX (표준치료) | 813 | 6.1 | 20.2 h | 11.2 |
| 내성균(MIC 4) · CEF+VAN, DEX 없음 | 1539 | 8.7 | 73.0 h | 43.0 |
| 내성균 · CEF+VAN**+DEX** | 1310 | 7.0 | 160.2 h | 36.1 |
| 내성균 · CEF+VAN+RIF+DEX | 789 | 6.1 | 39.9 h | 12.0 |
| 내성균 · CEF 단독 | — | 8.7 | 실패 | 78.6 |
읽는 법:
- 세팔로스포린 내성균(MIC 4 mg/L)에서 덱사메타손이 장벽을 조이면 반코마이신 CSF 노출이 15 % 줄고 멸균이 73 h → 160 h 로 2.2배 늦어집니다. 같은 조작이 감수성균의 세프트리악손에는 아무 해가 없습니다 — C/MIC 여유가 200배이기 때문입니다. 같은 약, 같은 용량, 반대 부호.
- 그리고 그 해는 친유성 약을 더하면 사라집니다: 리팜핀은 침투가 Pb 에 거의 무관하므로(a = 0.1), 스테로이드가 문을 닫아도 들어갑니다. 멸균 39.9 h, 사망 36.1 % → 12.0 %.
- 내성균에서 세프트리악손 단독은 살균 역치(4×MIC = 16 mg/L)를 336 시간 중 75.8 h 만 넘고 멸균에 실패합니다. 반코마이신을 병용하는 임상 근거가 가정이 아니라 출력으로 나옵니다.
β-락탐 살균은 용해성(lytic)이라 한 마리당 화물 수확량 Y 가 최대(1.00)이고,
리팜핀은 비용해성이라 0.15 입니다. 리팜핀을 2 시간 선행 투여하면:
| 멸균 시간 | AUC_lysis | 최고 CW | 최고 PLY | 최고 TNF | 사망 % | |
|---|---|---|---|---|---|---|
| CEF 단독 | 28.5 h | 22.3 | 17.4 | 5.58 | 737 | 26.9 |
| RIF 2 h 선행 → CEF | 17.5 h | 16.6 | 13.4 | 4.75 | 699 | 25.6 |
더 빨리 멸균되면서 화물은 26 % 적습니다. 죽이는 양은 같은데 수확량이 다르기 때문입니다.
| AUC_TNF | 최고 PMN | Pb 최고 | 청력 dB | 사망 % | |
|---|---|---|---|---|---|
| DEX 없음 | 39,800 | 3,364 | 8.4 | 0.3 | 26.9 |
| DEX 20분 전 | 19,400 | 1,847 | 6.1 | 0.1 | 11.9 |
| DEX 동시 | 18,700 | 1,849 | 6.1 | 0.1 | 11.3 |
| DEX +4 h | 20,100 | 2,327 | 6.4 | 0.1 | 11.7 |
| DEX +12 h | 22,100 | 2,960 | 6.9 | 0.1 | 12.9 |
| DEX 1일만 (4일 대신) | 24,700 | 1,868 | 6.5 | 0.2 | 16.8 |
주는 것 자체의 효과는 큽니다 (26.9 % → 11.9 %, AUC_TNF −51 %). 투여 기간도 분명합니다 (4일 11.9 % vs 1일 16.8 %). 그런데 투여 시점의 효과는 이 모델에서 거의 없습니다. 중증 환자 (균량 10⁸ · 항생제 6 h 지연)에서 다시 확인해도 −20분 21.8 % vs +4 h 21.8 % vs +24 h 24.1 % 입니다.
이유는 구조적입니다. 사이토카인 폭발은 2–6 시간 폭이지만 청력·인지·피질 손상 적분은 수십~수백 시간 규모이므로, 봉우리가 14일 적분에서 차지하는 비중이 작습니다. 이것은 편리한 결론이 아니라 반증 가능한 예측입니다. 빠진 기전의 가장 유력한 후보는 손상항의 문턱입니다: 이 모델의 모든 손상항은 구동자에 1차이므로 "3 h 동안 2배 높은 봉우리"와 "100 h 동안 1.06배 높은 수준"이 같은 손상을 냅니다. 실제 유모세포·뉴런 사멸에 문턱이 있다면 봉우리가 불균형하게 비싸지고, 그때 비로소 "항생제 전"이 산술적으로 중요해집니다. 문턱을 넣어 권고를 재현하는 것은 쉽지만 그것은 결론에 맞춰 구조를 고르는 일이므로, 여기서는 하지 않고 예측으로 남겨 둡니다.
CSF 저포도당(hypoglycorrhachia)은 보통 "세균이 포도당을 먹어서"로 설명됩니다. 이 모델에서는 두 소비항을 각각 독립적으로 계산합니다:
- 세균: 폐렴균 건조중량 0.3 pg/세포, 균질발효 수율 ≈ 25 g/mol 에서 → 10⁷ CFU/mL 당 2.0 mg/dL/h
- 호중구: 보고된 해당작용 속도에서 → 4,000 /µL 당 20 mg/dL/h
호중구가 열 배 더 먹습니다. 그래서 CSF 가 배양음성이 된 뒤에도 저포도당이 남고 (S24: 멸균 20.2 h, 포도당은 48 h 에 44.5 mg/dL, 120 h 에 51.9), CSF 포도당은 조기 멸균 지표가 될 수 없습니다. 이것은 가정이 아니라 두 계수의 비율에서 나오는 결과입니다.
그리고 GLUT1 은 양방향 운반체이므로 CSF 포도당이 떨어질수록 순유입이 커집니다. 그래서 치료된 환자에서 포도당이 0 에 닿지 않습니다.
위험함수가 만성적으로 적립되지 않는지 확인하려면 회복 궤적이 실제로 정상으로 돌아가야 합니다:
t[h] log10Nc PMN/µL 단백mg/dL 포도당 락테이트 TNF Mg Pb R_out ICP CPP MAP SOFA
0 6.97 1283 180 35.4 4.8 530 0.69 6.0 0.50 16.9 69.0 85.9 2.15
4 5.51 2017 180 37.6 3.7 682 0.95 6.3 0.65 19.7 54.0 73.7 3.04
12 2.37 2374 185 33.2 4.0 239 0.91 6.3 0.87 23.8 39.9 63.7 2.80
24 0.00 1691 173 39.0 3.5 164 0.86 5.6 0.97 25.7 37.8 63.6 1.65
48 0.00 1204 146 44.5 2.9 144 0.69 4.8 0.98 25.8 39.0 64.7 1.09
120 0.00 536 109 51.9 2.2 50 0.07 3.6 0.92 24.7 46.5 71.2 0.53
240 0.00 1 50 59.6 1.7 0 0.00 1.6 0.71 20.9 66.4 87.3 0.01
336 0.00 0 31 60.0 1.7 0 0.00 1.0 0.58 18.1 69.9 88.0 0.00
균은 24 시간에 사라지고, 장벽은 6.0 → 1.0 으로 닫히고, ICP 는 25.7 → 18.1 mmHg, CPP 는 37.8 → 69.9 mmHg, SOFA 는 0 으로 갑니다. 이 궤적이 맞지 않으면 아래의 어떤 엔드포인트도 의미가 없습니다.
Haz = h0·T + h_icp·I_icp + h_cpp·I_cpp + h_sofa·I_sofa + h_isch·I_isch
+ h_bact·I_bact + h_ncsf·I_ncsf (급성) + h_cort·(1 − N_cort) (구조)
| 시나리오 | ICP | CPP | SOFA | 허혈 | 균혈증 | CSF감염 | 구조 | 합계 | 사망 % |
|---|---|---|---|---|---|---|---|---|---|
| 무치료 | 0.006 | 0.644 | 0.119 | 0.010 | 0.106 | 0.457 | 0.325 | 1.671 | 81.2 |
| CEF 단독 | 0.002 | 0.168 | 0.015 | 0.002 | 0.001 | 0.015 | 0.108 | 0.314 | 26.9 |
| CEF+DEX | 0.000 | 0.062 | 0.011 | 0.000 | 0.001 | 0.016 | 0.035 | 0.128 | 11.9 |
| 12 h 지연 | 0.001 | 0.149 | 0.021 | 0.000 | 0.005 | 0.042 | 0.068 | 0.290 | 25.2 |
| 내성균+VAN+DEX | 0.001 | 0.190 | 0.035 | 0.001 | 0.011 | 0.115 | 0.091 | 0.448 | 36.1 |
| 표준치료 | 0.000 | 0.060 | 0.010 | 0.000 | 0.001 | 0.011 | 0.035 | 0.121 | 11.2 |
| 표준치료+EVD | 0.000 | 0.004 | 0.010 | 0.000 | 0.001 | 0.011 | 0.031 | 0.061 | 6.2 |
관류(CPP)가 모든 시나리오에서 가장 큰 항이고, 그 다음이 구조적 피질 손상입니다. 수막염은 혈액이 아니라 뇌를 통해 죽입니다. 이 분해를 상태변수로 만들어 두었기 때문에, 위험함수 계수를 눈으로 맞추지 않고 임상시험 목표치에 대해 산술로 풀 수 있었습니다 (아래).
| 환자 | 균량 log₁₀ | 지연 h | 숙주방어 | 사망%(DEX−) | 사망%(DEX+) | 청력dB(−) | 청력dB(+) |
|---|---|---|---|---|---|---|---|
| 1 | 6.0 | 0 | 1.00 | 15.9 | 6.8 | 0.1 | 0.1 |
| 2 | 6.5 | 2 | 1.00 | 27.0 | 11.4 | 0.3 | 0.1 |
| 3 | 7.0 | 1 | 0.90 | 29.4 | 12.6 | 2.3 | 0.2 |
| 4 | 7.5 | 4 | 0.85 | 38.4 | 19.5 | 21.7 | 8.5 |
| 5 | 8.0 | 6 | 0.70 | 41.0 | 22.3 | 27.4 | 15.5 |
| 6 | 7.0 | 3 | 0.55 | 36.1 | 17.5 | 16.7 | 1.7 |
| 7 | 6.7 | 1 | 1.00 | 26.2 | 10.8 | 0.3 | 0.1 |
| 8 | 7.3 | 8 | 0.60 | 41.5 | 23.0 | 28.1 | 16.2 |
| 9 | 7.9 | 2 | 0.80 | 37.8 | 18.9 | 20.8 | 7.3 |
| 10 | 7.5 | 24 | 0.40 | 48.7 | 32.5 | 43.4 | 33.6 |
| 평균 | 34.2 | 17.5 | 16.1 | 8.3 |
| 지표 | 모델 | 문헌 (폐렴균 성인) |
|---|---|---|
| 사망률, 스테로이드 없음 | 34.2 % | 34 % (de Gans & van de Beek, EDS 2002) |
| 사망률, 스테로이드 | 17.5 % | 14 % |
| 불량결과, 없음 → 있음 | 50.9 % → 26.1 % | 52 % → 26 % |
| 어떤 난청 (>25 dB) | 30 % → 10 % | 20–30 % |
위험함수 계수는 이 대조로부터 산술적으로 결정했습니다: 성분별 적분을
상태변수로 가지고 있으므로 Haz = h0·T + Σ hᵢ·Iᵢ 로 사후 분해가 되고,
코호트 평균에 대해 무-DEX 1.256·s = 0.416, DEX 0.538·s = 0.151 을 풀면
s = 0.331 / 0.281 이 나옵니다. 모델의 위험비(2.33)가 목표(2.76)보다 약간
작아 두 값이 갈리므로 타협값 s = 0.30 을 모든 계수에 적용했습니다.
맞지 않는 것도 적어 둡니다: 무치료 자연경과가 14일 사망 81.2 % 로 나오는데 실제로는 거의 100 % 입니다. 모델에 뇌간 사망·감입(herniation)을 흡수하는 상태가 없기 때문이고, 이것은 명시적 한계입니다. 중증 난청(>60 dB)도 코호트에서 한 명(늦게 온 환자)만 나와 보고치보다 적습니다.
모델을 돌리기 전에 맞아야 하는 것들 — 적합이 아니라 항등식입니다:
| 검산 | 모델 | 문헌 |
|---|---|---|
| ICP = P_ss + Q_f·R_out (정상 R_out) | 9.5 mmHg | 정상 7–15 |
| 같은 식, R_out 7배 | 31.2 mmHg | 중증 수막염 25–40 |
| 세프트리악손 유리분율 (30 → 250 mg/L) | 0.076 → 0.167 | 포화성 결합, 0.05 → 0.20 |
| CEF 총농도 CSF 침투 (Pb 1 → 9) | 0.021 → 0.094 | 0.015–0.10 |
| VAN 총농도 CSF 침투 (Pb 1 → 9) | 0.022 → 0.144 | 0.01–0.15 |
| RIF 총농도 CSF 침투 | 0.118 → 0.144 | 0.10–0.20, 염증 비의존 |
| DEX 총농도 CSF 침투 | 0.193 → 0.235 | 0.15–0.30 |
| CSF 포도당 60 mg/dL 균형 | 유입 9.3 = bulk 8.4 + 뇌 1.0 | — |
| RK4 안정조건 (ICP 고유값 33.1/h) | dt < 0.084 h, 사용 0.04 | — |
적분 수렴 검산(dt 절반)에서 최고 ICP·청력·사망확률·AUC_lysis 의 상대차는 모두 0.54 % 미만입니다.
| 파일 | 내용 |
|---|---|
abm_qsp_model.dot · .svg · .png |
기계론적 지도 — 164 노드 · 21 클러스터 · 146 간선 |
abm_mrgsolve_model.R |
63-ODE mrgsolve 모델 · 26 치료 시나리오 · 가상 코호트 |
abm_reference_python.py |
같은 63개 방정식의 의존성 없는 Python RK4 독립 구현 (검증용) |
abm_reference_output.txt |
위 실행 로그 전문 — README 의 모든 숫자의 출처 |
abm_shiny_app.R |
12 탭 인터랙티브 대시보드 |
abm_references.md |
문헌 650편 (25 섹션) — PubMed API 응답만 옮김 |
fetch_refs.py |
위 문헌 목록 생성 스크립트 (NCBI E-utilities) |
1 숙주 위험인자 · 2 집락형성→침습→CSF 진입 · 3 병독인자(방출될 화물) · 4 아집단과 살균항(곱이 만들어지는 자리) · 5 선천면역 인식 · 6 사이토카인 · 7 호중구와 단백분해 · 8 혈액-CSF 장벽(양방향 문) · 9 CSF 역학과 ICP · 10 관류·자동조절·허혈 · 11 CSF 대사 · 12 신경 손상(두 갈래) · 13 항생제 PK(전달 단계) · 14 항생제 PD(여유·용해·내성 절벽) · 15 덱사메타손(방패와 닫히는 문) · 16 보조요법 · 17 와우 손상 · 18 뇌혈관 합병증 · 19 전신 패혈증 · 20 진단과 엔드포인트 · 21 범례
| 구획 | 상태변수 |
|---|---|
| 1–3 | 세프트리악손: 중앙 · 말초 · CSF |
| 4–6 | 반코마이신: 중앙 · 말초 · CSF |
| 7–8 | 리팜핀: 중앙 · CSF |
| 9–11 | 덱사메타손: 혈장 · CSF · 비대칭 전사효과 구획 |
| 12–15 | 삼투요법: 만니톨 · 글리세롤(장·혈장) · 뇌내 축적 삼투질 |
| 16–20 | 세균: CSF 유리 · 부착/격리 · 혈중 · 세포벽 · 용소 |
| 21–30 | 대식세포활성 · TNF · IL-1β · IL-6 · IL-10 · CXCL8 · 보체 · 호중구 · MMP-9 · ROS |
| 31–40 | 장벽 Pb · CSF 알부민 · 단백 · 포도당 · 락테이트 · CSF 축적 · 유출저항 · ICP · 뇌수분 · 자동조절 |
| 41–44 | MAP · 체온 · SOFA · 용적상태 |
| 45–50 | 피질 뉴런 · 해마 치상핵 · 와우 유모세포 · 미로 골화 · 경련 부담 · 급성 위험 |
| 51–57 | 노출 추적자 (AUC_cef, AUC_van, AUC_lysis, AUC_TNF, T>4×MIC ×2, AUC_ICP) |
| 58–63 | 위험함수 성분별 적분 (I_icp, I_cpp, I_sofa, I_isch, I_bact, I_ncsf) |
핵심 기전 몇 가지:
- 비대칭 전사효과 구획: 글루코코르티코이드 억제는 빨리 켜지고(t½ 1.5 h) 느리게 꺼집니다(t½ 17 h). 이 비대칭이 4일 요법이 충분한 이유이면서, 동시에 (①에 의하면) 타이밍 창을 만들어야 하는 구조입니다.
- 부착·격리 아집단: 살균효과의 35 % 만 받으므로 언제나 늦게 사라집니다. CSF 유리균 멸균 20.2 h vs 부착 아집단 제거 30.2 h — 이 간극이 완전 치료기간이 필요한 이유입니다.
- 지수적 두개내 순응도: dP/dV = K·P (Marmarou PVI 25 mL). 항정상태에서 Davson 관계 ICP = P_ss + Q_f·R_out 로 환원되며, 그것으로 R_out 을 보정했습니다.
- 삼투요법 효력 ∝ 1/Pb: 부종을 만든 장벽 손상이 삼투제도 새게 합니다.
이 환경에는 R 런타임이 없었으므로, 63개 방정식을 의존성 없는 Python RK4로
독립 구현해서 실제로 적분했습니다. 그렇게 하지 않았으면 발견되지 않았을
결함이 33개 나왔습니다. 전체 목록은
abm_reference_python.py 머리말의 [F1]–[F33] 에
있고, 각각 고친 자리에 주석으로 남겼습니다. 구조를 바꾼 것들:
| # | 결함 | 증상 |
|---|---|---|
| F1 | 포도당 소비항에 기질 의존성이 없었다 | CSF 포도당이 0 이 된 뒤에도 소비가 계속되어 락테이트 243 mmol/L (실측 6–12). 락테이트는 포도당에서만 나온다 |
| F14 | 혈중 세균에 수용력이 없었다 | N_b 가 14일에 2×10¹⁰ CFU/mL 로 발산하고 그것이 CSF 를 역오염 |
| F15 | 삼투압 환산에 분자량 182 대신 18.2 | 만니톨 0.5 g/kg 이 혈장 삼투압을 113 mOsm/kg 올렸다 (실제 11.3). 10배 단위오류 |
| F16 | 위험함수가 영구손상을 시간에 곱했다 | 회복한 환자도 14일 내내 위험을 적립 → 26개 시나리오 중 22개가 사망확률 99 % 이상, 모든 치료 비교가 무의미 |
| F17 | 자동조절 공식이 자동조절을 해로운 것으로 만들었다 | CPP 38 에서 자동조절이 남은 환자가 압력수동 환자보다 더 허혈 (0.54 vs 0.07) |
| F25 | 피질 손실이 ROS 항으로 만성 진행 | 최적 치료 환자도 피질의 56 % 를 잃었다 |
| F28 | ICP·CPP 위험항이 선형이었다 | "CPP 45 를 100 h" 와 "CPP 5 를 10 h" 가 같은 위험 → 경증의 지속 이탈이 위험함수를 지배 |
| F29 | 소멸 하한을 도함수 블록 안에 쓰면 안 된다 | 불연속 스위치가 적분기를 멈춘다. mrgsolve 판은 상태를 건드리지 않고 보고 임계값으로만 판정 |
| F30 | 위험함수 계수를 눈으로 맞추고 있었다 | 성분별 적분을 상태변수로 추가해 임상시험 목표치에 대해 산술로 풀었다 (단일 배율 0.30) |
| F33 | t < 0 에 예정된 투여가 조용히 사라졌다 | "덱사메타손 20분 전" 팔의 실제 첫 투여가 +5.67 h 였다 — 즉 "먼저 준 팔"이 "+4 h 팔"보다 늦게 받고 있었고, 이 모델의 중심 주장이 검증 불가능한 상태였다. 증거는 출력에 있었다: "리팜핀 2 h 선행" 시나리오의 최고 세포벽·용소·TNF 가 리팜핀 없는 시나리오와 소수점 두 자리까지 같았다 (17.4 / 5.58 / 737) — 리팜핀을 전혀 받지 않았기 때문 |
F33 이 이 목록에서 가장 중요합니다. 그것을 고치기 전에는 "덱사메타손 타이밍은 중요하지 않다"가 모델의 결론처럼 보였지만, 그것은 생리가 아니라 일정표 버그였습니다. 고친 뒤에도 결론은 (그 근거는 완전히 달라진 채로) 유지되었고, 위에 적은 것처럼 반증 가능한 예측으로 남겨 두었습니다.
source("abm_mrgsolve_model.R")
structural_checks() # 손으로 검산되는 항등식들
claim1_burst() # 곱의 봉우리 (첫 24 h)
claim2_timing() # 덱사메타손 타이밍 스윕
claim3_signflip() # 부호가 갈리는 네 팔
claim4_delay() # 항생제 지연 스윕
hazard_decomposition() # 무엇이 실제로 죽이는가
run_all_scenarios() # 26 시나리오
run_cohort() # 가상 코호트 10명 × 2
shiny::runApp("abm_shiny_app.R") # 12 탭 대시보드R 이 없는 환경에서는 참조 구현을 그대로 돌릴 수 있습니다:
python3 abm_reference_python.py # 의존성 없음. abm_reference_output.txt 재생성교육·연구 목적의 반정량적 QSP 모델입니다. 공개 문헌과 임상시험 데이터로 구성했으나 독립적으로 검증·인증되지 않았으며, 실제 임상 의사결정·처방·규제 제출에 사용해서는 안 됩니다. 위에 적은 한계(무치료 사망률 과소예측, 중증 난청 꼬리 부족, 손상항의 문턱 부재)는 알려진 것만 적은 것이고, 알려지지 않은 것이 더 있을 것입니다.