Because secular trends in the use of cannabis are unlikely to be linear, we re-ran our models with linear and quadratic age adjustments on the PGD and CU thresholds. We found an identical pattern of results; the uni-directional PGD→CU model was rejected while the CU→PGD causal model provided the best fit to the data.