pith. sign in

arxiv: 2512.03760 · v2 · submitted 2025-12-03 · 📊 stat.AP · stat.ME

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

classification 📊 stat.AP stat.ME
keywords spatio-temporal modelmass drug administrationneglected tropical diseasesprevalence surveysdecay adjustmentgeostatistical modelsintervention evaluationdisease mapping
0
0 comments X

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.

The paper proposes a decay-adjusted spatio-temporal model, called DAST, to explicitly incorporate the changing effects of mass drug administration on the prevalence of neglected tropical diseases. This allows for more accurate estimation of intervention impacts from limited survey data collected at different times and places. Standard geostatistical models often overlook these time-varying effects, which can distort assessments of program effectiveness. The authors illustrate the approach with examples from soil-transmitted helminths and lymphatic filariasis, showing its utility for short-term forecasting in public health programs. They emphasize using parsimonious models due to data limitations and identifiability issues.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2512.03760 by Claudio Fronterre, Emanuele Giorgi, Peter J. Diggle.

Figure 1
Figure 1. Figure 1: Schematic representation of a spatially confounded survey design and its effects. Panel A: spatial [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Results from a simulated data-set using the model outlined in Section 2.2. Estimated and true decay [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Illustration of the expected effect of one round of mass drug administration on disease prevalence [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Mean standardised bias sBiast across simulation replicates, shown for each sampling scenario, time point, and penalisation setting. The points represent the mean across replicates (defined in the main text as sBiast), vertical bars indicate the 25th–75th percentile range. The shaded background marks the forecast time point t = 4. Unpenalised fits are shown in blue and penalised fits in red [PITH_FULL_IMAG… view at source ↗
Figure 5
Figure 5. Figure 5: Standardised mean square error sRMSEt across simulation replicates for each sampling scenario, time point, and penalisation setting. Shaded panels highlight the forecast time point (t = 4). Unpenalised fits are shown in blue and penalised fits in red. infections are endemic in many low-resource settings and are associated with anemia, growth retardation, and impaired cognitive development, particularly amo… view at source ↗
Figure 6
Figure 6. Figure 6: Spatio-temporal distribution of soil-transmitted helminth (STH) survey locations in Kenya, grouped [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Left panel: Timeline of surveys and mass drug administration (MDA) activities by subcounty, with [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Boxplots of the empirical prevalence of three soil-transmitted helminths species by cumulative rounds [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Estimated MDA impact functions for each of the three STH species, based on the fitted exponential [PITH_FULL_IMAGE:figures/full_fig_p020_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Average non-randomised PIT (AnPIT) curves for the held-out samples for the soil transmitted [PITH_FULL_IMAGE:figures/full_fig_p033_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Spatio-temporal distribution of lymphatic filariasis (LF) survey locations in Madagascar, grouped [PITH_FULL_IMAGE:figures/full_fig_p034_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Left panel: Timeline of surveys and mass drug administration (MDA) activities by implementation [PITH_FULL_IMAGE:figures/full_fig_p035_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Boxplots of the empirical prevalence of lymphatic filariasis by cumulative rounds of mass-drug [PITH_FULL_IMAGE:figures/full_fig_p035_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Average non-randomised PIT (AnPIT) curves for the held-out samples for the soil transmitted [PITH_FULL_IMAGE:figures/full_fig_p036_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Predicted number of additional annual mass drug administration (MDA) rounds required to reduce [PITH_FULL_IMAGE:figures/full_fig_p037_15.png] view at source ↗
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.

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 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)
  1. [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.
  2. [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)
  1. [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.
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

1 free parameters · 0 axioms · 0 invented entities

Based solely on the abstract, the model introduces a decay adjustment for time-varying MDA effects, which implies fitted parameters for decay rates and spatio-temporal components; no explicit axioms or invented entities are detailed. Identifiability challenges are noted, suggesting reliance on data-driven choices for model complexity.

free parameters (1)
  • decay adjustment parameters
    Parameters controlling the time-varying impact of MDA that are estimated from survey data to adjust prevalence trends.

pith-pipeline@v0.9.0 · 5434 in / 1141 out tokens · 40793 ms · 2026-05-17T02:07:29.082899+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

31 extracted references · 31 canonical work pages

  1. [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. [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. [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. [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. [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. [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. [7]

    doi: 10.1371/journal.pntd.0006780. C. J. Geyer. Likelihood and exponential families.Department of Statistics, University of Minnesota,

  8. [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. [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. [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. [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. [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. [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. [14]

    URLhttps://doi.org/10

    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. [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. [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. [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. [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. [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. [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. [21]

    32 Giorgi, E., Fronterre, C

    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. [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. [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. [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. [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. [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. [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. [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. [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. [30]

    Monitoring and Epidemiological Assessment of Mass Drug Administration in the Global Programme to Eliminate Lymphatic Filariasis: A Manual for National Elimination Programmes

    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. [31]

    10: Average non-randomised PIT (AnPIT) curves for the held-out samples for the soil transmitted helminths geostatistical analysis

    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...