The aim of this study was to develop a multistate model for overall survival analysis, based on parametric hazard functions and combined with an investigation of predictors derived from a longitudinal tumor size model on the transition hazards. Different states -stable disease, tumor response, progression, second-line treatment and death following docetaxel treatment initiation (stable state) in HER2-negative breast cancer patients (n=183) were used in model building. Past changes in tumor size prospectively predicts the probability of state changes. The hazard of death after progression was lower for subjects who had longer treatment response, i.e. longer time-to-progression. Young age increased the probability of receiving second-line treatment. The developed multistate model adequately described the transitions between different states and jointly the overall event and survival data. The multistate model allows for simultaneous estimation of transition rates along with their tumor model derived metrics. The metrics were evaluated in a prospective manner so not to cause immortal time bias. Investigation of predictors and characterization of the time to develop response, the duration of response, the progression-free survival and the OS can be performed in a single multistate modeling exercise. This modeling approach can be applied to other cancer types and therapies to provide a better understanding of efficacy of drug and characterizing different states, thereby facilitating early clinical interventions to improve anticancer therapy.