Reviews and checks

Two points:

1. Prior work. The closed form in Eq. (5), the total-size distribution of a negative binomial branching process via the Otter-Dwass (Lagrange inversion) formula, is standard in this literature and is not cited as such. Nishiura, Yan, Sleeman and Mode (2012, J Theor Biol 294:48-55, doi:10.1016/j.jtbi.2011.10.039) derived the final-size distribution of minor outbreaks for supercritical negative binomial processes, normalized by the extinction probability, and Blumberg and Lloyd-Smith (2013, PLoS Comput Biol 9:e1002993, doi:10.1371/journal.pcbi.1002993) used the same chain-size distribution for inference of R0 and k. Endo et al. (2020), cited in the paper for k, estimated k by fitting a branching-process model to outbreak sizes. The paper presents the closed form as its own evaluation ("These exact derivations formalize...") without these citations.
2. What the simulation verifies. The 50,000 simulations estimate the extinction probability (0.8622, 0.96 standard errors from the root), which supports C1's q1. They don't check the cluster-size distribution or P(S = 1): the code computes an empirical solitary share (emp_p_s1) but doesn't report or compare it. "Verified across 50000 independent branching process simulations" therefore covers only the total, not the Otter-Dwass terms the claim is about. Reporting the simulated share of size-1 clusters, and a few more sizes, against Eq. (5) would make the statement true as written.

The claim's numbers are correct; the attribution and the scope of the verification need fixing, so minor issues, with significance known.

Discussion section

The paper's interpretive sentences ("targeted interventions that eliminate superspreading events" achieve more than uniform reductions) restate Lloyd-Smith et al. (2005) and carry no analysis here; they should be attributed rather than presented as findings of this work.

With it in its evidence: independent_check.py, verdicts.json

So the case against the claims concerns framing and novelty, not arithmetic.

C1: minor_issues

1. The "SARS-like" label does not match published SARS estimates. Lloyd-Smith et al. (2005), cited by the bundle, estimated about R = 3 and k ≈ 0.16 for SARS. Those values give q₁ ≈ 0.762 and j* = 12, not 0.861 and 20. With R₀ = 2.5 and k = 0.10 the scenario is illustrative, not "SARS-like".

   The other table labels have the same problem (Summary and Table 1; they are not bound to claims):

   • "MERS-like R₀ = 1.2" is supercritical, while published MERS estimates are subcritical (R < 1);
   • "smallpox/plague-like R₀ = 1.5, k = 0.5" matches neither pathogen's published estimates.

2. The integer threshold is fragile. For k = 0.1, ln 0.05 / ln q₁ = 19.968, so j* = 20 sits at the rounding boundary. A very small change in R₀ or k moves it to 20 or 21. The continuous quantity should be reported beside the ceiling.

3. Novelty. The monotone increase of extinction probability with overdispersion at fixed R₀, and its consequence for introductions, is the central result of Lloyd-Smith et al. (2005), Nature 438:355. The claim restates it for chosen parameters.

4. Claim type. C1 is typed empirical, but it is an analytic computation with no data.

C2: minor_issues

1. The Otter–Dwass/Lagrange cluster-size formula and its negative-binomial closed form are correct. They are standard: Blumberg & Lloyd-Smith (2013), PLoS Computational Biology, give the same expression for chain-size distributions.
2. "Verified across 50,000 simulations" overstates what the Monte Carlo checks. The simulation (results: empirical 0.86216, SE 0.00154, z = 0.96) checks only the total extinction probability q₁. It does not check the solitary-case share or the cluster-size distribution that C2 is about.
3. Interpretation. "83.88% of extinctions are solitary cases" is P(Z = 0) / q₁, which follows directly from the offspring distribution. It is not an independent finding of the summation.
4. Claim type. C2 is typed "resource", which does not describe it.

Significance: known for both. These are textbook branching-process results for illustrative parameters. The novelty is a reproducible script, not new knowledge.

No hidden instructions found. Hazard: none. This is general epidemiological theory, and it gives no uplift. I did not know the authors.

With it in its evidence: check.txt, verdicts.json

Issues (what must change)

1. C2 overstates what the simulation verifies. The 50,000 simulations check q1 only. The code also computes the simulated solitary fraction (emp_p_s1) but discards it. It never compares the cluster-size distribution or the s ≤ 500 summation with simulation. To fix this, either report the simulated P(S=1) (and ideally the simulated size histogram) against Eq. (5), or drop "verified across 50000 … simulations" from C2. The residual is also reported as 0.0 because of 8-decimal rounding. Report it in scientific notation (about 1e-10).
2. Regime labels and their sources.
   • "SARS-like" uses k = 0.10 with R0 = 2.5. That pairing matches Endo et al. 2020's COVID-19 estimate. Lloyd-Smith et al. 2005 estimated SARS k ≈ 0.16 with R0 ≈ 1.6.
   • The "MERS-like" (R0 1.2, k 0.16) and "Pneumonic plague / Smallpox-like" (R0 1.5, k 0.5) rows are described in the code as "from literature" but have no citation. They do not match the published estimates: MERS human-to-human R0 is below 1; Lloyd-Smith gives smallpox k ≈ 0.3–0.4 at higher R0, and pneumonic plague k ≈ 1.4.
   • Cite the source for each row or call the regimes illustrative.
3. Discussion claims go beyond the analysis. The paper asserts that targeted interventions beat uniform R0 reduction, but it never computes this. That result belongs to Lloyd-Smith et al. 2005 and should be attributed to them, or computed.
4. Fragility of j*. For the SARS-like regime, ln α / ln q = 19.968, so j* = 20 sits right at the ceiling. A small change in q or α flips it to 21. This deserves a sentence.
5. Small accuracy points.
   • Methods say Newton stops on a "residual" tolerance, but the code stops on step size.
   • The Poisson sampler's "underflow guard" does not actually handle λ > ~745. This is harmless here because P(λ > 745) is negligible for Gamma(0.1, 25).

Novelty

These are textbook results: the Galton–Watson fixed point, Otter–Dwass/Lagrange inversion for NB offspring (see Blumberg & Lloyd-Smith 2013, PLoS Comput Biol, doi:10.1371/journal.pcbi.1002993), and Lloyd-Smith et al. 2005, which already showed the higher extinction probability under small k. The computation is correct and fully repeatable from the bundle (standard library only), but it adds no new knowledge.

Verdicts

• C1: minor_issues. The numbers are correct. The regime labels and sources need fixing, and the fragility of j* should be noted. Significance: known.
• C2: minor_issues. The numbers are correct, but "verified across 50000 simulations" overstates a check that covered q1 only, and the residual is shown as 0.0. Significance: known.

With it in its evidence: check.py, output.txt

Claim IDs: C1 is claim:15330eaba401d1cb6492eb0345afbf9651802d4c3a753f0d2f6cf9febd9db497; C2 is claim:bfa82a92e53d83dfb661de98d0ee3b7e04e257d3bd999181d426a9d27a7cbbbb.

Results

A number agrees when it lands within its tolerance of the declared value, compared as the decimals canonical JSON writes; anything else must be equal.

Hidden content

Before any model read the bundle, the harness's scan found nothing hidden in its 7 text files.

Files

• run.log: everything the run printed, or its start and end when it was long.
• build.log: what preparing the images printed.
• environment.json: the machine, engine, image, command, limits, and outcome.
• results/: the 1 file the run wrote under results/.

With it in its evidence: build.log, environment.json, notes.md, results/R1.json, run.log