Commit 74af5547 authored by Peter Jentsch's avatar Peter Jentsch
Browse files

notes from meeting

parent db80d49b
#here are the functions for
#do this only at the start
"""
Called in the inner loop of `generate_contact_vectors!`, this function removes the degree corresponding to kth entry of index_lists, replaces it with another randomly sampled degree, and returns `csum` adjusted to the new samples.
......
......@@ -91,6 +91,9 @@ function time_dep_mixing_graphs(len,base_network,demographics,index_vectors,ws_m
rest_weekly_edges = WeightedGraph(demographics,index_vectors,rest_matrix_tuple.twice_a_week,contact_time_distributions.rest)
rest_justonce_edges = WeightedGraph(demographics,index_vectors,rest_matrix_tuple.otherwise,contact_time_distributions.rest)
#weekly multiply durations by 1/5
#daily multiply by 1/2
inf_network_list = [home_static_edges,rest_static_edges,ws_justonce_edges,rest_justonce_edges]
soc_network_list = [home_static_edges,rest_static_edges,ws_static_edges]
......
......@@ -21,7 +21,7 @@
κ = 0.0,
ω = 0.0061,
ω_en = 0.05,
Γ = 1/7,
Γ = 1/7, #make this .906
ξ = 5.0,
notification_parameter = 0.0005,
vaccinator_prob = 0.6,
......
......@@ -4,7 +4,7 @@
using Printf
const ts_colors = cgrad(:PuBu_9)
#plot for each day some metric of repeated notifications
function plot_model(varname,univariate_series, output_list::Vector{T},p) where T
@unpack sim_length, immunization_intervals,infection_introduction_day = p
immunization_begin = first(minimum(immunization_intervals))
......
Markdown is supported
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