powered by
Plot seroprevalence estimates on top of the serosurvey
plot_seroprev_estimates( seromodel, serosurvey, alpha = 0.05, size_text = 11, bin_serosurvey = FALSE, bin_step = 5 )
ggplot object with seroprevalence estimates and serosurveys plots
stan_fit object obtained from sampling a model with fit_seromodel
survey_year
Year in which the survey took place (only needed to plot time models)
age_min
Floor value of the average between age_min and age_max
age_max
The size of the sample
n_sample
Number of samples for each age group
n_seropositive
Number of positive samples for each age group
1 - alpha indicates the credibility level to be used
Size of text for plotting (base_size in ggplot2)
base_size
If TRUE, serodata is binned by means of prepare_bin_serosurvey. Otherwise, age groups are kept as originally input.
TRUE
serodata
prepare_bin_serosurvey
Integer specifying the age groups bin size to be used when bin_serosurvey is set to TRUE.
bin_serosurvey
data(veev2012) seromodel <- fit_seromodel(veev2012, iter = 100) plot_seroprev_estimates(seromodel, veev2012)
Run the code above in your browser using DataLab