# 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)) NMB_indp = EconEval.INMB_Indp("Testing independent NMB", costs_new=cost_intervention, effects_new=effect_intervention, costs_base=cost_base, effects_base=effect_base) print('Independent NMB (confidence and prediction interval): ', NMB_indp.get_INMB(wtp=10000), NMB_indp.get_CI(wtp=10000, alpha=.05), NMB_indp.get_PI(wtp=10000, alpha=.05)) print('')
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{}'. format(NMB_paired.get_INMB(wtp=10000), NMB_paired.get_CI(wtp=10000, alpha=.05), NMB_paired.get_PI(wtp=10000, alpha=.05))) NMB_indp = EconEval.INMB_Indp("Testing independent NMB", cost_intervention, effect_intervention, cost_base, effect_base, health_measure='d') print( 'Independent NMB (confidence and prediction interval): \n\t{}\n\t{}\n\t{}'. format(NMB_indp.get_INMB(wtp=10000), NMB_indp.get_CI(wtp=10000, alpha=.05), NMB_indp.get_PI(wtp=10000, alpha=.05))) print('')