A decay-adjusted spatio-temporal model to account for the impact of mass drug administration on neglected tropical disease prevalence
Pith reviewed 2026-05-17 02:07 UTC · model grok-4.3
The pith
A decay-adjusted spatio-temporal model accounts for the time-varying effects of mass drug administration on neglected tropical disease prevalence from sparse surveys.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors propose the decay-adjusted spatio-temporal (DAST) model that explicitly accounts for the time-varying impact of MDA on NTD prevalence. Using case studies on soil-transmitted helminths and lymphatic filariasis, they demonstrate that DAST offers a practical alternative to standard geostatistical models for quantifying MDA impact and supporting short-term programmatic forecasting. The model is flexible and interpretable for estimating intervention effects from sparse survey data, while discussing extensions and identifiability challenges and advocating data-driven parsimony.
What carries the argument
The decay-adjusted spatio-temporal (DAST) model, which modifies a standard geostatistical framework by adding a time-varying decay component that represents the diminishing impact of mass drug administration on disease prevalence.
If this is right
- The DAST model enables quantification of the impact of mass drug administration on neglected tropical disease prevalence.
- It supports short-term programmatic forecasting to guide control decisions.
- It provides a more interpretable alternative to standard models when survey data are sparse.
- The framework highlights the need for parsimonious models in data-limited settings to avoid over-parameterization.
Where Pith is reading between the lines
- The same decay-adjustment idea could apply to monitoring other interventions with temporary effects, such as vaccination campaigns.
- Adding environmental or demographic covariates to the decay component might further improve estimates in future applications.
- Longer sequences of repeated surveys across multiple sites would provide a direct test of whether the model generalizes beyond the presented case studies.
Load-bearing premise
The decay adjustment for mass drug administration effects can be reliably identified and estimated from the available sparse prevalence survey data without major identifiability problems.
What would settle it
New prevalence surveys conducted after further rounds of mass drug administration could be compared to the model's short-term forecasts to check whether the predicted decay in disease rates matches the observed changes.
Figures
read the original abstract
Prevalence surveys are routinely used to monitor the effectiveness of mass drug administration (MDA) programmes for controlling neglected tropical diseases (NTDs). We propose a decay-adjusted spatio-temporal (DAST) model that explicitly accounts for the time-varying impact of MDA on NTD prevalence, providing a flexible and interpretable framework for estimating intervention effects from sparse survey data. Using case studies on soil-transmitted helminths and lymphatic filariasis, we show that DAST offers a practical alternative to standard geostatistical models when the objective includes quantifying MDA impact and supporting short-term programmatic forecasting. We also discuss extensions and identifiability challenges, advocating for data-driven parsimony over complexity in settings where the available data are too sparse to support the estimation of highly parameterised models.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a decay-adjusted spatio-temporal (DAST) model to explicitly incorporate the time-varying effects of mass drug administration (MDA) on neglected tropical disease (NTD) prevalence. It presents the model as a flexible, interpretable alternative to standard geostatistical approaches for estimating intervention impacts and supporting short-term forecasting from sparse prevalence surveys. Case studies on soil-transmitted helminths and lymphatic filariasis are used to illustrate performance, while the discussion addresses identifiability challenges and advocates data-driven parsimony over complex parameterizations.
Significance. If the decay adjustment can be shown to be reliably estimable, the framework would provide a useful tool for NTD control programs by separating MDA effects from underlying spatio-temporal prevalence trends. The explicit focus on parsimony when data are sparse is a constructive contribution, but the significance is limited by the need to demonstrate that the added decay component yields identifiable, non-degenerate estimates rather than collapsing into a standard smoother.
major comments (2)
- [Abstract] Abstract: The claim that DAST 'offers a practical alternative' for quantifying MDA impact rests on the case studies, yet the abstract itself flags identifiability challenges without indicating how the decay parameters were shown to remain estimable (with quantified uncertainty) at the observation densities typical of NTD monitoring programs.
- [Case studies] Case studies section: To support the central claim that the decay term can be separated from the spatio-temporal field, the results should include explicit diagnostics (e.g., posterior intervals for decay parameters, or fits under data subsampling) demonstrating that the adjustment does not reduce to a constant or become unidentified when survey points are thinned to realistic sparsity levels.
minor comments (2)
- [Methods] Clarify the exact functional form of the decay adjustment (e.g., exponential, power-law) and its interaction with the spatio-temporal random effects in the model equations.
- [Results] Add a table comparing key metrics (e.g., WAIC, predictive scores) between DAST and the baseline geostatistical model for both case studies.
Simulated Author's Rebuttal
We thank the referee for their detailed and constructive report. Their comments highlight important aspects of model identifiability that we address below. We have revised the manuscript to strengthen the presentation of estimability evidence while preserving the focus on parsimony for sparse NTD data.
read point-by-point responses
-
Referee: [Abstract] Abstract: The claim that DAST 'offers a practical alternative' for quantifying MDA impact rests on the case studies, yet the abstract itself flags identifiability challenges without indicating how the decay parameters were shown to remain estimable (with quantified uncertainty) at the observation densities typical of NTD monitoring programs.
Authors: We agree that the abstract would benefit from a concise indication of how estimability was assessed. In the revised manuscript we have updated the abstract to note that posterior credible intervals for the decay parameters remained well-defined and excluded degeneracy in both case studies at the observed survey densities, consistent with the data-driven parsimony approach discussed in the main text. revision: yes
-
Referee: [Case studies] Case studies section: To support the central claim that the decay term can be separated from the spatio-temporal field, the results should include explicit diagnostics (e.g., posterior intervals for decay parameters, or fits under data subsampling) demonstrating that the adjustment does not reduce to a constant or become unidentified when survey points are thinned to realistic sparsity levels.
Authors: We concur that additional explicit diagnostics would reinforce the separation of the decay term. The original manuscript already reports posterior intervals for the decay parameters in the case-study results; to directly address the referee's suggestion we have added a new supplementary analysis that applies progressive data thinning to realistic NTD survey sparsity levels. This analysis shows that the posterior means and intervals for the decay rate remain stable and do not collapse to a constant or become unidentified. revision: yes
Circularity Check
No significant circularity in derivation chain
full rationale
The paper proposes a decay-adjusted spatio-temporal model as an extension of standard geostatistical approaches for NTD prevalence, with explicit discussion of identifiability challenges and advocacy for parsimony when data are sparse. No load-bearing step reduces a claimed prediction or result to its own inputs by construction, self-citation, or fitted parameter renaming. The central framework is defined independently via the added decay term and remains falsifiable against external survey data without requiring the target estimates as inputs. This is the expected honest non-finding for a model-building paper that flags its own limitations.
Axiom & Free-Parameter Ledger
free parameters (1)
- decay adjustment parameters
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
f(v) = α exp(−v/γ), v > 0 (eq. 10); P(x,t) = P*(x) ∏ [1−f(t−uj)] I(x,uj) (eq. 8); penalised logL with g_pen(α) = −(λ1 log α + λ2 log(1−α)) (eq. 17)
-
IndisputableMonolith/Foundation/DimensionForcing.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
separable double-exponential covariance; Matérn with exponential limit; MCML + Laplace sampling for latent GP
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
doi: 10.1186/s40249-018-0420-9. O. Christensen, G. Roberts, and M. Sk¨ old. Robust markov chain monte carlo methods for spatial generalized linear mixed models.Journal of Computational and Graphical Statistics, 15(1):1–17,
-
[2]
doi: 10.1093/cid/ciw506. P. J. Diggle and E. Giorgi. Model-based geostatistics for prevalence mapping in low-resource settings.Journal of the American Statistical Association, 111(515):1096–1120,
-
[3]
doi: 10.1080/01621459.2015.1123158. P. J. Diggle, J. A. Tawn, and R. A. Moyeed. Model-based geostatistics.Journal of the Royal Statistical Society: Series C (Applied Statistics), 47(3):299–350,
-
[4]
doi: 10.1111/1467-9876.00113. B. Efron. Bootstrap methods: Another look at the jackknife.The Annals of Statistics, 7(1):1–26,
-
[5]
The Annals of Statistics , author =
doi: 10.1214/aos/1176344552. URLhttps://projecteuclid.org/journals/annals-of-statistics/volume-7/issue-1/ Bootstrap-methods-another-look-at-the-jackknife/10.1214/aos/1176344552.full. D. Firth. Bias reduction of maximum likelihood estimates.Biometrika, 80(1):27–38,
-
[6]
doi: 10.1093/biomet/ 80.1.27. A. Garchitorena, E. M. Raza-Fanomezanjanahary, S. A. Mioramalala, C. B. Chesnais, C. A. Ratsimbasoa, H. Ramarosata, M. H. Bonds, and H. Rabenantoandro. Towards elimination of lymphatic filariasis in southeastern madagascar: Successes and challenges for interrupting transmission.PLOS Neglected Tropical Diseases, 12(9):e0006780,
-
[7]
doi: 10.1371/journal.pntd.0006780. C. J. Geyer. Likelihood and exponential families.Department of Statistics, University of Minnesota,
-
[8]
doi: https://doi.org/10.1111/insr.12268. E. Giorgi, C. Fronterr` e, P. M. Macharia, V. A. Alegana, R. W. Snow, and P. J. Diggle. Model building and assessment of the impact of covariates for disease prevalence mapping in low-resource settings: to explain and to predict.Journal of the Royal Society Interface, 18(178):20210104,
-
[9]
doi: 10.1098/rsif.2021.0104. T. Gneiting. Nonseparable, stationary covariance functions for space–time data.Journal of the American Statistical Association, 97(458):590–600,
-
[10]
doi: 10.1371/journal.pmed.1000324. G. S. A. Gunawardena, N. D. Karunaweera, and M. M. Ismail. Soil-transmitted helminth infections among plantation sector schoolchildren in sri lanka: prevalence after ten years of preventive chemotherapy.PLOS Neglected Tropical Diseases, 5(9):e1341,
-
[11]
doi: 10.1371/journal.pntd.0001341. S. I. Hay. Global mapping of infectious disease.Philosophical Transactions of the Royal Society B: Biological Sciences, 368(1623):20120250,
-
[12]
doi: 10.1098/rstb.2012.0250. T. D. Hollingsworth, R. L. Pullan, D. Gurarie, J. E. Truscott, S. J. Brooker, and R. M. Anderson. Elimination of transmission of soil-transmitted helminths: a revised breakpoint analysis accounting for density dependence and stochasticity.Journal of Theoretical Biology, 487:110102,
-
[13]
doi: 10.1016/j.jtbi.2019.110102. D. A. Karagiannis-Voules, M. B¨ urgi, N. Schur, T. B¨ arnighausen, J. Utzinger, and P. Vounatsou. Mapping the global distribution of trachoma: the global trachoma mapping project.PLOS Neglected Tropical Diseases, 15 (2):e0009233,
-
[14]
doi: 10.1371/journal.pntd.0012782. URLhttps://doi.org/10. 1371/journal.pntd.0012782. C. Linard, M. Gilbert, R. W. Snow, A. M. Noor, and A. J. Tatem. Population distribution, settlement patterns and accessibility across africa in 2010.PLOS ONE, 7(2):e31743,
-
[15]
doi: 10.1371/journal.pone.0031743. A. Loukouri, A. M´ e¨ ıt´ e, B. G. Koudou, C. W. Goss, D. Lew, G. J. Weil, E. K. N’Goran, and P. U. Fischer. Impact of annual and semi-annual mass drug administration for lymphatic filariasis and onchocerciasis on hookworm infection in cˆ ote d’ivoire.PLOS Neglected Tropical Diseases, 15(3):e0009250,
-
[16]
doi: 10.1371/journal.pntd. 0009250. B. Mat´ ern.Spatial variation: Stochastic models and their application to some problems in forest surveys and other sampling investigations, volume 49 ofMeddelanden fran Statens Skogsforskningsinstitut. Almqvist & Wiksell, Stockholm,
-
[17]
doi: 10.1371/journal.pone.0139779. K. Mohammed, O. Debrah, D. Batsa, and et al. Feasibility and safety of mass drug coadministration with ivermectin, albendazole, and azithromycin for the integrated treatment of neglected tropical diseases: a single- arm intervention trial.The Lancet Global Health, 6(10):e1132–e1138,
-
[18]
URLhttps: //www.bmj.com/content/358/bmj.j4307
doi: 10.1136/bmj.j4307. URLhttps: //www.bmj.com/content/358/bmj.j4307. C. Mwandawiro, C. Okoyo, J. Kihara, E. Simiyu, S. Kepha, S. J. Campbell, M. C. Freeman, S. J. Brooker, and S. M. Njenga. Results of a national school-based deworming programme on soil-transmitted helminths infections and schistosomiasis in kenya: 2012–2017.Parasites & Vectors, 12(76),
-
[19]
URLhttps://parasitesandvectors.biomedcentral.com/articles/10.1186/s13071-019-3322-1
doi: 10.1186/s13071-019-3322-1. URLhttps://parasitesandvectors.biomedcentral.com/articles/10.1186/s13071-019-3322-1. B. M. Psaty, T. D. Koepsell, D. Lin, N. S. Weiss, D. S. Siscovick, F. R. Rosendaal, M. Pahor, and C. D. Furberg. Assessment and control for confounding by indication in observational studies.Journal of the American Geriatrics Society, 47(6)...
-
[20]
doi: 10.1111/j.1532-5415.1999.tb01603.x. R. L. Pullan, J. L. Smith, R. Jasrasaria, and S. J. Brooker. Global numbers of infection and disease burden of soil transmitted helminth infections in 2010.Parasites & vectors, 7(1):37,
-
[21]
doi: 10.1214/ss/1015346320. 32 Giorgi, E., Fronterre, C. and Diggle, P. J. M. Sasanami, B. Amoah, A. Amza, A. S. Y. Souley, A. N. Diori, A. Bakhtiari, B. Kadri, C. L. Szwarcwald, D. V. F. Gomez, I. Almou, M. F. C. Lopes, M. P. Masika, N. Beidou, S. Boyd, E. M. Harding-Esch, A. W. Solomon, and E. Giorgi. Using model-based geostatistics for assessing the el...
-
[22]
doi: 10.1371/journal.pntd.0011476. C. A. Schmidt, E. A. Cromwell, and et al. Prevalence of onchocerciasis in africa 2000–2018: A geospatial modelling study.The Lancet Microbe, 3(7):e431–e440,
-
[23]
doi: 10.1016/S2666-5247(22)00101-8. J. Senecal, A. Nordin, and B. Vinner˚ as. Fate of Ascaris at various ph, temperature and moisture levels.Journal of Water and Health, 18(3):375–382,
-
[24]
doi: 10.2166/wh.2020.264. A. Srivathsan, A. Abdou, T. Al-Khatib, S.-C. Apadinuwe, M. D. Badiane, V. Bucumi, T. Chisenga, G. Kabona, M. Kabore, S. K. Kanyi, et al. District-level forecast of achieving trachoma elimination as a public health problem by 2030: An ensemble modelling approach.Clinical Infectious Diseases, 78(Supplement 2):S101–S109,
-
[25]
doi: 10.1093/cid/ciad071. W. A. Stolk, Q. A. ten Bosch, S. J. de Vlas, P. U. Fischer, G. J. Weil, and A. S. Goldman. Modeling the impact and costs of semiannual mass drug administration for accelerated elimination of lymphatic filariasis.PLOS Neglected Tropical Diseases, 7(1):e1984,
-
[26]
doi: 10.1371/journal.pntd.0001984. W. A. Stolk, S. J. de Vlas, and J. D. F. Habbema. Are we on track to eliminate lymphatic filariasis?Parasitology, 142(8):1123–1136,
-
[27]
doi: 10.1016/bs.apar.2016.08.002. J. E. Truscott, M. Werkman, J. E. Wright, S. H. Farrell, R. Sarkar, J. Toor, and R. M. Anderson. Measuring the impact of deworming programmes by community-wide treatment coverage and worm intensity.Parasites & Vectors, 10(1):199,
-
[28]
doi: 10.1186/s13071-017-2122-y. J. Wakefield, T. Okonek, and J. Pedersen. Small area estimation for disease prevalence mapping.International Statistical Review, 88(2):398–418,
-
[29]
doi: 10.1111/insr.12400. M. Werkman, J. E. Truscott, J. Toor, J. E. Wright, and R. M. Anderson. Testing for soil-transmitted helminth transmission elimination: analysis of a new model-based framework for evaluating interruption of transmission through mass drug administration.PLOS Neglected Tropical Diseases, 12(5):e0006331,
-
[30]
doi: 10.1016/S1473-3099(19)30633-2. Monitoring and Epidemiological Assessment of Mass Drug Administration in the Global Programme to Eliminate Lymphatic Filariasis: A Manual for National Elimination Programmes. WHO, Geneva,
-
[31]
A decay-adjusted spatio-temporal model for prevalence mapping33 AscarisHookwormTrichuris 0.00 0.25 0.50 0.75 1.00 0.00 0.25 0.50 0.75 1.00 0.00 0.25 0.50 0.75 1.00 0.00 0.25 0.50 0.75 1.00 u Average non−Randomized Probability Integral Tranform (AnPIT) DAST S ST Fig. 10: Average non-randomised PIT (AnPIT) curves for the held-out samples for the soil transm...
work page 2004
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.