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,
EconEval.HealthMeasure.DISUTILITY)
print('Paired NMB (confidence and prediction interval): ',
NMB_paired.get_NMB(wtp=10000), NMB_paired.get_CI(wtp=10000, alpha=.05),
NMB_paired.get_PI(wtp=10000, alpha=.05))
NMB_indp = EconEval.NMB_indp("Testing independent NMB", cost_intervention,
effect_intervention, cost_base, effect_base,
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)
print('Paired NMB (confidence and prediction interval): ',
NMB_paired.get_NMB(wtp=10000),
NMB_paired.get_CI(wtp=10000, alpha=.05),
NMB_paired.get_PI(wtp=10000, alpha=.05))
NMB_indp = EconEval.NMB_indp("Testing independent NMB",
cost_intervention, effect_intervention, cost_base, effect_base)
