Predicted median eGFR trajectories from the multivariable linear mixed model for eGFR levels (with baseline eGFR as part of the dependent variable) obtained by AIC-based backward elimination (solid line). Shaded areas indicate the interquartile range of predictions. Superimposed box plots show the observed values summarized at each yearly FU visit.