Author: Camilo Vargas
- Uruguay was the first country to legalize cannabis for recreational consumption in 2013. Its effects have been a public health concern in the last decade across the continent as more countries examine similar strategies.
- This paper aims to measure the post-legalization impact on the prevalence of cannabis use disorder among Uruguayan consumers by constructing a synthetic control method covering data from 1990 to 2017 through 19 Latin American countries.
- The findings suggest that the prevalence of cannabis use disorder decreased in Uruguay after the legalization. However, these results are not unique since placebo tests indicate that other countries would have experienced similar outcomes.
- Thus, it cannot be confirmed that the decreasing effect in the prevalence of cannabis disorders is only due to the legalization in 2013.
- The added value of this paper is to add empirical econometric research to public health studies in Latin America. For further analysis and discussion, the full printed version of this project is available upon request.
- Abuse of illegal substances is a global struggle that significantly impacts individuals and society.
- Roth et al. (2018) reported that worldwide smoking, alcohol, and illicit drug use kill directly and indirectly 11.8 million people annually, exceeding the total number of deaths caused by all types of cancer.
- Direct deaths from alcohol and illicit drug use disorders estimate that over 350,000 people die each year globally, with over half of these deaths occurring in males younger than 50.
- Substance abuse also contributes to the global disease burden, with 1.5% resulting from alcohol and illicit drug addiction; in some countries like the U.S., it is over 5%. Additionally, over 2% of the world population has an alcohol or illicit drug addiction.
Panel data across 27 years from 1990 to 2017 for a total sample of 19 Latin American countries.
- Data on the prevalence of drug use disorders is drawn from the Institute for Health Metrics and Evaluation (IHME), Global Burden of Disease (GBD 2019). These data measure the number of persons of both sexes and all ages per country and year that have reported the prevalence use disorder of the following substances: amphetamine, cannabis, cocaine, opioid and other drugs.
- Public healthcare indicators per country and year are drawn from the World Health Organization Global Health Expenditure database.
- Macroeconomic and demographic variables per country and year are drawn from the World Bank – World Development Indicators.
The descriptive statistics of the variables are summarized as follows.
This paper uses a counterfactual analytical approach to evaluate what would have been the outcome in the prevalence of cannabis use disorders if Uruguay had not legalized cannabis for recreational use in 2013 by implementing the synthetic control method (SCM) and following the guidelines from Abadie (2021), creating a combination of weights during the pre-legalization period from 18 non-treated Latin American countries, which are considered the donor pool or the control group because these countries did not approve of cannabis legalization during the period of study. This weighted combination, the synthetic control, closely estimates a counterfactual scenario without legalization. Then the synthetic trend is forecasted in the post-legalization period to be compared with treated Uruguay.
Figure below illustrates the prevalence of drug use disorders detailed by substances in some of the most populous Latin American countries measured as a per capita share of their population. From there, it can be assumed that Uruguay’s prevalence of cannabis disorder has remained constant at 0.47%; on average, 15823 people reported this condition annually. The figure also highlights that cannabis disorders are the most prevalent among the observed countries.
#delimit;
use ${RAW}DataCannabisLatam.dta, clear ;
keep if Code == "URY" | Code == "ARG" | Code == "COL" | Code == "CHL" | Code == "BRA" | Code == "MEX" ;
twoway line
PrevalenceCannabisShare PrevalenceCocaineShare PrevalenceAmphetamineShare PrevalenceOpioidShare PrevalenceOtherDrugShare PrevalenceAlcoholusedisorde
Year ,
sort lwidth(thick)
by(Entity)
scheme(s1color) ytitle(Share of Pop. (%))
legend(off rows(2) forces size(small)
lab(1 "Cannabis") lab(2 "Cocaine") lab(3 "Amphetamine") lab(4 "Opioid") lab(5 "Other Drug") lab(6 "Alcohol"))
xlabel(, labsize(small) )
;
#delimit cr
In the SCM, the control group and the experimental unit are equally balanced pre-intervention on several predictors that could predict the dependent variable, creating a quasi-experimental setting. In this current study, the period before intervention was from 1990 to 2012, and the posterior intervention was from 2013 to 2017. The predictors are the prevalence use disorder of illegal substances like amphetamine, cocaine, opioid, alcohol and other public health indicators described in the previous section of the text. The primary outcome of concern has been stated to be the prevalence of cannabis use disorder.
* Install the synthetic control method package:
ssc install synth, all
*net install synth_runner, from(https://raw.github.com/bquistorff/synth_runner/master/) replace
* Set the data as time series, using tsset
use ${RAW}DataCannabisLatam.dta, clear
tsset CodeCountry Year
* Perform synt control method:
#delimit;
synth PrevalenceCannabisShare // depvar
PrevalenceCannabisShare(1990) PrevalenceCannabisShare(2000)
PrevalenceCannabisShare(2012) PrevalenceOtherDrugShare(2012)
PrevalenceAmphetamineShare(2012) PrevalenceCocaineShare(2012)
PrevalenceOpioidShare(2012) PrevalenceAlcoholusedisorde(2012)
Currenthealthexpenditureof DomGenGovHeaExpGDP(2012)
Domesticprivatehealthexpendit(2012) Outofpocketexpenditureof(2012)
lnGDPpercapita PopulationtotalSPPOPTOTL , // predictors
trunit(858) // trunit(858) = Uruguay
trperiod(2013) unitnames(Entity)
mspeperiod(1990(1)2013) resultsperiod(1990(1)2017)
fig
nested keep(${RAW}synth_results_data.dta) replace
;
#delimit cr
The table below shows the relative contribution of each of the 18 countries and their respective weights to the synthetic control of Uruguay. The synthetic version is a weighted average of mainly Panama and Suriname, followed by Belize, Chile and Colombia, with weights in decreasing order.
The figure below displays the results of the synthetic control method examining the effects of legalizing cannabis in Uruguay on the prevalence of cannabis disorder. Synthetic Uruguay, composed of the five countries from the donor pool indicated previously, traced well the trend beside actual Uruguay’s outcome throughout the pre-legalization period (1990 – 2012). However, after the policy was enacted (2013), the true prevalence of cannabis use disorder fell remarkably below 0.47%, whereas the synthetic estimates that the prevalence will increase at pre-legalization levels above 0.48%.
From the previous results, the difference between treatment and counterfactual can be computed as shown next. Therefore, the evidence suggests that the treatment effect post-legalization of cannabis in Uruguay caused a decrease in the share of the population that reported the prevalence of cannabis disorder compared to the synthetic control group.
use ${RAW}synth_results_data.dta, clear
drop _Co_Number _W_Weight // Drops the columns that store the donor state weights
gen GapUry = _Y_treated - _Y_synthetic // The counterfactual is _Y_synthetic
twoway line (GapUry _time), scheme(s1color) lcolor(green) lwidth(thick) xline(2013) xmlabel(2013 "Legalization") yline(0)
xtitle(Year) ytitle(Gap in Prevalence Cannabis Prediction Error) legend(on) // title("Difference between treatment and counterfactual")
To measure the change in the outcome of interest, i.e., the prevalence of cannabis disorder, it can be inferred by comparing the pre and post-treatment performance of the predictors balance with that of the synthetic control. Table below summarizes the considered predictors' results, which can be interpreted as the treatment effects.
* Predictors – Pre Treatment "2012" are taken from previous steps.
use ${RAW}DataCannabisLatam.dta, clear
tsset CodeCountry Year
* Predictors – Post Treatment. "2017"
#delimit;
synth PrevalenceCannabisShare
PrevalenceCannabisShare(2017) PrevalenceOtherDrugShare(2017)
PrevalenceAmphetamineShare(2017) PrevalenceCocaineShare(2017)
PrevalenceOpioidShare(2017) PrevalenceAlcoholusedisorde(2017)
Currenthealthexpenditureof DomGenGovHeaExpGDP(2017)
Domesticprivatehealthexpendit(2017) Outofpocketexpenditureof(2017)
lnGDPpercapita PopulationtotalSPPOPTOTL
, trunit(858) // trunit(858) = Uruguay
trperiod(2013) unitnames(Entity)
mspeperiod(1990(1)2015) resultsperiod(1990(1)2017) //
// fig
// nested keep(${RAW}synth_results_data.dta) replace
;
#delimit cr
Finally, to test whether the cannabis legalization policy had a statistically significant effect on Uruguay’s consumers with this drug disorder, several placebo estimates will be performed by assigning the same treatment period to all 18 potential donors, recalculating the model’s coefficients, and collecting them into a distribution which is then used for the analysis.
Table below displays a set of RMSPE values for the pre and post-legalization period across all countries as a test statistic used for inference. From there, it can be argued that Peru and Costa Rica have statistically significant p-values at a 10% confidence level, implying that Uruguay’s legalization effect was not unique.
* Inference 2: Estimate the pre- and post-RMSPE and calculate the ratio of post-pre RMSPE
foreach i of global Countries {
use ${RAW}synth_gap_`i', clear
gen gap3 = gap`i'* gap`i'
egen postmean = mean(gap3) if year > 2013
egen premean = mean(gap3) if year <= 2013
gen rmspe = sqrt(premean) if year <= 2013
replace rmspe = sqrt(postmean) if year > 2013
gen ratio = rmspe/rmspe[_n-1] if year == 2014
gen rmspe_post = sqrt(postmean) if year > 2013
gen rmspe_pre = rmspe[_n-1] if year == 2014
mkmat rmspe_pre rmspe_post ratio if year == 2014, matrix (country`i')
}
* Show pre/post expansion RMSPE ratio for all countries
* Table 3: Pre - Post expansion RMSPE ratio.
foreach i of global Countries {
matrix rownames country`i' = `i'
matlist country`i', names(row)
}
The placebo gaps across 18 control countries and Uruguay are illustrated below, assuming they were exposed to the same legalization of cannabis in 2013. Agreeing that there are no meaningful differences in the reported prevalence of cannabis use disorders post-legalization in Uruguay compared to the estimates from placebo tests, it could be claimed that the legalization in Uruguay had no significant effect according to this test.
This study found that legalizing cannabis in Uruguay was not associated with an increase in the prevalence of cannabis use disorders among its consumers, as was expected. In contrast, evidence suggests that the actual outcome decreased after the legalization period in 2013. However, this result is insignificant since other countries would have experienced the same results if exposed to the treatment. As such, it is unlikely that cannabis legalization alone was the cause of the reduction in the prevalence of cannabis use disorders in Uruguay.
- Abadie, Diamond, & Hainmueller. (2010). Synthetic Control Methods for Comparative Case Studies: Estimating the Effect of California’s Tobacco Control Program. Journal of the American Statistical Association.
- Abadie, A., Diamond, A., Hainmueller, J., (2011). "SYNTH: Stata module to implement Synthetic Control Methods for Comparative Case Studies," Statistical Software Components S457334, Boston College Department of Economics revised 09 May 2020.
- Abadie, Diamond, & Hainmueller. (2015). Comparative politics and the synthetic control method. American Journal of Political Science, 59(2), 495-510.
- Abadie, A. (2021). Using synthetic controls: Feasibility, data requirements, and methodological aspects. Journal of Economic Literature, 59(2), 391–425.
- Cunningham, S. (2021). Causal inference: The mixtape. Yale university press.
- Roth, G. A., Abate, D., Abate, K. H., Abay, S. M., Abbafati, C., Abbasi, N., … & Abdollahpour, I. (2018). Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980–2017: a systematic analysis for the Global Burden of Disease Study 2017. The Lancet, 392(10159), 1736-1788.
To construct the SCM for this paper, the synth package in Stata developed by Abadie et al. (2011) was implemented following the guidelines of Cunningham (2021).
For further analysis and discussion, the full printed version of this project is available upon request.