pith. sign in

arxiv: 2605.03608 · v1 · submitted 2026-05-05 · 📊 stat.ME

Bayesian copula-based modelling for multi-type spatio-temporal epidemic data

Matthew Adeoye , Simon E.F. Spencer , Xavier Didelot This is my paper

Pith reviewed 2026-05-07 13:56 UTC · model grok-4.3

classification 📊 stat.ME
keywords Bayesian modelingcopulaspatio-temporalepidemic modelingmulti-type pathogensMCMCinfectious diseasemodel comparison
0
0 comments X

The pith

A Bayesian model uses copulas to capture how different strains of a pathogen interact across space and time.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a new way to model epidemics caused by multiple strains of the same pathogen that spread over geography and time. It combines a joint state-space representation of all strain epidemics with copula functions to describe their statistical dependencies. This allows researchers to infer how strains influence each other biologically while accounting for spatial and temporal patterns. The approach is tested on simulated data where it recovers the true interactions and parameters, and then applied to real data on meningococcal disease in Europe. If successful, it provides a flexible framework for understanding complex multi-type outbreaks without assuming independence between strains.

Core claim

The authors develop a multi-type spatio-temporal model that formulates a joint state-space for epidemics of different strains and uses copula models to uncover the dependence structure between them, supported by biologically informed interaction terms. They provide an efficient MCMC sampling scheme and robust model comparison via bridge and importance sampling. On simulated data, the models correctly identify epidemics and infer parameters, and they are fitted to monthly incidence data for invasive meningococcal disease across 26 European countries.

What carries the argument

Copula models integrated into a joint state-space representation of multi-strain epidemic dynamics, with biologically informed interaction terms, enabling inference via MCMC.

Load-bearing premise

The copula functions, together with the chosen interaction terms, accurately capture the real dependence structure among strains without distorting the inference.

What would settle it

If the model is fitted to data from a known multi-strain outbreak where strain interactions are independently measured and the inferred dependencies do not match or lead to poor predictive performance on held-out data.

Figures

Figures reproduced from arXiv: 2605.03608 by Matthew Adeoye, Simon E.F. Spencer, Xavier Didelot.

Figure 1
Figure 1. Figure 1: 2.2.1 Spatio-temporal model for the sporadic cases The spatio-temporal model for the sporadic cases excludes the epidemic component in Equation 4, where 𝑎𝑘 remains the strain-specific intercept. The trend, seasonal and spatial components are given priors from the family of intrinsic Gaussian Markov random fields [31], designed to induce smoothness in either space or time. The trend component 𝑟𝑡 is assumed … view at source ↗
Figure 1
Figure 1. Figure 1: Directed acyclic graph (DAG) illustrating the relationships between all model components and the observed view at source ↗
Figure 2
Figure 2. Figure 2: Multi-type model simulations with 𝐾 = 5 strains from (A) No epidemic model (B) Independent 1 model (C) Independent 2 model (D) Frank copula 1 model (F) Frank copula 2 model (G) Gaussian factor copula 1 model (H) Gaussian factor copula 2 model, and (H) General-dependent model. For each of these eight simulated datasets, we wanted to perform inference under the same model as was used for simulation. We tried… view at source ↗
Figure 3
Figure 3. Figure 3: Posterior means, 95% credible intervals and the true values of the temporal trend and seasonal components view at source ↗
Figure 4
Figure 4. Figure 4: Posterior densities and the true values of the spatial components, intercepts, epidemic effects, transition view at source ↗
Figure 5
Figure 5. Figure 5: Heatmap showing simulated epidemics per strain and inferred probability of the epidemic state for strains 1 to view at source ↗
Figure 6
Figure 6. Figure 6: Posterior means for the meningococcal disease application of the trend components, and the seasonal view at source ↗
Figure 7
Figure 7. Figure 7: Posterior median relative risks (compared to the geometric mean risk, the mean of the log of the risk) for the view at source ↗
Figure 8
Figure 8. Figure 8: Heat-maps showing for the meningococcal disease application the posterior probability of the epidemic state view at source ↗
read the original abstract

