Pith. sign in

REVIEW 3 major objections 6 minor 40 references

Doubly robust Methods for Recurrent Event Outcomes: Causal Effects of Blood Pressure Medications on Acute Kidney Injuries

T0 review · 3 major / 6 minor · reviewed 2026-07-07 · glm-5.2

Pith's one-line read Doubly robust method ties intensive BP therapy to higher AKI risk

desk verdict Solid integration of longitudinal TMLE with semi-competing risks for recurrent events; the PP estimates rest on a thin time-varying covariate set and need a sensitivity analysis for unmeasured confounding. read the letter →

arxiv 2607.05293 v1 pith:TE3OV7JB submitted 2026-07-06 stat.ME

classification stat.ME
keywords bloodcausaleffectsacuteaveragedataeventkidney
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper extends longitudinal targeted maximum likelihood estimation (TMLE) to recurrent event outcomes under time-varying treatment and semi-competing risks. The central methodological contribution is an efficient influence function for the cumulative count of recurrent events that yields a doubly robust estimator: it remains consistent if either the treatment-adherence model or the outcome model is correctly specified, but not necessarily both. Applied to the SPRINT trial, the method estimates that intensive blood-pressure-lowering therapy causes approximately 0.019 additional acute kidney injury (AKI) events per participant over four years compared to standard therapy (95% CI: [0.011, 0.026]) under intent-to-treat, and 0.017 under per-protocol analysis (95% CI: [0.007, 0.027]). The authors handle death as a semi-competing risk that truncates the recurrent outcome process rather than as a censoring event to be eliminated, arguing this captures the real-world total effect of treatment. The near-equality of ITT and PP estimates is explained by negligible adherence effects: the population intervention effect of full versus observed adherence is essentially zero in both arms.

What carries the argument

Longitudinal TMLE built on the efficient influence function in Eq. (3), which sums inverse-propensity-weighted residuals over time intervals with backward-recursive outcome regression updates. Death enters as a semi-competing risk via the time-varying covariate vector L*_k, with a deterministic link (D_k=1 implies Y_k=0 thereafter). The causal estimand is identified through an iteratively conditional expectation (ICE) form derived from the longitudinal g-formula, and the TMLE updating procedure (Algorithm 1) performs propensity-weighted fluctuation updates to reduce residual bias.

What would settle it

If either the propensity score model or the outcome model is misspecified in a way that biases the weighted residuals, and the other model is also wrong, the TMLE estimator is inconsistent. More specifically for the SPRINT application: if unmeasured confounders of the adherence-AKI relationship exist (e.g., medication side effects influencing both adherence and kidney outcomes), the conditional exchangeability assumption fails and the PP effect estimate of 0.017 is invalid.

Watch

Extended reading notes

Core claim

The efficient influence function for the cumulative recurrent event count, given in Equation (3), decomposes the estimation problem into a sequence of weighted residuals across time intervals. Each term combines an inverse-propensity-weighted outcome residual with backward-recursive updates of conditional outcome expectations. This structure guarantees double robustness: consistency holds if either the propensity score (treatment and adherence mechanism) or the outcome regression is correctly specified. When applied to SPRINT, the estimator produces tight, machine-learning-compatible estimates showing that intensive blood pressure therapy increases AKI recurrences by about 1.7 to 1.9 events,

Load-bearing premise

The per-protocol estimates rely on conditional exchangeability of adherence: that measured covariate history (essentially just mean arterial pressure) captures all common causes of medication adherence and future AKI outcomes. If unmeasured factors like side effects, patient motivation, or access to care drive both adherence and kidney injury risk, the PP estimates are biased. This assumption cannot be tested from the data.

Editorial extensions

