haiti1_joint: Build pomp object for Model 1 - Joint

View source: R/haiti1_joint.R

haiti1_jointR Documentation

Build pomp object for Model 1 - Joint

Description

Generate a class ‘pomp’ object for fitting to epidemic/endemic Haiti cholera data jointly.

Usage

haiti1_joint(
  vacscen = "id0",
  breakpoint = 232,
  rho_flag = FALSE,
  tau_flag = TRUE,
  sig_sq_flag = TRUE,
  beta_flag = FALSE,
  nu_flag = FALSE
)

Arguments

vacscen

ID code for vaccination scenario

breakpoint

the week number at which to start re-estimation

rho_flag

boolean indicating whether to re-estimate \rho at the breakpoint

tau_flag

boolean indicating whether to re-estimate \tau at the breakpoint

sig_sq_flag

boolean indicating whether to re-estimate \sigma^2 at the breakpoint

beta_flag

boolean indicating whether to re-estimate \beta's at the breakpoint

nu_flag

boolean indicating whether to re-estimate \nu at the breakpoint

Value

An object of class ‘pomp’

Examples

## Not run: 
m1 <- haiti1_joint(
  vacscen = "id0", breakpoint = 232, rho_flag = T,
  tau_flag = F, sig_sq_flag = T, beta_flag = F, nu_flag = F
  )

## End(Not run)

zjiang2/haitipkg documentation built on March 1, 2024, 11:34 p.m.