- The
drjacoby
package can now be used to fit mcmc for deterministic model drjacoby_mcmc
function analogous topmcmc
for fitting modelsconvert_log_likelihood_func_for_drjacoby
helps convert log likelihood functions, e.g.calc_loglikelihood
into suitable likelihood functions that can be correctly passed aslog_likelihood
todrjacoby_mcmc
.convert_log_prior_func_for_drjacoby
does the same for priors- The above changes are because
drjacoby
requires initial parameters and ranges to be passed as adata.frame
internally, which means dates are internally handled as numerics, e.g. "2020-01-01" is 18262. Consequently, the sampled parameters fromdrjacoby_mcmc
will havestart_date
values as numerics. Seesample_drjacoby
function for how these are handled correctly to align with converted log likelihood function draws.
parse_durations
bug fix. Handles dur_E, dur_IMild and dur_ICase correctly now
sample_pmcmc
bug fix dor negative posterior checks incorrectly.
atol
andrtol
can be passed topmcmc
fitting as named arguments inpars_obs
to be used inrun_deterministic_comparison
pmcmc
can now be used to fit to PCR prevalence data (deterministic model only) by passingpcr_df
andpcr_det
as named list elements ofpars_obs
run_deterministic_comparison
patches for sero fitting and argument checks
pmcmc
can now be used to fit to serology data (deterministic model only) by passingsero_df
andsero_det
as names list elements ofpars_obs
projections
can be used now fornimue
models
pmcmc
now accepts dates of change for vaccine efficacy fornimue
pmcmc
changes to adapt fornimue
class models and vaccines
beta_est
acceptsnimue
models
- Typo in
dur_not_get_mv_survive
pmcmc(..., scaling_factor)
allows scaling factors to be passed to pmcmc- Correction to
deterministic_model
model class
- Waning added to squire
- Updated to use odin 1.1.0
odin_is_discrete
new unexported function
- Github actions included for automatic release tagging
pmcmc(..., gibbs_sampling = TRUE)
implements Gibbs sampler
population$n
are all integer class
beta_est
accepts apothecary models- Updated json odin files that weren't compiled correctly
default_durations
contains all durationsdefault_durations
anddefault_probs
now exported for easerun_deterministic_SEIR_model
has odin pars exported now
Parameters
vignette updated to include new parametersrun_deterministic_comparison
bug fix for when transmission is eliminated
projections
bug fixed for setting time varying durations- Taiwan healthcare capacity added
default_durations
updated for incorrect hospital durations
- Parameter Updates. See parameters vignette for update.
triggering
function for doingprojections
with model state triggered R0 changes- Time-varying durations of ICU and hospital stay
projections
allows for different args to be passed to each replicateprojections
allows for additional user model args to be passed
pmcmc
argumenttreated_deaths_only
will use deaths from treatment in likelihood
format_output
acceptsdeaths_treatment
anddeaths_no_treatment
intervention_dates_for_odin
fills dates earlier thandate_R0_change
- Correction to random walk parameters to more accurately continue the walk
evaluate_Rt_pmcmc
takesRt_rw_n
for random walk esque shifts
format_output
backwards compatible for incidence
calibrate
patched asdt
argument was not being passed throughformat_output
corrected for hospital/ICU incidence
evaluate_Rt_pmcmc
patch for non specifiedRt_shift
in 3 parameter version
pmcmc
more flexible for varying which parameters are fitted.evaluate_Rt_pmcmc
altered for increased flexibility
format_output
now returnshospital_incidence
andICU_incidence
evaluate_Rt_pmcmc
patch forMeff_switch_date
earlier thandate_R0_change[1]
Meff_pl
can be scaled or not in relation to Meff
pmcmc
no longer takes Rt_func as hard coded given Meff_pl complexity
- Johnstone Change approach implemented for pmcmc RM
evaluate_Rt
patch forroll
of Meff
evaluate_Rt
allows fordate_Meff_change
to be later/earlier thandate_R0_change
ends
- Large overhaul to make
projections
and plotting work with deterministic models
calc_loglikelihood
fix for reduced date vectors
pars_init
forpmcmc
can be provided as a list for where each chain starts
pmcmc
method for 4 parameter estimation.
projections
now matches by column names when including projection outputs- Patch to
t0_variables
for new Rt_func approach
- replace chained binomials with mutivariate hypergeometric sampling in stochastic model
- stochastic model to output number newly requiring ICU beds at each timepoint
calibrate
now takes anRt_func
argument for specifying Rt calculation
- Explicit Meff inclusion in
R0_date_particle_filter
andprojections
fix
- patching
projections
when projecting further than original calibrate run
calibrate
new argument forR0_prior
calibrate
now allows intervetnion dates to be before epidemic start date
calibrate
argumentreporting_fraction
implemented now in likelihood
calibrate
can do 3D grid scans investigating Meffplot.squire_scan
will collapse n-dimensional scans usingshow
argument
- Patch for
particle_filter_data
when only 1 row of data is provided
- Add functionality for running the deterministic model in javascript
- Patch for
calibrate
to ensure correct dt values provided
calibrate
can now takesquire:::deterministic_model()
as model type.- Output of above does not work with plotting or projections yet.
- Patch for
plot.scan_results
where log likelihood was being plotted for both
- Patch for passing in user defined population and matrices
calibrate
now uses particle filter. See README for examples (#12)- Dates are to be provided for timings of interventions before today (#44, #17)
- Custom populations, contact matrices and bed capacity allowed (#73, #80, #28, #16)
- Increased modularity for new models (#62, #61, #49)
- Returns lieklihood of different baseline R0s and start dates (#14)
get_mixing_matrix
takes iso3c (#84)- Number of plotting enhancement (#72, #69)
- Parameter references vignette and website update (#13, #1, #5, #78)
projection
function for doing scenario projections
format_output
and calibrate orderly plotting functioanlity in
calibrate
now allows different seeding events (#50) and undereporting (#11)- New plotting functions for
plot.squire_calibration
- iso3c country codes in (#52)
extract_all_outputs
andextract_specific_output
(#37)
- New minor version for iterative PR development cycle from now
- initial version of the model