R/childhood.R

In pomp: Statistical Inference for Partially Observed Markov Processes

##' Historical childhood disease incidence data
##'
##' @description
##' \code{LondonYorke} is a data frame containing the monthly number of reported cases of chickenpox, measles, and mumps from two American cities (Baltimore and New York) in the mid-20th century (1928--1972).
##'
##' \code{ewmeas} and \code{ewcitmeas} are data frames containing weekly reported cases of measles in England and Wales.
##' \code{ewmeas} records the total measles reports for the whole country, 1948--1966.
##' One questionable data point has been replaced with an NA.
##' \code{ewcitmeas} records the incidence in seven English cities 1948--1987.
##' These data were kindly provided by Ben Bolker, who writes:
##' \dQuote{Most of these data have been manually entered from published records by various people, and are prone to errors at several levels.
##' All data are provided as is; use at your own risk.}
##'
##' @name childhood_disease_data
##' @rdname childhood
##' @aliases LondonYorke ewmeas ewcitmeas
##' @docType data
##' @seealso \link[=compartmental_models]{compartmental models}, \code{\link{bsflu}}
##' @references
##'
##' \London1973
##' 
##' @family pomp datasets
##' @family pomp examples
##'
##' @examples
##'
##' plot(cases~time,data=LondonYorke,subset=disease=="measles",type='n',main="measles",bty='l')
##' lines(cases~time,data=LondonYorke,subset=disease=="measles"&town=="Baltimore",col="red")
##' lines(cases~time,data=LondonYorke,subset=disease=="measles"&town=="New York",col="blue")
##' legend("topright",legend=c("Baltimore","New York"),lty=1,col=c("red","blue"),bty='n')
##'
##' plot(
##'      cases~time,
##'      data=LondonYorke,
##'      subset=disease=="chickenpox"&town=="New York",
##'      type='l',col="blue",main="chickenpox, New York",
##'      bty='l'
##'     )
##'
##' plot(
##'      cases~time,
##'      data=LondonYorke,
##'      subset=disease=="mumps"&town=="New York",
##'      type='l',col="blue",main="mumps, New York",
##'      bty='l'
##'     )
##'
##' plot(reports~time,data=ewmeas,type='l')
##'
##' plot(reports~date,data=ewcitmeas,subset=city=="Liverpool",type='l')
##'
NULL