Time series analysis (TSA) of human invasive listeriosis trends, 2008–2015
Creators
- EFSA Panel on Biological Hazards (BIOHAZ)
- Ricci, Antonia1
- Allende, Ana2
- Bolton, Declan3
- Chemaly, Marianne4
- Davies, Robert5
- Fernández Escámez, Pablo Salvador6
- Girones, Rosina7
- Herman, Lieve8
- Koutsoumanis, Konstantinos9
- Nørrung, Birgit10
- Robertson, Lucy11
- Ru, Giuseppe12
- Sanaa, Moez4
- Simmons, Marion13
- Skandamis, Panagiotis14
- Snary, Emma15
- Speybroeck, Niko16
- Ter Kuile, Benno17
- Threlfall, John13
- Wahlström, Helene18
- Takkinen, Johanna19
- Wagner, Martin20
- Arcella, Davide21
- Da Silva Felicio, Maria Teresa21
- Georgiadis, Marios21
- Messens, Winy21
- Lindqvist, Roland22
- 1. Istituto Zooprofilattico Sperimentale delle Venezie, Italy
- 2. CEBAS-CSIC, Spain: Departamento de Ciencia y Tecnología de los Alimentos, Grupo de Calidad, Seguridad y Bioactividad de Alimentos Vegetables, Centro de Edafología y Biología aplicada del Segura (CEBAS) Consejo Superior de Investigaciones Científicas (CSIC), Spain
- 3. Teagasc, Agriculture and Food Development Authority, Ireland
- 4. French Agency for Food, Environmental and Occupational Health & Safety, France
- 5. Animal and Plant Health Agency, Bacteriology Department, United Kingdom
- 6. Technical University of Cartagena, Spain
- 7. University of Barcelona, Spain
- 8. Research Institute for Agriculture, Fisheries and Food, Belgium
- 9. Aristotle University of Thessaloniki, Greece
- 10. University of Copenhagen, Denmark
- 11. Norwegian University of Life Sciences, Norway
- 12. Istituto Zooprofilattico Sperimentale del Piemonte, Liguria e Valle d'Aosta, Italy
- 13. No current affiliation
- 14. Agricultural University of Athens, Greece
- 15. Animal and Plant Health Agency, Department of Epidemiological Sciences, United Kingdom
- 16. Université Catholique de Louvain, Belgium
- 17. Netherlands Food and Consumer Product Safety Authority and University of Amsterdam, The Netherlands
- 18. National Veterinary Institute, Sweden
- 19. European Centre for Disease Prevention and Control (ECDC), Sweden
- 20. University of Veterinary Medicine Vienna, Austria
- 21. European Food Safety Authority, Italy
- 22. National Food Agency, Sweden
Description
The R codes were developed and applied by the EFSA Working Group on Listeria monocytogenes contamination of ready-to-eat foods during the preparatory work on the Scientific Opinion ‘Listeria monocytogenes contamination of ready-to-eat foods and the risk for human health in the EU’ (see 10.2903/j.efsa.2018.5134).
The codes were used to perform time-series analyses of the number of confirmed invasive listeriosis cases in the EU/EEA for the period 2008–2015. Data from The European Surveillance System – TESSy, provided by Austria, Belgium, Cyprus, Czech Republic, Denmark, Estonia, Finland, France, Germany, Greece, Hungary, Iceland, Ireland, Italy, Latvia, Luxembourg, Malta, Netherlands, Norway, Poland, Romania, Slovakia, Slovenia, Spain, Sweden, United Kingdom, and released by ECDC.
Two R codes are provided, one has been used for an analysis using the entire dataset of 14,002 confirmed cases (aggregated analysis) and another one has been used for the analysis using each time a subset of the population (14 subgroups defined by gender and age – disaggregated analysis). The following 14 gender–age group combinations were used: 1–4, 5–14, 15–24, 25–44, 45–64, 65–74, ≥ 75 years old, for both males and females.
The analyses are implemented in R. Two codes are provided, one for the aggregated data analysis and one for the disaggregated data analysis. The working directory needs to be defined in the code (currently there is a placeholder named 'WORKING DIRECTORY'). The data are provided as .csv files ('totals.csv' for the aggregated analysis; 'merged_eu_wide.csv' and 'merged_eu.csv' for the disaggregated analysis) and can be read by the R code if placed in the same working directory. R packages that are required to be installed for the aggregated analysis are 'dlm' and 'strucchange' together with its dependencies, and for the disaggregated analysis the package 'epitools'. Also the script 'pests.R' needs to be downloaded from http://www.utdallas.edu/~pbrandt/code/pests.r and saved in the respective working directory, for both analyses.
The analyses are explained step by step in the R script, which is intended to be run line by line and not entirely in one run. Indeed, the R code currently does not allow for opening new graph windows for new plots, which each time overwrite the previously created ones. It is recommended that the user ignores the specific warning messages that are produced during the disaggregated PAR analysis for the subgroups Female1524, male2544, male6574 and male75, since these have to do with the starting values used for the models but the results are still correct.
Notes
Files
merged_eu.csv
Files
(109.4 kB)
Name | Size | Download all |
---|---|---|
md5:5d8bfbb7a9f6abec6e41f20f7d28f4ec
|
63.7 kB | Preview Download |
md5:30e22b5892c3875c6746b33380ad3d36
|
17.0 kB | Preview Download |
md5:5927ee7c66320032b76bb4736306cbab
|
2.2 kB | Preview Download |
md5:deab5e4942e5e03dda88fed0e35d2dd3
|
7.2 kB | Download |
md5:686f45dca7e520b31cbe1ad5d723f876
|
19.3 kB | Download |
Additional details
Related works
- Is supplement to
- 10.2903/j.efsa.2018.5134 (DOI)
References
- Brandt P and Williams JT, 2001. A linear Poisson autoregressive model: the Poisson AR(p) model. Political Analysis, 9, 164-184.
- Petris G, 2010. An R Package for Dynamic Linear Models. Journal of Statistical Software, 36, 1-16.
- Petris G, Petrone S and Campagnoli P, 2009. Dynamic linear models. In: Eds Gentleman R, Hornik K and Parmigiani G. Dynamic Linear Models with R. Springer, Milano, Italy, pp.31-84.
Subjects
- Risk assessment
- http://id.agrisemantics.org/gacs/C1470
- Risk analysis
- http://id.agrisemantics.org/gacs/C21254
- Listeria monocytogenes
- http://id.agrisemantics.org/gacs/C3027
- Listeriosis
- http://id.agrisemantics.org/gacs/C19651
- Epidemiology
- http://id.agrisemantics.org/gacs/C155
- Disease incidence
- http://id.agrisemantics.org/gacs/C2711
- Time series analysis
- http://id.agrisemantics.org/gacs/C6489
- Trends
- http://id.agrisemantics.org/gacs/C2594
- Autocorrelation
- http://id.agrisemantics.org/gacs/C16575
- Seasonality
- http://id.agrisemantics.org/gacs/C5567
- Forecasting
- http://id.agrisemantics.org/gacs/C420