skill-pharmacoeconomic-evaluation（技能药物经济学评价）

Pharmacoeconomic Evaluation Skill

Overview

This skill provides comprehensive guidance for conducting pharmacoeconomic evaluations, including cost-effectiveness analysis, cost-utility analysis, cost-benefit analysis, sensitivity analysis, and model construction. Following the Chinese Pharmacoeconomic Evaluation Guidelines (up-to-date Edition), it provides complete workflows, calculation tools, and reference materials for economic evaluation of healthcare interventions.

Evaluation Types

Choose the appropriate evaluation type based on research objectives and data characteristics:

• Cost-Effectiveness Analysis (CEA): Use when intervention effect can be measured by a single clinical indicator (e.g., life years, survival rate)
• Cost-Utility Analysis (CUA): Use when both quality and quantity of life need to be considered; outcome measure is QALYs
• Cost-Benefit Analysis (CBA): Use when intervention effects can be expressed in monetary terms
• Cost-Minimization Analysis (CMA): Use when two interventions have proven equivalent efficacy; only compare costs

Core Workflow

Step 1: Define Research Framework

1. Define research question

- Identify target disease and population - Determine intervention and comparator - Set perspective (recommended: societal)

1. Select evaluation type

- Choose CEA, CUA, CBA, or CMA based on outcome measure - Consider discounting for long-term studies (both costs and outcomes) - Recommended discount rate: 3.5%

1. Determine time horizon

- Chronic diseases: lifetime or sufficiently long (>95% patients have reached the death situation) - Acute diseases: short-term follow-up (1-3 years)

Step 2: Identify and Measure Costs

Identify costs following Chinese Pharmacoeconomic Evaluation Guidelines:

Direct Medical Costs

• Medication costs
• Outpatient costs
• Inpatient costs
• Diagnostic and test costs
• Surgical treatment costs
• Adverse event treatment costs

Direct Non-Medical Costs

• Transportation
• Accommodation
• Nutritional support
• Unprofessional caregiving

Indirect Costs

• Productivity loss (premature death or sick leave)
• Caregiver burden

Intangible Costs

• Pain, anxiety, quality of life reduction not included in monetary costs; considered in utility analysis

Cost Data Sources:

• Hospital Information Systems
• Insurance Databases
• Epidemiological Studies
• Literature Review
• Questionnaire Surveys

Step 3: Measure Effects/Utilities

Effect Measure Selection

• Survival indicators: Life Years (LY), Survival Rate
• Disease-specific indicators: Event-free survival, Symptom improvement
• Others: Complication rate, Hospitalization frequency

Utility Measurement (Recommended Indirect Methods)

• EQ-5D (EuroQol Five-Dimensional Questionnaire)
• SF-6D (Based on SF-36)
• QWB (Quality of Well-Being Index)

Utility Value Source Priority:

1. Primary data from target population (best)
2. Published Chinese population utility values
3. Data from other countries (requires adjustment)

Step 4: Build Decision Analytic Models

Select appropriate model type based on research characteristics:

Decision Tree Model

• Scenarios: Short-term, single decision, clear event sequence
• Advantages: Intuitive, easy to understand, suitable for analyzing decision processes
• Steps:

1. Define decision nodes, chance nodes, terminal nodes 2. Assign probabilities to each chance node (sum to 1) 3. Assign costs and effects to each terminal node 4. Roll back to calculate expected values 5. Compare decision options

Markov Model

• Usage scenario: Chronic diseases, long-term follow-up, recurrent events
• Advantages: Can handle cyclical state transitions, clear structure
• Steps:

1. Define health states (e.g., healthy, mild, moderate, severe, death) 2. Build transition matrix (describe state-to-state transition probabilities) 3. Estimate transition probabilities (from incidence, survival curves, or literature) 4. Assign cycle costs and utilities to each state 5. Set cycle length (typically 1 year) and model time horizon 6. Run Markov simulation

Discrete Event Simulation (DES)

• Scenarios: Large individual variation, irregular event timing, resource constraints
• Advantages: Most flexible, can simulate individual paths, precise time-dependent modeling
• 步骤:

1. Define entities (patients) and their attributes 2. Define possible event types 3. Establish event scheduling mechanism 4. Run simulation 5. Aggregate results

Partitioned Survival Model (PSM)

• Scenarios: Oncology research, based on survival curves
• Advantages: Directly based on survival data, reasonable extrapolation
• Steps:

1. Obtain PFS and OS survival curves 2. Fit parametric distributions (exponential, Weibull, etc.) 3. Extrapolate to model time horizon 4. Calculate population distribution across partitions 5. Accumulate costs and utilities

