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))

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))
# 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{}'.
      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{}'.