plot_posterior_density: Plot posterior density

In jameshay218/virosolver: Infer epidemic dynamics using viral load data

Plot posterior density

Description

Plot posterior density from model output

Usage

plot_posterior_density(
  chain,
  var_name,
  parTab,
  prior_mean,
  prior_sd,
  prior_dist = "normal",
  real_data = FALSE
)

Arguments

chain
var_name	Character string denoting variable name from parameter table.
parTab	Dataframe specifying prior distribution for model parameter.
prior_mean	Mean of the prior distribution for a particular parameter.
prior_sd	Standard deviation of the prior distribution for a particular parameter.
real_data	Plot true parameter values to simulate the data. Defaults to FALSE.

Value

ggplot

Author(s)

James Hay, [e-mail]

Examples

data(example_seir_partab)
chains <- load_mcmc_chains(location="mcmc_chains/readme_multiple_cross_section/",
parTab=par_tab,
burnin=mcmc_pars["adaptive_period"],
chainNo=FALSE,
unfixed=FALSE,
multi=FALSE)

var_name <- "tshift"
prior_mean <- subset(example_seir_partab, names == "tshift")$values 
prior_sd <- 0.2 
plot_posterior_density <- function(chain, var_name, parTab, prior_mean, prior_sd, real_data=FALSE)

jameshay218/virosolver documentation built on April 17, 2025, 2:57 p.m.