pith. sign in

arxiv: 2604.12740 · v1 · submitted 2026-04-14 · 📊 stat.AP

Bayesian Joint Modelling of Longitudinal Creatinine Trajectories in Children with Auto-Immune Disorders to Predict Paediatric Kidney Disease Risk in a Single Centre Study

Pith reviewed 2026-05-10 14:20 UTC · model grok-4.3

classification 📊 stat.AP
keywords joint modelingcreatinine trajectoriespediatric kidney diseaseautoimmune disordersBayesian analysisdynamic risk predictionlongitudinal datatime-to-event analysis
0
0 comments X

The pith

Joint Bayesian models link rising creatinine trajectories to higher risk of kidney events in children with autoimmune disorders.

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

This paper applies a joint modeling framework to repeated serum creatinine measurements and time-to-event data in children with autoimmune disorders at one hospital. It establishes that evolving creatinine profiles are strongly associated with the composite outcome of death, acute kidney injury, or chronic kidney disease. Certain treatments raise that risk while others lower it, and creatinine levels themselves shift with age and body-mass index. The model produces updated risk estimates for each patient as new creatinine readings arrive, aiming to support real-time clinical decisions in pediatric nephrology.

Core claim

The paper claims a strong association between evolving creatinine profiles and the risk of the composite event. Treatment with corticosteroids and calcium channel blockers was associated with an increased event risk, whereas immunosuppressive therapy was associated with a reduced risk. Creatinine trajectories were significantly influenced by age and BMI z-score. Dynamic risk predictions generated from patients' observed creatinine trajectories illustrate the framework's utility for personalized risk assessment.

What carries the argument

The joint modelling framework that simultaneously models repeated creatinine measurements and the time to the composite event, letting the biomarker trajectory directly inform the hazard of the event.

If this is right

  • Corticosteroid and calcium channel blocker use increases the risk of the composite kidney event.
  • Immunosuppressive therapy decreases the risk of death, acute kidney injury, or chronic kidney disease.
  • Creatinine trajectories vary significantly with patient age and BMI z-score.
  • Dynamic risk predictions can be updated for individual patients using their ongoing creatinine measurements.

Where Pith is reading between the lines

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

  • If validated, the approach could allow earlier treatment changes based on a child's creatinine path before an event occurs.
  • Extending the same joint-model structure to other repeated biomarkers could improve outcome forecasting in related pediatric conditions.
  • Multi-center data would be required to test whether the observed treatment associations generalize beyond the original hospital population.

Load-bearing premise

The joint model accurately captures the true relationship between creatinine changes and event risk without unmeasured confounding, model misspecification, or selection bias from single-center data.

What would settle it

In an external group of similar children, the model's predicted risks from their creatinine trajectories fail to match the actual rates of death, acute kidney injury, or chronic kidney disease.

Figures

Figures reproduced from arXiv: 2604.12740 by Christiana Charalambous, John Booth, Qendresa Selimi, Stephen D Marks, Taban Baghfalaki.

