{"id":"cb86581c-931c-4fad-96f4-9f92cab53844","arxiv_id":"2507.08895","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":8,"one_line_summary":"The paper presents an optimal control analysis of a rabies model with environmental transmission, but the effective reproduction number omits the environment and contradicts the parameter table.","lead":"The paper adds four time-dependent control measures to the authors' earlier rabies model for humans, free-ranging dogs, domestic dogs, and an environmental virus compartment, then uses optimal control to argue that vaccination and post-exposure treatment can eliminate rabies within five years. It is a standard optimal control extension, but the central reproduction number is miscomputed and numerically inconsistent with the reported parameters.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"With the paper's own fitted parameters, the free-ranging dog subsystem has R21 ≈ 1.39 > 1 and no control acts on it, so the claimed five-year elimination is inconsistent with the model unless different parameters or a different R_e formula were used.","rationale":"The most consequential check on the central claim is not the technical gap in the uniqueness proof or even the missing environmental terms in the F matrix, although both are real problems. It is whether the model can eliminate rabies at all when the fitted parameters are used. The free-ranging dog compartment is the only one with no direct control, and its isolated reproduction number R21 is > 1 under the fitted values in Table 2. Since R21 is a nonnegative additive component of R_e in (16), this also explains why the reported R_e < 1 contour plots cannot have been generated from the stated formula and fitted parameters. A simulation of the closed free-range subsystem would settle the point: if IF persists, then the five-year elimination curves in Figures 9–11 are not solutions of (1) with Table 2 parameters. This is a concrete, falsifiable check against the paper's own equations and data. I agree with the reader's weakest assumption and with the REJECT verdict; the concern is load-bearing because it targets the very mechanism by which the paper claims elimination, and it uses the paper's own fitted parameter values rather than an external modeling assumption.","tokens_in":28323,"tokens_out":6002,"duration_ms":70729,"concrete_test":"Set ID = ED = RD = M = 0 and integrate only the free-range subsystem of (1) with the fitted Table 2 parameters, a small positive IF(0), and controls set to u1 = u2 = u3 = u4 = 1 (which do not appear in these equations). Run for 20 years and check whether IF(t) decays to zero or approaches a positive endemic level. If IF(t) persists, the five-year elimination claim is contradicted by the model's own equations and fitted parameters. Equivalently, computing R21 = κ1θ2γ/[μ2(μ2+γ)(σ2+μ2)] from Table 2 gives approximately 1.39, and no control variable enters this expression.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the integrated strategy (u1–u4) eliminates rabies within five years. In model (1), controls appear only in the human and domestic-dog equations; the free-ranging dog equations for SF, EF, IF contain no uj. Consequently, IF can only be suppressed indirectly, through reduced ID and M. But the free-range subsystem alone has the reproduction number R21 = κ1θ2γ/[μ2(μ2+γ)(σ2+μ2)], defined in §3.2. Substituting the fitted values from Table 2 (θ2 = 1004.12, κ1 = 2.0×10^-5, γ = 0.1664, μ2 = 0.0663, σ2 = 0.0896) gives R21 ≈ 1.39 > 1. Since R_e in equation (16) is R21 plus nonnegative terms, the paper's own threshold also exceeds 1 under these parameters, yet Figures 5–7 plot R_e < 1. If R21 > 1, the free-range dog infection is self-sustaining: even with ID, ED, M and all human/domestic compartments at zero, a small IF(0) grows back to a positive endemic level and reseeds domestic dogs and humans through τ1, ψ1, and ν2. No control in the model acts on that subsystem. Thus the five-year elimination shown in Figures 9–11 cannot follow from system (1) with the parameters used for fitting; either the simulations used different parameters, an altered R_e formula, or the claimed DFE stability threshold is not the one being simulated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a deterministic 12-compartment ODE model for rabies transmission among humans, free-ranging dogs, domestic dogs, and an environmental virus compartment, with four time-dependent controls (health practice, domestic-dog vaccination, community education, and post-exposure prophylaxis). The authors derive a disease-free equilibrium, define an effective reproduction number R_e, state global stability results for the disease-free and endemic equilibria, perform a PRCC sensitivity analysis, solve an optimal control problem with Pontryagin's maximum principle, fit the model to reported human rabies incidence in Tanzania, and simulate four control strategies. The central claim is that an integrated strategy, especially vaccination of domestic dogs and PEP, can eliminate rabies within five years.","tokens_in":28696,"tokens_out":11188,"duration_ms":113555,"significance":"If the model and its quantitative predictions were correct, the paper would provide useful evidence for combining dog vaccination, education, and PEP in rabies control programs, and the inclusion of an environmental virus compartment addresses a genuinely underexplored transmission route. The authors also provide a standard optimal-control formulation and a data-fitting exercise. However, the central quantitative claims are undermined by internal inconsistencies: the reproduction number omits the environmental route that is a focus of the paper, the analytic R_e formulas contain incorrect pathway compositions, and the reported figures are inconsistent with the paper's own fitted parameters. The five-year elimination claim is not supported by the stated model equations.","major_comments":[{"comment":"The next-generation matrix F in Eq. (14) has zero entries in every column corresponding to the environmental compartment M, despite the fact that the incidence terms in system (1) contain τ3λ(M), κ3λ(M), and ψ3λ(M)/(1+ρ3). Since λ'(0)=1/C ≠ 0, the partial derivatives ∂F_i/∂M at the disease-free equilibrium are nonzero. The environmental transmission route is therefore absent from R_e, which invalidates the threshold condition R_e<1 as a statement about interruption of transmission through the environment emphasized in the title and conclusions.","section":"Sec. 3.2, Eq. (14)"},{"comment":"The formulas for R23, R31, and R33 in Eq. (16) are incorrectly composed. The path from infected domestic dogs to infected free-ranging dogs is ID→SF→EF→IF, so R23 should carry the factor γ/[(μ2+γ)(σ2+μ2)], not the domestic-dog factor a3=γ/[(μ3+γ1+γ2+u4)(σ3+μ3)]. Similarly, the paths from IF and ID to ID/H-domestic infection are IF→SD→ED→ID and ID→SD→ED→ID, so R31 and R33 should carry γ1/[(μ3+γ1+γ2+u4)(σ3+μ3)], not the free-range progression rate γ. As written, Eq. (16) is not the spectral radius of F V^{-1} even for the dog-only subsystem.","section":"Sec. 3.2, Eq. (16)"},{"comment":"Using the fitted values in Table 2 (κ1=2.0×10^-5, θ2=1004.12, γ=0.1664, μ2=0.0663, σ2=0.0896) gives R21≈1.39>1. Because Eq. (16) expresses R_e as R21 plus nonnegative terms, the model with Table 2 parameters implies R_e>1 at every control level, and u2 and u4 do not appear in R21 at all. Yet Figures 5–7 plot R_e between approximately 0.72 and 0.9. These figures therefore cannot have been generated from Eq. (16) with the parameters stated in Table 2; the paper does not disclose which parameters or formula were actually used.","section":"Sec. 4.1, Table 2 and Sec. 4.2, Figs. 5–7"},{"comment":"The claim that the integrated strategy eliminates rabies within five years is inconsistent with system (1) and Table 2. Controls u1–u4 appear only in the human and domestic-dog equations; the equations for S_F, E_F, and I_F contain no control term. With R21≈1.39>1, the free-ranging dog subsystem is self-sustaining: even if domestic dogs, humans, and the environment are all set to zero, a small positive I_F(0) grows back to a positive endemic level and reseeds the other compartments through τ1, ψ1, and ν2. No control in the model acts on that subsystem, so the elimination trajectories in Figures 9–11 cannot be solutions of (1) with the fitted parameters unless the simulations used different parameters or a modified model.","section":"Sec. 4.4, Strategy A, Figs. 9–11"},{"comment":"The parameter estimates that drive all later simulations are obtained by fitting the single equation (30) for I_H to data labelled 'Cumulative incidence of rabies' in Figure 4, but I_H in system (1) is an infected-state variable, not a cumulative incidence. No goodness-of-fit statistic, residual analysis, or uncertainty interval is reported, and the full model (1) is not fitted. The fitted parameters therefore provide a weak empirical basis for the quantitative predictions in Section 4.4.","section":"Sec. 4.1, Eqs. (30)–(31) and Fig. 4"}],"minor_comments":[{"comment":"The uniqueness proof is invalid: δu_j = \\bar u_j - \\tilde u_j may be negative, so the lower bound 0 ≤ 2A_j \\tilde u_j δu_j in Eq. (27) is false; the expansion of J also has the wrong coefficients (2A_j and A_j instead of A_j and A_j/2). The uniqueness claim is therefore not established.","section":"Sec. 3.7.1, Theorem 6"},{"comment":"The proof of Lemma 1 establishes Lipschitz continuity of the vector field, which gives local existence and uniqueness, but it does not prove positivity of solutions; a standard invariant-region argument for the nonnegative orthant is missing.","section":"Sec. 3.1, Lemma 1"},{"comment":"In the Hamiltonian, the cost term for u4 is written as (1/2)u4^2 A3; it should be (1/2)u4^2 A4. Also, the V matrix in Eq. (15) uses δ3 where the model and Table 2 use σ3.","section":"Sec. 3.7.2, Eq. (28)"},{"comment":"The objective functional contains the term -K6 S_D, and the text states that K6 is the weight for vaccinated domestic dogs, but S_D is the susceptible domestic-dog compartment. The modeling intention should be clarified.","section":"Sec. 3.6, Eq. (22)"},{"comment":"The expression for E_H^* in Eq. (18) is independent of the transmission rates λ_H, τ1, τ2, and τ3, which is not plausible for an endemic equilibrium; the formula also appears dimensionally inconsistent. This part of the endemic-equilibrium analysis needs correction.","section":"Sec. 3.4, Eq. (18)"},{"comment":"The data availability statement says the data came from existing literature, but Figure 4 presents Tanzania incidence data from 1990–2018 without a specific source or citation; please identify the dataset.","section":"Data Availability"}],"recommendation":"reject","confidential_remarks":"The central inconsistency between Table 2 (which implies R21≈1.39 and therefore R_e>1) and Figures 5–7 (which plot R_e<1) is easily checkable and is not a presentation issue: it invalidates the paper's main quantitative claims. The five-year elimination result in Figures 9–11 similarly cannot be reproduced from the stated model with the stated parameters. These are load-bearing errors, not local gaps, so I recommend rejection. If the authors supply the simulation code and a corrected R_e derivation, a resubmission could be considered."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — quick take on 2507.08895. The novelty is modest: it takes the authors' earlier rabies model [3], adds four time-dependent controls (health practice, domestic dog vaccination, education, PEP), applies Pontryagin's principle, and runs LHS/PRCC sensitivity plus Matlab simulations. That is a legitimate exercise; the conclusion that vaccination plus PEP matters is not new, but a calibrated optimal-control version for Tanzania could be useful.\n\nWhere it falls apart: Eq. (16) for Re has no terms involving the environmental transmission rates tau3, kappa3, psi3 or the shedding rates nu_i. The F matrix in Eq. (14) has zeros in the M column, so the 'environmental effects' in the title do not enter the threshold at all. That is a serious omission for a model whose stated novelty is environmental transmission. Worse, the stress-test note is correct. With the fitted values in Table 2, R21 is about 1.39, so Re > 1, yet Figures 5–7 plot Re < 1. The baseline model simulations show persistent infection, which is consistent with Re>1 but not with the claimed stable disease-free equilibrium. Since no control acts directly on the free-ranging dog subsystem and R21>1, the five-year elimination shown in Figures 9–11 cannot follow from the model with the estimated parameters. Either the simulations used a different parameter set or a different Re formula, and the paper does not say which.\n\nOther soft spots: Theorem 6's uniqueness proof is invalid—it bounds a nonnegative integral by an integral of absolute values and concludes zero, ignoring that delta u_j can be negative. The global stability proof for the endemic equilibrium in Appendix A has variable typos, hand-wavy 'set to zero' steps, and does not actually establish that dH/dt <= 0. The model is also fitted to the same Tanzania incidence series used to motivate the study, so the numerical 'elimination in five years' is partly a restatement of the fitting objective; that circularity is real but secondary.\n\nCredit where it is earned: positivity and boundedness are handled cleanly, the parameter estimation is plausibly described, and the optimal control characterization follows the standard Pontryagin template. The problem is not the framework; it is that the paper's own numbers contradict its headline threshold and elimination claim.\n\nWho this is for: someone interested in the structure of an optimal-control rabies model might read the formulation, but the quantitative claims should not be relied on. I would not cite it for the five-year result. It is worth sending to a serious referee—fixing Re and rerunning the simulations is feasible—but as it stands this needs major revision and re-verification before any policy conclusion is drawn.","headline":"The optimal-control framework is standard and the model builds on their earlier work, but the five-year elimination claim is contradicted by their own fitted parameters and the reproduction number drops the environmental route entirely.","tokens_in":29232,"tokens_out":3878,"would_cite":false,"duration_ms":46683,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["92D30","49K15","34D23"],"pacs":[],"model":"deepseek-v4-flash","headline":"Vaccinating domestic dogs and giving bite victims post-exposure prophylaxis can eliminate rabies within five years, according to an optimal control model fitted to Tanzanian data.","keywords":["rabies","optimal control theory","effective reproduction number","dog vaccination","post-exposure prophylaxis","free-ranging dogs","environmental transmission","Tanzania"],"falsifier":"Compute the reproduction number of the free-ranging dog subsystem alone, $R_{21}$, from the fitted parameters; if $R_{21}>1$, run the optimal-control simulation with only domestic-dog vaccination and human PEP active. If infected free-ranging dogs and environmental virus do not decay to zero, the five-year elimination claim is contradicted.","tokens_in":28139,"feed_emoji":"🐕","tokens_out":8982,"duration_ms":93795,"temperature":0.7,"pith_summary":"This paper tries to show that a coordinated package of time-dependent interventions can break the rabies transmission cycle among humans, domestic dogs, free-ranging dogs, and the environment. The authors extend a previously published deterministic model by adding four controls: health promotion and surveillance ($u_1$), domestic-dog vaccination ($u_2$), community education ($u_3$), and post-exposure prophylaxis (PEP) for exposed individuals ($u_4$). After fitting the uncontrolled system to reported human rabies incidence in Tanzania from 1990 to 2018, they compute an effective reproduction number $R_e$ and show numerically that increased control lowers it. The central result is that the integrated optimal strategy, with vaccination and PEP held near their upper limits in the early years, drives all infected human and dog classes to zero within five years. If the model is right, a practicable combination of dog vaccination, education, and prompt treatment could meet the 2030 goal of eliminating dog-mediated rabies.","feed_headline":"Dog vaccines plus bite-victim treatment can end rabies in five years","feed_subtitle":"An optimal-control model fitted to Tanzanian data shows the combined strategy pushes the reproduction number below 1.","key_machinery":"The load-bearing object is a 12-compartment deterministic ODE system: susceptible, exposed, infected, and recovered classes for humans, free-ranging dogs, and domestic dogs, plus a state $M$ for rabies virus in the environment with saturating contamination $\\lambda(M)=M/(M+C)$. Controls enter multiplicatively: $(1-u_1-u_3)$ scales human exposure, $(1-u_1-u_2)$ scales domestic-dog exposure, and PEP ($u_4$) moves exposed humans and dogs into recovered classes. Two tools carry the argument: the next-generation matrix gives the reproduction number $R_e$ that sets the stability threshold, and Pontryagin's maximum principle gives explicit characterizations of the four optimal controls in terms of adjoint variables, solved numerically by a forward-backward fourth-order Runge–Kutta scheme.","core_discovery":"On the paper's own terms, the discovery is that domestic-dog vaccination and post-exposure prophylaxis are the two controls that carry the elimination, and that they work without a direct control on free-ranging dogs. The effective reproduction number takes the form $R_e = (R_{33}+R_{21}+\\sqrt{R_{21}^2-2R_{33}R_{21}+4R_{31}R_{23}+R_{33}^2})/2$, where $R_{21}$ is free-ranging dog transmission, $R_{33}$ is domestic-dog transmission, and $R_{31}, R_{23}$ are cross-transmission terms; the dog terms are multiplied by $(1-u_1-u_2)$ and the human terms by $(1-u_1-u_3)$, so controls enter the threshold directly. The paper proves that the rabies-free equilibrium is globally stable when $R_e<1$ and that the endemic equilibrium is globally stable when $R_e>1$. In simulations with the fitted Tanzanian parameters, the optimal control problem keeps $u_2$ and $u_4$ at 100% for the first several to sixteen years and gradually reduces them, and all exposed and infected compartments, including free-ranging dogs and environmental virus, fall to near zero by year five.","pith_inferences":["If the paper's fitted parameters leave the free-ranging dog subsystem with $R_{21}>1$, the five-year result rests on the cross-species and environmental links being interrupted indirectly; a closed free-ranging-population variant of the model is the natural stress test.","Because the paper reports a real mass-vaccination campaign at $1.61 per dog, its quadratic cost weights could be replaced by unit costs to turn the elimination claim into a budget projection.","The same control framework could be applied to imperfect PEP uptake, which the paper notes leaves 20–30% of bite victims untreated, to see how much the five-year horizon stretches.","Setting the environmental shedding rates $\\nu_1=\\nu_2=\\nu_3=0$ would isolate how much of the result depends on the virus reservoir $M$ rather than on direct bites."],"forward_implications":["If $R_e<1$ is achieved by the optimal controls, the disease-free equilibrium is globally stable, so once rabies is cleared it does not rebound while the controls are sustained.","Vaccination of domestic dogs ($u_2$) and post-exposure prophylaxis ($u_4$) are the controls held at 100% the longest, making them the main drivers of the five-year elimination.","The integrated strategy also clears the free-ranging dog and environmental compartments, even though no term in the model directly reduces those infections.","Higher dog-to-dog contact rates ($\\psi_1$, $\\psi_2$) raise $R_e$, so the same package must run longer or harder in high-contact settings.","Deterrence factors ($\\rho_1$, $\\rho_2$, $\\rho_3$) lower the exposed and infected domestic-dog counts, so discouraging contact complements vaccination."],"supporting_citations":[{"why":"Supplies the base three-host model that this paper extends with time-dependent controls.","marker":"[3]"},{"why":"Supplies the Metzler-matrix method used to prove global stability of the disease-free equilibrium.","marker":"[18]"},{"why":"Supply the existence and uniqueness results for the optimal control problem.","marker":"[19, 20]"},{"why":"Supplies Pontryagin's maximum principle used to characterize the four optimal controls.","marker":"[24]"},{"why":"Supplies the forward-backward numerical algorithm used to solve the optimality system.","marker":"[25]"},{"why":"Supplies baseline transmission-rate values used in the Tanzania data fitting.","marker":"[26]"},{"why":"Supplies parameter values and the SEIR-style rabies model structure the fitted model draws on.","marker":"[27]"}],"fun_headline_variants":["Vaccinate dogs, treat bites: rabies gone in 5 years","Optimal control: dog vaccines and PEP eliminate rabies by year 5","Rabies elimination: vaccines and post-exposure therapy do it in 5 years","Key to rabies: combine dog vaccination with bite treatment","Model shows dog shots plus PEP eradicate rabies in 5 years"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The elimination result assumes that free-ranging dogs can be cleared indirectly by controlling domestic dogs and humans, even though no control acts directly on the free-ranging dog transmission terms.","fun_headline_variants_meta":{"raw":{"variants":["Vaccinate dogs, treat bites: rabies gone in 5 years","Optimal control: dog vaccines and PEP eliminate rabies by year 5","Rabies elimination: vaccines and post-exposure therapy do it in 5 years","Key to rabies: combine dog vaccination with bite treatment","Model shows dog shots plus PEP eradicate rabies in 5 years"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000194,"raw_usage":{"total_tokens":1342,"prompt_tokens":923,"completion_tokens":419,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":539,"completion_tokens_details":{"reasoning_tokens":327}},"tokens_in":539,"tokens_out":419,"duration_ms":4518,"temperature":1.0,"reasoning_tokens":327,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:26:39.966889+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the reproduction number of the free-ranging dog subsystem alone, $R_{21}$, from the fitted parameters; if $R_{21}>1$, run the optimal-control simulation with only domestic-dog vaccination and human PEP active. If infected free-ranging dogs and environmental virus do not decay to zero, the five-year elimination claim is contradicted.","supporting_citations":[{"cited_title":"Charles, V","cited_arxiv_id":null,"evidence_quote":"Supplies the base three-host model that this paper extends with time-dependent controls."},{"cited_title":"Castillo-Chavez, Z","cited_arxiv_id":null,"evidence_quote":"Supplies the Metzler-matrix method used to prove global stability of the disease-free equilibrium."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies Pontryagin's maximum principle used to characterize the four optimal controls."},{"cited_title":"Campos, C","cited_arxiv_id":null,"evidence_quote":"Supplies the forward-backward numerical algorithm used to solve the optimality system."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies baseline transmission-rate values used in the Tanzania data fitting."},{"cited_title":"Zhang, Z","cited_arxiv_id":null,"evidence_quote":"Supplies parameter values and the SEIR-style rabies model structure the fitted model draws on."}],"review_version":1}