If this is right

  • Clinicians weighing intensive blood pressure targets must account for a small but real increase in recurrent AKI risk, estimated at roughly 0.02 additional episodes per patient over four years.
  • The near-equivalence of ITT and PP estimates suggests that adherence interventions alone would not meaningfully change the AKI risk profile of intensive therapy in trial-like populations.
  • Treating death as a semi-competing risk rather than censoring provides a more honest total-effect estimate in populations where mortality differs between treatment arms, since censoring-based approaches would miss the protective effect of death against later AKI events.
  • The doubly robust framework can be extended to other recurrent event settings with time-varying treatments and competing terminal events, such as repeated hospitalizations in heart failure or recurrent strokes in atrial fibrillation patients.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. This manuscript develops a longitudinal TMLE estimator for the average causal effect of treatment strategies on recurrent event outcomes, handling time-varying treatment, semi-competing risks, and model misspecification via double robustness. The method is applied to the SPRINT trial to estimate the effect of intensive versus standard blood pressure therapy on recurrent AKI episodes. The identification results (Section 2, Eqs. 1–2) and the efficient influence function (Eq. 3) follow standard semiparametric theory. The TMLE algorithm is described for K=2 (Algorithm 1) with generalization deferred to Appendix A.3. The SPRINT analysis reports consistent estimates across GLM, XGB, and GAM nuisance models, with total ITT effect 0.019 (95% CI [0.011, 0.026]) and total PP effect 0.017 (95% CI [0.007, 0.027]).

Significance. The paper addresses a practically important problem: causal inference for recurrent events under time-varying treatment with semi-competing risks. The application to SPRINT is well-motivated and the consideration of both ITT and PP estimands in the presence of death as a semi-competing risk is a genuine contribution. The authors provide a transparent TMLE implementation and report consistent estimates across multiple nuisance model specifications (GLM, XGB, GAM), which strengthens the applied contribution. The simulation studies in Appendix A.4 and the time-discretization sensitivity analysis (Section 4.3.3) provide additional support. The work bridges established longitudinal TMLE theory with a concrete clinical application, which is valuable for practitioners.

