measles {pomp} | R Documentation

## Historical childhood disease incidence data

### Description

`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).

`ewmeas` and `ewcitmeas` are data frames containing weekly reported cases of measles in England and Wales. `ewmeas` records the total measles reports for the whole country, 1948–1966. One questionable data point has been replaced with an NA. `ewcitmeas` records the incidence in seven English cities 1948–1987. These data were kindly provided by Ben Bolker, who writes: “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.”

### References

W. P. London and J. A. Yorke, Recurrent Outbreaks of Measles, Chickenpox and Mumps: I. Seasonal Variation in Contact Rates, American Journal of Epidemiology, 98:453–468, 1973.

### See Also

Other datasets: `blowflies`, `bsflu`, `dacca`, `ebola`, `parus`

Other pomp examples: `blowflies`, `dacca`, `ebola`, `gompertz`, `ou2`, `ricker`, `rw2`, `sir_models`, `verhulst`

### 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')

[Package pomp version 2.3 Index]