Figure 1
Figure 1. Figure 1: Kaplan–Meier survival curve estimating the probability of remaining free from acute kidney [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Distribution of serum creatinine and its log-transformed values in the GPKC dataset. [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Log-transformed serum creatinine trajectories for a randomly selected subset of patients from [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Dynamic predictions of survival probability for three patients (38, 313, and 415), evaluated at [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
read the original abstract

This study investigates the relationship between longitudinal serum creatinine measurements and the risk of adverse kidney outcomes in paediatric patients with auto-immune disorders at Great Ormond Street Hospital for Children NHS Foundation Trust, London. To jointly analyse repeated biomarker measurements and time-to-event outcomes, we employed a joint modelling framework that combines the creatinine trajectories with the time to death or diagnosis of acute kidney injury or chronic kidney disease. Covariates considered in analysis included demographic and clinical characteristics. The results demonstrate a strong association between evolving creatinine profiles and the risk of the composite event. Specifically, treatment with corticosteroids and calcium channel blockers was associated with an increased event risk, whereas immunosuppressive therapy was associated with a reduced risk. The longitudinal component showed that creatinine trajectories were significantly influenced by age and BMI z-score. To demonstrate the practical utility of the proposed framework, dynamic risk predictions were generated using patients' observed creatinine trajectories. Model performance was compared using model selection criteria, alongside area under the curve and Brier score to evaluate the accuracy of dynamic risk predictions. These predictions illustrate the potential of joint models to support personalised medicine and clinical decision making in paediatric nephrology through real-time risk assessment.

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

3 major / 2 minor

Summary. The paper presents a Bayesian joint model combining longitudinal serum creatinine trajectories with a time-to-event submodel for a composite outcome (death, AKI or CKD) in a single-center cohort of children with auto-immune disorders. It reports that age and BMI z-score influence creatinine trajectories, that corticosteroids and calcium-channel blockers are associated with increased event risk while immunosuppressive therapy is associated with reduced risk, and that dynamic risk predictions derived from observed trajectories achieve useful AUC and Brier-score performance.

Significance. If the joint-model assumptions hold and treatment associations are robust to confounding, the work offers a practical demonstration of how longitudinal biomarker data can be used for real-time, individualized risk assessment in pediatric nephrology. The Bayesian framework and emphasis on dynamic predictions are appropriate strengths for this setting.

major comments (3)
  1. [Methods] Methods section: the joint-model specification (longitudinal submodel form, random-effects structure, association parameter linking creatinine process to hazard, priors, and software implementation) is not described in sufficient detail to allow reproduction or evaluation of the reported associations and predictions.
  2. [Results] Results / Discussion: the directional treatment associations (corticosteroids and calcium-channel blockers increasing composite-event risk; immunosuppressives decreasing it) are presented without addressing confounding by indication, time-varying treatment effects, or sensitivity analyses, which is load-bearing for the highlighted clinical claims in an observational single-center study.
  3. [Abstract/Methods] Abstract and Methods: no information is given on missing-data mechanisms for the longitudinal creatinine measurements, model diagnostics, or any form of internal or external validation of the dynamic predictions beyond AUC/Brier scores on the same data.
minor comments (2)
  1. [Abstract] Abstract: the phrase 'model selection criteria' is used without naming the specific criteria (DIC, WAIC, etc.) employed.
  2. [Discussion] The single-center design and lack of external validation should be more explicitly acknowledged when discussing implications for 'personalised medicine'.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed review of our manuscript. We have addressed each major comment below and agree that the suggested additions will improve clarity, reproducibility, and appropriate interpretation of the findings. We plan to submit a revised version incorporating these changes.

read point-by-point responses
  1. Referee: [Methods] Methods section: the joint-model specification (longitudinal submodel form, random-effects structure, association parameter linking creatinine process to hazard, priors, and software implementation) is not described in sufficient detail to allow reproduction or evaluation of the reported associations and predictions.

    Authors: We agree that the Methods section requires greater detail to support reproducibility. In the revised manuscript we will expand this section to specify: the longitudinal submodel as a linear mixed-effects model on log-creatinine with fixed effects for age and BMI z-score plus random intercepts and slopes; the random-effects covariance structure; the association structure linking the current value of the longitudinal process to the hazard; the prior distributions employed in the Bayesian estimation; and the software implementation (including package and version). These additions will allow readers to fully evaluate the model and reported results. revision: yes

  2. Referee: [Results] Results / Discussion: the directional treatment associations (corticosteroids and calcium-channel blockers increasing composite-event risk; immunosuppressives decreasing it) are presented without addressing confounding by indication, time-varying treatment effects, or sensitivity analyses, which is load-bearing for the highlighted clinical claims in an observational single-center study.

    Authors: We acknowledge that confounding by indication is a valid concern in this observational single-center study and that the reported treatment associations should not be interpreted as causal. We will revise the Results and Discussion to explicitly note the observational nature of the data, discuss the potential for confounding by indication (e.g., more severely affected patients receiving corticosteroids), and highlight the absence of granular time-varying treatment information that would permit sensitivity analyses such as propensity-score methods. These clarifications will temper the clinical interpretation while retaining the reported associations as adjusted model outputs. revision: partial

  3. Referee: [Abstract/Methods] Abstract and Methods: no information is given on missing-data mechanisms for the longitudinal creatinine measurements, model diagnostics, or any form of internal or external validation of the dynamic predictions beyond AUC/Brier scores on the same data.

    Authors: We will add the requested information to the Methods section of the revised manuscript. We will state that the joint model assumes a missing-at-random mechanism and naturally incorporates all observed creatinine measurements; describe model diagnostics including MCMC convergence diagnostics and posterior predictive checks; and clarify that the AUC and Brier scores were obtained via an internal landmarking approach on the full dataset. We will also note as a limitation that external validation is not feasible with the current single-center cohort. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the joint modeling derivation

full rationale

The paper applies a standard Bayesian joint model to link longitudinal creatinine trajectories with a composite time-to-event outcome, estimating associations and generating dynamic predictions from the fitted parameters and observed patient trajectories. This is conventional practice in joint modeling and does not reduce any claimed result to a tautology or self-definition by construction. No equations or steps in the provided description show fitted inputs relabeled as independent predictions, self-citation load-bearing premises, uniqueness theorems imported from the authors' prior work, or ansatzes smuggled via citation. The reported treatment associations and trajectory effects are direct outputs of the model fit on the single-center data; while external validity may be limited by design, the derivation chain itself remains independent of its inputs and self-contained.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard joint modeling assumptions and data-specific parameter estimates from a single-center cohort; no new entities are postulated.

free parameters (2)
  • association parameter linking longitudinal creatinine process to event hazard
    Estimated from data to quantify how creatinine trajectories influence event risk.
  • random effects variances and covariance structure
    Fitted to capture individual variability in creatinine trajectories.
axioms (2)
  • domain assumption The joint model correctly specifies the dependence between the longitudinal biomarker process and the time-to-event process.
    Core assumption required for the association and dynamic prediction claims.
  • domain assumption The single-center cohort is representative of the broader pediatric autoimmune population.
    Invoked implicitly for generalizability of risk predictions.

pith-pipeline@v0.9.0 · 5530 in / 1376 out tokens · 43611 ms · 2026-05-10T14:20:35.373548+00:00 · methodology

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

49 extracted references · 49 canonical work pages

  1. [1]

    On estimating the relationship between longitudinal measurements and time-to-event data using a simple two-stage procedure.Biometrics, 66(3):983–987, 2010

    Paul S Albert and Joanna H Shih. On estimating the relationship between longitudinal measurements and time-to-event data using a simple two-stage procedure.Biometrics, 66(3):983–987, 2010

  2. [2]

    Armero, A

    C. Armero, A. Forte, H. Perpi˜ n´ an, M. J. Sanahuja, and S. Agust´ ı. Bayesian joint modeling for as- sessing the progression of chronic kidney disease in children.Statistical Methods in Medical Research, 27(1):298–311, 2018

  3. [3]

    The ckid study: overview and summary of findings related to kidney disease progression.Pediatric Nephrology, 36(3):527–538, 2021

    Meredith A Atkinson, Derek K Ng, Bradley A Warady, Susan L Furth, and Joseph T Flynn. The ckid study: overview and summary of findings related to kidney disease progression.Pediatric Nephrology, 36(3):527–538, 2021

  4. [4]

    A bayesian approach for joint modeling of skew-normal longitudinal measurements and time to event data.REVSTAT-Statistical Journal, 13(2):169–191, 2015

    Taban Baghfalaki and Mojtaba Ganjali. A bayesian approach for joint modeling of skew-normal longitudinal measurements and time to event data.REVSTAT-Statistical Journal, 13(2):169–191, 2015

  5. [5]

    Bayesian joint modeling of longitudinal mea- surements and time-to-event data using robust distributions.Journal of biopharmaceutical statistics, 24(4):834–855, 2014

    Taban Baghfalaki, Mojtaba Ganjali, and Reza Hashemi. Bayesian joint modeling of longitudinal mea- surements and time-to-event data using robust distributions.Journal of biopharmaceutical statistics, 24(4):834–855, 2014

  6. [6]

    Samson Belay, Dessie Melese, and Kasim Muhammed. Joint modeling on serum creatinine and time to end stage of renal disease for chronic kidney disease patients under treatment at the university of gondar referral hospital.Health Science Reports, 6(9):e1563, 2023

  7. [7]

    Samson Belay, Dessie Melese, and Kasim Muhammed. Joint modeling on serum creatinine and time to end stage of renal disease for chronic kidney disease patients under treatment at the University of Gondar Referral Hospital.Health Science Reports, page e1563, 2023

  8. [8]

    Global, regional, and national burden of chronic kidney disease, 1990–2017: a systematic analysis for the global burden of disease study 2017.The lancet, 395(10225):709–733, 2020

    Boris Bikbov, Caroline A Purcell, Andrew S Levey, Mari Smith, Amir Abdoli, Molla Abebe, Oladimeji M Adebayo, Mohsen Afarideh, Sanjay Kumar Agarwal, Marcela Agudelo-Botero, et al. Global, regional, and national burden of chronic kidney disease, 1990–2017: a systematic analysis for the global burden of disease study 2017.The lancet, 395(10225):709–733, 2020

  9. [9]

    Paul Blanche, Cecile Proust-Lima, Lucie Loubere, Claudine Berr, Jean-Francois Dartigues, and Helene Jacqmin-Gadda. Quantifying and comparing dynamic predictive accuracy of joint models for longitudinal marker and time-to-event in presence of censoring and competing risks.Biometrics, 71:102–113, 2015

  10. [10]

    British 1990 growth reference centiles for weight, height, body mass index and head circumference fitted by maximum penalized likelihood

    Tim J Cole, Jenny V Freeman, and Michael A Preece. British 1990 growth reference centiles for weight, height, body mass index and head circumference fitted by maximum penalized likelihood. Statistics in medicine, 17(4):407–429, 1998

  11. [11]

    D. R. Cox. Regression models and life-tables.Journal of the Royal Statistical Society. Series B (Methodological), 34(2):187–220, 1972

  12. [12]

    Daniels and Joseph W

    Michael J. Daniels and Joseph W. Hogan.Missing Data in Longitudinal Studies: Strategies for Bayesian Modeling and Sensitivity Analysis. CRC Press, 2008

  13. [13]

    Predicting kidney graft failure using time-dependent renal function covariates.Journal of clinical epidemiology, 56(5):448–455, 2003

    Mattheus HJ de Bruijne, Yvo WJ Sijpkens, Leendert C Paul, Rudi GJ Westendorp, Hans C van Houwelingen, and Aeilko H Zwinderman. Predicting kidney graft failure using time-dependent renal function covariates.Journal of clinical epidemiology, 56(5):448–455, 2003. 16

  14. [14]

    Paul, Rudi G.J Westendorp, Hans C

    Mattheus H.J de Bruijne, Yvo W.J Sijpkens, Leendert C. Paul, Rudi G.J Westendorp, Hans C. van Houwelingen, and Aeilko H. Zwinderman. Predicting kidney graft failure using time-dependent renal function covariates.Journal of Clinical Epidemiology, 56(5):448–455, 2003

  15. [15]

    The benefit of early treatment with immunosuppressive agents in lupus nephritis.The Journal of Rheumatology, 21(11):2046–2051, 1994

    JM Esdaile, L Joseph, T MacKenzie, M Kashgarian, and JP Hayslett. The benefit of early treatment with immunosuppressive agents in lupus nephritis.The Journal of Rheumatology, 21(11):2046–2051, 1994

  16. [16]

    Immunosuppressive treatment of the glomerulonephritis of systemic lupus.Clinical nephrology, 44(6):367–375, 1995

    R Faedda, D Palomba, A Satta, Mario Pirisi, F Tanda, and E Bartoli. Immunosuppressive treatment of the glomerulonephritis of systemic lupus.Clinical nephrology, 44(6):367–375, 1995

  17. [17]

    Glucocorticoids in renal disease: Theoretical basis, consequences and efficacy of use in the pediatric patient.American journal of nephrology, 2(6):330– 341, 1982

    Aaron L Friedman and Russell W Chesney. Glucocorticoids in renal disease: Theoretical basis, consequences and efficacy of use in the pediatric patient.American journal of nephrology, 2(6):330– 341, 1982

  18. [18]

    Model determination using predictive distributions with imple- mentation via sampling-based methods.Bayesian Statistics, 4:147–167, 1994

    Alan E Gelfand and Dipak K Dey. Model determination using predictive distributions with imple- mentation via sampling-based methods.Bayesian Statistics, 4:147–167, 1994

  19. [19]

    Andrew Gelman and Donald B. Rubin. Inference from iterative simulation using multiple sequences. Statistical Science, 7(4):457–472, 1992

  20. [20]

    Separate and joint modeling of longitudinal and event time data using standard computer packages.The american statistician, 58(1):16–24, 2004

    Xu Guo and Bradley P Carlin. Separate and joint modeling of longitudinal and event time data using standard computer packages.The american statistician, 58(1):16–24, 2004

  21. [21]

    Hypertension in children with chronic kidney disease: pathophysiology and management.Pediatric Nephrology, 23(3):363–371, 2008

    Charlotte Hadtstein and Franz Schaefer. Hypertension in children with chronic kidney disease: pathophysiology and management.Pediatric Nephrology, 23(3):363–371, 2008

  22. [22]

    Jager, Csaba Kovesdy, Robyn Langham, Mark Rosenberg, Vivekanand Jha, and Carmine Zoccali

    Kitty J. Jager, Csaba Kovesdy, Robyn Langham, Mark Rosenberg, Vivekanand Jha, and Carmine Zoccali. A single number for advocacy and communication—worldwide more than 850 million indi- viduals have kidney diseases.Kidney International, 96(5):1048–1050, 2019

  23. [23]

    Csaba P. Kovesdy. Epidemiology of chronic kidney disease: an update 2022.Kidney International Supplements, 12(1):7–11, 2022. Aldosterone and the Mineralocorticoid Receptor: Emerging Thera- peutic Paradigms for Cardiorenal Protection

  24. [24]

    Children with CKD Are Not Little Adults with CKD: Pediatric Considerations for the Advancing American Kidney Health Initiative.CJASN, 3:470–472, March 2021

    Michael J.G.2; on behalf of the American Society of Pediatric Kula, Alexander J.1; Somers. Children with CKD Are Not Little Adults with CKD: Pediatric Considerations for the Advancing American Kidney Health Initiative.CJASN, 3:470–472, March 2021

  25. [25]

    Random-effects models for longitudinal data.Biometrics, 38(4):963–974, 1982

    Nan M Laird and James H Ware. Random-effects models for longitudinal data.Biometrics, 38(4):963–974, 1982

  26. [26]

    An approach for personalized dynamic assessment of chronic kidney disease progression using joint model.Biomedicines, 12(3):622, 2024

    Chen-Mao Liao, Yi-Wei Kao, Yi-Ping Chang, and Chih-Ming Lin. An approach for personalized dynamic assessment of chronic kidney disease progression using joint model.Biomedicines, 12(3):622, 2024

  27. [27]

    John Wiley & Sons, 2019

    Roderick JA Little and Donald B Rubin.Statistical analysis with missing data. John Wiley & Sons, 2019

  28. [28]

    Incidence of initial renal replacement therapy over the course of kidney disease in children.American journal of epidemiology, 188(12):2156–2164, 2019

    Derek K Ng, Matthew B Matheson, Bradley A Warady, Susan R Mendley, Susan L Furth, and Alvaro Mu˜ noz. Incidence of initial renal replacement therapy over the course of kidney disease in children.American journal of epidemiology, 188(12):2156–2164, 2019

  29. [29]

    An overview of joint modeling of time-to-event and longitudinal outcomes.Annual Review of Statistics and Its Application, 6:223–240, 2019

    Grigorios Papageorgiou, Katya Mauff, Anirudh Tomer, and Dimitris Rizopoulos. An overview of joint modeling of time-to-event and longitudinal outcomes.Annual Review of Statistics and Its Application, 6:223–240, 2019

  30. [30]

    Glucocorticoids in the treatment of glomerular diseases: pitfalls and pearls.Clinical Journal of the American Society of Nephrology, 13(5):815–822, 2018

    Claudio Ponticelli and Francesco Locatelli. Glucocorticoids in the treatment of glomerular diseases: pitfalls and pearls.Clinical Journal of the American Society of Nephrology, 13(5):815–822, 2018. 17

  31. [31]

    Dynamic predictions and prospective accuracy in joint models for longitudinal and time-to-event data.Biometrics, 67(3):819–829, 2011

    Dimitris Rizopoulos. Dynamic predictions and prospective accuracy in joint models for longitudinal and time-to-event data.Biometrics, 67(3):819–829, 2011

  32. [32]

    CRC Press, Boca Raton, FL, 2012

    Dimitris Rizopoulos.Joint Models for Longitudinal and Time-to-Event Data: With Applications in R. CRC Press, Boca Raton, FL, 2012

  33. [33]

    The r package jmbayes for fitting joint models for longitudinal and time-to-event data using mcmc.Journal of statistical software, 72:1–46, 2016

    Dimitris Rizopoulos. The r package jmbayes for fitting joint models for longitudinal and time-to-event data using mcmc.Journal of statistical software, 72:1–46, 2016

  34. [34]

    Hatfield, Bradley P

    Dimitris Rizopoulos, Laura A. Hatfield, Bradley P. Carlin, and Johanna J. M. Takkenberg. Com- bining dynamic predictions from joint models for longitudinal and time-to-event data using bayesian model averaging.Journal of the American Statistical Association, 109(508):1385–1397, 2014

  35. [35]

    R package version 0.6-0

    Dimitris Rizopoulos, Pedro Miranda-Afonso, and Grigorios Papageorgiou.JMbayes2: Extended Joint Models for Longitudinal and Time-to-Event Data, 2026. R package version 0.6-0

  36. [36]

    Dynamic predictions with time-dependent covariates in survival analysis using joint modeling and landmarking.Biometrical Journal, 59(6):1261–1276, 2017

    Dimitris Rizopoulos, Geert Molenberghs, and Emmanuel MEH Lesaffre. Dynamic predictions with time-dependent covariates in survival analysis using joint modeling and landmarking.Biometrical Journal, 59(6):1261–1276, 2017

  37. [37]

    John Wiley & Sons, New York, NY, 1987

    Donald B Rubin.Multiple Imputation for Nonresponse in Surveys. John Wiley & Sons, New York, NY, 1987

  38. [38]

    Abbott, and Jennifer Bragg-Gresham et al

    Rajiv Saran, Bruce Robinson, Kevin C. Abbott, and Jennifer Bragg-Gresham et al. Us renal data system 2019 annual data report: Epidemiology of kidney disease in the united states.American Journal of Kidney Diseases, 75(1, Supplement 1):A6–A7, 2020. US Renal Data System 2019 Annual Data Report

  39. [39]

    Hojjat Sayyadi, Farid Zayeri, Ahmad Reza Baghestani, Taban Baghfalaki, Ali Taghizadeh Afshari, Mohsen Mohammadrahimi, Javid Fereidoni, and Khadijeh Makhdoomi. Assessing risk indicators of allograft survival of renal transplant: An application of joint modeling of longitudinal and time-to- event analysis.Iran Red Crescent Med J, 19(3):e40583, 2017

  40. [40]

    A prognostic tool for individualized prediction of graft failure risk within ten years after kidney transplantation

    D Stamenic, A Rousseau, M Essig, P Gatault, M Buchler, M Filloux, P Marquet, and A Pr´ emaud. A prognostic tool for individualized prediction of graft failure risk within ten years after kidney transplantation. j transplant. 2019; 2019: 7245142

  41. [41]

    Taylor, Y

    J.M.G. Taylor, Y. Park, et al. Real-time individual predictions of prostate cancer recurrence using joint models.Biometrics, 69(1):206–214, 2013

  42. [42]

    A. A. Tsiatis and M. Davidian. Joint modeling of longitudinal and time-to-event data: An overview. Statistica Sinica, 14:809–834, 2004

  43. [43]

    Joint modeling of longitudinal and time-to-event data: an overview.Statistica Sinica, pages 809–834, 2004

    Anastasios A Tsiatis and Marie Davidian. Joint modeling of longitudinal and time-to-event data: an overview.Statistica Sinica, pages 809–834, 2004

  44. [44]

    Verbeke and G

    G. Verbeke and G. Molenberghs.Linear Mixed Models for Longitudinal Data. Springer, 2000

  45. [45]

    R package version 0.9.11

    Mandy Vogel.childsds: Data and Methods Around Reference Values in Pediatrics, 2025. R package version 0.9.11

  46. [46]

    Asymptotic equivalence of bayes cross validation and widely applicable information criterion in singular learning theory.Journal of Machine Learning Research, 11:3571–3594, 2010

    Sumio Watanabe. Asymptotic equivalence of bayes cross validation and widely applicable information criterion in singular learning theory.Journal of Machine Learning Research, 11:3571–3594, 2010

  47. [47]

    Plasma soluble urokinase plasminogen activator receptor (supar) and ckd progression in children

    Darcy K Weidemann, Alison G Abraham, Jennifer L Roem, Susan L Furth, and Bradley A Warady. Plasma soluble urokinase plasminogen activator receptor (supar) and ckd progression in children. American Journal of Kidney Diseases, 76(2):194–202, 2020

  48. [48]

    Are all antihypertensive drugs renoprotective?Herz, 29(3):248, 2004

    Sabine Wolf and Teut Risler. Are all antihypertensive drugs renoprotective?Herz, 29(3):248, 2004

  49. [49]

    M. S. Wulfsohn and A. A. Tsiatis. A joint model for survival and longitudinal data measured with error.Biometrics, 53(1):330–339, 1997. 18 Appendix A Table A. 1: Classification of medications commonly used by study participants, grouped by therapeutic target and subdivided by mechanism of action. A binary covariate was created for each subgroup for modell...