Finally, multivariable regression was used to examine changes in sources of prescription opioids used for non-medical reasons over time within the POUD group, modeling source as a function of time while adjusting for covariates. Models were adjusted for past-year non-heroin illicit drug use and past-year heroin use in supplementary analyses.