The study of infectious disease epidemiology for multi-type disease pathogens requires modelling techniques that account for the complex interactions existing between strains across geography and time. In this paper, we propose a novel multi-type spatio-temporal infectious disease model to better support the understanding of these pathogens. We formulate a joint state-space for all epidemics arising for a given multi-type pathogen as well as biologically informed representations of how these epidemic states may interact. We introduce the use of several copula models to uncover the dependence structure of epidemics between strains. We develop a computationally efficient Markov chain Monte Carlo (MCMC) sampling scheme for all proposed models. We also provide robust model comparison techniques using bridge sampling and importance sampling to evaluate model evidence in high-dimensional space. We demonstrate the performance of our proposed models using simulated datasets, where simulated epidemics were successfully identified and associated parameters correctly inferred. The proposed models were also fitted to monthly multi-type incidence data on invasive meningococcal disease from 26 European countries. The accompanying software is freely available as a R package at https://github.com/Matthewadeoye/MultiOutbreaks.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper proposes a novel Bayesian multi-type spatio-temporal model for infectious disease epidemics that uses a joint state-space formulation with biologically informed interaction terms and several copula families to capture dependence structures among pathogen strains. It develops an efficient MCMC sampler, provides bridge and importance sampling for model evidence, validates parameter recovery on simulated data, and applies the framework to monthly invasive meningococcal disease incidence across 26 European countries, with accompanying open-source R software.

Significance. If the central claims hold, the work offers a flexible, computationally tractable approach to joint modeling of interacting strains that could improve understanding of multi-type pathogen dynamics; the provision of reproducible code and MCMC/bridge-sampling machinery are concrete strengths that lower the barrier for adoption.

major comments (2)
  1. [Abstract / Simulation study] Abstract and simulation results: the claim that 'simulated epidemics were successfully identified and associated parameters correctly inferred' lacks any quantitative recovery metrics (bias, coverage, RMSE, or comparison to independent baselines), which is load-bearing for establishing that the copula-plus-interaction parameterization recovers the target quantities rather than fitting by construction.
  2. [Application to real data] Real-data application: the meningococcal incidence fit is presented without reported posterior predictive checks, sensitivity to copula choice, or external validation against known cross-strain dependence patterns; this leaves open the skeptic concern that the chosen copulas may introduce artifacts when the true process contains unmodeled spatio-temporal heterogeneity or reporting delays.
minor comments (2)
  1. [Abstract] The abstract lists 'several copula models' but does not name the families (Gaussian, Clayton, etc.) or the specific interaction terms; adding these would improve immediate clarity.
  2. [Software availability] The GitHub repository is cited but the manuscript contains no usage example or function signatures; a short code snippet or table of exported functions would aid reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their detailed and constructive review of our manuscript. We address each of the major comments below, indicating the revisions we plan to make to strengthen the paper.

read point-by-point responses

  1. Referee: [Abstract / Simulation study] Abstract and simulation results: the claim that 'simulated epidemics were successfully identified and associated parameters correctly inferred' lacks any quantitative recovery metrics (bias, coverage, RMSE, or comparison to independent baselines), which is load-bearing for establishing that the copula-plus-interaction parameterization recovers the target quantities rather than fitting by construction.

    Authors: We agree that providing quantitative metrics would better substantiate the simulation results. In the revised manuscript, we will augment the simulation study section with tables reporting bias, RMSE, and 95% credible interval coverage for key parameters across 100 simulation replicates. We will also include a comparison to an independent baseline model (e.g., separate univariate fits) to highlight the benefits of the joint copula-based approach. These changes will clarify that the parameterization recovers the target quantities effectively. revision: yes

  2. Referee: [Application to real data] Real-data application: the meningococcal incidence fit is presented without reported posterior predictive checks, sensitivity to copula choice, or external validation against known cross-strain dependence patterns; this leaves open the skeptic concern that the chosen copulas may introduce artifacts when the true process contains unmodeled spatio-temporal heterogeneity or reporting delays.

    Authors: We acknowledge the value of these additional checks for the real-data application. We will add posterior predictive checks, including plots of replicated data versus observed incidence, to the revised paper. A sensitivity analysis to different copula families (e.g., Gaussian, Clayton, Gumbel) will be included, with comparisons of posterior inferences and model evidence via bridge sampling. For external validation, we will reference and compare our estimated dependence structures to published studies on meningococcal strain interactions in Europe. While unmodeled factors like reporting delays are a general challenge in epidemic modeling, the spatio-temporal structure and copula flexibility help mitigate this; we will discuss this limitation explicitly. revision: yes

