print(path1.stat.get_mean())
# second sample path with initial size = 1
path2 = Path.IncidenceSamplePath('Path 2', delta_t=1, sim_rep=1)
# record the observations
path2.record_increment(0.5, 4)
path2.record_increment(1.8, 2)
path2.record_increment(5.5, 1)
path2.close(6)
# stats
print(path2.stat.get_mean())
# plot path 1 only
Path.plot_sample_path(
sample_path=path1,
title='Plotting a single sample path',
x_label='Time Period',
y_label='observed value',
legend='Path 1',
color_code='r',
connect='line')
Path.plot_sample_path(
sample_path=path2,
title='Plotting a single sample path',
x_label='Time Period',
y_label='observed value',
legend='Path 2',
color_code='r',
connect='line')
import SimPy.Plots.Histogram as Hist
import SimPy.Plots.SamplePaths as Path
# selected therapy
therapy = P.Therapies.COMBO
# create a cohort
myCohort = Cls.Cohort(id=1,
pop_size=D.POP_SIZE,
parameters=P.Parameters(therapy=therapy))
# simulate the cohort over the specified time steps
myCohort.simulate(sim_length=D.SIM_LENGTH)
# plot the sample path (survival curve)
Path.plot_sample_path(sample_path=myCohort.cohortOutcomes.nLivingPatients,
title='Survival Curve',
x_label='Time-Step (Year)',
y_label='Number Survived')
# plot the histogram of survival times
Hist.plot_histogram(data=myCohort.cohortOutcomes.survivalTimes,
title='Histogram of Patient Survival Time',
x_label='Survival Time (Year)',
y_label='Count',
bin_width=1)
# print the outcomes of this simulated cohort
Support.print_outcomes(sim_outcomes=myCohort.cohortOutcomes,
therapy_name=therapy)
increments=[2, -1, 0],
sim_rep=1)
path4 = Path.PrevalencePathBatchUpdate('Path 4',
0,
times_of_changes=[0.5, 4, 4.5],
increments=[2.5, 1, 1],
sim_rep=1)
# stats
print(path1.stat.get_mean())
# plot path 1 only
Path.plot_sample_path(sample_path=path1,
title='Plotting a single sample path',
x_label='time',
y_label='observed value',
legend='Path 1',
color_code='r',
connect='step')
# plot path 3 only
Path.plot_sample_path(
sample_path=path3,
title='Plotting a single sample path that is updated in batch',
x_label='time',
y_label='observed value',
legend='Path 3',
color_code='r')
# plot both paths
Path.plot_sample_paths(sample_paths=[path1, path2],
