costs_new=cost_intervention,
effects_new=effect_intervention,
costs_base=cost_base,
effects_base=effect_base)
print('Paired ICER:'
'\n\tICER: {} '
'\n\tCI (boostrap): {} '
'\n\tCI (Bayesian): {} '
'\n\tPI: {}'.format(ICER_paired.get_ICER(), ICER_paired.get_CI(0.05),
ICER_paired.get_CI(0.05, method='Bayesian'),
ICER_paired.get_PI(0.05)))
# ICER calculation assuming independent observations
ICER_indp = EconEval.ICER_Indp('Testing independent ICER',
costs_new=cost_intervention,
effects_new=effect_intervention,
costs_base=cost_base,
effects_base=effect_base)
print('Independent ICER (confidence and prediction interval): ',
ICER_indp.get_ICER(), ICER_indp.get_CI(0.05), ICER_indp.get_PI(0.05, ))
# NMB
NMB_paired = EconEval.INMB_Paired("Testing paired NMB",
costs_new=cost_intervention,
effects_new=effect_intervention,
costs_base=cost_base,
effects_base=effect_base)
print('Paired NMB (confidence and prediction interval): ',
NMB_paired.get_INMB(wtp=10000), NMB_paired.get_CI(wtp=10000, alpha=.05),
NMB_paired.get_PI(wtp=10000, alpha=.05))
ICER_paired = EconEval.ICER_Paired('Testing paired ICER',
cost_intervention,
effect_intervention,
cost_base,
effect_base,
health_measure='d')
print(
'Paired ICER:\n\tICER: {} \n\tCI (boostrap): {} \n\tCI (Bayesian): {} \n\tPI: '
.format(ICER_paired.get_ICER(), ICER_paired.get_CI(0.05, 1000),
ICER_paired.get_CI(0.05, 1000, method='Bayesian'),
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,
health_measure='d')
print(
'Independent ICER (confidence and prediction interval): \n\t{}\n\t{}\n\t{}'
.format(ICER_indp.get_ICER(), ICER_indp.get_CI(0.05, 1000),
ICER_indp.get_PI(0.05)))
# try NMB
NMB_paired = EconEval.INMB_Paired("Testing paired NMB",
cost_intervention,
effect_intervention,
cost_base,
effect_base,
health_measure='d')
print('Paired NMB (confidence and prediction interval): \n\t{}\n\t{}\n\t{}'.