Circularity Check

0 steps flagged

No circularity: standard Bayesian modeling with copulas and MCMC applied to epidemic data

full rationale

The paper formulates a joint state-space model with biologically informed interaction terms, introduces copula functions for dependence, and develops an MCMC sampler with bridge sampling for evidence. Simulation studies recover parameters from data generated under the model, which is standard validation rather than a reduction of outputs to inputs by construction. The real-data application to meningococcal incidence is presented without any self-referential fitting that renames parameters as predictions. No self-citation load-bearing steps, uniqueness theorems, or ansatz smuggling appear in the derivation chain. The framework remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The framework rests on standard Bayesian modeling assumptions plus domain-specific choices about how strains interact; multiple parameters are estimated from data.

free parameters (2)
  • Copula dependence parameters
    Parameters that control the strength and form of dependence between strain-specific epidemic states; fitted via MCMC.
  • State-space transition rates
    Parameters governing how epidemic states evolve and interact across strains, space, and time; estimated from incidence data.
axioms (2)
  • domain assumption A joint state-space can represent the epidemic processes of multiple strains with biologically informed interactions.
    Invoked in the model formulation section of the abstract.
  • domain assumption Copula functions can capture the dependence structure among epidemic states while leaving marginal distributions flexible.
    Core modeling choice stated in the abstract.

pith-pipeline@v0.9.0 · 5490 in / 1367 out tokens · 83477 ms · 2026-05-07T13:56:05.303903+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

