pith. sign in

arxiv: 2606.08391 · v1 · pith:G5ICKV3Xnew · submitted 2026-06-07 · 🧬 q-bio.PE · q-bio.QM

Cruise Ship-Associated Andes Virus Cluster aboard MV Hondius, 2026: A Stochastic Scenario Analysis

Pith reviewed 2026-06-27 17:59 UTC · model grok-4.3

classification 🧬 q-bio.PE q-bio.QM
keywords Andes viruscruise ship outbreakstochastic epidemic modelhantaviruslatent infectionsreproductive numberapproximate Bayesian computation
0
0 comments X

The pith

Stochastic modeling indicates two latent infections at embarkation best explain the 2026 Andes virus cluster on the MV Hondius.

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

The paper applies a stochastic epidemic model to four embarkation scenarios for the Andes hantavirus outbreak that produced 13 cases on the MV Hondius cruise ship. It finds that the scenario with two latent infected persons boarding is most consistent with the observed final size and takeoff probability. The analysis anchors reproductive numbers to published ANDV estimates and uses approximate Bayesian computation to assess latent infection counts. Findings underscore the value of exposure-history checks, early surveillance, case isolation, and post-disembarkation monitoring for travelers from endemic areas.

Core claim

Scenario D with two latent infected persons at embarkation was most consistent with the observed outbreak, yielding P(final size >= 13) = 11.6% and P(takeoff) = 58.5% at R0 = 2.12. Approximate Bayesian computation provided complementary support for multiple latent infections at embarkation, especially E1(0)=1 and E3(0)=2, though R0 remained weakly identifiable. A day-35 transmission reduction changed takeoff probability little.

What carries the argument

Stochastic epidemic model run on four embarkation scenarios with reproductive numbers taken from published ANDV estimates, used to compute consistency probabilities for observed outbreak size.

Load-bearing premise

Reproductive numbers from published ANDV estimates accurately represent transmission dynamics aboard the ship.

What would settle it

A new outbreak with known initial latent count and final size that ranks a different embarkation scenario highest would contradict the consistency ranking of Scenario D.

Figures

Figures reproduced from arXiv: 2606.08391 by Gerardo Chowell, Hamed Karami, Kaustubh Wagh, Kenji Mizumoto, Raj Kumar Subedi.

