Using a Bayesian model of the joint distribution of pain and time on medication to decide on pain medication for neuropathy.

The PAIN-CONTRoLS trial compared four medications in treating Cryptogenic sensory polyneuropathy. The primary outcome was a utility function that combined two outcomes, patients' pain score reduction and patients' quit rate. However, additional analysis of the individual outcomes could also be leveraged to inform selecting an optimal medication for future patients. We demonstrate how joint modeling of longitudinal and time-to-event data from PAIN-CONTRoLS can be used to predict the effects of medication in a patient-specific manner and helps to make patient-focused decisions. A joint model was used to evaluate the two outcomes while accounting for the association between the longitudinal process and the time-to-event processes. Results suggested no significant association between the patients' pain scores and time to the medication quit in the PAIN-CONTRoLS study, but the joint model still provided robust estimates and a better model fit. Using the model estimates, given patients' baseline characteristics, a drug profile on both the pain reduction and medication time could be obtained for each drug, providing information on how likely they would quit and how much pain reduction they should expect. Our analysis suggested that drugs viable for one patient may not be beneficial for others.