Combined effects of T_eff, t̄_★, [Fe/H], and direction of motion on the predicted survival time (in this case, the predicted A(Li)), accounting for all covariates simultaneously. Predictions are shown for stars with detected A(Li) measurements (orange), upper limits (purple), and missing A(Li) values (black), with colours consistent across all figures. The curves represent the model’s output from the parametric survival regression (survreg with logistic distribution), integrating the contributions of all covariates.