major comments (3)
  1. Section 4.1: The per-protocol (PP) estimates depend on the sequential exchangeability assumption for adherence (Section 2, condition 1). In the SPRINT application, the time-varying confounder set L*_k consists solely of mean arterial pressure (MAP). Plausible unmeasured confounders of the adherence–AKI relationship—such as kidney function (eGFR), medication side effects, or changes in comorbidity status—are not included. The authors acknowledge this limitation in the Discussion (Section 5, final paragraph), but no formal sensitivity analysis for unmeasured confounding of the adherence–outcome relationship is provided. The time-discretization sensitivity analysis (Section 4.3.3) addresses a distinct concern (coarsening) and does not speak to exchangeability. Since the PP effect is the methodologically novel estimand, a sensitivity analysis (e.g., bounds or tipping-point analysis for unmea
  2. Algorithm 1 (Section 3.3): The TMLE algorithm is presented only for K=2, with the general case deferred to Appendix A.3. For a methods paper whose central contribution is the longitudinal TMLE estimator, the main-text algorithm should either handle general K or clearly state why the K=2 case is sufficient to convey the general procedure. As written, a reader cannot verify from the main text alone how the algorithm scales, how the iterative backward updating proceeds for arbitrary K, or where the K=2 presentation breaks down. This is a load-bearing presentation gap for the central methodological claim.
  3. Section 4.2 and Table 4: The PP effect estimate (0.017, 95% CI [0.007, 0.027]) is nearly identical to the ITT estimate (0.019, 95% CI [0.011, 0.026]). The authors attribute this to negligible adherence effects (Section 4.3.1, PIE results). However, this near-equivalence also raises the question of whether the PP estimator is meaningfully identifying a different estimand or whether the adherence intervention has near-zero practical variation in this population. The manuscript should more explicitly discuss whether the PP and ITT estimands are nearly identified by the same functional of the data distribution in this setting (given high adherence rates), and what the implications are for interpreting the PP result as evidence of the method's value versus an artifact of the data.
minor comments (6)
  1. Table 3: The MAP values (85.0 intensive, 95.0 standard) appear inconsistent with the expected direction of intensive therapy lowering blood pressure. Please verify these values and clarify whether they reflect post-treatment measurements or another definition.
  2. Section 3.2, Eq. (3): The EIF expression uses notation Q^{z,ā}_{k,h} that is defined in Section 3.1 but the connection between the two sections could be made more explicit for readability.
  3. Table 4 caption: 'repectively' should be 'respectively'; 'CDE' is listed in the caption abbreviation list but the controlled direct effect results are discussed in Section 4.3.2 and Appendix A.5.4, not shown in Table 4 itself.
  4. Section 4.1: The text states '9,361 subjects' in the SPRINT cohort and then '9,322 subjects with complete baseline variables.' The 39 excluded subjects should be briefly characterized (e.g., which baseline variables had missingness).
  5. Section 2: The notation R_{K+1} = sum of Y_k uses index k in the sum but the sum runs over j=1 to K+1. Consistent indexing would improve clarity.
  6. Appendix references (A.1–A.4) are cited throughout the main text but not included in the provided manuscript. The referee assumes these appendices contain the identification proofs, general TMLE algorithm, simulation results, and additional SPRINT analyses referenced. These should be included for review if not already submitted.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for a careful and constructive report. The referee raises three major comments concerning (1) the absence of a formal sensitivity analysis for unmeasured confounding of the adherence–outcome relationship in the PP analysis, (2) the presentation of Algorithm 1 only for K=2 with generalization deferred to the appendix, and (3) the near-equivalence of ITT and PP estimates and its implications for interpreting the PP result. We address each comment below. We agree with comments 1 and 2 and will revise the manuscript accordingly; on comment 3 we agree the discussion should be expanded and will partially revise, while noting that the near-equivalence is a substantive finding about the data rather than a methodological limitation.

read point-by-point responses
  1. Referee: Section 4.1: The per-protocol (PP) estimates depend on the sequential exchangeability assumption for adherence. In the SPRINT application, the time-varying confounder set L*_k consists solely of MAP. Plausible unmeasured confounders of the adherence–AKI relationship—such as eGFR, medication side effects, or changes in comorbidity status—are not included. No formal sensitivity analysis for unmeasured confounding of the adherence–outcome relationship is provided. The time-discretization sensitivity analysis addresses a distinct concern and does not speak to exchangeability. Since the PP effect is the methodologically novel estimand, a sensitivity analysis (e.g., bounds or tipping-point analysis) for unmeasured confounding is needed.

    Authors: The referee is correct that the PP estimand relies on sequential exchangeability for adherence, that the confounder set is limited to MAP, and that no formal sensitivity analysis for unmeasured confounding of the adherence–AKI relationship is currently provided. We agree this is an important gap, particularly given that the PP estimand is the methodologically novel contribution. We will add a formal sensitivity analysis for unmeasured confounding in the revised manuscript. Specifically, we plan to implement a tipping-point analysis examining how large an unmeasured confounder would need to be (in terms of its association with both adherence and AKI recurrence) to shift the PP estimate meaningfully—for example, to include the null or to equal the ITT estimate. We will also expand the confounder set where the SPRINT data permit: eGFR was measured in SPRINT and can be incorporated as an additional time-varying covariate, which directly addresses one of the referee's specific concerns. We acknowledge that medication side effects and certain comorbidity changes may not be available in the data and will state this explicitly. The revised manuscript will include this sensitivity analysis as a new subsection in Section 4.3 and will revise the Discussion to reference its results. revision: yes

  2. Referee: Algorithm 1 (Section 3.3): The TMLE algorithm is presented only for K=2, with the general case deferred to Appendix A.3. For a methods paper whose central contribution is the longitudinal TMLE estimator, the main-text algorithm should either handle general K or clearly state why the K=2 case is sufficient to convey the general procedure. A reader cannot verify from the main text alone how the algorithm scales, how the iterative backward updating proceeds for arbitrary K, or where the K=2 presentation breaks down. This is a load-bearing presentation gap.

    Authors: The referee is correct that presenting only the K=2 case in the main text is a presentation gap for a methods paper. We chose K=2 for conciseness, but we agree that a reader should be able to understand the general procedure from the main text. In the revision, we will expand Algorithm 1 (or provide a companion general-K algorithm in the main text) that makes the iterative backward updating structure explicit for arbitrary K. Specifically, we will show how the targeting step at each time point k (updating Q_{k,h} for h = k-1, ..., 0 using inverse propensity score weighting) proceeds recursively backward, and how the K=2 case generalizes. We will also add a brief remark in the main text clarifying that the K=2 presentation illustrates the essential features—nested outcome regression, sequential targeting, and propensity-score-weighted fluctuation—that carry over directly to general K, with the full general-K algorithm retained in the appendix for reference. This addresses the referee's concern that a reader cannot currently verify the scaling from the main text alone. revision: yes

  3. Referee: Section 4.2 and Table 4: The PP effect estimate (0.017, 95% CI [0.007, 0.027]) is nearly identical to the ITT estimate (0.019, 95% CI [0.011, 0.026]). The authors attribute this to negligible adherence effects. However, this near-equivalence also raises the question of whether the PP estimator is meaningfully identifying a different estimand or whether the adherence intervention has near-zero practical variation. The manuscript should more explicitly discuss whether the PP and ITT estimands are nearly identified by the same functional of the data distribution in this setting (given high adherence rates), and what the implications are for interpreting the PP result as evidence of the method's value versus an artifact of the data.

    Authors: The referee raises a substantive point that deserves more explicit discussion. We agree that the near-equivalence of ITT and PP estimates should be interpreted carefully. In the SPRINT setting, adherence rates were high in both arms, which means the PP estimand (full adherence to assigned therapy) differs only slightly from the ITT estimand (assigned therapy under observed adherence) in terms of the target population's outcome distribution. This is a substantive feature of the data, not an artifact of the estimator. The PIE results in Section 4.3.1 directly support this interpretation: the population intervention effect of adherence is estimated at approximately zero in both arms, indicating that shifting from observed to full adherence would not materially change outcomes. We agree, however, that the manuscript should more explicitly discuss the formal relationship between the ITT and PP functionals when adherence is near-uniform. Specifically, when Pr(A_k = 1 | history) is close to 1 for all strata, the inverse probability weights in the PP estimator are close to 1, and the PP and ITT estimating functionals become nearly identical. We will add this discussion to Section 4.3.1 or 4.4. We want to be clear, however, that the value of the PP analysis lies not in producing a different point estimate but in (a) formally defining and identifying the per-protocol estimand under explicit causal assumptions, (b) providing a valid confidence interval for that estimand, and (c) demonstrating that the ITT–PP gap is negligible in this population—which is itself a meaningful clinical finding. The methodological contribution is the framework and estimator, not the magnitude of the ITT–PP difference in this particular application. We will revise the manuscript to make this explicit. revision: partial

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found. The derivation chain builds entirely on established external theory (Robins g-formula, van der Laan TMLE, standard EIF theory), and the SPRINT application uses external trial data with no fitted constants relabeled as predictions.

full rationale

The paper's derivation chain is self-contained against external benchmarks and free of circularity. (1) Causal identification (Eqs. 1–2) follows the standard longitudinal g-formula from Robins [17] and counterfactual framework from Hernán & Robins [13], with explicitly stated exchangeability, positivity, and consistency assumptions. No step reduces to its inputs by construction. (2) The EIF (Eq. 3) is derived via pathwise derivatives using standard semiparametric theory cited to Tsiatis [23], Levy [24], Hines et al. [25], and Hampel [28] — none of which are self-citations (the authors Zhang, Cotton, Wen do not appear in the reference list as cited authors). (3) The TMLE algorithm (Algorithm 1) follows the established targeted learning framework of van der Laan and colleagues [18–22], with no novel ansatz introduced via self-citation. (4) The SPRINT application estimates effects (ITT: 0.019, PP: 0.017) from external trial data using the proposed estimator; no parameter is fitted to a subset and then 'predicted' on a closely related quantity. (5) Simulations in Appendix A.4 validate finite-sample performance independently. The skeptic's valid concern about sequential exchangeability for adherence (with only MAP as time-varying confounder) is a correctness risk about untestable assumptions, not a circularity issue — the estimand is identified from observed data under stated assumptions, not defined in terms of the estimation output.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new entities, particles, or postulated objects. It uses standard causal inference machinery (counterfactual framework, g-formula, TMLE, EIF) applied to a recurrent event outcome. The free parameters are modeling choices (time discretization, nuisance model type) rather than fitted constants. The axioms are standard causal identification assumptions (exchangeability, positivity, consistency) plus the rare-event-per-interval assumption specific to the recurrent event setting.

free parameters (2)
  • Number of time intervals K = 2 (bi-yearly coarsening)
    Chosen for the SPRINT analysis; sensitivity analysis with annual coarsening also conducted. Not a fitted parameter per se but a modeling choice that affects results.
  • Nuisance model specifications (GLM/XGB/GAM) = Multiple specifications compared
    The choice of nuisance model affects estimates; the paper compares parametric GLMs, XGB, and GAMs and reports similar results across specifications.
assumptions (5)
  • domain assumption Conditional exchangeability for adherence: no unmeasured confounders of adherence and outcome given observed history
    Invoked in Section 2 for identification of per-protocol effects. Not testable; plausibility depends on richness of measured covariates, which in SPRINT is primarily MAP.
  • domain assumption Positivity: every adherence level has positive probability within observed covariate strata
    Stated in Section 2; assessed empirically by examining adherence distribution across strata.
  • domain assumption Consistency: adherence indicators are well-defined under interventions
    Stated in Section 2; relies on clarity of the causal question. The authors note this is not empirically testable.
  • domain assumption At most one recurrent event per time interval (rare event assumption)
    Stated in Section 2 and verified in Section 4.1 (less than 0.24% of subjects had multiple AKIs per interval). Justifies binary Y_k.
  • domain assumption Missing completely at random for baseline covariates
    Invoked in Section 4.1; approximately 0.4% missingness in baseline covariates.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Doubly robust Methods for Recurrent Event Outcomes: Causal Effects of Blood Pressure Medications on Acute Kidney Injuries." pith.science (2026). https://pith.science/paper/TE3OV7JB

@misc{pith2026260705293,
  author       = {Pith},
  title        = {Pith review of: Doubly robust Methods for Recurrent Event Outcomes: Causal Effects of Blood Pressure Medications on Acute Kidney Injuries},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TE3OV7JB}},
  note         = {Machine review of arXiv:2607.05293}
}
read the original abstract

