R code produced using the 'dlnm' package. Intended for reproduction of main results