See references/model_methods.md for detailed modeling methods.

Step 5: Calculate Key Metrics

Use calculation tools in scripts/:

Incremental Cost-Effectiveness Ratio (ICER)

Use calculate_icere() from scripts/cost_effectiveness_analysis.py:

from cost_effectiveness_analysis import calculate_icere

result = calculate_icere(
    cost_intervention,  # Intervention group cost
    effect_intervention,  # Intervention group effect (e.g., QALYs)
    cost_control,  # Control group cost
    effect_control,  # Control group effect
    threshold=30000  # Threshold (30KUSD for US & UK, and close to 2x GDP per QALY of China)
)

ICER Formula: \[ ICER = \frac{C_A - C_B}{E_A - E_B} = \frac{\Delta C}{\Delta E} \]

Quality-Adjusted Life Years (QALYs)

Use calculate_qaly() from scripts/cost_effectiveness_analysis.py:

from cost_effectiveness_analysis import calculate_qaly

qalys = calculate_qaly(
    life_years=10,  # Life years
    utility_scores=np.array([...]),  # Utility scores for each period
    discount_rate=0.03  # Discount rate
)

QALY Formula: \[ QALY = \sum_{t=1}^{T} U_t \ imes \frac{1}{(1+r)^{t-1}} \]

Net Benefit

\[ NB = \lambda \ imes E - C \]

Where:

• NB = Net Benefit
• λ = Willingness-to-pay threshold
• E = Effect
• C = Cost

Step 6: Conduct Sensitivity Analysis

One-Way Sensitivity Analysis

Use deterministic_sensitivity_analysis() from scripts/cost_effectiveness_analysis.py:

from cost_effectiveness_analysis import deterministic_sensitivity_analysis

# Define parameter ranges
param_ranges = {
    'drug_cost': (10000, 20000),
    'hospital_cost': (5000, 15000),
    'effectiveness': (0.8, 1.2)
}

# Run sensitivity analysis
results_df = deterministic_sensitivity_analysis(
    base_params=base_parameters,
    param_ranges=param_ranges,
    outcome_func=outcome_function
)

Tornado Plot Data: Use tornado_plot_data() function

Probabilistic Sensitivity Analysis (PSA)

Use MonteCarloSimulator from scripts/monte_carlo_simulation.py:

from monte_carlo_simulation import MonteCarloSimulator

# Create simulator
simulator = MonteCarloSimulator(n_simulations=10000, seed=42)

# Define parameter distributions
parameters = {
    'cost': {
        'distribution': 'gamma',
        'params': (2, 15000),  # shape, scale
        'min_value': 0
    },
    'effect': {
        'distribution': 'beta',
        'params': (5, 3),  # alpha, beta
        'min_value': 0,
        'max_value': 10
    }
}

# Run PSA
results_df = simulator.probabilistic_sensitivity_analysis(
    parameters=parameters,
    outcome_func=outcome_function,
    threshold=120000
)

Generate CEAC: Use generate_ceac() function

Value of Information (VOI): Use value_of_information_analysis() function

Step 7: Interpret and Report Results

Willingness-to-Pay Threshold (Reference)

• UK: 25000~30000 GBP
• US: 50000~100000 USD

Interpretation

• ICER ≤ Threshold: Cost-effective
• ICER > Threshold: Not cost-effective
• Strict Dominance: Lower cost and better effect
• Strict Disadvantage: Higher cost and worse effect

Reporting Requirements

Follow CHEERS 2022 and Chinese Pharmacoeconomic Evaluation Guidelines:

1. Clearly describe research design and methods
2. Report baseline analysis results
3. Provide sensitivity analysis results (one-way and probabilistic)
4. Report confidence intervals
5. Discuss limitations and generalizability
6. Clearly state funding sources and potential conflicts of interest

See references/guidelines.md for detailed guidelines.

Scripts Guide

cost_effectiveness_analysis.py

Core functions: Cost-effectiveness analysis, ICER calculation, QALY calculation, deterministic sensitivity analysis

Main functions:

• calculate_icere(): Calculate ICER
• calculate_qaly(): Calculate QALYs
• calculate_ceac(): Calculate Cost-Effectiveness Acceptability Curve
• deterministic_sensitivity_analysis(): One-way sensitivity analysis
• tornado_plot_data(): Prepare tornado plot data
• markov_model_transition(): Markov model simulation
• discount_costs(): Cost discounting

monte_carlo_simulation.py