Evaluating the average causal effects of treatment strategies on recurrent event outcomes, such as heart attacks or renal failure, is important in clinical and medical research. However, the analysis becomes increasingly complex as multiple interacting factors are considered within a longitudinal setting. In this paper, we use advanced methodologies to estimate the average causal effects of standard versus intensive blood pressure-lowering therapies on acute kidney injury recurrences. We address time-varying treatment and confounding, and model misspecification during the identification and estimation processes for the effect estimands. We analyze the Systolic Blood Pressure Intervention Trial data set using our proposed method, accounting for medication adherence and the semi-competing risk of death observed in the data.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

40 extracted references · 40 canonical work pages

  1. [1]

    Seventh report of the joint national committee on pre- vention, detection, evaluation, and treatment of high blood pressure.hypertension, 42(6):1206–1252, 2003

    Aram V Chobanian, George L Bakris, Henry R Black, William C Cushman, Lee A Green, Joseph L Izzo Jr, Daniel W Jones, Barry J Materson, Suzanne Oparil, Jack- son T Wright Jr, et al. Seventh report of the joint national committee on pre- vention, detection, evaluation, and treatment of high blood pressure.hypertension, 42(6):1206–1252, 2003

  2. [2]

    A randomized trial of intensive versus standard blood-pressure control.The New England journal of medicine, 374(23):2294–2294, 2016

    Jackson T Wright Jr, Paul K Whelton, and David M Reboussin. A randomized trial of intensive versus standard blood-pressure control.The New England journal of medicine, 374(23):2294–2294, 2016

  3. [3]

    Effects of intensive bp control in ckd.Journal of the American Society of Nephrology, 28(9):2812–2823, 2017

    Alfred K Cheung, Mahboob Rahman, David M Reboussin, Timothy E Craven, Tom Greene, Paul L Kimmel, William C Cushman, Amret T Hawfield, Karen C Johnson, Cora E Lewis, et al. Effects of intensive bp control in ckd.Journal of the American Society of Nephrology, 28(9):2812–2823, 2017

  4. [4]

    Adam P Bress, Holly Kramer, Rasha Khatib, Srinivasan Beddhu, Alfred K Cheung, Rachel Hess, Vinod K Bansal, Guichan Cao, Jerry Yee, Andrew E Moran, et al. Po- tential deaths averted and serious adverse events incurred from adoption of the sprint (systolic blood pressure intervention trial) intensive blood pressure regimen in the united states: projections ...

  5. [5]

    Hypertension in chronic kidney disease.Hypertension: from basic research to clinical practice, pages 307–325, 2017

    Seyed Mehrdad Hamrahian and Bonita Falkner. Hypertension in chronic kidney disease.Hypertension: from basic research to clinical practice, pages 307–325, 2017

  6. [6]

    Blood pressure targets in adults with hypertension.Cochrane Database of Systematic Reviews, (12), 2020

    Jose Agustin Arguedas, Viriam Leiva, and James M Wright. Blood pressure targets in adults with hypertension.Cochrane Database of Systematic Reviews, (12), 2020

  7. [7]

    Modelling recurrent events: A tutorial for analysis in epidemiology.International Journal of Epidemiology, 44:324–333, 2015

    Leila D A F Amorim and Jianwen Cai. Modelling recurrent events: A tutorial for analysis in epidemiology.International Journal of Epidemiology, 44:324–333, 2015

  8. [8]

    Causal inference for recurrent events data with all-or- none compliance.Communications in Statistics-Theory and Methods, 45(24):7306– 7325, 2016

    Xiang Gao and Ming Zheng. Causal inference for recurrent events data with all-or- none compliance.Communications in Statistics-Theory and Methods, 45(24):7306– 7325, 2016

