We characterized the association between tumor size kinetics and survival in patients with advanced urothelial carcinoma treated with atezolizumab (anti‐programmed death‐ligand 1, Tecentriq) using a joint model. The model, developed on data from 309 patients of a phase II clinical trial, identified the time‐to‐tumor growth and the instantaneous changes in tumor size as the best on‐treatment predictors of survival. On the validation dataset containing data from 457 patients from a phase III study, the model predicted individual survival probability using 3‐month or 6‐month tumor size follow‐up data with an area under the receptor‐occupancy curve between 0.75 and 0.84, as compared with values comprised between 0.62 and 0.75 when the model included only information available at treatment initiation. Including tumor size kinetics in a relevant statistical framework improves the prediction of survival probability during immunotherapy treatment and may be useful to identify most‐at‐risk patients in “real‐time.”

In metastatic castration-resistant prostate cancer (mCRPC) clinical trials, the assessment of treatment efficacy essentially relies on the time to death and the kinetics of prostate-specific antigen (PSA). Joint modeling has been increasingly used to characterize the relationship between a time to event and a biomarker kinetics, but numerical difficulties often limit this approach to linear models. Here, we evaluated by simulation the capability of a new feature of the Stochastic Approximation Expectation-Maximization algorithm in Monolix to estimate the parameters of a joint model where PSA kinetics was defined by a mechanistic nonlinear mixed-effect model. The design of the study and the parameter values were inspired from one arm of a clinical trial. Increasingly high levels of association between PSA and survival were considered, and results were compared with those found using two simplified alternatives to joint model, a two-stage and a joint sequential model. We found that joint model allowed for a precise estimation of all longitudinal and survival parameters. In particular, the effect of PSA kinetics on survival could be precisely estimated, regardless of the strength of the association. In contrast, both simplified approaches led to bias on longitudinal parameters, and two-stage model systematically underestimated the effect of PSA kinetics on survival. In summary, we showed that joint model can be used to characterize the relationship between a nonlinear kinetics and survival. This opens the way for the use of more complex and physiological models to improve treatment evaluation and prediction in oncology.