adaptiveSIS_PseudoMarginalMCMC: Adaptive Pseudo-Marginal MCMC for SIS Epidemic Panel Data

Description Usage Arguments Value See Also

View source: R/AdaptiveSIS_PM-MCMC.R

Description

Adapts proposal parameters for with a view of optimal of a target using Pseudo-Marginal MCMC scheme.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
adaptiveSIS_PseudoMarginalMCMC(
  obsTransData,
  I_0,
  obsTimes,
  N,
  beta0,
  gamma0,
  lambda0 = 2.38/sqrt(2),
  V0 = diag(c(1/N, 1)),
  noSims,
  noIts,
  burnIn,
  lagMax = NA,
  thinningFactor = 1,
  parallel = FALSE,
  noCores,
  delta = 0.05
)

Arguments

obsTransData

Interpanel transition data.

I_0

Initial number of infectives in the population.

obsTimes

Times at which epidemic cohort were followed up.

N

Population size.

beta0

Starting value for infectious process parameter.

gamma0

Starting value for removal/recovery process parameter.

lambda0

Starting value for RWM proposal parameter which is to be adapted.

V0

Starting state for RWM proposal Covariance matrix which is to be adapted.

noSims

Number of epidemic simulations to be used per iteration for estimation of likelihood.

noIts

Number of MCMC iterations.

lagMax

Plotting parameter for acf() function.

thinningFactor

Controls the factor by which MCMC samples are thinned, to reduce dependency.

parallel

Are epidemic simulations run in parallel?

noCores

If epidemic simulations are run in parallel, this is the number of cores utilised.

delta

Probability that a proposal is made based on the starting proposal conditions.

Value

Proposal parameters which can be used for more optimal exploration of target distribution (plus MCMC summary). Create Storage Matrix Proposal Acceptance Counter Propose new beta and gamma using Multiplicative RW propsal Store State Calculating Summary Statistics for samples

See Also

Other Panel Data MCMC: SIR_fsMCMC_blockedIS(), adaptiveSIR_PseudoMarginalMCMC(), adaptiveSISfsMCMC()


JMacDonaldPhD/Epidemics documentation built on Jan. 10, 2020, 2:48 a.m.