Show all 40 references
  1. [9]

    Doubly robust estimation and causal inference for recurrent event data.Statistics in Medicine, 39:2324–2338, 2020

    Chien-Lin Su, Russell Steele, and Ian Shrier. Doubly robust estimation and causal inference for recurrent event data.Statistics in Medicine, 39:2324–2338, 2020

  2. [10]

    Causal inference for re- current event data using pseudo-observations.Biostatistics, 23(1):189–206, 2022

    Chien-Lin Su, Robert W Platt, and Jean-Fran¸ cois Plante. Causal inference for re- current event data using pseudo-observations.Biostatistics, 23(1):189–206, 2022. 12

  3. [11]

    Per-protocol analyses of pragmatic trials

    Miguel A Hern´ an, James M Robins, et al. Per-protocol analyses of pragmatic trials. N Engl J Med, 377(14):1391–1398, 2017

  4. [12]

    Discordance between reported intention-to-treat and per protocol analyses.Journal of clinical epidemiology, 60(7):663–669, 2007

    N´ uria Porta, Catalina Bonet, and Erik Cobo. Discordance between reported intention-to-treat and per protocol analyses.Journal of clinical epidemiology, 60(7):663–669, 2007

  5. [13]

    Robins.Causal Inference: What If

    Miguel A Hern´ an and James M. Robins.Causal Inference: What If. Boca Raton: Chapman & Hall/CRC, 2020

  6. [14]

    Estimating causal effects of treatment in randomized and nonran- domized studies.Journal of Educational Psychology, 66:688–701, 1974

    Donald B Rubin. Estimating causal effects of treatment in randomized and nonran- domized studies.Journal of Educational Psychology, 66:688–701, 1974

  7. [15]

    Statistics and causal inference.Journal of the American Statistical Association, 81:945–960, 1986

    Paul W Holland. Statistics and causal inference.Journal of the American Statistical Association, 81:945–960, 1986

  8. [16]

    Causal inference in statistics: An overview.Statistics Surveys, 3:96–146, 2009

    Judea Pearl. Causal inference in statistics: An overview.Statistics Surveys, 3:96–146, 2009

  9. [17]

    James Robins. A new approach to causal inference in mortality studies with a sustained exposure period—application to control of the healthy worker survivor effect.Mathematical modelling, 7(9-12):1393–1512, 1986

  10. [18]

    Targeted maximum likelihood learning

    Mark J Van Der Laan and Daniel Rubin. Targeted maximum likelihood learning. The international journal of biostatistics, 2(1), 2006

  11. [19]

    Springer, 2011

    Mark J Van der Laan, Sherri Rose, et al.Targeted learning: causal inference for observational and experimental data, volume 4. Springer, 2011

  12. [20]

    ltmle: an r package implementing targeted minimum loss-based estimation for longitudinal data.Journal of Statistical Software, 81:1–21, 2017

    Samuel D Lendle, Joshua Schwab, Maya L Petersen, and Mark J van der Laan. ltmle: an r package implementing targeted minimum loss-based estimation for longitudinal data.Journal of Statistical Software, 81:1–21, 2017

  13. [21]

    Using longitudinal targeted maximum likelihood estimation in complex settings with dynamic interventions.Statistics in medicine, 38(24):4888–4911, 2019

    Michael Schomaker, Miguel Angel Luque-Fernandez, Valeriane Leroy, and Mary- Ann Davies. Using longitudinal targeted maximum likelihood estimation in complex settings with dynamic interventions.Statistics in medicine, 38(24):4888–4911, 2019

  14. [22]

    Studying continuous, time-varying, and/or complex exposures using longitudinal modified treatment policies.Epidemiology, 35(5):667–675, 2024

    Katherine L Hoffman, Diego Salazar-Barreto, Nicholas T Williams, Kara E Rudolph, and Iv´ an D´ ıaz. Studying continuous, time-varying, and/or complex exposures using longitudinal modified treatment policies.Epidemiology, 35(5):667–675, 2024

  15. [23]

    Semiparametric theory and missing data

    Anastasios A Tsiatis. Semiparametric theory and missing data. 2006

  16. [24]

    Tutorial: Deriving the efficient influence curve for large models, 2019

    Jonathan Levy. Tutorial: Deriving the efficient influence curve for large models, 2019

  17. [25]

    Demystifying statistical learning based on efficient influence functions.The American Statistician, 76(3):292–304, 2022

    Oliver Hines, Oliver Dukes, Karla Diaz-Ordaz, and Stijn Vansteelandt. Demystifying statistical learning based on efficient influence functions.The American Statistician, 76(3):292–304, 2022

  18. [26]

    Does water kill? a call for less casual causal inferences.Annals of epidemiology, 26(10):674–680, 2016

    Miguel A Hern´ an. Does water kill? a call for less casual causal inferences.Annals of epidemiology, 26(10):674–680, 2016

  19. [27]

    Story-led causal inference.Epidemiology, 35(3):289–294, 2024

    Jessica G Young. Story-led causal inference.Epidemiology, 35(3):289–294, 2024

  20. [28]

    The influence curve and its role in robust estimation.Journal of the american statistical association, 69(346):383–393, 1974

    Frank R Hampel. The influence curve and its role in robust estimation.Journal of the american statistical association, 69(346):383–393, 1974. 13

  21. [29]

    Brandon K Bellows, Yiyi Zhang, Zugui Zhang, Donald M Lloyd-Jones, Adam P Bress, Jordan B King, Paul Kolm, William C Cushman, Karen C Johnson, Leonardo Tamariz, et al. Estimating systolic blood pressure intervention trial participant posttrial survival using pooled epidemiologi...

  22. [30]

    Formula and nomogram for the sphygmomanometric calcula- tion of the mean arterial pressure.Heart, 84(1):64–64, 2000

    Eduardo Meaney, Felix Alva, Rafael Moguel, Alejandra Meaney, JUAN ALVA, and RICHARD WEBEL. Formula and nomogram for the sphygmomanometric calcula- tion of the mean arterial pressure.Heart, 84(1):64–64, 2000

  23. [31]

    Physiology, mean arterial pressure

    Daniel DeMers and Daliah Wachs. Physiology, mean arterial pressure. InStatPearls [Internet]. StatPearls Publishing, 2021

  24. [32]

    Young, Mats J

    Jessica G. Young, Mats J. Stensrud, Eric J. Tchetgen Tchetgen, and Miguel A. Hern´ an. A causal framework for classical statistical estimands in failure-time settings with competing events.Statistics in Medicine, 39:1199–1236, 2020

  25. [33]

    Causal in- ference with recurrent and competing events.Lifetime data analysis, 30(1):59–118, 2024

    Matias Janvin, Jessica G Young, P˚ al C Ryalen, and Mats J Stensrud. Causal in- ference with recurrent and competing events.Lifetime data analysis, 30(1):59–118, 2024

  26. [34]

    Springer Science & Business Media, 2008

    Odd Aalen, Ornulf Borgan, and Hakon Gjessing.Survival and event history analysis: a process point of view. Springer Science & Business Media, 2008

  27. [35]

    Sebastien Haneuse and Kyu Ha Lee. Semi-competing risks data analysis: accounting for death as a competing risk when the outcome of interest is nonterminal.Circula- tion: Cardiovascular Quality and Outcomes, 9(3):322–331, 2016

  28. [36]

    glm2: Fitting generalized linear models with convergence problems

    Ian C Marschner. glm2: Fitting generalized linear models with convergence problems. The R Journal, 3(2), 2011

  29. [37]

    Xgboost: extreme gradient boosting.R package version 0.4-2, 1(4):1–4, 2015

    Tianqi Chen, Tong He, Michael Benesty, Vadim Khotilovich, Yuan Tang, Hyunsu Cho, Kailong Chen, Rory Mitchell, Ignacio Cano, Tianyi Zhou, et al. Xgboost: extreme gradient boosting.R package version 0.4-2, 1(4):1–4, 2015

  30. [38]

    Wood.Generalized Additive Models: An Introduction with R

    Simon N. Wood.Generalized Additive Models: An Introduction with R. Chapman and Hall/CRC, 2 edition, 2017

  31. [39]

    From patients to policy: population intervention effects in epi- demiology.Epidemiology, 28(4):525–528, 2017

    Daniel Westreich. From patients to policy: population intervention effects in epi- demiology.Epidemiology, 28(4):525–528, 2017

  32. [40]

    Population intervention effects in observational studies to emulate target trial results: reconciling the effects of improved sanitation on child growth

    Elizabeth T Rogawski McQuade, Jade Benjamin-Chung, Daniel Westreich, and Ben- jamin F Arnold. Population intervention effects in observational studies to emulate target trial results: reconciling the effects of improved sanitation on child growth. International journal of epid...

Pith tools

Reviewed July 7, 2026 · model on record in the stance chip above.