Methodology, Parameters, and Calculations
health economics methodology, clinical trial cost analysis, medical research ROI, cost-benefit analysis healthcare, sensitivity analysis, Monte Carlo simulation, DALY calculation, pragmatic clinical trials
Version note. This website is the maintained interactive version of this paper. Citable archival snapshot: https://doi.org/10.5281/zenodo.18356231.
Overview
This appendix documents all 55 parameters used in the analysis, organized by type:
- External sources (peer-reviewed): 27
- Calculated values: 23
- Core definitions: 5
Calculated Values
Parameters derived from mathematical formulas and economic models.
Annual Chronic Disease Patients Treated: 982 million people
Estimated unique patients receiving chronic disease treatment annually. Derived from IQVIA days of therapy (1.28T) divided by 365 days divided by 2.5 average medications per patient times 70% post-1962 drugs.
Inputs:
- Annual Days of Chronic Disease Therapy 📊: 1.28 trillion days (95% CI: 1 trillion days - 1.5 trillion days)
\[ \begin{gathered} N_{treated} \\ = DOT_{chronic} \times 0.000767 \\ = 1.28T \times 0.000767 \\ = 982M \end{gathered} \]
Methodology:43
? Low confidence
Sensitivity Analysis
Sensitivity Indices for Annual Chronic Disease Patients Treated
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Annual Days of Chronic Disease Therapy (days) | 1.0000 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Annual Chronic Disease Patients Treated
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 982 million |
| Mean (expected value) | 981 million |
| Median (50th percentile) | 976 million |
| Standard Deviation | 98.4 million |
| 90% Range (5th-95th percentile) | [827 million, 1.15 billion] |
The histogram shows the distribution of Annual Chronic Disease Patients Treated across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Annual Chronic Disease Patients Treated will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Total DALYs Lost from Disease Eradication Delay: 7.94 billion DALYs
Total Disability-Adjusted Life Years lost from disease eradication delay (PRIMARY estimate)
Inputs:
- Years of Life Lost from Disease Eradication Delay 🔢: 7.07 billion years
- Years Lived with Disability During Disease Eradication Delay 🔢: 873 million years
\[ DALYs_{lag} = YLL_{lag} + YLD_{lag} = 7.07B + 873M = 7.94B \] where: \[ \begin{gathered} YLL_{lag} \\ = Deaths_{lag} \times (LE_{global} - Age_{death,delay}) \\ = 416M \times (79 - 62) \\ = 7.07B \end{gathered} \] where: \[ \begin{gathered} Deaths_{lag} \\ = T_{lag} \times Deaths_{disease,daily} \times 338 \\ = 8.2 \times 150{,}000 \times 338 \\ = 416M \end{gathered} \] where: \[ \begin{gathered} YLD_{lag} \\ = Deaths_{lag} \times T_{suffering} \times DW_{chronic} \\ = 416M \times 6 \times 0.35 \\ = 873M \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Total DALYs Lost from Disease Eradication Delay
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Years of Life Lost from Disease Eradication Delay (years) | 0.7043 | Strong driver |
| Years Lived with Disability During Disease Eradication Delay (years) | 0.3107 | Moderate driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Total DALYs Lost from Disease Eradication Delay
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 7.94 billion |
| Mean (expected value) | 8.05 billion |
| Median (50th percentile) | 7.89 billion |
| Standard Deviation | 2.31 billion |
| 90% Range (5th-95th percentile) | [4.43 billion, 12.1 billion] |
The histogram shows the distribution of Total DALYs Lost from Disease Eradication Delay across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Total DALYs Lost from Disease Eradication Delay will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Total Deaths from Disease Eradication Delay: 416 million deaths
Total eventually avoidable deaths from delaying disease eradication by 8.2 years (PRIMARY estimate, conservative). Excludes fundamentally unavoidable deaths (primarily accidents ~7.9%).
Inputs:
- Regulatory Delay for Efficacy Testing Post-Safety Verification 📊: 8.2 years (SE: ±2 years)
- Global Daily Deaths from Disease and Aging 📊: 150 thousand deaths/day (SE: ±7.5 thousand deaths/day)
\[ \begin{gathered} Deaths_{lag} \\ = T_{lag} \times Deaths_{disease,daily} \times 338 \\ = 8.2 \times 150{,}000 \times 338 \\ = 416M \end{gathered} \]
~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Total Deaths from Disease Eradication Delay
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Regulatory Delay for Efficacy Testing Post-Safety Verification (years) | 1.1404 | Strong driver |
| Global Daily Deaths from Disease and Aging (deaths/day) | -0.1422 | Weak driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Total Deaths from Disease Eradication Delay
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 416 million |
| Mean (expected value) | 420 million |
| Median (50th percentile) | 414 million |
| Standard Deviation | 122 million |
| 90% Range (5th-95th percentile) | [225 million, 630 million] |
The histogram shows the distribution of Total Deaths from Disease Eradication Delay across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Total Deaths from Disease Eradication Delay will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Total Economic Loss from Disease Eradication Delay: $1.19 quadrillion
Total economic loss from delaying disease eradication by 8.2 years (PRIMARY estimate, 2024 USD). Values global DALYs at standardized US/International normative rate ($150k) rather than local ability-to-pay, representing the full human capital loss.
Inputs:
- Total DALYs Lost from Disease Eradication Delay 🔢: 7.94 billion DALYs
- Standard Economic Value per QALY 📊: $150K (SE: ±$30K)
\[ \begin{gathered} Value_{lag} \\ = DALYs_{lag} \times Value_{QALY} \\ = 7.94B \times \$150K \\ = \$1190T \end{gathered} \] where: \[ DALYs_{lag} = YLL_{lag} + YLD_{lag} = 7.07B + 873M = 7.94B \] where: \[ \begin{gathered} YLL_{lag} \\ = Deaths_{lag} \times (LE_{global} - Age_{death,delay}) \\ = 416M \times (79 - 62) \\ = 7.07B \end{gathered} \] where: \[ \begin{gathered} Deaths_{lag} \\ = T_{lag} \times Deaths_{disease,daily} \times 338 \\ = 8.2 \times 150{,}000 \times 338 \\ = 416M \end{gathered} \] where: \[ \begin{gathered} YLD_{lag} \\ = Deaths_{lag} \times T_{suffering} \times DW_{chronic} \\ = 416M \times 6 \times 0.35 \\ = 873M \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Total Economic Loss from Disease Eradication Delay
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Total DALYs Lost from Disease Eradication Delay (DALYs) | 1.0671 | Strong driver |
| Standard Economic Value per QALY (USD/QALY) | -0.0733 | Minimal effect |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Total Economic Loss from Disease Eradication Delay
| Statistic | Value |
|---|---|
| Baseline (deterministic) | $1.19 quadrillion |
| Mean (expected value) | $1.27 quadrillion |
| Median (50th percentile) | $1.18 quadrillion |
| Standard Deviation | $581T |
| 90% Range (5th-95th percentile) | [$443T, $2.41 quadrillion] |
The histogram shows the distribution of Total Economic Loss from Disease Eradication Delay across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Total Economic Loss from Disease Eradication Delay will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Years Lived with Disability During Disease Eradication Delay: 873 million years
Years Lived with Disability during disease eradication delay (PRIMARY estimate)
Inputs:
- Total Deaths from Disease Eradication Delay 🔢: 416 million deaths
- Pre-Death Suffering Period During Post-Safety Efficacy Delay 📊: 6 years (95% CI: 4 years - 9 years)
- Disability Weight for Untreated Chronic Conditions 📊: 0.35 weight (SE: ±0.07 weight)
\[ \begin{gathered} YLD_{lag} \\ = Deaths_{lag} \times T_{suffering} \times DW_{chronic} \\ = 416M \times 6 \times 0.35 \\ = 873M \end{gathered} \] where: \[ \begin{gathered} Deaths_{lag} \\ = T_{lag} \times Deaths_{disease,daily} \times 338 \\ = 8.2 \times 150{,}000 \times 338 \\ = 416M \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Years Lived with Disability During Disease Eradication Delay
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Pre-Death Suffering Period During Post-Safety Efficacy Delay (years) | 2.0883 | Strong driver |
| Disability Weight for Untreated Chronic Conditions (weight) | -0.9003 | Strong driver |
| Total Deaths from Disease Eradication Delay (deaths) | -0.2255 | Weak driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Years Lived with Disability During Disease Eradication Delay
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 873 million |
| Mean (expected value) | 1.02 billion |
| Median (50th percentile) | 846 million |
| Standard Deviation | 716 million |
| 90% Range (5th-95th percentile) | [217 million, 2.43 billion] |
The histogram shows the distribution of Years Lived with Disability During Disease Eradication Delay across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Years Lived with Disability During Disease Eradication Delay will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Years of Life Lost from Disease Eradication Delay: 7.07 billion years
Years of Life Lost from disease eradication delay deaths (PRIMARY estimate)
Inputs:
- Total Deaths from Disease Eradication Delay 🔢: 416 million deaths
- Global Life Expectancy (2024) 📊: 79 years (SE: ±2 years)
- Mean Age of Preventable Death from Post-Safety Efficacy Delay 📊: 62 years (SE: ±3 years)
\[ \begin{gathered} YLL_{lag} \\ = Deaths_{lag} \times (LE_{global} - Age_{death,delay}) \\ = 416M \times (79 - 62) \\ = 7.07B \end{gathered} \] where: \[ \begin{gathered} Deaths_{lag} \\ = T_{lag} \times Deaths_{disease,daily} \times 338 \\ = 8.2 \times 150{,}000 \times 338 \\ = 416M \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Years of Life Lost from Disease Eradication Delay
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Global Life Expectancy (2024) (years) | 2.0066 | Strong driver |
| Mean Age of Preventable Death from Post-Safety Efficacy Delay (years) | -1.3852 | Strong driver |
| Total Deaths from Disease Eradication Delay (deaths) | 0.3779 | Moderate driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Years of Life Lost from Disease Eradication Delay
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 7.07 billion |
| Mean (expected value) | 7.03 billion |
| Median (50th percentile) | 7.05 billion |
| Standard Deviation | 1.62 billion |
| 90% Range (5th-95th percentile) | [4.21 billion, 9.68 billion] |
The histogram shows the distribution of Years of Life Lost from Disease Eradication Delay across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Years of Life Lost from Disease Eradication Delay will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
dFDA Trial Cost Reduction Percentage: 97.7%
Trial cost reduction percentage: 1 - (dFDA pragmatic cost / traditional Phase 3 cost)
Inputs:
- dFDA Pragmatic Trial Cost per Patient 📊: $929 (95% CI: $97 - $3K)
- Phase 3 Cost per Patient 📊: $41K (95% CI: $20K - $120K)
\[ \begin{gathered} Reduce_{pct} \\ = 1 - \frac{Cost_{pragmatic,pt}}{Cost_{P3,pt}} \\ = 1 - \frac{\$929}{\$41K} \\ = 97.7\% \end{gathered} \]
✓ High confidence
Sensitivity Analysis
Sensitivity Indices for dFDA Trial Cost Reduction Percentage
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| dFDA Pragmatic Trial Cost per Patient (USD/patient) | -6.4207 | Strong driver |
| Phase 3 Cost per Patient (USD/patient) | 5.6539 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: dFDA Trial Cost Reduction Percentage
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 97.7% |
| Mean (expected value) | 98% |
| Median (50th percentile) | 97.9% |
| Standard Deviation | 0.401% |
| 90% Range (5th-95th percentile) | [97.5%, 98.9%] |
The histogram shows the distribution of dFDA Trial Cost Reduction Percentage across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that dFDA Trial Cost Reduction Percentage will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Diseases Without Effective Treatment: 6.65 thousand diseases
Number of diseases without effective treatment. 95% of 7,000 rare diseases lack FDA-approved treatment (per Orphanet 2024). This represents the therapeutic search space that remains unexplored.
Inputs:
- Total Number of Rare Diseases Globally 📊: 7 thousand diseases (95% CI: 6 thousand diseases - 10 thousand diseases)
\[ \begin{gathered} N_{untreated} \\ = N_{rare} \times 0.95 \\ = 7{,}000 \times 0.95 \\ = 6{,}650 \end{gathered} \]
Methodology:137
~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Diseases Without Effective Treatment
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Total Number of Rare Diseases Globally (diseases) | 1.0000 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Diseases Without Effective Treatment
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 6.65 thousand |
| Mean (expected value) | 6.73 thousand |
| Median (50th percentile) | 6.64 thousand |
| Standard Deviation | 835 |
| 90% Range (5th-95th percentile) | [5.7 thousand, 8.24 thousand] |
The histogram shows the distribution of Diseases Without Effective Treatment across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Diseases Without Effective Treatment will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Total Drugs Approved Since 1962: 3.1 thousand drugs
Estimated total drugs approved globally since 1962 (62 years × average approval rate). Conservative: uses current rate, actual historical rate was lower in 1960s-80s.
Inputs:
- Average Annual New Drug Approvals Globally 📊: 50 drugs/year (95% CI: 45 drugs/year - 60 drugs/year)
\[ \begin{gathered} N_{drugs,62} \\ = Drugs_{ann,curr} \times 62 \\ = 50 \times 62 \\ = 3{,}100 \end{gathered} \]
Methodology:14
~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Total Drugs Approved Since 1962
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Average Annual New Drug Approvals Globally (drugs/year) | 1.0000 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Total Drugs Approved Since 1962
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 3.1 thousand |
| Mean (expected value) | 3.11 thousand |
| Median (50th percentile) | 3.09 thousand |
| Standard Deviation | 220 |
| 90% Range (5th-95th percentile) | [2.79 thousand, 3.5 thousand] |
The histogram shows the distribution of Total Drugs Approved Since 1962 across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Total Drugs Approved Since 1962 will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Drug Cost Increase: 1980s to Current: 13.4x
Drug development cost increase from 1980s to current
Inputs:
- Drug Development Cost (1980s) 📊: $194M (95% CI: $146M - $242M)
- Pharma Drug Development Cost (Current System) 📊: $2.6B (95% CI: $1.5B - $4B)
\[ \begin{gathered} k_{cost,80s} \\ = \frac{Cost_{dev,curr}}{Cost_{dev,80s}} \\ = \frac{\$2.6B}{\$194M} \\ = 13.4 \end{gathered} \]
Methodology:22
✓ High confidence
Sensitivity Analysis
Sensitivity Indices for Drug Cost Increase: 1980s to Current
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Pharma Drug Development Cost (Current System) (USD) | 1.6909 | Strong driver |
| Drug Development Cost (1980s) (USD) | -0.7048 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Drug Cost Increase: 1980s to Current
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 13.4x |
| Mean (expected value) | 13.3x |
| Median (50th percentile) | 13.3x |
| Standard Deviation | 0.915x |
| 90% Range (5th-95th percentile) | [11.9x, 14.7x] |
The histogram shows the distribution of Drug Cost Increase: 1980s to Current across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Drug Cost Increase: 1980s to Current will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Cumulative Efficacy Testing Cost (1962-2024): $4.84T
Cumulative Phase 2/3 efficacy testing cost since 1962. Uses direct Phase 2/3 cost per drug - this is a LOWER BOUND because it excludes opportunity cost of delays, compounds abandoned due to cost barrier, and regulatory overhead.
Inputs:
- Pharma Phase 2/3 Cost Barrier Per Drug: $1.56B (SE: ±$200M)
- Total Drugs Approved Since 1962 🔢: 3.1 thousand drugs
\[ \begin{gathered} Cost_{eff,cumul} \\ = Cost_{P2+P3} \times N_{drugs,62} \\ = \$1.56B \times 3{,}100 \\ = \$4.84T \end{gathered} \] where: \[ \begin{gathered} N_{drugs,62} \\ = Drugs_{ann,curr} \times 62 \\ = 50 \times 62 \\ = 3{,}100 \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Cumulative Efficacy Testing Cost (1962-2024)
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Total Drugs Approved Since 1962 (drugs) | 0.5385 | Strong driver |
| Pharma Phase 2/3 Cost Barrier Per Drug (USD) | 0.4652 | Moderate driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Cumulative Efficacy Testing Cost (1962-2024)
| Statistic | Value |
|---|---|
| Baseline (deterministic) | $4.84T |
| Mean (expected value) | $4.88T |
| Median (50th percentile) | $4.81T |
| Standard Deviation | $977B |
| 90% Range (5th-95th percentile) | [$3.42T, $6.62T] |
The histogram shows the distribution of Cumulative Efficacy Testing Cost (1962-2024) across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Cumulative Efficacy Testing Cost (1962-2024) will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Efficacy Lag Deaths (9/11 Equivalents): 34.1 thousand 9/11s
Total deaths from efficacy lag expressed in 9/11 equivalents. Makes the mortality cost viscerally understandable: how many September 11ths worth of deaths did the 1962 efficacy requirements cause?
Inputs:
- Total Deaths from Historical Progress Delays 🔢: 102 million deaths
- September 11 Deaths 📊: 2.98 thousand people
\[ \begin{gathered} N_{9/11,equiv} \\ = \frac{Deaths_{lag,total}}{N_{9/11}} \\ = \frac{102M}{2{,}980} \\ = 34{,}100 \end{gathered} \] where: \[ \begin{gathered} Deaths_{lag,total} \\ = Lives_{saved,annual} \times T_{lag} \\ = 12.4M \times 8.2 \\ = 102M \end{gathered} \] where: \[ \begin{gathered} Lives_{saved,annual} \\ = \frac{LY_{saved,annual}}{T_{ext}} \\ = \frac{149M}{12} \\ = 12.4M \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Efficacy Lag Deaths (9/11 Equivalents)
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Total Deaths from Historical Progress Delays (deaths) | 1.0000 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Efficacy Lag Deaths (9/11 Equivalents)
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 34.1 thousand |
| Mean (expected value) | 36 thousand |
| Median (50th percentile) | 32.7 thousand |
| Standard Deviation | 17.8 thousand |
| 90% Range (5th-95th percentile) | [12.4 thousand, 71.8 thousand] |
The histogram shows the distribution of Efficacy Lag Deaths (9/11 Equivalents) across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Efficacy Lag Deaths (9/11 Equivalents) will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Treatment Delay YLD - Annual: 2.01 billion DALYs
Annual YLD from treatment delay: patients receiving chronic disease treatment would have collectively avoided this disability if treatments were available 8.2 years earlier. Represents morbidity burden for treatment beneficiaries (distinct from mortality burden).
Inputs:
- Annual Chronic Disease Patients Treated 🔢: 982 million people
- Regulatory Delay for Efficacy Testing Post-Safety Verification 📊: 8.2 years (SE: ±2 years)
- Treatment Disability Reduction 📊: 0.25 weight (95% CI: 0.15 weight - 0.35 weight)
\[ \begin{gathered} YLD_{treat\_delay} \\ = N_{treated} \times T_{lag} \times \Delta DW_{treat} \\ = 982M \times 8.2 \times 0.25 \\ = 2.01B \end{gathered} \] where: \[ \begin{gathered} N_{treated} \\ = DOT_{chronic} \times 0.000767 \\ = 1.28T \times 0.000767 \\ = 982M \end{gathered} \] ? Low confidence
Sensitivity Analysis
Sensitivity Indices for Treatment Delay YLD - Annual
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Annual Chronic Disease Patients Treated (people) | 3.0959 | Strong driver |
| Treatment Disability Reduction (weight) | -2.4506 | Strong driver |
| Regulatory Delay for Efficacy Testing Post-Safety Verification (years) | 0.3319 | Moderate driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Treatment Delay YLD - Annual
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 2.01 billion |
| Mean (expected value) | 2.2 billion |
| Median (50th percentile) | 1.99 billion |
| Standard Deviation | 1.18 billion |
| 90% Range (5th-95th percentile) | [661 million, 4.41 billion] |
The histogram shows the distribution of Treatment Delay YLD - Annual across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Treatment Delay YLD - Annual will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Total Deaths from Historical Progress Delays: 102 million deaths
Total deaths from delaying existing drugs over 8.2-year efficacy lag. One-time impact of eliminating Phase 2-4 testing delay for drugs already approved 1962-2024. Based on Lichtenberg (2019) estimate of 12M lives saved annually × 8.2 years efficacy lag. Excludes innovation acceleration effects.
Inputs:
- Annual Lives Saved by Pharmaceuticals 🔢: 12.4 million deaths
- Regulatory Delay for Efficacy Testing Post-Safety Verification 📊: 8.2 years (SE: ±2 years)
\[ \begin{gathered} Deaths_{lag,total} \\ = Lives_{saved,annual} \times T_{lag} \\ = 12.4M \times 8.2 \\ = 102M \end{gathered} \] where: \[ \begin{gathered} Lives_{saved,annual} \\ = \frac{LY_{saved,annual}}{T_{ext}} \\ = \frac{149M}{12} \\ = 12.4M \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Total Deaths from Historical Progress Delays
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Annual Lives Saved by Pharmaceuticals (deaths) | 1.2721 | Strong driver |
| Regulatory Delay for Efficacy Testing Post-Safety Verification (years) | -0.2811 | Weak driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Total Deaths from Historical Progress Delays
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 102 million |
| Mean (expected value) | 107 million |
| Median (50th percentile) | 97.3 million |
| Standard Deviation | 53 million |
| 90% Range (5th-95th percentile) | [36.9 million, 214 million] |
The histogram shows the distribution of Total Deaths from Historical Progress Delays across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Total Deaths from Historical Progress Delays will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Annual Lives Saved by Pharmaceuticals: 12.4 million deaths
Annual lives saved by pharmaceutical interventions globally. Derived from Lichtenberg (2019) finding of 148.7M life-years saved, divided by assumed 12-year average life extension per beneficiary. Note: Life-years is the primary metric; lives is an approximation for intuitive communication.
Inputs:
- Annual Life-Years Saved by Pharmaceuticals 📊: 149 million life-years (95% CI: 79.4 million life-years - 240 million life-years)
- Average Life Extension per Beneficiary: 12 years (95% CI: 8 years - 18 years)
\[ \begin{gathered} Lives_{saved,annual} \\ = \frac{LY_{saved,annual}}{T_{ext}} \\ = \frac{149M}{12} \\ = 12.4M \end{gathered} \]
Methodology:74
? Low confidence
Sensitivity Analysis
Sensitivity Indices for Annual Lives Saved by Pharmaceuticals
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Annual Life-Years Saved by Pharmaceuticals (life-years) | 1.0000 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Annual Lives Saved by Pharmaceuticals
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 12.4 million |
| Mean (expected value) | 12.3 million |
| Median (50th percentile) | 11.9 million |
| Standard Deviation | 3.2 million |
| 90% Range (5th-95th percentile) | [7.6 million, 18.6 million] |
The histogram shows the distribution of Annual Lives Saved by Pharmaceuticals across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Annual Lives Saved by Pharmaceuticals will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Thalidomide DALYs Per Event: 41.8 thousand DALYs
Total DALYs per US-scale thalidomide event (YLL + YLD)
Inputs:
- Thalidomide YLD Per Event 🔢: 13 thousand years
- Thalidomide YLL Per Event 🔢: 28.8 thousand years
\[ \begin{gathered} DALY_{thal} \\ = YLD_{thal} + YLL_{thal} \\ = 13{,}000 + 28{,}800 \\ = 41{,}800 \end{gathered} \] where: \[ \begin{gathered} YLD_{thal} \\ = DW_{thal} \times N_{thal,survive} \times LE_{thal} \\ = 0.4 \times 540 \times 60 \\ = 13{,}000 \end{gathered} \] where: \[ \begin{gathered} N_{thal,survive} \\ = N_{thal,US,prevent} \times (1 - Rate_{thal,mort}) \\ = 900 \times (1 - 40\%) \\ = 540 \end{gathered} \] where: \[ \begin{gathered} N_{thal,US,prevent} \\ = N_{thal,global} \times Pct_{US,1960} \\ = 15{,}000 \times 6\% \\ = 900 \end{gathered} \] where: \[ \begin{gathered} YLL_{thal} \\ = Deaths_{thal} \times 80 \\ = 360 \times 80 \\ = 28{,}800 \end{gathered} \] where: \[ \begin{gathered} Deaths_{thal} \\ = Rate_{thal,mort} \times N_{thal,US,prevent} \\ = 40\% \times 900 \\ = 360 \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Thalidomide DALYs Per Event
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Thalidomide YLL Per Event (years) | 0.6300 | Strong driver |
| Thalidomide YLD Per Event (years) | 0.3701 | Moderate driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Thalidomide DALYs Per Event
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 41.8 thousand |
| Mean (expected value) | 42.5 thousand |
| Median (50th percentile) | 40.8 thousand |
| Standard Deviation | 12.2 thousand |
| 90% Range (5th-95th percentile) | [24.8 thousand, 67.1 thousand] |
The histogram shows the distribution of Thalidomide DALYs Per Event across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Thalidomide DALYs Per Event will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Thalidomide Deaths Per Event: 360 deaths
Deaths per US-scale thalidomide event
Inputs:
- Thalidomide Mortality Rate 📊: 40% (95% CI: 35% - 45%)
- Thalidomide US Cases Prevented 🔢: 900 cases
\[ \begin{gathered} Deaths_{thal} \\ = Rate_{thal,mort} \times N_{thal,US,prevent} \\ = 40\% \times 900 \\ = 360 \end{gathered} \] where: \[ \begin{gathered} N_{thal,US,prevent} \\ = N_{thal,global} \times Pct_{US,1960} \\ = 15{,}000 \times 6\% \\ = 900 \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Thalidomide Deaths Per Event
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Thalidomide US Cases Prevented (cases) | 1.5027 | Strong driver |
| Thalidomide Mortality Rate (percentage) | -0.5048 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Thalidomide Deaths Per Event
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 360 |
| Mean (expected value) | 364 |
| Median (50th percentile) | 353 |
| Standard Deviation | 95.8 |
| 90% Range (5th-95th percentile) | [223, 556] |
The histogram shows the distribution of Thalidomide Deaths Per Event across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Thalidomide Deaths Per Event will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Thalidomide Survivors Per Event: 540 cases
Survivors per US-scale thalidomide event
Inputs:
- Thalidomide Mortality Rate 📊: 40% (95% CI: 35% - 45%)
- Thalidomide US Cases Prevented 🔢: 900 cases
\[ \begin{gathered} N_{thal,survive} \\ = N_{thal,US,prevent} \times (1 - Rate_{thal,mort}) \\ = 900 \times (1 - 40\%) \\ = 540 \end{gathered} \] where: \[ \begin{gathered} N_{thal,US,prevent} \\ = N_{thal,global} \times Pct_{US,1960} \\ = 15{,}000 \times 6\% \\ = 900 \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Thalidomide Survivors Per Event
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Thalidomide Mortality Rate (percentage) | 0.5607 | Strong driver |
| Thalidomide US Cases Prevented (cases) | 0.4398 | Moderate driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Thalidomide Survivors Per Event
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 540 |
| Mean (expected value) | 537 |
| Median (50th percentile) | 531 |
| Standard Deviation | 86.3 |
| 90% Range (5th-95th percentile) | [399, 698] |
The histogram shows the distribution of Thalidomide Survivors Per Event across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Thalidomide Survivors Per Event will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Thalidomide US Cases Prevented: 900 cases
Estimated US thalidomide cases prevented by FDA rejection
Inputs:
- Thalidomide Cases Worldwide 📊: 15 thousand cases (95% CI: 10 thousand cases - 20 thousand cases)
- US Population Share 1960 📊: 6% (95% CI: 5.5% - 6.5%)
\[ \begin{gathered} N_{thal,US,prevent} \\ = N_{thal,global} \times Pct_{US,1960} \\ = 15{,}000 \times 6\% \\ = 900 \end{gathered} \]
~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Thalidomide US Cases Prevented
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Thalidomide Cases Worldwide (cases) | 1.3746 | Strong driver |
| US Population Share 1960 (percentage) | -0.3756 | Moderate driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Thalidomide US Cases Prevented
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 900 |
| Mean (expected value) | 901 |
| Median (50th percentile) | 884 |
| Standard Deviation | 182 |
| 90% Range (5th-95th percentile) | [622, 1.25 thousand] |
The histogram shows the distribution of Thalidomide US Cases Prevented across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Thalidomide US Cases Prevented will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Thalidomide YLD Per Event: 13 thousand years
Years Lived with Disability per thalidomide event
Inputs:
- Thalidomide Disability Weight 📊: 0.4:1 (95% CI: 0.32:1 - 0.48:1)
- Thalidomide Survivors Per Event 🔢: 540 cases
- Thalidomide Survivor Lifespan 📊: 60 years (95% CI: 50 years - 70 years)
\[ \begin{gathered} YLD_{thal} \\ = DW_{thal} \times N_{thal,survive} \times LE_{thal} \\ = 0.4 \times 540 \times 60 \\ = 13{,}000 \end{gathered} \] where: \[ \begin{gathered} N_{thal,survive} \\ = N_{thal,US,prevent} \times (1 - Rate_{thal,mort}) \\ = 900 \times (1 - 40\%) \\ = 540 \end{gathered} \] where: \[ \begin{gathered} N_{thal,US,prevent} \\ = N_{thal,global} \times Pct_{US,1960} \\ = 15{,}000 \times 6\% \\ = 900 \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Thalidomide YLD Per Event
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Thalidomide Disability Weight (ratio) | 28.4785 | Strong driver |
| Thalidomide Survivor Lifespan (years) | -23.4440 | Strong driver |
| Thalidomide Survivors Per Event (cases) | -4.0444 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Thalidomide YLD Per Event
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 13 thousand |
| Mean (expected value) | 13.3 thousand |
| Median (50th percentile) | 12.6 thousand |
| Standard Deviation | 4.5 thousand |
| 90% Range (5th-95th percentile) | [6.94 thousand, 22.6 thousand] |
The histogram shows the distribution of Thalidomide YLD Per Event across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Thalidomide YLD Per Event will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Thalidomide YLL Per Event: 28.8 thousand years
Years of Life Lost per thalidomide event (infant deaths)
Inputs:
- Thalidomide Deaths Per Event 🔢: 360 deaths
\[ \begin{gathered} YLL_{thal} \\ = Deaths_{thal} \times 80 \\ = 360 \times 80 \\ = 28{,}800 \end{gathered} \] where: \[ \begin{gathered} Deaths_{thal} \\ = Rate_{thal,mort} \times N_{thal,US,prevent} \\ = 40\% \times 900 \\ = 360 \end{gathered} \] where: \[ \begin{gathered} N_{thal,US,prevent} \\ = N_{thal,global} \times Pct_{US,1960} \\ = 15{,}000 \times 6\% \\ = 900 \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Thalidomide YLL Per Event
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Thalidomide Deaths Per Event (deaths) | 1.0000 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Thalidomide YLL Per Event
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 28.8 thousand |
| Mean (expected value) | 29.2 thousand |
| Median (50th percentile) | 28.2 thousand |
| Standard Deviation | 7.67 thousand |
| 90% Range (5th-95th percentile) | [17.9 thousand, 44.5 thousand] |
The histogram shows the distribution of Thalidomide YLL Per Event across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Thalidomide YLL Per Event will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Ratio of Type II Error Cost to Type I Error Benefit: 3.07k:1
Ratio of Type II error cost to Type I error benefit (harm from delay vs. harm prevented)
Inputs:
- Total DALYs Lost from Disease Eradication Delay 🔢: 7.94 billion DALYs
- Maximum DALYs Saved by FDA Preventing Unsafe Drugs (1962-2024) 🔢: 2.59 million DALYs
\[ \begin{gathered} Ratio_{TypeII} \\ = \frac{DALYs_{lag}}{DALY_{TypeI}} \\ = \frac{7.94B}{2.59M} \\ = 3{,}070 \end{gathered} \] where: \[ DALYs_{lag} = YLL_{lag} + YLD_{lag} = 7.07B + 873M = 7.94B \] where: \[ \begin{gathered} YLL_{lag} \\ = Deaths_{lag} \times (LE_{global} - Age_{death,delay}) \\ = 416M \times (79 - 62) \\ = 7.07B \end{gathered} \] where: \[ \begin{gathered} Deaths_{lag} \\ = T_{lag} \times Deaths_{disease,daily} \times 338 \\ = 8.2 \times 150{,}000 \times 338 \\ = 416M \end{gathered} \] where: \[ \begin{gathered} YLD_{lag} \\ = Deaths_{lag} \times T_{suffering} \times DW_{chronic} \\ = 416M \times 6 \times 0.35 \\ = 873M \end{gathered} \] where: \[ \begin{gathered} DALY_{TypeI} \\ = DALY_{thal} \times 62 \\ = 41{,}800 \times 62 \\ = 2.59M \end{gathered} \] where: \[ \begin{gathered} DALY_{thal} \\ = YLD_{thal} + YLL_{thal} \\ = 13{,}000 + 28{,}800 \\ = 41{,}800 \end{gathered} \] where: \[ \begin{gathered} YLD_{thal} \\ = DW_{thal} \times N_{thal,survive} \times LE_{thal} \\ = 0.4 \times 540 \times 60 \\ = 13{,}000 \end{gathered} \] where: \[ \begin{gathered} N_{thal,survive} \\ = N_{thal,US,prevent} \times (1 - Rate_{thal,mort}) \\ = 900 \times (1 - 40\%) \\ = 540 \end{gathered} \] where: \[ \begin{gathered} N_{thal,US,prevent} \\ = N_{thal,global} \times Pct_{US,1960} \\ = 15{,}000 \times 6\% \\ = 900 \end{gathered} \] where: \[ \begin{gathered} YLL_{thal} \\ = Deaths_{thal} \times 80 \\ = 360 \times 80 \\ = 28{,}800 \end{gathered} \] where: \[ \begin{gathered} Deaths_{thal} \\ = Rate_{thal,mort} \times N_{thal,US,prevent} \\ = 40\% \times 900 \\ = 360 \end{gathered} \] ~ Medium confidence
Sensitivity Analysis
Sensitivity Indices for Ratio of Type II Error Cost to Type I Error Benefit
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Total DALYs Lost from Disease Eradication Delay (DALYs) | 7.2872 | Strong driver |
| Maximum DALYs Saved by FDA Preventing Unsafe Drugs (1962-2024) (DALYs) | -7.1207 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Ratio of Type II Error Cost to Type I Error Benefit
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 3.07k:1 |
| Mean (expected value) | 3.05k:1 |
| Median (50th percentile) | 3.09k:1 |
| Standard Deviation | 101:1 |
| 90% Range (5th-95th percentile) | [2.88k:1, 3.12k:1] |
The histogram shows the distribution of Ratio of Type II Error Cost to Type I Error Benefit across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Ratio of Type II Error Cost to Type I Error Benefit will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
Maximum DALYs Saved by FDA Preventing Unsafe Drugs (1962-2024): 2.59 million DALYs
Maximum DALYs saved by FDA preventing unsafe drugs over 62-year period 1962-2024 (extreme overestimate: one Thalidomide-scale event per year)
Inputs:
- Thalidomide DALYs Per Event 🔢: 41.8 thousand DALYs
\[ \begin{gathered} DALY_{TypeI} \\ = DALY_{thal} \times 62 \\ = 41{,}800 \times 62 \\ = 2.59M \end{gathered} \] where: \[ \begin{gathered} DALY_{thal} \\ = YLD_{thal} + YLL_{thal} \\ = 13{,}000 + 28{,}800 \\ = 41{,}800 \end{gathered} \] where: \[ \begin{gathered} YLD_{thal} \\ = DW_{thal} \times N_{thal,survive} \times LE_{thal} \\ = 0.4 \times 540 \times 60 \\ = 13{,}000 \end{gathered} \] where: \[ \begin{gathered} N_{thal,survive} \\ = N_{thal,US,prevent} \times (1 - Rate_{thal,mort}) \\ = 900 \times (1 - 40\%) \\ = 540 \end{gathered} \] where: \[ \begin{gathered} N_{thal,US,prevent} \\ = N_{thal,global} \times Pct_{US,1960} \\ = 15{,}000 \times 6\% \\ = 900 \end{gathered} \] where: \[ \begin{gathered} YLL_{thal} \\ = Deaths_{thal} \times 80 \\ = 360 \times 80 \\ = 28{,}800 \end{gathered} \] where: \[ \begin{gathered} Deaths_{thal} \\ = Rate_{thal,mort} \times N_{thal,US,prevent} \\ = 40\% \times 900 \\ = 360 \end{gathered} \] ? Low confidence
Sensitivity Analysis
Sensitivity Indices for Maximum DALYs Saved by FDA Preventing Unsafe Drugs (1962-2024)
Regression-based sensitivity showing which inputs explain the most variance in the output.
| Input Parameter | Sensitivity Coefficient | Interpretation |
|---|---|---|
| Thalidomide DALYs Per Event (DALYs) | 1.0000 | Strong driver |
Interpretation: Standardized coefficients show the change in output (in SD units) per 1 SD change in input. Values near ±1 indicate strong influence; values exceeding ±1 may occur with correlated inputs.
Monte Carlo Distribution
Simulation Results Summary: Maximum DALYs Saved by FDA Preventing Unsafe Drugs (1962-2024)
| Statistic | Value |
|---|---|
| Baseline (deterministic) | 2.59 million |
| Mean (expected value) | 2.63 million |
| Median (50th percentile) | 2.53 million |
| Standard Deviation | 754 thousand |
| 90% Range (5th-95th percentile) | [1.54 million, 4.16 million] |
The histogram shows the distribution of Maximum DALYs Saved by FDA Preventing Unsafe Drugs (1962-2024) across 10,000 Monte Carlo simulations. The CDF (right) shows the probability of the outcome exceeding any given value, which is useful for risk assessment.
Exceedance Probability
This exceedance probability chart shows the likelihood that Maximum DALYs Saved by FDA Preventing Unsafe Drugs (1962-2024) will exceed any given threshold. Higher curves indicate more favorable outcomes with greater certainty.
External Data Sources
Parameters sourced from peer-reviewed publications, institutional databases, and authoritative reports.
Disability Weight for Untreated Chronic Conditions: 0.35 weight
Disability weight for untreated chronic conditions (WHO Global Burden of Disease)
Source:4
Uncertainty Range
Technical: Distribution: Normal (SE: 0.07 weight)
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence • 📊 Peer-reviewed
Average Annual New Drug Approvals Globally: 50 drugs/year
Average annual new drug approvals globally
Source:14
Uncertainty Range
Technical: 95% CI: [45 drugs/year, 60 drugs/year] • Distribution: Lognormal
What this means: This estimate has moderate uncertainty. The true value likely falls between 45 drugs/year and 60 drugs/year (±15%). This represents a reasonable range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
✓ High confidence
dFDA Pragmatic Trial Cost per Patient: $929
dFDA pragmatic trial cost per patient. Uses ADAPTABLE trial ($929) as DELIBERATELY CONSERVATIVE central estimate. Ramsberg & Platt (2018) reviewed 108 embedded pragmatic trials; 64 with cost data had median of only $97/patient - our estimate may overstate costs by 10x. Confidence interval spans meta-analysis median to complex chronic disease trials.
Source:1
Uncertainty Range
Technical: 95% CI: [$97, $3K] • Distribution: Lognormal
What this means: This estimate is highly uncertain. The true value likely falls between $97 and $3K (±156%). This represents a very wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence
Drug Development Cost (1980s): $194M
Drug development cost in 1980s (compounded to approval, 1990 dollars)
Source:22
Uncertainty Range
Technical: 95% CI: [$146M, $242M] • Distribution: Lognormal
What this means: This estimate has moderate uncertainty. The true value likely falls between $146M and $242M (±25%). This represents a reasonable range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
✓ High confidence
Regulatory Delay for Efficacy Testing Post-Safety Verification: 8.2 years
Regulatory delay for efficacy testing (Phase II/III) post-safety verification. Based on BIO 2021 industry survey. Note: This is for drugs that COMPLETE the pipeline - survivor bias means actual delay for any given disease may be longer if candidates fail and must restart.
Source:23
Uncertainty Range
Technical: Distribution: Normal (SE: 2 years)
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence • 📊 Peer-reviewed • Updated 2021
FDA Phase 1 to Approval Timeline: 10.5 years
FDA timeline from Phase 1 start to approval. Derived from BIO 2021 industry survey: Phase 1 (2.3 years) + efficacy lag (8.2 years) = 10.5 years. Consistent with PMC meta-analysis finding 9.1 years median (95% CI: 8.2-10.0).
Source:23
Uncertainty Range
Technical: 95% CI: [6 years, 12 years] • Distribution: Gamma (SE: 2 years)
What this means: There’s significant uncertainty here. The true value likely falls between 6 years and 12 years (±29%). This represents a wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The gamma distribution means values follow a specific statistical pattern.
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
✓ High confidence
Annual Days of Chronic Disease Therapy: 1.28 trillion days
Annual days of therapy for chronic conditions globally (diabetes, CVD, respiratory, cancer). IQVIA reports 1.8 trillion total days of therapy in 2019, with 71% for chronic conditions.
Source:43
Uncertainty Range
Technical: 95% CI: [1 trillion days, 1.5 trillion days] • Distribution: Lognormal
What this means: This estimate has moderate uncertainty. The true value likely falls between 1 trillion days and 1.5 trillion days (±20%). This represents a reasonable range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence
Global Daily Deaths from Disease and Aging: 150 thousand deaths/day
Total global deaths per day from all disease and aging (WHO Global Burden of Disease 2024)
Source:4
Uncertainty Range
Technical: Distribution: Normal (SE: 7.5 thousand deaths/day)
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
✓ High confidence • 📊 Peer-reviewed
Global Life Expectancy (2024): 79 years
Global life expectancy (2024)
Source:4
Uncertainty Range
Technical: Distribution: Normal (SE: 2 years)
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
✓ High confidence • 📊 Peer-reviewed • Updated 2024
Pharma Drug Development Cost (Current System): $2.6B
Average cost to develop one drug in current system
Source:72
Uncertainty Range
Technical: 95% CI: [$1.5B, $4B] • Distribution: Lognormal (SE: $500M)
What this means: There’s significant uncertainty here. The true value likely falls between $1.5B and $4B (±48%). This represents a wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
✓ High confidence • 📊 Peer-reviewed
Annual Life-Years Saved by Pharmaceuticals: 149 million life-years
Annual life-years saved by pharmaceutical innovations globally. Lichtenberg (2019, NBER WP 25483) found that drugs launched after 1981 saved 148.7M life-years in 2013 across 22 countries using 3-way fixed-effects regression (disease-country-year). 95% CI [79.4M, 239.8M] propagated from Table 2 regression standard errors (β₀₋₁₁=-0.031±0.008, β₁₂₊=-0.057±0.013).
Source:74
Uncertainty Range
Technical: 95% CI: [79.4 million life-years, 240 million life-years] • Distribution: Lognormal
What this means: This estimate is highly uncertain. The true value likely falls between 79.4 million life-years and 240 million life-years (±54%). This represents a very wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence
Pharma Drug Success Rate (Current System): 10%
Percentage of drugs that reach market in current system
Source:76
✓ High confidence • 📊 Peer-reviewed
Phase I Safety Trial Duration: 2.3 years
Post-1962 Drug Approval Reduction: 70%
Reduction in new drug approvals after 1962 Kefauver-Harris Amendment (70% drop from 43→17 drugs/year)
Source:82
✓ High confidence • Updated 1962-1970
Pre-1962 Physician Count (Unverified): 144 thousand physicians
Estimated physicians conducting real-world efficacy trials pre-1962 (unverified estimate)
Source:84
? Low confidence
Total Number of Rare Diseases Globally: 7 thousand diseases
Total number of rare diseases globally
Source:85
Uncertainty Range
Technical: 95% CI: [6 thousand diseases, 10 thousand diseases] • Distribution: Normal
What this means: There’s significant uncertainty here. The true value likely falls between 6 thousand diseases and 10 thousand diseases (±29%). This represents a wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The normal distribution means values cluster around the center with equal chances of being higher or lower.
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
✓ High confidence
Mean Age of Preventable Death from Post-Safety Efficacy Delay: 62 years
Mean age of preventable death from post-safety efficacy testing regulatory delay (Phase 2-4)
Source:4
Uncertainty Range
Technical: Distribution: Normal (SE: 3 years)
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence • 📊 Peer-reviewed
Pre-Death Suffering Period During Post-Safety Efficacy Delay: 6 years
Pre-death suffering period during post-safety efficacy testing delay (average years lived with untreated condition while awaiting Phase 2-4 completion)
Source:4
Uncertainty Range
Technical: 95% CI: [4 years, 9 years] • Distribution: Lognormal
What this means: There’s significant uncertainty here. The true value likely falls between 4 years and 9 years (±42%). This represents a wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence • 📊 Peer-reviewed
September 11 Deaths: 2.98 thousand people
Total deaths in the September 11, 2001 attacks. 2,977 victims (excluding 19 hijackers). Used as a reference point for scale comparisons.
Source:88
Uncertainty Range
Technical: Distribution: Fixed
✓ High confidence
Standard Economic Value per QALY: $150K
Standard economic value per QALY
Source:94
Uncertainty Range
Technical: Distribution: Normal (SE: $30K)
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
✓ High confidence
Thalidomide Cases Worldwide: 15 thousand cases
Total thalidomide birth defect cases worldwide (1957-1962)
Source:101
Uncertainty Range
Technical: 95% CI: [10 thousand cases, 20 thousand cases] • Distribution: Lognormal
What this means: There’s significant uncertainty here. The true value likely falls between 10 thousand cases and 20 thousand cases (±33%). This represents a wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence
Thalidomide Disability Weight: 0.4:1
Disability weight for thalidomide survivors (limb deformities, organ damage)
Source:102
Uncertainty Range
Technical: 95% CI: [0.32:1, 0.48:1] • Distribution: Lognormal
What this means: This estimate has moderate uncertainty. The true value likely falls between 0.32:1 and 0.48:1 (±20%). This represents a reasonable range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence
Thalidomide Mortality Rate: 40%
Mortality rate for thalidomide-affected infants (died within first year)
Source:101
Uncertainty Range
Technical: 95% CI: [35%, 45%] • Distribution: Lognormal
What this means: This estimate has moderate uncertainty. The true value likely falls between 35% and 45% (±13%). This represents a reasonable range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
✓ High confidence
Thalidomide Survivor Lifespan: 60 years
Average lifespan for thalidomide survivors
Source:102
Uncertainty Range
Technical: 95% CI: [50 years, 70 years] • Distribution: Lognormal
What this means: This estimate has moderate uncertainty. The true value likely falls between 50 years and 70 years (±17%). This represents a reasonable range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence
Phase 3 Cost per Patient: $41K
Phase 3 cost per patient (median from FDA study)
Source:104
Uncertainty Range
Technical: 95% CI: [$20K, $120K] • Distribution: Lognormal
What this means: This estimate is highly uncertain. The true value likely falls between $20K and $120K (±122%). This represents a very wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The lognormal distribution means values can’t go negative and have a longer tail toward higher values (common for costs and populations).
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
✓ High confidence
Treatment Disability Reduction: 0.25 weight
Average disability weight reduction from pharmaceutical treatment. Untreated chronic disease averages 0.35 disability weight, treated disease averages 0.10, difference is 0.25.
Source:105
Uncertainty Range
Technical: 95% CI: [0.15 weight, 0.35 weight] • Distribution: Normal
What this means: There’s significant uncertainty here. The true value likely falls between 0.15 weight and 0.35 weight (±40%). This represents a wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The normal distribution means values cluster around the center with equal chances of being higher or lower.
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
~ Medium confidence • 📊 Peer-reviewed
Core Definitions
Fundamental parameters and constants used throughout the analysis.
Average Life Extension per Beneficiary: 12 years
Average years of life extension per person saved by pharmaceutical interventions. Assumption used to convert life-years saved to approximate lives saved. Based on Lichtenberg’s methodology where life-years are calculated from Years of Life Lost (YLL) reductions.
Uncertainty Range
Technical: 95% CI: [8 years, 18 years] • Distribution: Triangular
What this means: There’s significant uncertainty here. The true value likely falls between 8 years and 18 years (±42%). This represents a wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The triangular distribution means values cluster around a most-likely point but can range higher or lower.
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
Core definition
Eventually Avoidable Death Percentage: 92.6%
Percentage of deaths that are eventually avoidable with sufficient biomedical research and technological advancement. Central estimate ~92% based on ~7.9% fundamentally unavoidable (primarily accidents). Wide uncertainty reflects debate over: (1) aging as addressable vs. fundamental, (2) asymptotic difficulty of last diseases, (3) multifactorial disease complexity.
Uncertainty Range
Technical: 95% CI: [50%, 98%] • Distribution: Beta
What this means: There’s significant uncertainty here. The true value likely falls between 50% and 98% (±26%). This represents a wide range that our Monte Carlo simulations account for when calculating overall uncertainty in the results.
The beta distribution means values are bounded and can skew toward one end.
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
Core definition
Standard Discount Rate for NPV Analysis: 3%
Standard discount rate for NPV analysis (3% annual, social discount rate)
Uncertainty Range
Technical: Distribution: Fixed
Core definition
Pharma Phase 2/3 Cost Barrier Per Drug: $1.56B
Average Phase 2/3 efficacy testing cost per drug that pharma must fund (~60% of total drug development cost)
Uncertainty Range
Technical: Distribution: Normal (SE: $200M)
Input Distribution
This chart shows the assumed probability distribution for this parameter. The shaded region represents the 95% confidence interval where we expect the true value to fall.
Core definition





























































































