import SimPy.EconEvalClasses as EconEval import numpy as np np.random.seed(573) cost_base = np.random.normal(loc=10000, scale=100, size=1000) effect_base = np.random.normal(loc=2, scale=.1, size=1000) cost_intervention = np.random.normal(loc=20000, scale=200, size=1000) effect_intervention = np.random.normal(loc=1, scale=.2, size=1000) print('') # ICER calculation assuming paired observations ICER_paired = EconEval.ICER_paired('Testing paired ICER', cost_intervention, effect_intervention, cost_base, effect_base, EconEval.HealthMeasure.DISUTILITY) print('Paired ICER (confidence and prediction interval): ', ICER_paired.get_ICER(), ICER_paired.get_CI(0.05, 1000), ICER_paired.get_PI(0.05, )) # ICER calculation assuming independent observations ICER_indp = EconEval.ICER_indp('Testing independent ICER', cost_intervention, effect_intervention, cost_base, effect_base, EconEval.HealthMeasure.DISUTILITY) print('Independent ICER (confidence and prediction interval): ', ICER_indp.get_ICER(), ICER_indp.get_CI(0.05, 1000), ICER_indp.get_PI(0.05, )) # try NMB NMB_paired = EconEval.NMB_paired("Testing paired NMB", cost_intervention, effect_intervention, cost_base, effect_base,
import SimPy.EconEvalClasses as EconEval import numpy as np np.random.seed(573) cost_base = np.random.normal(loc=10000, scale=100, size=1000) effect_base = np.random.normal(loc=1, scale=.1, size=1000) cost_intervention = np.random.normal(loc=20000, scale=200, size=1000) effect_intervention = np.random.normal(loc=2, scale=.2, size=1000) print('') # ICER calculation assuming paired observations ICER_paired = EconEval.ICER_paired('Testing paired ICER', cost_intervention, effect_intervention, cost_base, effect_base) print('Paired ICER (confidence and prediction interval): ', ICER_paired.get_ICER(), ICER_paired.get_CI(0.05, 1000), ICER_paired.get_PI(0.05, )) # ICER calculation assuming independent observations ICER_indp = EconEval.ICER_indp('Testing independent ICER', cost_intervention, effect_intervention, cost_base, effect_base) print('Independent ICER (confidence and prediction interval): ', ICER_indp.get_ICER(), ICER_indp.get_CI(0.05, 1000), ICER_indp.get_PI(0.05, )) # try NMB NMB_paired = EconEval.NMB_paired("Testing paired NMB", cost_intervention, effect_intervention, cost_base, effect_base)