The remaining model parameters are fitted to smoking and vaping prevalence data; we assume the remaining parameters may differ between adult and youth populations due to both physiological and social differences. Best fit parameters were estimated by minimizing the error between the model and the data using the Nelder-Mead algorithm [1]}. To reduce the risk of overfitting, we fit the model to each dataset in two steps.