Core functions: Monte Carlo simulation, probabilistic sensitivity analysis, value of information analysis

Main classes and methods:

• MonteCarloSimulator: Monte Carlo simulator

- generate_samples(): Generate samples from specified distribution - probabilistic_sensitivity_analysis(): Run PSA - generate_ceac(): Generate CEAC - value_of_information_analysis(): VOI analysis - scatter_plot_data(): Prepare cost-effectiveness scatter plot data

References Guide

guidelines.md

Summary of key content from ISPOR Good Practices, including:

• Evaluation framework and perspective
• Cost identification and measurement
• Effect/utility measurement
• Model construction methods
• Discounting principles
• Sensitivity analysis requirements
• Result presentation and reporting standards
• Common calculation formulas

Use case: Query specific requirements, standards, and methods for Chinese pharmacoeconomic evaluation

model_methods.md

Detailed decision analytic model construction methods, including:

• Markov model (basic concepts, transition matrix, probability estimation)
• Decision tree model (structure, probability assignment, rollback calculation)
• Discrete event simulation (core elements, advantages/disadvantages)
• Partitioned survival model (survival curve fitting)
• Model comparison and selection
• Modeling best practices

Use case: Learn specific modeling methods, build decision analytic models

Common Task Scenarios

Scenario 1: Conduct Cost-Effectiveness Analysis for New Drug

1. Determine research perspective (societal)
2. Identify direct medical and non-medical costs
3. Collect clinical trial data for effect measures (survival, QALYs)
4. Build Markov model to simulate disease progression
5. Calculate ICER and compare with threshold
6. Conduct one-way and probabilistic sensitivity analysis
7. Write report following CHEERS standards

Scenario 2: Model Building and Validation

1. Select model type based on disease characteristics
2. Learn modeling methods from references/model_methods.md
3. Estimate model parameters from literature or clinical trials
4. Validate model (internal and external validation)
5. Run baseline analysis
6. Conduct sensitivity analysis to verify model stability

Scenario 3: Probabilistic Sensitivity Analysis

1. Specify probability distributions for each key parameter
2. Run 10,000+ simulations using MonteCarloSimulator
3. Generate cost-effectiveness scatter plot
4. Generate Cost-Effectiveness Acceptability Curve (CEAC)
5. Conduct Value of Information (VOI) analysis
6. Report cost-effectiveness probability and confidence intervals

Parameter Management Best Practices

Parameter Organization

Organize parameters by category:

• Research Framework Parameters: Perspective, time horizon, discount rate, threshold
• Model Structure Parameters: Health states, initial distribution
• Transition Probability Parameters: State-to-state transition probabilities
• Cost Parameters: Annual costs by state
• Utility Parameters: Utility values by state
• Sensitivity Analysis Parameters: Parameter ranges and probability distributions
• Simulation Parameters: Number of simulations, random seed, etc.

Parameter Source Documentation

Each parameter value must have a clear data source:

• Literature Citation: Author, journal, year, pages
• Database: Database name, version, access date
• Guidelines/Standards: Guideline name, version, issuing organization
• Expert Opinion: Expert source and judgment basis
• Research Assumption: Rationale for assumption

Example Code Format

# ========== Parameter Category Title ==========
PARAMETER_NAME = {
    'parameter_key': value,  # Source: Detailed source description
    'another_key': value,    # Source: Reference [Author, Journal, Year]
}

See scripts/example.py for complete parameter organization format.

Important Notes

1. Follow Chinese Guidelines: Ensure research methods meet requirements of Chinese Pharmacoeconomic Evaluation Guidelines (2023)

1. Transparency: Clearly describe all assumptions, data sources, and calculation methods

1. Parameter Source Documentation: All parameter values must cite sources for traceability and verification

1. Discounting: Both costs and outcomes need discounting; recommended rate is 3.5%

1. Sensitivity Analysis: Conduct sufficient sensitivity analysis to evaluate uncertainty

1. Model Validation: Validate model internally; conduct external validation if possible

1. Reporting Standards: Follow CHEERS 2022 reporting standards

1. Threshold: Clearly state the threshold used and its basis (Reference: 1-3x GDP/QALY)

1. Time Horizon: Select sufficiently long time horizon to capture all relevant costs and outcomes

1. Cost Measurement: Avoid using payment prices (reimbursed prices); use actual costs or standardized charges

1. Utility Measurement: Prioritize Chinese population utility values; note applicability of measurement tools

1. Parameter Organization: Reference format in scripts/example.py, organize parameters neatly and document sources in detail