Figure 1
Figure 1. Figure 1: Epidemic curve by reported symptom-onset date and case status, MV Hondius Andes virus outbreak, 2026. Bars show the 9 confirmed/probable cases with publicly reported onset dates; four additional confirmed/probable cases lacked reported onset dates and are not displayed. Vertical dashed line indicates cruise embarkation from Ushuaia on April 1, 2026. Source: WHO/ECDC public updates, data as of May 27, 2026 … view at source ↗
Figure 2
Figure 2. Figure 2: Stochastic Poisson SE1E2E3IR trajectories conditioned on epidemic takeoff (final size ≥ 5) under the two-beta model. Solid lines show median simulated cumulative infections; darker bands show the interquartile range and lighter bands show the 2.5th￾97.5th percentile simulation envelope. Dashed line shows the observed cumulative onset curve among cases with known onset dates (n = 9). Dotted vertical line in… view at source ↗
Figure 3
Figure 3. Figure 3: R0 sensitivity: P(final size ≥ 13) across R0 ∈ {0.96, 1.19, 2.12} for Scenarios B, C, and D, two-beta model. R0 values from Martínez et al. 2020. Post-intervention β1 was set as R1 × γ, with R1 = 0.96 from day 35 (May 6) [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Approximate Bayesian computation posterior support over initial conditions E1(0) × E3(0) from ABC simulation-based inference. ABC used 100,000 simulations; the top 1% by weekly incidence RMSE were retained as the approximate posterior. The prior was R0 ~ U(0.5, 3.0). The most supported pair was E1(0) = 1 and E3(0) = 2, with posterior probability 33.8%. The U(0.5, 4.0) sensitivity gave the same posterior mo… view at source ↗
Figure 10
Figure 10. Figure 10 [PITH_FULL_IMAGE:figures/full_fig_p028_10.png] view at source ↗
Figure 11
Figure 11. Figure 11 [PITH_FULL_IMAGE:figures/full_fig_p028_11.png] view at source ↗
Figure 12
Figure 12. Figure 12 [PITH_FULL_IMAGE:figures/full_fig_p029_12.png] view at source ↗
read the original abstract

In April 2026, the MV Hondius expedition cruise ship became the site of the first documented cruise ship-associated Andes hantavirus (ANDV) cluster, with 13 confirmed and probable cases and 3 deaths among 149 passengers and crew. We applied a stochastic epidemic model to evaluate four embarkation scenarios under reproductive numbers anchored to published ANDV estimates. Scenario D, involving two latent infected persons at embarkation, was most consistent with the observed outbreak, yielding P(final size >= 13) = 11.6% and P(takeoff) = 58.5% at R0 = 2.12. Approximate Bayesian computation provided complementary support for multiple latent infections at embarkation, especially E1(0)=1 and E3(0)=2, but R0 remained weakly identifiable. A day-35 transmission reduction changed takeoff probability little in this counterfactual model. Findings support exposure-history assessment, early onboard surveillance, rapid isolation of symptomatic cases, and postdisembarkation monitoring for travelers from ANDV-endemic regions.

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 / 1 minor

Summary. The paper applies a stochastic epidemic model to four embarkation scenarios for an Andes hantavirus (ANDV) outbreak aboard the MV Hondius (149 passengers/crew, 13 cases). It concludes that Scenario D (two latent infections at embarkation) is most consistent with the data, yielding P(final size >=13)=11.6% and P(takeoff)=58.5% at R0=2.12 anchored to published ANDV estimates; approximate Bayesian computation (ABC) supports multiple latent infections at embarkation but finds R0 weakly identifiable. A day-35 transmission reduction is examined as a counterfactual.

Significance. If the anchored R0 and latent-period assumptions hold for shipboard mixing, the scenario probabilities and ABC results offer a quantitative basis for prioritizing exposure-history assessment, onboard surveillance, and post-disembarkation monitoring. The use of ABC to evaluate initial conditions is a methodological strength that allows direct comparison of embarkation hypotheses.

major comments (2)
  1. [Abstract] Abstract: The headline ranking of Scenario D as most consistent rests on probabilities computed at R0=2.12 taken from external ANDV estimates, yet the same abstract states that ABC left R0 weakly identifiable from the outbreak data. This means the reported P(final size >=13) and P(takeoff) values are not constrained by the observed cluster and inherit any mismatch between published household/rodent-contact R0 and the cruise-ship contact matrix (shared cabins, dining, excursions).
  2. [Abstract] Abstract: No model equations, state variables, transition rates, or contact assumptions are supplied, nor are any validation checks, sensitivity analyses, or raw data shown. Without these, it is impossible to verify whether the stochastic structure reproduces the observed final size or whether the latent-period distribution matches the ship setting.
minor comments (1)
  1. [Abstract] Abstract: The day-35 transmission-reduction counterfactual is mentioned without reporting the resulting change in takeoff probability or the precise implementation (e.g., which compartments are affected).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We respond point-by-point to the major comments below and will make revisions to improve clarity and completeness.

read point-by-point responses
  1. Referee: [Abstract] Abstract: The headline ranking of Scenario D as most consistent rests on probabilities computed at R0=2.12 taken from external ANDV estimates, yet the same abstract states that ABC left R0 weakly identifiable from the outbreak data. This means the reported P(final size >=13) and P(takeoff) values are not constrained by the observed cluster and inherit any mismatch between published household/rodent-contact R0 and the cruise-ship contact matrix (shared cabins, dining, excursions).

    Authors: We agree that the reported probabilities for Scenario D are computed under a fixed R0=2.12 drawn from external ANDV literature and are not directly constrained by the ship outbreak data, as the ABC analysis shows R0 to be weakly identifiable. The abstract already notes the ABC result on R0, but we will revise the abstract to more explicitly distinguish the anchored-R0 analysis from the ABC results and to note the possibility that shipboard contact patterns differ from household or rodent-contact settings used in the external R0 estimates. We will also add discussion of this limitation. revision: yes

  2. Referee: [Abstract] Abstract: No model equations, state variables, transition rates, or contact assumptions are supplied, nor are any validation checks, sensitivity analyses, or raw data shown. Without these, it is impossible to verify whether the stochastic structure reproduces the observed final size or whether the latent-period distribution matches the ship setting.

    Authors: The methods section of the manuscript describes the stochastic SEIR model, including state variables, transition rates, and contact assumptions based on ship layout. However, to make the abstract self-contained and address the concern, we will add a concise description of the model structure and key assumptions to the abstract. We will also incorporate validation checks, sensitivity analyses, and a supplementary table of raw case timing data in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity; external R0 anchoring is independent

full rationale

The paper applies a stochastic SEIR-type model to four embarkation scenarios, anchoring R0 to published ANDV estimates from the literature and using ABC to assess consistency with the observed 13-case cluster. R0 is explicitly stated as weakly identifiable from the ship data, so the reported probabilities (e.g., P(final size >=13)=11.6% for Scenario D at R0=2.12) are conditional outputs given those external parameters rather than a self-referential fit or prediction. No equation reduces to its own inputs by construction, no self-citation chain is load-bearing for the central claim, and the model structure is applied to external benchmarks. This is standard parameter use from prior independent studies.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

Abstract-only review limits extraction; key anchored value is R0 drawn from prior literature and initial infection counts inferred via ABC.

free parameters (2)
  • R0 = 2.12
    Anchored to published ANDV estimates and used across scenarios; value 2.12 appears in reported probabilities.
  • initial latent infections
    E1(0) and E3(0) values tested and supported by ABC; central to scenario ranking.
axioms (2)
  • domain assumption Published ANDV reproductive number estimates apply directly to the ship transmission setting.
    Used to anchor all four scenarios.
  • domain assumption Stochastic model structure and contact assumptions match the cruise ship environment.
    Required for the probability calculations to be meaningful.

pith-pipeline@v0.9.1-grok · 5740 in / 1455 out tokens · 20687 ms · 2026-06-27T17:59:30.254724+00:00 · methodology

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

22 extracted references

  1. [1]

    Evidence for human-to-human transmission of hantavirus: a systematic review

    Toledo J, Haby MM, Reveiz L, Sosa Leon L, Angerami R, Aldighieri S. Evidence for human-to-human transmission of hantavirus: a systematic review. J Infect Dis. 2022;226(8):1362-71

  2. [6]

    Transmission potential of the novel coronavirus (COVID-

    Mizumoto K, Chowell G. Transmission potential of the novel coronavirus (COVID-

  3. [7]

    Infect Dis Model

    onboard the diamond Princess Cruises Ship, 2020. Infect Dis Model. 2020;5:264-70

  4. [8]

    Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020

    Mizumoto K, Kagaya K, Zarebski A, Chowell G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Euro Surveill. 2020;25(10):2000180

  5. [9]

    Counting everyone onboard is not enough: modelling lessons from the MV Hondius Andes virus outbreak

    Chowell G, Mizumoto K. Counting everyone onboard is not enough: modelling lessons from the MV Hondius Andes virus outbreak. J Travel Med. 2026 [In Press]

  6. [10]

    Cruise Ships as Amplifiers of Emerging Zoonotic Threats: Implications of the Suspected Andes Hantavirus Cluster aboard MV Hondius

    Thaledi S, Pal S, Juyal D. Cruise Ships as Amplifiers of Emerging Zoonotic Threats: Implications of the Suspected Andes Hantavirus Cluster aboard MV Hondius. J Travel Med. 2026:taag034

  7. [12]

    Hantavirus-associated cluster of illness on a cruise ship: ECDC assessment and recommendations

    ECDC. Hantavirus-associated cluster of illness on a cruise ship: ECDC assessment and recommendations. 2026 [cited 2026 May 17]; Available from: https://www.ecdc.europa.eu/sites/default/files/documents/TAB-hantavirus-06052026.pdf

  8. [13]

    Andes hantavirus outbreak in cruise ship

    European Centre for Disease Prevention and Control. Andes hantavirus outbreak in cruise ship. 2026 [cited 5/27/2026]; Available from: https://www.ecdc.europa.eu/en/infectious-disease-topics/hantavirus- infection/surveillance-and-updates/andes-hantavirus-outbreak

  9. [15]

    Approximate accelerated stochastic simulation of chemically reacting systems

    Gillespie DT. Approximate accelerated stochastic simulation of chemically reacting systems. J Chem Phys. 2001;115(4):1716-33

  10. [22]

    Boundary hit

    R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://wwwR-projectorg/; 2024. Author biographical sketch: Raj Kumar Subedi is a doctoral candidate and a 2CI Fellow in the Department of Population Health Sciences, School of Public Health, Georgia State University, Atlanta, Geor...

  11. [23]

    Disease outbreak news: hantavirus cluster linked to cruise ship travel, multi-country

    World Health Organization. Disease outbreak news: hantavirus cluster linked to cruise ship travel, multi-country. Geneva: World health organization; 2026

  12. [24]

    Viral shedding and viraemia of Andes virus during acute hantavirus infection: a prospective study

    Ferrés M, Martínez-Valdebenito C, Henriquez C, Marco C, Angulo J, Barrera A, et al. Viral shedding and viraemia of Andes virus during acute hantavirus infection: a prospective study. Lancet Infect Dis. 2024;24(7):775-82

  13. [25]

    Equivalence of the Erlang-distributed SEIR epidemic model and the renewal equation

    Champredon D, Dushoff J, Earn DJ. Equivalence of the Erlang-distributed SEIR epidemic model and the renewal equation. SIAM J Appl Math. 2018;78(6):3258-78

  14. [26]

    Hantavirus in humans: a review of clinical aspects and management

    Vial PA, Ferrés M, Vial C, Klingström J, Ahlm C, López R, et al. Hantavirus in humans: a review of clinical aspects and management. Lancet Infect Dis. 2023;23(9):e371-e82

  15. [27]

    Super-spreaders

    Martínez VP, Di Paola N, Alonso DO, Pérez-Sautu U, Bellomo CM, Iglesias AA, et al. “Super-spreaders” and person-to-person transmission of Andes virus in Argentina. N Engl J Med. 2020;383(23):2230-41

  16. [28]

    Harnessing testing strategies and public health measures to avert COVID-19 outbreaks during ocean cruises

    Chowell G, Dahal S, Bono R, Mizumoto K. Harnessing testing strategies and public health measures to avert COVID-19 outbreaks during ocean cruises. Sci Rep. 2021;11(1):15482

  17. [29]

    Epydemix: An open-source Python package for epidemic modeling with integrated approximate Bayesian calibration

    Gozzi N, Chinazzi M, Davis JT, Gioannini C, Rossi L, Ajelli M, et al. Epydemix: An open-source Python package for epidemic modeling with integrated approximate Bayesian calibration. PLOS Comp Bio. 2025;21(11):e1013735

  18. [30]

    Approximate Bayesian computation in population genetics

    Beaumont MA, Zhang W, Balding DJ. Approximate Bayesian computation in population genetics. Genetics. 2002;162(4):2025-35

  19. [31]

    Approximate Bayesian computation (ABC) in practice

    Csilléry K, Blum MG, Gaggiotti OE, François O. Approximate Bayesian computation (ABC) in practice. Trends Ecol Evol. 2010;25(7):410-8

  20. [32]

    Population growth of human Y chromosomes: a study of Y chromosome microsatellites

    Pritchard JK, Seielstad MT, Perez-Lezaun A, Feldman MW. Population growth of human Y chromosomes: a study of Y chromosome microsatellites. Mol Biol Evol. 1999;16(12):1791-8

  21. [33]

    Differential elimination for dynamical models via projections with applications to structural identifiability

    Dong R, Goodbrake C, Harrington HA, Pogudin G. Differential elimination for dynamical models via projections with applications to structural identifiability. SIAM Journal on Applied Algebra and Geometry. 2023;7(1):194-235

  22. [34]

    A tutorial on structural identifiability of epidemic models using StructuralIdentifiability

    Liyanage YR, Saucedo O, Tuncer N, Chowell G. A tutorial on structural identifiability of epidemic models using StructuralIdentifiability. jl. arXiv preprint arXiv:250510517. 2025