Commit 0b817375 authored by Peter Jentsch's avatar Peter Jentsch
Browse files

count assortativity at pre-season

parent c1fec531
......@@ -140,10 +140,10 @@ function plot_parameter_plane(input_fname)
"wdg_change_O.pdf",
"final_size_change.pdf",
"final_size_standard_dev.pdf",
"assortativity.pdf",
"join_11_counts.pdf",
"join_01_counts.pdf",
"join_00_counts.pdf",
"assortativity_preseason_count.pdf",
"join_11_counts_preseason_count.pdf",
"join_01_counts_preseason_count.pdf",
"join_00_counts_preseason_count.pdf",
]
titles = [
"Change in w. deg. of vaccinators with EN",
......
......@@ -199,7 +199,7 @@ function agents_step!(t,modelsol)
else
π_base_vec = @SVector [π_base_y*ζ,π_base_m*ζ,π_base_o*ζ]
end
if t == modelsol.sim_length
if t == last(minimum(immunization_intervals)) + 14
weighted_assortativity, join_counts = assortativity(t,soc_network,u_inf)
fit!(modelsol.output_data.assortativity_at_end,weighted_assortativity)
fit!(modelsol.output_data.vac_join_11_counts_at_end,join_counts[1])
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment