Pith. sign in

REVIEW 2 major objections 6 minor 47 references

For trials where patients die before follow-up ends, the paper argues that pairing the survivor average causal effect with restricted mean survival time gives an interpretable treatment-effect picture without ranking death against disease s

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-02 15:15 UTC pith:JHUS6MZN

load-bearing objection Useful taxonomy and a defensible SACE+RMST recommendation, but the SACE estimator rests on an unstated conditional-independence assumption and the simulation priors are partly circular—still deserves refereeing. the 2 major comments →

arxiv 2604.26410 v3 pith:JHUS6MZN submitted 2026-04-29 stat.ME stat.AP

Longitudinal Outcomes Truncated by Death: Causal Estimands and Bayesian Estimators

classification stat.ME stat.AP MSC 62D2062F1562N0162P10
keywords causal estimandstruncation by deathprincipal stratificationsurvivor average causal effectrestricted mean survival timeBayesian estimationcomposite outcomesamyotrophic lateral sclerosis
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

When patients die before the end of follow-up, their longitudinal outcome is no longer a real number, and the treatment effect on that outcome is no longer well defined for everyone. The paper shows that existing estimands split into two classes: those that assume an ordering between death and the outcome scale—such as composite scores that treat death as the worst outcome—and those that avoid such an ordering, notably the survivor average causal effect (SACE) and restricted mean survival time (RMST). Because ordering assumptions encode either a biological claim or a patient preference that is rarely justified in chronic disease, the paper argues that SACE and RMST should be reported together: SACE captures the effect among patients who would survive under either treatment, and RMST captures the effect on survival itself. The authors support this by comparing Bayesian estimators in simulations and in an amyotrophic lateral sclerosis trial, finding that composite estimands can obscure an unfavorable survival effect when the treatment helps function but harms survival. If the paper is right, trials with death truncation should report these two estimands jointly rather than a single composite score.

Core claim

The paper's central claim is that estimands for longitudinal outcomes truncated by death fall into two classes: those that assume an ordering between death and the outcome scale (pairwise comparison, median of the composite outcome), and those that do not (SACE, RMST). Composite estimands encode a preference—commonly that death is worst—which is often inappropriate in chronic disease, where survival and function are driven by the same biological process and patients may value them differently. SACE avoids the ordering by restricting to always-survivors, whose longitudinal outcomes are real numbers under both treatments; RMST captures the survival effect over a fixed horizon. The paper theref

What carries the argument

The central machinery is the principal-stratification decomposition of the population by joint survival under both treatments: always-survivors (LL_t, live under either treatment), protected (DL_t), harmed (LD_t), and never-survivors (DD_t). SACE is the average treatment effect on the longitudinal outcome within LL_t; RMST is the integrated difference of survival curves. The classification of estimands turns on whether they order the post-death symbol * with the real-valued outcome scale. The Bayesian estimators combine a piecewise-constant proportional-hazards/Poisson model for survival (from which the always-survivor probability Ŝ^mis_i,t is computed) with a normal linear model for longitu

Load-bearing premise

The load-bearing premise is that, given baseline covariates, a patient's potential survival time under one treatment carries no information about their potential survival time under the other; if an unmeasured frailty drives both (so the diagonal-Poisson independence is false), the estimated always-survivor weights and the resulting SACE are systematically biased.

What would settle it

Simulate trials with a shared patient-level frailty that multiplies both potential hazards, so T(0) and T(1) are dependent given covariates; run the paper's Bayesian SACE estimator and compare it to the true always-survivor effect. If the estimated SACE shows bias that grows with the frailty variance, the conditional-independence assumption—not just finite-sample noise—is what the estimator depends on.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If SACE and RMST are reported together, a trial can separately quantify a treatment's effect on function among always-survivors and its effect on survival, without collapsing them into one number that presupposes how death compares to disability.
  • Composite estimands such as pairwise comparison and median of composite outcome are only interpretable if the ordering assumption (e.g., death as worst outcome) is accepted; in chronic diseases that assumption can mask a harmful survival effect when the treatment benefits function.
  • The simulation shows that when treatment helps function but harms survival, composite estimands may flip from beneficial to harmful over time, while SACE alone only shows the functional benefit; RMST is needed to reveal the survival cost.
  • The MCO estimand can take infinite values in finite samples, limiting its practicality, whereas SACE and RMST always produce finite estimates.
  • Bayesian estimation with weakly informative priors gives posterior intervals for SACE and RMST, and the paper's application to the ALS trial shows credible intervals that include the null for both, suggesting no clear treatment effect on either dimension.

Where Pith is reading between the lines

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

  • The ordering assumption could be relaxed by eliciting patient-specific utilities for survival-with-disability versus death, making composite estimands preference-sensitive rather than assuming a universal order; the paper's classification suggests a template for such preference-weighted estimands.
  • The SACE estimator's reliance on the estimated always-survivor probability Ŝ^mis_i,t inherits the survival model's assumption that potential survival times under the two treatments are conditionally independent given covariates; a sensitivity analysis introducing a shared frailty would test how robust SACE is to unmeasured common causes of death under both arms.
  • The two-step weighting structure suggests a diagnostic: when the estimated always-survivor probability is near zero for many observed survivors, SACE estimates will be unstable, and reporting the distribution of Ŝ^mis_i,t alongside SACE would make that fragility visible.
  • The estimand classification generalizes beyond ALS to any chronic disease with mortality and patient-reported or functional outcomes, providing a decision rule for choosing estimands before trial analysis.

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

2 major / 6 minor

Summary. The paper reviews causal estimands for longitudinal outcomes truncated by death, proposes a classification into estimands that require an ordering assumption placing death on the real-valued outcome scale (composite estimands such as PC and MCO) and estimands that do not (SACE and RMST). It develops Bayesian estimators for missing potential survival and longitudinal outcomes under a piecewise-constant proportional hazards model and a normal longitudinal model, evaluates them in simulations with four scenarios, and illustrates them on the olesoxime ALS RCT. The paper's central recommendation is that, in the presence of truncation by death, SACE and RMST should be reported together to give an interpretable, multidimensional characterization of treatment effects.

Significance. If the central claim holds, the paper provides a useful practical synthesis for trialists in chronic, fatal diseases: it clarifies that composite estimands encode value judgments about death versus disability that may not be clinically appropriate, while SACE and RMST avoid such ordering assumptions. The paper is commendable for its explicit finite-sample formulation, its reproducible public code, its honest tabulation of bias and coverage, and its application to a real ALS trial. However, the SACE estimator relies on a conditional independence assumption for potential survival times that is not stated, and the simulation study calibrates priors to the true generating values, so the empirical support for the recommendation is partly circular. The conceptual taxonomy is defensible, but the paper needs revision before the practical guidance can be accepted.

major comments (2)
  1. [§4.1.1 and §4.2.2] The SACE estimator weights observed survivors by S_hat_mis_i,t, the model-based probability that the counterfactual survival time exceeds t. In §4.1.1, the joint model for (d_i,j(0), d_i,j(1)) is a diagonal bivariate Poisson, which implies T(0) ⟂ T(1) | X. This conditional independence is never stated or justified in the identification sections (§3.3.1, §4.1.1). Under unmeasured shared frailty or common genetic risk, observing T_obs > t under the assigned arm is informative about T_mis; the correct weight is Pr(T_mis > t | T_obs > t, X), not Pr(T_mis > t | X). Because the paper's central recommendation—report SACE alongside RMST—depends on SACE being reliably estimable, this missing assumption is load-bearing. Please state the assumption explicitly, justify it in the ALS context, or implement a sensitivity analysis with correlated potential survival times.
  2. [§5.1 and Tables 8–9] The simulation sets the longitudinal slope prior β_u,0 ~ N(−2,3) to the simulated control slope and the hazard prior λ_u,j ~ Γ(0.035,0.1) to the true hazard range. These are not weakly informative priors in the usual sense; they are calibrated to the data-generating truth. Favorable coverage results (e.g., SACE 65–99%) may therefore reflect prior tuning rather than estimator performance. Appendix C checks only one alternative slope prior in the beneficial scenario and does not vary the hazard prior. Please add a systematic prior sensitivity analysis—varying both means and variances, including priors centered away from the truth—and discuss how bias and coverage change. This is load-bearing for the empirical claim that the proposed estimators are reliable.
minor comments (6)
  1. [Table 3] The PC row appears to contain a duplicated term: the expression for τ_PC includes two identical '1/2 1(Z_i,t(1)=Z_i,t(0))' terms. Please correct.
  2. [References] Reference [31] is incomplete: it contains placeholder text '<NOTEBOOK AUTHORS,S.A. <notebook title>'. This must be fixed.
  3. [§5.1] In the sentence listing the estimators, 'MOC' should be 'MCO' (median of the composite outcome).
  4. [§3.3.1] Typo: 'Conditionning' should be 'Conditioning'.
  5. [§5.2 and Table 8] The PC estimator has 0% coverage in the beneficial and mixed scenarios, which is striking. The text notes prior sensitivity but does not interpret this coverage failure; please add a sentence explaining whether this is an artifact of the symmetric no-effect prior or a more general limitation of the PC estimator.
  6. [§4.2.2] The formula for the RMST estimator is typeset awkwardly: the integral of S_hat_mis appears without a clear differential placement. Please reformat to avoid ambiguity.

Circularity Check

1 steps flagged

Simulation priors encode true data-generating values, making coverage validation partially circular; central SACE+RMST recommendation remains independent.

specific steps
  1. other [Section 5.1 (Simulation Study); Supplementary A Tables 7-9]
    "We assume prior knowledge on the longitudinal slope, and the priors are set to the simulated value β u,0_t ∼ N(−2,3). The simulated median survival lies between 18 and 22 months under a Weibull distribution (Supplementary A). Assuming a constant hazard over time and with equation 7, this corresponds to λ u,j ∈ [0.032,0.038]. Accordingly, we specify λ u,j ∼ Γ(0.035,0.1)."

    The simulation study's validation targets—coverage and bias (Supplementary A Tables 8–9)—are evaluated using priors centered on the data-generating truth: the longitudinal slope prior N(−2,3) is centered on the true slope a0_0 = −2, and the hazard prior Γ(0.035,0.1) is centered on the true median survival range (18–22 months). Thus the reported coverage values are partly a consequence of the prior matching the truth, not an independent test of the estimators' calibration. This is a circular element in the simulation evidence, although it does not affect the paper's conceptual recommendation to pair SACE with RMST, which rests on the definitions and interpretation of the estimands rather than on these simulation results.

full rationale

The paper's central contribution is a conceptual review and recommendation: estimands for death-truncated longitudinal outcomes fall into two classes, and pairing SACE with RMST gives an interpretable characterization. This conclusion follows from the definitions of the estimands (SACE conditions on always-survivors; RMST summarizes survival) and from the stated implausibility of ordering death on the outcome scale. It is not derived by fitting a parameter and then predicting the same quantity, and there are no load-bearing self-citations: the cited works (Frangakis & Rubin, Imbens & Rubin, Mattei et al., Grossi et al.) are external, and the modeling choices are stated as assumptions, not imported as unverified author-specific results. The one genuine circular element is the simulation study: the priors for the longitudinal slope and hazard are explicitly set to the simulated truth, so the favorable coverage results are partly in-sample. The paper does include a prior-sensitivity check in Appendix C for the slope prior, which mitigates but does not eliminate this concern. The skeptical point about unstated conditional independence of potential survival times (diagonal Poisson in §4.1.1) is a substantive identification/correctness concern about bias under unmeasured frailty, not a circularity, because the model explicitly posits conditional independence and the paper acknowledges strong ignorability assumptions in the Discussion. Overall, the circularity is limited to a validation artifact and does not undermine the central argument's independence.

Axiom & Free-Parameter Ledger

4 free parameters · 8 axioms · 0 invented entities

The central recommendation rests on standard causal assumptions that are plausible in RCTs, plus two stronger untestable assumptions: conditional independence of potential outcomes and, for composite estimands, a value judgment ordering death as worst. The simulation evaluation uses priors centered at the true generating values, which is a free-parameter concern.

free parameters (4)
  • Longitudinal slope prior mean β_{u,0} = N(-2,3) in simulation; N(-1,3) in ALS application
    Hand-chosen: in the simulation it is set to the true simulated slope; in the application it is set to the literature progression rate. It affects SACE, PC, and MCO estimates.
  • Hazard prior baseline λ_{u,j} = Γ(0.035,0.1)
    Chosen to match an expected median survival of 18-30 months; in the simulation this equals the true generating hazard range, making the prior unusually well calibrated.
  • Survival covariate effects α_u and piecewise hazards λ_{u,j} posterior = posterior samples
    Estimated from data under the piecewise-constant proportional hazards model; standard model parameters.
  • Longitudinal intercepts/slopes β_{u,k,t} and residual SD σ_t = posterior samples
    Estimated from data via the Normal linear model per treatment arm; weakly informative priors for covariates and residual SD.
axioms (8)
  • domain assumption SUTVA and no interference for survival and longitudinal outcomes
    Stated in §4.1.1 and §4.1.2; standard and reasonable in an RCT.
  • domain assumption Strong ignorable treatment assignment given covariates for survival
    Stated in §4.1.1; holds by randomization.
  • domain assumption Strong ignorable treatment assignment for longitudinal outcome conditional on principal strata
    Stated in §4.1.2 and §7; acknowledged as strong and requiring substantive knowledge.
  • ad hoc to paper Conditional independence of potential survival times T(0) and T(1) given X
    Implied by the diagonal Poisson model in §4.1.1 but never stated; load-bearing for the SACE weights.
  • ad hoc to paper Uncorrelated potential longitudinal outcomes given X and LL stratum
    Stated in §4.1.2 as an assumption to avoid unsupported dependence; unidentifiable from data but does not affect the SACE mean parameter.
  • domain assumption Death as worst possible outcome ordering for PC and MCO estimands
    Assumed in §3.2; a value judgment that may not reflect individual patient preferences in chronic disease.
  • domain assumption Missing at random for longitudinal outcomes conditional on covariates
    Stated in §6.1 for the ALS application.
  • domain assumption Normal approximation for ALSFRS-R without floor or ceiling effects
    Stated in §6.1; reasonable given baseline scores but an approximation.

pith-pipeline@v1.3.0-alltime-deepseek · 21131 in / 16729 out tokens · 177507 ms · 2026-08-02T15:15:42.753085+00:00 · methodology

0 comments
read the original abstract

In randomized controlled trials with longitudinal outcomes, death before the end of follow-up poses a fundamental challenge: after death, the outcome is no longer a real-valued measurement. This complicates the definition and interpretation of causal estimands, particularly when treatment may affect both survival and longitudinal outcomes. We review existing estimands for longitudinal outcomes truncated by death and clarify the assumptions required for their identification and estimation. We show that these estimands fall into two broad classes, distinguished by whether they require additional assumptions to compare longitudinal outcomes beyond death. Such assumptions may be inappropriate in chronic diseases, either because i) death and longitudinal outcomes are driven by the same underlying biological process or ii) the relative desirability of survival with poor function versus death may depend on individual preferences. We compare the behavior of the estimands in a simulation study using Bayesian estimators and illustrate their use with data from a randomized controlled trial in amyotrophic lateral sclerosis. We argue that, in the presence of death truncation, pairing the survivor average causal effect with the restricted mean survival time estimand provides an interpretable characterization of treatment effects on longitudinal and survival outcomes.

discussion (0)

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

Reference graph

Works this paper leans on

47 extracted references · 21 canonical work pages

  1. [1]

    J., KOCHUROV, M., KUMAR, R., LAO, J., LUHMANN, C

    ABRIL-PLA, O., ANDREANI, V., CARROLL, C., DONG, L., FONNESBECK, C. J., KOCHUROV, M., KUMAR, R., LAO, J., LUHMANN, C. C., MARTIN, O. A., OSTHEGE, M., VIEIRA, R., WIECKI, T. and ZINKOV, R. (2023). PyMC: A Modern and Comprehensive Probabilistic Program- ming Framework in Python.PeerJ Computer Science9. https://doi.org/10.7717/peerj-cs.1516

  2. [2]

    R., CAIN, K

    ANDERSON, J. R., CAIN, K. C. and GELBER, R. D. (1983). Analysis of survival by tumor response.Journal of Clinical On- cology: Official Journal of the American Society of Clinical On- cology1710–719. https://doi.org/10.1200/JCO.1983.1.11.710

  3. [3]

    A., BRUIJN, L

    ANDREWS, J. A., BRUIJN, L. I. and SHEFNER, J. M. (2019). ALS drug development guidances and trial guide- lines.Neurology9366–71. https://doi.org/10.1212/WNL. 0000000000007695

  4. [5]

    A., MALASPINA, A., ROGERS, M.-L., HORNSTEIN, E., LOMBARDI, V., REN- FREY, D., SHEPHEARD, S., MAGEN, I., COHEN, Y., GRANIT, V., STATLAND, J

    BENATAR, M., MACKLIN, E. A., MALASPINA, A., ROGERS, M.-L., HORNSTEIN, E., LOMBARDI, V., REN- FREY, D., SHEPHEARD, S., MAGEN, I., COHEN, Y., GRANIT, V., STATLAND, J. M., HECKMANN, J. M., RADE- MAKERS, R., MCHUTCHISON, C. A., PETRUCELLI, L., MCMILLAN, C. T. and WUU, J. (2024). Prognostic Clinical and Biological Markers for Amyotrophic Lateral Sclerosis Dise...

  5. [6]

    D., MILLER, R., MOORE, D

    BERRY, J. D., MILLER, R., MOORE, D. H., CUDKOW- ICZ, M. E., VANDENBERG, L. H., KERR, D. A., DONG, Y., INGERSOLL, E. W. and ARCHIBALD, D. (2013). The Com- bined Assessment of Function and Survival (CAFS): A new end- point for ALS clinical trials.Amyotrophic Lateral Sclerosis and Frontotemporal Degeneration14162–168. Publisher: Taylor & Francis _eprint: htt...

  6. [7]

    CARLIN, J. B. and MORENO-BETANCUR, M. (2024). On the uses and abuses of regression models: a call for reform of sta- tistical practice and teaching. arXiv:2309.06668 [stat]. https: //doi.org/10.48550/arXiv.2309.06668

  7. [8]

    L., SIMPSON, E., SHEFNER, J

    CASTRILLO-VIGUERA, C., GRASSO, D. L., SIMPSON, E., SHEFNER, J. and ET., A. (2010). Clinical significance in the change of decline in ALSFRS-R.Amyotrophic Lateral Sclerosis. https://doi.org/10.3109/17482960903093710

  8. [9]

    and YANG, S

    COLNET, B., MAYER, I., CHEN, G., DIENG, A., LI, R., VARO- QUAUX, G., VERT, J.-P., JOSSE, J. and YANG, S. (2023). Causal inference methods for combining randomized trials and observational studies: a review. arXiv:2011.08047 [stat]. https: //doi.org/10.48550/arXiv.2011.08047

  9. [10]

    C., SEIFU, Y., MÜTZE, T

    DONG, G., HUANG, B., VERBEECK, J., CUI, Y., SONG, J., GAMALO-SIEBERS, M., WANG, D., HOAGLIN, D. C., SEIFU, Y., MÜTZE, T. and KOLASSA, J. (2023). Win statis- tics (win ratio, win odds, and net benefit) can complement one another to show the strength of the treatment effect on time- to-event outcomes.Pharmaceutical Statistics2220–33. https: //doi.org/10.100...

  10. [11]

    (EMA), E. M. A. (2020). E9 (R1) Step 5 addendum on esti- mands and Sensitivity Analysis in Clinical Trials to the guideline on statistical principles for clinical trials. LONGITUDINAL OUTCOMES TRUNCATED BY DEATH: CAUSAL ESTIMANDS AND BAYESIAN ESTIMATORS13

  11. [12]

    and G, T

    F, T., G, D., S, B. and G, T. (2020). Genetics and Sex in the Pathogenesis of Amyotrophic Lateral Sclerosis (ALS): Is There a Link?International journal of molecular sciences21. https: //doi.org/10.3390/ijms21103647

  12. [13]

    P., BRITTAIN, E

    FAY, M. P., BRITTAIN, E. H., SHIH, J. H., FOLLMANN, D. A. and GABRIEL, E. E. (2018). Causal estimands and confidence intervals associated with Wilcoxon-Mann-Whitney tests in ran- domized experiments.Statistics in Medicine372923–2937. https://doi.org/10.1002/sim.7799

  13. [14]

    FAY, M. P. and LI, F. (2024). Causal Interpretation of the Haz- ard Ratio in Randomized Clinical Trials.Clinical trials (London, England)21623. https://doi.org/10.1177/17407745241243308

  14. [15]

    Adjusting for Covariates in Randomized Clinical Trials for Drugs and Biological Products.FDA guidance docu- ment

    FDA (2023). Adjusting for Covariates in Randomized Clinical Trials for Drugs and Biological Products.FDA guidance docu- ment

  15. [16]

    FRANGAKIS, C. E. and RUBIN, D. B. (2002). Principal stratifi- cation in causal inference.Biometrics5821–29. https://doi.org/ 10.1111/j.0006-341x.2002.00021.x

  16. [17]

    H., CHENG, B., SALACHAS, F., PRADAT, P.- F., BRUNETEAU, G., CORCIA, P., LACOMBLEZ, L

    GORDON, P. H., CHENG, B., SALACHAS, F., PRADAT, P.- F., BRUNETEAU, G., CORCIA, P., LACOMBLEZ, L. and MEININGER, V. (2010). Progression in ALS is not linear but is curvilinear.Journal of Neurology2571713–1717. https://doi. org/10.1007/s00415-010-5609-1

  17. [18]

    and CHIÒ, A

    GRASSANO, M., MOGLIA, C., PALUMBO, F., KOUMAN- TAKIS, E., CUGNASCO, P., CALLEGARO, S., CANOSA, A., MANERA, U., VASTA, R., DEMATTEI, F., MAT- TEONI, E., FUDA, G., SALAMONE, P., MARCHESE, G., CASALE, F., DEMARCHI, F., MAZZINI, L., MORA, G., CALVO, A. and CHIÒ, A. Sex Differences in Amyotrophic Lateral Sclerosis Survival and Progression: A Multidi- mensional...

  18. [19]

    and MEALLI, F

    GROSSI, G., MARIANI, M., MATTEI, A. and MEALLI, F. (2025). Bayesian Principal Stratification with Longitudinal Data and Truncation by Death.Econometrics and Statistics. https: //doi.org/10.1016/j.ecosta.2025.01.003

  19. [20]

    and JUNG, I

    HAN, K. and JUNG, I. (2022). Restricted Mean Survival Time for Survival Analysis: A Quick Guide for Clinical Researchers. Korean Journal of Radiology23495–499. https://doi.org/10. 3348/kjr.2022.0061

  20. [21]

    HAND, D. J. (1992). On Comparing Two Treatments.The Amer- ican Statistician46190–192. https://doi.org/10.1080/00031305. 1992.10475881

  21. [22]

    A., HERNÁNDEZ-DÍAZ, S

    HERNÁN, M. A., HERNÁNDEZ-DÍAZ, S. and ROBINS, J. M. (2004). A Structural Approach to Selection Bias.Epidemiology 15615. https://doi.org/10.1097/01.ede.0000135174.63482.43

  22. [23]

    HOLFORD, T. R. (1980). The Analysis of Rates and of Survivor- ship Using Log-Linear Models.Biometrics36299–305. Pub- lisher: [Wiley, International Biometric Society]. https://doi.org/ 10.2307/2529982

  23. [24]

    HOFFMAN, M. D. and GELMAN, A. (2014). The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo.Journal of Machine Learning Research

  24. [25]

    IMBENS, G. W. and RUBIN, D. B. (2015).Causal Inference for Statistics, Social, and Biomedical Sciences: An Introduc- tion, 1 ed. Cambridge University Press. https://doi.org/10.1017/ CBO9781139025751

  25. [26]

    KARRISON, T. G. (1997). Use of Irwin’s restricted mean as an index for comparing survival in different treatment groups– interpretation and power considerations.Controlled Clinical Trials18151–167. https://doi.org/10.1016/s0197-2456(96) 00089-x

  26. [27]

    and OLIVIER, D

    LAIRD, N. and OLIVIER, D. (1981). Covariance Analysis of Censored Survival Data Using Log-Linear Analysis Techniques. Journal of the American Statistical Association76231–240. Publisher: [American Statistical Association, Taylor & Francis, Ltd.]. https://doi.org/10.2307/2287816

  27. [28]

    L., LU- DOLPH, A., MORA, J

    LENGLET, T., LACOMBLEZ, L., ABITBOL, J. L., LU- DOLPH, A., MORA, J. S., ROBBERECHT, W., SHAW, P. J., PRUSS, R. M., CUVIER, V., MEININGER, V. andTHEMITO- TARGET STUDY GROUP(2014). A Phase II-III Trial of Ole- soxime in Subjects with Amyotrophic Lateral Sclerosis.Euro- pean Journal of Neurology21529–536. https://doi.org/10.1111/ ene.12344

  28. [29]

    and MEALLI, F

    MATTEI, A., DING, P., BALLERINI, V. and MEALLI, F. (2025). Assessing Causal Effects in the Presence of Treatment Switching Through Principal Stratification.Bayesian Analysis20. https: //doi.org/10.1214/24-BA1425

  29. [31]

    <notebook title>

    <NOTEBOOK AUTHORS,S.A. <notebook title>. InPyMC exam- ples(P. Team, ed.). https://doi.org/10.5281/zenodo.5654871

  30. [32]

    OLARTEPARRA, C., DANIEL, R. M. and BARTLETT, J. W. (2023). Hypothetical Estimands in Clinical Trials: A Unification of Causal Inference and Missing Data Methods.Statistics in Bio- pharmaceutical Research15421–432. https://doi.org/10.1080/ 19466315.2022.2081599

  31. [33]

    M., WRIGHT, D

    OLARTEPARRA, C., DANIEL, R. M., WRIGHT, D. and BARTLETT, J. W. (2025). Estimating hypothetical estimands with causal inference and missing data estimators in a diabetes trial case study.Biometrics81ujae167. https://doi.org/10.1093/ biomtc/ujae167

  32. [34]

    J., ARITI, C

    POCOCK, S. J., ARITI, C. A., COLLIER, T. J. and WANG, D. (2012). The win ratio: a new approach to the analysis of com- posite endpoints in clinical trials based on clinical priorities. European Heart Journal33176–182. https://doi.org/10.1093/ eurheartj/ehr352

  33. [35]

    R., VESTRUCCI, M., DE- TRY, M

    QUINTANA, M., SAVILLE, B. R., VESTRUCCI, M., DE- TRY, M. A., CHIBNIK, L., SHEFNER, J., BERRY, J. D., CHASE, M., ANDREWS, J., SHERMAN, A. V., YU, H., DRAKE, K., CUDKOWICZ, M., PAGANONI, S., MACK- LIN, E. A. and HEALEY ALS PLATFORMTRIALSTUDY GROUP(2023). Design and Statistical Innovations in a Platform Trial for Amyotrophic Lateral Sclerosis.Annals of Neuro...

  34. [36]

    S., EVANS, R

    RICHARDSON, T. S., EVANS, R. J. and ROBINS, J. M. (2011). Transparent Parametrizations of Models for Potential Outcomes. InBayesian Statistics 9(J. M. Bernardo, M. J. Ba- yarri, J. O. Berger, A. P. Dawid, D. Heckerman, A. F. M. Smith and M. West, eds.) 569–610. Oxford University Press. https: //doi.org/10.1093/acprof:oso/9780199694587.003.0019

  35. [37]

    (2012).Joint Models for Longitudinal and Time-to-Event Data: With Applications in R

    RIZOPOULOS, D. (2012).Joint Models for Longitudinal and Time-to-Event Data: With Applications in R. Chapman and Hall/CRC, New York. https://doi.org/10.1201/b12208

  36. [38]

    ROBINS, J. M. and GREENLAND, S. (1992). Identifiability and exchangeability for direct and indirect effects.Epidemiol- ogy (Cambridge, Mass.)3143–155. https://doi.org/10.1097/ 00001648-199203000-00013

  37. [39]

    RODRÍGUEZ, G. (2007). Lecture Notes on Generalized Linear Models. https://grodri.github.io/glms/notes/

  38. [40]

    and HARDIMAN, O

    ROONEY, J., BURKE, T., VAJDA, A., HEVERIN, M. and HARDIMAN, O. (2017). What does the ALSFRS-R really mea- sure? A longitudinal and survival analysis of functional dimen- sion subscores in amyotrophic lateral sclerosis.Journal of Neu- 14 rology, Neurosurgery, and Psychiatry88381–385. https://doi. org/10.1136/jnnp-2016-314661

  39. [41]

    ROSENBAUM, P. R. and RUBIN, D. B. (1983). The central role of the propensity score in observational studies for causal ef- fects.Biometrika7041–55. https://doi.org/10.1093/biomet/70. 1.41

  40. [42]

    and PARMAR, M

    ROYSTON, P. and PARMAR, M. K. B. (2011). The use of restricted mean survival time to estimate the treatment effect in randomized clinical trials when the proportional hazards assump- tion is in doubt.Statistics in Medicine302409–2421. _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1002/sim.4274. https://doi.org/10.1002/sim.4274

  41. [43]

    Censoring

    RUBIN, D. B. (2006). Causal Inference Through Potential Out- comes and Principal Stratification: Application to Studies with “Censoring” Due to Death.Statistical Science21. https://doi. org/10.1214/088342306000000114

  42. [44]

    O., MALEK, A

    TALBOTT, E. O., MALEK, A. M. and LACOMIS, D. (2016). The epidemiology of amyotrophic lateral sclerosis.Handbook of Clinical Neurology138225–238. https://doi.org/10.1016/ B978-0-12-802973-2.00013-6

  43. [45]

    D., REYNDERS, D., MUSORO, J

    THOMASSEN, D., ROYCHOUDHURY, S., AMDAL, C. D., REYNDERS, D., MUSORO, J. Z., SAUERBREI, W., GOET- GHEBEUR, E.,LECESSIE, S., KAMALAKAR, R., SAIL, K., BASCH, E., REIJNEVELD, J., KEATING, K., SU, Y., ALANYA, A., RUMPOLD, G., SASSEVILLE, M., BLACK, J., MOLENBERGHS, G., RANTELL, K., SCHLICHTING, M., REG- NAULT, A., NESS, D.,TENSELDAM, S., REKSTEN, T. R., SCHIEL...

  44. [46]

    VANEIJK, R. P. A., VANDENBERG, L. H. and LU, Y. (2022). Composite endpoint for ALS clinical trials based on pa- tient preference: Patient-Ranked Order of Function (PROOF). Journal of Neurology, Neurosurgery & Psychiatry93539–546. https://doi.org/10.1136/jnnp-2021-328194

  45. [47]

    XIANG, Q., BOSCH, R. J. and LOK, J. J. (2023). The Survival- incorporated Median vs the Median in the Survivors or in the Always-survivors: What Are We Measuring? And Why?Statis- tics in Medicine425479–5490. https://doi.org/10.1002/sim. 9922

  46. [48]

    A., SOLOMON, S

    ZHAO, L., CLAGGETT, B., TIAN, L., UNO, H., PF- EFFER, M. A., SOLOMON, S. D., TRIPPA, L. and WEI, L. J. (2016). On the restricted mean survival time curve in survival analysis.Biometrics72215–221. _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1111/biom.12384. https://doi.org/10.1111/biom.12384

  47. [49]

    D., PFEF- FER, M

    ZHAO, L., TIAN, L., UNO, H., SOLOMON, S. D., PFEF- FER, M. A., SCHINDLER, J. S. and WEI, L. J. (2012). Utilizing the integrated difference of two survival functions to quantify the treatment contrast for designing, monitoring, and analyzing a comparative clinical study.Clinical Trials (London, England)9 570–577. https://doi.org/10.1177/1740774512455464 LO...