expand state vector by age classes and population distribution
Source:R/ageify.R
expand_state_age.Rdepidemiological state varies fast, age category varies slowly
expand_state_age(x, age_cat = mk_agecats(), Nvec = NULL)
ageify(ratemat, params, age_cat = mk_agecats())Arguments
- x
state vector
- age_cat
vector of age categories
- Nvec
population distribution (as counts)
- ratemat
rate matrix
- params
parameter vector
See also
Other classic_macpan:
add_d_log(),
add_updated_vaxrate(),
aggregate_agecats(),
calibrate_comb(),
calibrate(),
check_age_cat_compatibility(),
check_contact_rate_setting(),
col_multiply(),
condense_age(),
condense_params_vax(),
condense_state(),
condense_vax(),
dev_is_tikz(),
do_step(),
expand_params_age(),
expand_params_desc_age(),
expand_params_desc_variant(),
expand_params_desc_vax(),
expand_params_mistry(),
expand_params_variant(),
expand_params_vax(),
expand_state_vax(),
expand_stateval_testing(),
fix_pars(),
fix_stored(),
forecast_ensemble(),
forecast_sim(),
getData(),
get_GI_moments(),
get_Gbar(),
get_R0(),
get_doses_per_day(),
get_evec(),
get_kernel_moments(),
get_opt_pars(),
get_r(),
invlink_trans(),
make_betavec(),
make_beta(),
make_jac(),
make_ratemat(),
make_state(),
make_test_wtsvec(),
make_vaxrate(),
mk_Nvec(),
mk_agecats(),
mk_contact_rate_setting(),
mk_mistry_Nvec(),
mk_pmat(),
mk_vaxcats(),
mle_fun(),
non_expanded_states,
rExp(),
read_params(),
repair_names_age(),
restore(),
run_sim_ageify(),
run_sim_break(),
run_sim_loglin(),
run_sim_mobility(),
run_sim_range(),
run_sim(),
show_ratemat(),
testify(),
texify(),
trans_state_vars(),
update_contact_rate_setting(),
update_foi(),
update_params_mistry(),
vis_model(),
write_params()
Other classic_macpan:
add_d_log(),
add_updated_vaxrate(),
aggregate_agecats(),
calibrate_comb(),
calibrate(),
check_age_cat_compatibility(),
check_contact_rate_setting(),
col_multiply(),
condense_age(),
condense_params_vax(),
condense_state(),
condense_vax(),
dev_is_tikz(),
do_step(),
expand_params_age(),
expand_params_desc_age(),
expand_params_desc_variant(),
expand_params_desc_vax(),
expand_params_mistry(),
expand_params_variant(),
expand_params_vax(),
expand_state_vax(),
expand_stateval_testing(),
fix_pars(),
fix_stored(),
forecast_ensemble(),
forecast_sim(),
getData(),
get_GI_moments(),
get_Gbar(),
get_R0(),
get_doses_per_day(),
get_evec(),
get_kernel_moments(),
get_opt_pars(),
get_r(),
invlink_trans(),
make_betavec(),
make_beta(),
make_jac(),
make_ratemat(),
make_state(),
make_test_wtsvec(),
make_vaxrate(),
mk_Nvec(),
mk_agecats(),
mk_contact_rate_setting(),
mk_mistry_Nvec(),
mk_pmat(),
mk_vaxcats(),
mle_fun(),
non_expanded_states,
rExp(),
read_params(),
repair_names_age(),
restore(),
run_sim_ageify(),
run_sim_break(),
run_sim_loglin(),
run_sim_mobility(),
run_sim_range(),
run_sim(),
show_ratemat(),
testify(),
texify(),
trans_state_vars(),
update_contact_rate_setting(),
update_foi(),
update_params_mistry(),
vis_model(),
write_params()
Examples
params <- read_params("PHAC_testify.csv")
ss <- make_state(params=params)
#> Warning: CHECK: may not be working properly for testify?
# ss2 <- expand_state_age(ss)
params <- read_params("PHAC_testify.csv")
state <- make_state(params=params)
#> Warning: CHECK: may not be working properly for testify?
M <- make_ratemat(state,params, sparse=TRUE)
Ma <- ageify(M, params)
library(Matrix)
#>
#> Attaching package: ‘Matrix’
#> The following object is masked from ‘package:McMasterPandemic’:
#>
#> unpack
Matrix::image(Ma)
Mta <- ageify(testify(M,params),params)
Matrix::image(Mta)
