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
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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.
- [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)
- [Abstract] Abstract: the phrase 'model selection criteria' is used without naming the specific criteria (DIC, WAIC, etc.) employed.
- [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
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
-
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
-
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
-
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
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
free parameters (2)
- association parameter linking longitudinal creatinine process to event hazard
- random effects variances and covariance structure
axioms (2)
- domain assumption The joint model correctly specifies the dependence between the longitudinal biomarker process and the time-to-event process.
- domain assumption The single-center cohort is representative of the broader pediatric autoimmune population.
Reference graph
Works this paper leans on
-
[1]
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
work page 2010
- [2]
-
[3]
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
work page 2021
-
[4]
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
work page 2015
-
[5]
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
work page 2014
-
[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
work page 2023
-
[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
work page 2023
-
[8]
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
work page 1990
-
[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
work page 2015
-
[10]
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
work page 1990
-
[11]
D. R. Cox. Regression models and life-tables.Journal of the Royal Statistical Society. Series B (Methodological), 34(2):187–220, 1972
work page 1972
-
[12]
Michael J. Daniels and Joseph W. Hogan.Missing Data in Longitudinal Studies: Strategies for Bayesian Modeling and Sensitivity Analysis. CRC Press, 2008
work page 2008
-
[13]
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
work page 2003
-
[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
work page 2003
-
[15]
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
work page 2046
-
[16]
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
work page 1995
-
[17]
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
work page 1982
-
[18]
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
work page 1994
-
[19]
Andrew Gelman and Donald B. Rubin. Inference from iterative simulation using multiple sequences. Statistical Science, 7(4):457–472, 1992
work page 1992
-
[20]
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
work page 2004
-
[21]
Charlotte Hadtstein and Franz Schaefer. Hypertension in children with chronic kidney disease: pathophysiology and management.Pediatric Nephrology, 23(3):363–371, 2008
work page 2008
-
[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
work page 2019
-
[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
work page 2022
-
[24]
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
work page 2021
-
[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
work page 1982
-
[26]
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
work page 2024
-
[27]
Roderick JA Little and Donald B Rubin.Statistical analysis with missing data. John Wiley & Sons, 2019
work page 2019
-
[28]
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
work page 2019
-
[29]
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
work page 2019
-
[30]
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
work page 2018
-
[31]
Dimitris Rizopoulos. Dynamic predictions and prospective accuracy in joint models for longitudinal and time-to-event data.Biometrics, 67(3):819–829, 2011
work page 2011
-
[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
work page 2012
-
[33]
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
work page 2016
-
[34]
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
work page 2014
-
[35]
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
work page 2026
-
[36]
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
work page 2017
-
[37]
John Wiley & Sons, New York, NY, 1987
Donald B Rubin.Multiple Imputation for Nonresponse in Surveys. John Wiley & Sons, New York, NY, 1987
work page 1987
-
[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
work page 2019
-
[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
work page 2017
-
[40]
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
work page 2019
- [41]
-
[42]
A. A. Tsiatis and M. Davidian. Joint modeling of longitudinal and time-to-event data: An overview. Statistica Sinica, 14:809–834, 2004
work page 2004
-
[43]
Anastasios A Tsiatis and Marie Davidian. Joint modeling of longitudinal and time-to-event data: an overview.Statistica Sinica, pages 809–834, 2004
work page 2004
-
[44]
G. Verbeke and G. Molenberghs.Linear Mixed Models for Longitudinal Data. Springer, 2000
work page 2000
-
[45]
Mandy Vogel.childsds: Data and Methods Around Reference Values in Pediatrics, 2025. R package version 0.9.11
work page 2025
-
[46]
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
work page 2010
-
[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
work page 2020
-
[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
work page 2004
-
[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...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.