import SimParameters as P
from SimPy import EconEvalClasses as ce
import numpy as np
# run simulation of patients
year1 = P.YearofPatients(id=1)
simulation = year1.simulate()
alpha = 0.05
# summary statistics for all conditions
print('Summary Stats OS Cost', simulation.get_sumStat_OS_cost().get_mean())
print('95% Confidence Interval for OS Cost',
simulation.get_sumStat_OS_cost().get_t_CI(alpha))
print('Summary Stats No OS Cost',
simulation.get_sumStat_NoOS_cost().get_mean())
print('95% Confidence Interval for No OS Cost',
simulation.get_sumStat_NoOS_cost().get_t_CI(alpha))
print('Summary Stats OS Utility',
simulation.get_sumStat_OS_utility().get_mean())
print('95% Confidence Interval for OS Utility',
simulation.get_sumStat_OS_utility().get_t_CI(alpha))
print('Summary Stats No OS Utility',
simulation.get_sumStat_NoOS_utility().get_mean())
print('95% Confidence Interval for No OS Utility',
simulation.get_sumStat_NoOS_utility().get_t_CI(alpha))
print(" ")
print(" ")
# CEA plot
def __init__(self, params_list, init_r, init_rho):
self.params = SimParameters.Initialization(params_list, init_r,
init_rho)