70 extracted references

  1. [1]

    Ueber die agglutination der pneumokokken und über die theorieen der agglutination.Zeitschrift für Hygiene und Infektionskrankheiten, 40(1):54–72, 1902

    Fred Neufeld. Ueber die agglutination der pneumokokken und über die theorieen der agglutination.Zeitschrift für Hygiene und Infektionskrankheiten, 40(1):54–72, 1902

    1902

  2. [2]

    Streptococcus pneumoniae epidemiology, pathogenesis and control.Nature Reviews Microbiology, 23(4):256–271, 2025

    Ana Rita Narciso, Rebecca Dookie, Priyanka Nannapaneni, Staffan Normark, and Birgitta Henriques-Normark. Streptococcus pneumoniae epidemiology, pathogenesis and control.Nature Reviews Microbiology, 23(4):256–271, 2025

    2025

  3. [3]

    Neisseria meningitidis: biology, microbiology, and epidemiology

    Nadine G Rouphael and David S Stephens. Neisseria meningitidis: biology, microbiology, and epidemiology. Neisseria meningitidis: advanced methods and protocols, pages 1–20, 2011

    2011

  4. [4]

    Evolution of invasive meningococcal disease epidemiology in Europe, 2008 to 2017.Eurosurveillance, 27(3):2002075, 2022

    Charles Nuttens, Jamie Findlow, Paul Balmer, David L Swerdlow, and Myint Tin Tin Htar. Evolution of invasive meningococcal disease epidemiology in Europe, 2008 to 2017.Eurosurveillance, 27(3):2002075, 2022

    2008

  5. [5]

    Foot-and-mouth disease: past, present and future.Veterinary research, 44(1):116, 2013

    Syed M Jamal and Graham J Belsham. Foot-and-mouth disease: past, present and future.Veterinary research, 44(1):116, 2013

    2013

  6. [6]

    Multilocus sequence typing of bacteria.Annu

    Martin CJ Maiden. Multilocus sequence typing of bacteria.Annu. Rev. Microbiol., 60(1):561–588, 2006

    2006

  7. [7]

    Real-time PCR in virology.Nucleic acids research, 30(6):1292–1305, 2002

    Ian M Mackay, Katherine E Arden, and Andreas Nitsche. Real-time PCR in virology.Nucleic acids research, 30(6):1292–1305, 2002

    2002

  8. [8]

    Real-time reverse transcription PCR (qRT-PCR) and its potential use in clinical diagnosis.Clinical Science, 109(4):365–379, 2005

    Stephen A Bustin and Reinhold Mueller. Real-time reverse transcription PCR (qRT-PCR) and its potential use in clinical diagnosis.Clinical Science, 109(4):365–379, 2005

    2005

  9. [9]

    Compartmental models in epidemiology

    Fred Brauer. Compartmental models in epidemiology. InMathematical epidemiology, pages 19–79. Springer, 2008

    2008

  10. [10]

    A contribution to the mathematical theory of epidemics

    William Ogilvy Kermack and Anderson G McKendrick. A contribution to the mathematical theory of epidemics. Proceedings of the Royal Society of London. Series A, Containing papers of a mathematical and physical character, 115(772):700–721, 1927

    1927

  11. [11]

    Dynamics and selection of many-strain pathogens.Proceedings of the National Academy of Sciences, 99(26):17209–17214, 2002

    Julia R Gog and Bryan T Grenfell. Dynamics and selection of many-strain pathogens.Proceedings of the National Academy of Sciences, 99(26):17209–17214, 2002

    2002

  12. [12]

    Castillo-Chavez, H

    C. Castillo-Chavez, H. W. Hethcote, V . Andreasen, S. A. Levin, and W. M. Liu. Epidemiological models with age structure, proportionate mixing, and cross-immunity.Journal of Mathematical Biology, 27(3):233–258, May 1989

    1989

  13. [13]

    Cross immunity and vaccination against multiple microparasite strains

    LJ White, MJ Cox, and GF Medley. Cross immunity and vaccination against multiple microparasite strains. Mathematical Medicine and Biology: A Journal of the IMA, 15(3):211–233, 1998

    1998

  14. [14]

    Integrating life history and cross-immunity into the evolutionary dynamics of pathogens.Proceedings of the Royal Society B: Biological Sciences, 273(1585):409–416, February 2006

    Olivier Restif and Bryan T Grenfell. Integrating life history and cross-immunity into the evolutionary dynamics of pathogens.Proceedings of the Royal Society B: Biological Sciences, 273(1585):409–416, February 2006

    2006

  15. [15]

    Ted Cohen, Marc Lipsitch, Rochelle P Walensky, and Megan Murray. Beneficial and perverse effects of isoniazid preventive therapy for latent tuberculosis infection in HIV–tuberculosis coinfected populations.Proceedings of the National Academy of Sciences, 103(18):7042–7047, 2006

    2006

  16. [16]

    A modeling framework for the evolution and spread of antibiotic resistance: literature review and model categorization.American Journal of Epidemiology, 178(4):508–520, 2013

    Ian H Spicknall, Betsy Foxman, Carl F Marrs, and Joseph NS Eisenberg. A modeling framework for the evolution and spread of antibiotic resistance: literature review and model categorization.American Journal of Epidemiology, 178(4):508–520, 2013

    2013

  17. [17]

    Leonhard Knorr-Held and Sylvia Richardson. A hierarchical model for space–time surveillance data on meningo- coccal disease incidence.Journal of the Royal Statistical Society Series C: Applied Statistics, 52(2):169–183, 2003

    2003

  18. [18]

    spTimer: Spatio-temporal Bayesian modeling using R.Journal of Statistical Software, 63:1–32, 2015

    Khandoker Shuvo Bakar and Sujit K Sahu. spTimer: Spatio-temporal Bayesian modeling using R.Journal of Statistical Software, 63:1–32, 2015

    2015

  19. [19]

    Spatio-temporal analysis of epidemic phenomena using the R package surveillance.Journal of Statistical Software, 77:1–55, 2017

    Sebastian Meyer, Leonhard Held, and Michael Höhle. Spatio-temporal analysis of epidemic phenomena using the R package surveillance.Journal of Statistical Software, 77:1–55, 2017

    2017

  20. [20]

    Strategies for containing an emerging influenza pandemic in Southeast Asia.Nature, 437(7056):209–214, 2005

    Neil M Ferguson, Derek AT Cummings, Simon Cauchemez, Christophe Fraser, Steven Riley, Aronrag Meeyai, Sopon Iamsirithaworn, and Donald S Burke. Strategies for containing an emerging influenza pandemic in Southeast Asia.Nature, 437(7056):209–214, 2005

    2005

  21. [21]

    Springer Science & Business Media, 2012

    Håkan Andersson and Tom Britton.Stochastic epidemic models and their statistical analysis, volume 151. Springer Science & Business Media, 2012

    2012

  22. [22]

    Antibiotic-resistant Neisseria gonorrhoeae spread faster with more treatment, not more sexual partners.PLoS pathogens, 12(5):e1005611, 2016

    Stephanie M Fingerhuth, Sebastian Bonhoeffer, Nicola Low, and Christian L Althaus. Antibiotic-resistant Neisseria gonorrhoeae spread faster with more treatment, not more sexual partners.PLoS pathogens, 12(5):e1005611, 2016. 29 Bayesian copula-based modelling for multi-type spatio-temporal epidemic data

    2016

  23. [23]

    Estimating the fitness cost and benefit of cefixime resistance in Neisseria gonorrhoeae to inform prescription policy: a modelling study.PLoS medicine, 14(10):e1002416, 2017

    Lilith K Whittles, Peter J White, and Xavier Didelot. Estimating the fitness cost and benefit of cefixime resistance in Neisseria gonorrhoeae to inform prescription policy: a modelling study.PLoS medicine, 14(10):e1002416, 2017

    2017

  24. [24]

    The impact of COVID-19 and strategies for mitigation and suppression in low-and middle-income countries.Science, 369(6502):413–422, 2020

    Patrick GT Walker, Charles Whittaker, Oliver J Watson, Marc Baguelin, Peter Winskill, Arran Hamlet, Bimandra A Djafaara, Zulma Cucunubá, Daniela Olivera Mesa, Will Green, et al. The impact of COVID-19 and strategies for mitigation and suppression in low-and middle-income countries.Science, 369(6502):413–422, 2020

    2020

  25. [25]

    Possible future waves of SARS-CoV-2 infection generated by variants of concern with a range of characteristics.Nature communications, 12(1):5730, 2021

    Louise Dyson, Edward M Hill, Sam Moore, Jacob Curran-Sebastian, Michael J Tildesley, Katrina A Lythgoe, Thomas House, Lorenzo Pellis, and Matt J Keeling. Possible future waves of SARS-CoV-2 infection generated by variants of concern with a range of characteristics.Nature communications, 12(1):5730, 2021

    2021

  26. [26]

    The detection of spatially localised outbreaks in campylobacteriosis notification data.Spatial and spatio-temporal epidemiology, 2(3):173–183, 2011

    Simon EF Spencer, Jonathan Marshall, Ruth Pirie, Donald Campbell, and Nigel P French. The detection of spatially localised outbreaks in campylobacteriosis notification data.Spatial and spatio-temporal epidemiology, 2(3):173–183, 2011

    2011

  27. [27]

    Matthew Adeoye, Xavier Didelot, and Simon E.F. Spencer. Bayesian spatio-temporal modelling for infectious disease outbreak detection.Epidemics, 54:100879, 2026

    2026

  28. [28]

    Springer, 2006

    Roger B Nelsen.An introduction to copulas. Springer, 2006

    2006

  29. [29]

    Riemann manifold Langevin and Hamiltonian Monte Carlo methods.Journal of the Royal Statistical Society Series B: Statistical Methodology, 73(2):123–214, 2011

    Mark Girolami and Ben Calderhead. Riemann manifold Langevin and Hamiltonian Monte Carlo methods.Journal of the Royal Statistical Society Series B: Statistical Methodology, 73(2):123–214, 2011

    2011

  30. [30]

    ECDC ATLAS Database, 2024

    European Centre for Disease Prevention and Control (ECDC). ECDC ATLAS Database, 2024. Accessed: 2024-10-30

    2024

  31. [31]

    Chapman and Hall/CRC, 2005

    Havard Rue and Leonhard Held.Gaussian Markov random fields: theory and applications. Chapman and Hall/CRC, 2005

    2005

  32. [32]

    Fonctions de répartition à n dimensions et leurs marges

    M Sklar. Fonctions de répartition à n dimensions et leurs marges. InAnnales de l’ISUP, volume 8, pages 229–231, 1959

    1959

  33. [33]

    Springer, 2018

    Marius Hofert, Ivan Kojadinovic, Martin Mächler, and Jun Yan.Elements of copula modeling with R. Springer, 2018

    2018

  34. [34]

    CRC press, 2014

    Harry Joe.Dependence modeling with copulas. CRC press, 2014

    2014

  35. [35]

    Factor copula models for multivariate data.Journal of Multivariate Analysis, 120:85–101, 2013

    Pavel Krupskii and Harry Joe. Factor copula models for multivariate data.Journal of Multivariate Analysis, 120:85–101, 2013

    2013

  36. [36]

    Chapman and Hall/CRC, 2009

    Walter Zucchini and Iain L MacDonald.Hidden Markov models for time series: an introduction using R. Chapman and Hall/CRC, 2009

    2009

  37. [37]

    Intermittent missing observations in discrete-time hidden Markov models.Communications in Statistics-Simulation and Computation, 41(2):167–181, 2012

    Hung-Wen Yeh, Wenyaw Chan, and Elaine Symanski. Intermittent missing observations in discrete-time hidden Markov models.Communications in Statistics-Simulation and Computation, 41(2):167–181, 2012

    2012

  38. [38]

    Still ‘dairy farm fever’? A Bayesian model for leptospirosis notification data in New Zealand.Journal of the Royal Society Interface, 18(175), 2021

    Jackie Benschop, Shahista Nisa, and Simon EF Spencer. Still ‘dairy farm fever’? A Bayesian model for leptospirosis notification data in New Zealand.Journal of the Royal Society Interface, 18(175), 2021

    2021

  39. [39]

    On the use of non-local prior densities in Bayesian hypothesis tests.Journal of the Royal Statistical Society Series B: Statistical Methodology, 72(2):143–170, 2010

    Valen E Johnson and David Rossell. On the use of non-local prior densities in Bayesian hypothesis tests.Journal of the Royal Statistical Society Series B: Statistical Methodology, 72(2):143–170, 2010

    2010

  40. [40]

    Hybrid Monte Carlo.Physics letters B, 195(2):216–222, 1987

    Simon Duane, Anthony D Kennedy, Brian J Pendleton, and Duncan Roweth. Hybrid Monte Carlo.Physics letters B, 195(2):216–222, 1987

    1987

  41. [41]

    MCMC using Hamiltonian dynamics.Handbook of markov chain monte carlo, 2(11):2, 2011

    Radford M Neal et al. MCMC using Hamiltonian dynamics.Handbook of markov chain monte carlo, 2(11):2, 2011

    2011

  42. [42]

    Exponential convergence of Langevin distributions and their discrete approximations.Bernoulli, 2:341–363, 1996

    Gareth O Roberts and Richard L Tweedie. Exponential convergence of Langevin distributions and their discrete approximations.Bernoulli, 2:341–363, 1996

    1996

  43. [43]

    Optimal scaling of discrete approximations to Langevin diffusions

    Gareth O Roberts and Jeffrey S Rosenthal. Optimal scaling of discrete approximations to Langevin diffusions. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 60(1):255–268, 1998

    1998

  44. [44]

    Equation of state calculations by fast computing machines.The journal of chemical physics, 21(6):1087–1092, 1953

    Nicholas Metropolis, Arianna W Rosenbluth, Marshall N Rosenbluth, Augusta H Teller, and Edward Teller. Equation of state calculations by fast computing machines.The journal of chemical physics, 21(6):1087–1092, 1953

    1953

  45. [45]

    Monte Carlo sampling methods using Markov chains and their applications.Biometrika, 57:97–109, 1970

    W Keith Hastings. Monte Carlo sampling methods using Markov chains and their applications.Biometrika, 57:97–109, 1970

    1970

  46. [46]

    Efficient Bayesian structural equation modeling in Stan.Journal of Statistical Software, 100:1–22, 2021

    Edgar C Merkle, Ellen Fitzsimmons, James Uanhoro, and Ben Goodrich. Efficient Bayesian structural equation modeling in Stan.Journal of Statistical Software, 100:1–22, 2021. 30 Bayesian copula-based modelling for multi-type spatio-temporal epidemic data

    2021

  47. [47]

    A stochastic approximation method.The annals of mathematical statistics, pages 400–407, 1951

    Herbert Robbins and Sutton Monro. A stochastic approximation method.The annals of mathematical statistics, pages 400–407, 1951

    1951

  48. [48]

    Estimating HIV, HCV and HSV2 incidence from emergency department serosurvey.Gates Open Research, 5(116):116, 2021

    Simon EF Spencer, Oliver Laeyendecker, Louise Dyson, Yu-Hsiang Hsieh, Eshan U Patel, Richard E Rothman, Gabor D Kelen, Thomas C Quinn, and T Deirdre Hollingsworth. Estimating HIV, HCV and HSV2 incidence from emergency department serosurvey.Gates Open Research, 5(116):116, 2021

    2021

  49. [49]

    On the stability and ergodicity of adaptive scaling Metropolis algorithms.Stochastic processes and their applications, 121(12):2839–2860, 2011

    Matti Vihola. On the stability and ergodicity of adaptive scaling Metropolis algorithms.Stochastic processes and their applications, 121(12):2839–2860, 2011

    2011

  50. [50]

    Accelerating adaptation in the adaptive Metropolis–Hastings random walk algorithm.Aus- tralian & New Zealand Journal of Statistics, 63(3):468–484, 2021

    Simon EF Spencer. Accelerating adaptation in the adaptive Metropolis–Hastings random walk algorithm.Aus- tralian & New Zealand Journal of Statistics, 63(3):468–484, 2021

    2021

  51. [51]

    Stan: A probabilistic programming language.Journal of statistical software, 76, 2017

    Bob Carpenter, Andrew Gelman, Matthew D Hoffman, Daniel Lee, Ben Goodrich, Michael Betancourt, Marcus A Brubaker, Jiqiang Guo, Peter Li, and Allen Riddell. Stan: A probabilistic programming language.Journal of statistical software, 76, 2017

    2017

  52. [52]

    A direct optimization approach to hidden Markov modeling for single channel kinetics.Biophysical journal, 79(4):1915–1927, 2000

    Feng Qin, Anthony Auerbach, and Frederick Sachs. A direct optimization approach to hidden Markov modeling for single channel kinetics.Biophysical journal, 79(4):1915–1927, 2000

    1915

  53. [53]

    Exact computation of the observed information matrix for hidden Markov models.Journal of Computational and Graphical Statistics, 11(3):678–689, 2002

    Theodore C Lystig and James P Hughes. Exact computation of the observed information matrix for hidden Markov models.Journal of Computational and Graphical Statistics, 11(3):678–689, 2002

    2002

  54. [54]

    Direct maximization of the likelihood of a hidden Markov model.Computational Statistics & Data Analysis, 52(9):4147–4160, 2008

    Rolf Turner. Direct maximization of the likelihood of a hidden Markov model.Computational Statistics & Data Analysis, 52(9):4147–4160, 2008

    2008

  55. [55]

    Finding the observed information matrix when using the EM algorithm.Journal of the Royal Statistical Society Series B: Statistical Methodology, 44(2):226–233, 1982

    Thomas A Louis. Finding the observed information matrix when using the EM algorithm.Journal of the Royal Statistical Society Series B: Statistical Methodology, 44(2):226–233, 1982

    1982

  56. [56]

    CRC press, 2014

    Randal Douc, Eric Moulines, and David Stoffer.Nonlinear time series: Theory, methods and applications with R examples. CRC press, 2014

    2014

  57. [57]

    Monte Carlo on manifolds: sampling densities and integrating functions.Communications on Pure and Applied Mathematics, 71(12):2609–2647, 2018

    Emilio Zappa, Miranda Holmes-Cerfon, and Jonathan Goodman. Monte Carlo on manifolds: sampling densities and integrating functions.Communications on Pure and Applied Mathematics, 71(12):2609–2647, 2018

    2018

  58. [58]

    Bayes factors.Journal of the American Statistical Association, 90(430):773– 795, 1995

    Robert E Kass and Adrian E Raftery. Bayes factors.Journal of the American Statistical Association, 90(430):773– 795, 1995

    1995

  59. [59]

    OuP Oxford, 1998

    Harold Jeffreys.The theory of probability. OuP Oxford, 1998

    1998

  60. [60]

    A tutorial on bridge sampling.Journal of Mathematical Psychology, 81:80–97, 2017

    Quentin F Gronau, Alexandra Sarafoglou, Dora Matzke, Alexander Ly, Udo Boehm, Maarten Marsman, David S Leslie, Jonathan J Forster, Eric-Jan Wagenmakers, and Helen Steingroever. A tutorial on bridge sampling.Journal of Mathematical Psychology, 81:80–97, 2017

    2017

  61. [61]

    Default Bayesian model determination methods for generalised linear mixed models.Computational Statistics & Data Analysis, 54(12):3269–3288, 2010

    Antony M Overstall and Jonathan J Forster. Default Bayesian model determination methods for generalised linear mixed models.Computational Statistics & Data Analysis, 54(12):3269–3288, 2010

    2010

  62. [62]

    Panayiota Touloupou, Naif Alzahrani, Peter Neal, Simon E. F. Spencer, and Trevelyan J. McKinley. Efficient model comparison techniques for models requiring large scale data augmentation.Bayesian Analysis, 13(2):437 – 459, 2018

    2018

  63. [63]

    Gronau, Henrik Singmann, and Eric-Jan Wagenmakers

    Quentin F. Gronau, Henrik Singmann, and Eric-Jan Wagenmakers. bridgesampling: An R package for estimating normalizing constants.Journal of Statistical Software, 92(10):1–29, 2020

    2020

  64. [64]

    Eurostat Database, 2024

    Eurostat. Eurostat Database, 2024. Accessed: 2024-11-05

    2024

  65. [65]

    Mapping and measuring country shapes: The CShapes package.R J, 2:18–23, 2010

    K Gleditsch and NB Weidmann. Mapping and measuring country shapes: The CShapes package.R J, 2:18–23, 2010

    2010

  66. [66]

    CODA: convergence diagnosis and output analysis for MCMC.R news, 6(1):7–11, 2006

    Martyn Plummer, Nicky Best, Kate Cowles, Karen Vines, et al. CODA: convergence diagnosis and output analysis for MCMC.R news, 6(1):7–11, 2006

    2006

  67. [67]

    Shamez N Ladhani, Kazim Beebeejaun, Jay Lucidarme, Helen Campbell, Steve Gray, Ed Kaczmarski, Mary E Ramsay, and Ray Borrow. Increase in endemic Neisseria meningitidis capsular group W sequence type 11 complex associated with severe invasive disease in England and Wales.Clinical Infectious Diseases, 60(4):578–585, 2015

    2015

  68. [68]

    Increase of invasive meningococcal serogroup W disease in Europe, 2013 to 2017.Eurosurveillance, 24(14):1800245, 2019

    Manuel Krone, Steve Gray, Raquel Abad, Anna Skoczy ´nska, Paola Stefanelli, Arie van der Ende, Georgina Tzanakaki, Paula Mölling, Maria João Simões, Pavla Kˇrížová, et al. Increase of invasive meningococcal serogroup W disease in Europe, 2013 to 2017.Eurosurveillance, 24(14):1800245, 2019

    2013

  69. [69]

    Posterior predictive assessment of model fitness via realized discrepancies.Statistica sinica, pages 733–760, 1996

    Andrew Gelman, Xiao-Li Meng, and Hal Stern. Posterior predictive assessment of model fitness via realized discrepancies.Statistica sinica, pages 733–760, 1996

    1996

  70. [70]

    Bayesian image restoration, with two applications in spatial statistics.Annals of the institute of statistical mathematics, 43:1–20, 1991

    Julian Besag, Jeremy York, and Annie Mollié. Bayesian image restoration, with two applications in spatial statistics.Annals of the institute of statistical mathematics, 43:1–20, 1991. 31 Bayesian copula-based modelling for multi-type spatio-temporal epidemic data Supplementary Material u 1 u 4 u 7 u 2 u 3 u 5 u 6 u 8 u 9 Figure S1: Adjacency structure of ...

    1991