Globally, it is estimated that there may be as many as 15.9 million (range 11.0–21.2 million) people who inject drugs (PWID), with as many as 3 million (range 0.8–6.6) being HIV-positive ; this accounts for roughly 13% of all people with HIV/AIDS outside of sub-Saharan Africa [2,3]. The largest numbers of injectors were found in China, the USA and Russia, where estimates of HIV-prevalence among PWID were 12, 16 and 37%, respectively . Past studies have shown that HAART can reduce plasma HIV RNA thereby reducing morbidity and mortality and reducing HIV transmission [4–6], with a specific demonstration of the effect of diminishing plasma viral load on HIV incidence within a community of urban PWID in a Canadian setting . Unfortunately, PWID are known to have poor access and adherence to HAART. A recently published global survey suggests that only four out of every 100 HIV-positive PWID (range: 2–18) are accessing HAART, with many of those concentrated in Europe, and sparse reporting in many regions .
Opioid substitution treatment (OST) with methadone or buprenorphine is currently indicated for the treatment of opioid dependence in most developed-world countries . In British Columbia (BC) Canada, methadone is covered for reimbursement under standard provincial drug benefit formularies, although buprenorphine is available as a second-line option at the physician's discretion. There is clearly a potential for positive synergies between OST and HAART for HIV-positive, opioid-dependent individuals, yet OST remains highly restricted in many settings . Several studies have reported that methadone treatment was independently associated with decreased time to HAART initiation [8,9]. Engagement in OST has also been associated with decreases in HIV risk behaviour and adherence to HAART and resulting health gains [10–14]. Once HAART adherent, PWID have been noted to have a similar risk of mortality to noninjectors . Although OST has no direct effect on HIV disease progression or HIV-related mortality, it has the potential to improve HAART uptake and adherence, which impact these more distal outcomes.
Measuring the causal impact, or the true extent to which OST may influence HAART adherence is, however, complicated by competing hypotheses on the effect of OST on HAART, and HAART on OST. Specifically, OST is hypothesized to influence HAART adherence via the stability provided to its clients through reduced need to constantly address opioid withdrawal and craving . Linkage to the healthcare system afforded by stable access to OST may further influence HAART adherence; however, the same hypothesis can be made for the effect of HAART on OST . Marginal structural models are designed specifically to handle time-varying confounding effects such as that of HAART on OST in our context, allowing for causal interpretations of results.
Thus, our objective is to determine the causal effect of OST exposure, as well as a range of other individual-level factors, on HAART medication refill adherence among HIV-positive PWID using data from a linked population-level database with comprehensive OST and HAART dispensation records.
Materials and methods
This study was based on a provincial-level linkage of seven health administrative databases and disease registries, including the province-wide antiretroviral dispensation, virology and HIV-testing registries [antiretroviral dispensations, plasma viral load (pVL), CD4+ cell count tests, nominal HIV diagnoses], the Medical Services Plan database (physician billing records), the discharge abstract database (hospitalizations), the BC PharmaNet database (nonantiretroviral drug dispensations) and the BC Vital statistics database (deaths). Further details regarding the construction of the HIV-positive cohort and available databases were described elsewhere .
We selected all individuals identified as HIV-positive and either identified as having a history of injection drug use at diagnosis or having accessed OST before HAART eligibility, as indicated by methadone or buprenorphine dispensation records in the BC PharmaNet database, from 1 April 1996 to 31 March 2010. OST dispensations were identified by unique drug identification numbers (DINs). These codes encompass current and past DINs used for the medications in question and are specific to indications for opioid dependence (as opposed to pain). Previous studies have defined methods to clean and process PharmaNet data for accurate estimation of OST episodes . OST dispensation record capture was complete for the province of BC.
For each individual, we consider follow-up data from the point of initial HAART eligibility, as determined by IAS-USA criteria from 1996 to 2010, based on CD4+ cell count and pVL tests and AIDS-defining illness after the HIV diagnosis date, or by HAART initiation. Individuals were followed up until death, administrative loss to follow-up or censorship (ongoing care as of 30 March 2010). Administrative loss was defined as no records in any of the linked databases for a period of at least 18 months. The database is thus arranged longitudinally, with one row per month of observation following the date of HAART eligibility, until the last observed contact date within the database or the end of the study follow-up period.
The study cohort was followed in a unique environment wherein medical care (HIV and non-HIV related), antiretroviral treatment and laboratory monitoring is fully subsidized by the BC provincial government for people with HIV/AIDS. Antiretroviral drugs are centrally distributed according to provincial treatment guidelines that have remained consistent with those put forward by the International AIDS Society-USA (‘IAS-USA’) since the summer of 1996  and until the most recent guidelines released in 2012 . Antiretroviral medication and viral load data capture in the BC-CfE registries are complete for the population of individuals with HIV/AIDS in the province, and CD4+ cell count data capture have previously been estimated at 80% . This study received approval from the University of British Columbia/Providence Healthcare Research Ethics Board.
The dependent variable in the study was HAART adherence, defined as at least 95% pharmacy refill adherence in a given calendar month. Exposure to OST, defined by at least 95% days of OST dispensed in a calendar month, was the key independent variable considered. We considered HAART medication refill adherence as our primary outcome as opposed to more distal outcomes such as viral suppression or HIV-related mortality because it is the most proximal outcome in the causal pathway between OST and HIV-related outcomes, and because OST is not known to have a direct effect on HIV disease progression or HIV-related mortality.
We also considered a selection of other baseline and time-varying covariates hypothesized to be associated with the dependent and primary independent variable. Baseline covariates included age (30–39, 40–49, 50+), sex, aboriginal ethnicity, BC health authority of residence, CD4+ cell count (<200, 200–499, ≥500 cells/μl), calendar year (1996–1999, 2000–2003, 2004–2010) and AIDS status, all made available through the BC-CfE HIV registry. The Charlson comorbidity index (CCI) score  (≥1 vs. 0) and OST history (number of prior OST episodes) at the time of HAART eligibility were also incorporated, as well as time-varying measures of CD4+ cell count, health authority of residence, time since HAART eligibility, AIDS status and CCI.
Marginal structural models are designed to handle cases in which time-dependent variables are simultaneously confounders of the effect of interest and are predicted by previous treatment. Past HAART medication refill adherence can be considered a time-dependent confounder for the effect of OST on future HAART adherence, as it may be hypothesized to not only predict future HAART adherence but also subsequent initiation of OST, and as past OST history is an independent predictor of subsequent HAART adherence. The hypothesized relationship is characterized in the directed acyclic graph illustrated in Fig. 1. We hypothesized that exposure to OST at time t affects HAART medication refill adherence at time t + 1, which in turn influences OST exposure at time t + 2. Time-invariant and time-varying confounding at time t affects OST exposure, HAART medication refill adherence and confounders at time t + 1, and those confounders at t + 1 are also predicted by OST exposure at t. This relationship is analogous to the longitudinal relationship between zidovudine therapy and CD4+ cell count, which has served as an example of a repeated-measures application of marginal structural modelling .
To control for time-varying confounding in the exposure-outcome relationship, inverse probability of treatment weights (IPTWs) were estimated for each time point of the study. Under the four assumptions of consistency, exchangeability, positivity and no misspecification of the model used to estimate the weights, weighting intuitively creates a pseudo-population in which the exposure is independent of the measured confounders . The pseudo-population is the result of assigning each participant a weight that is inversely proportional to the participant's probability of receiving her own exposure history. In such a pseudo-population, one can regress the outcome on the exposure using a conventional regression model that does not include the measured confounders as covariates. Fitting a model in the pseudo-population is equivalent to fitting a weighted model in the study population. The parameters of such weighted regression models, which equal the parameters of marginal structural models , can be used to estimate the average causal effect of exposure in the original study population.
Weighted estimation of the parameters of marginal structural models requires fitting several models: the structural (i.e. weighted) model; the exposure model; and the censoring model. The exposure and censoring models were used to estimate the IPTW by a pooled logistic regression. The IPTWs were constructed according to their standard formulation , and the model specification was presented in a supplementary appendix, https://links.lww.com/QAD/A670. The final structural model was performed using PROC GENMOD [generalized estimating equation (GEE) regression model using independent working covariate matrix] for HAART medication refill adherence (≥95%), with independent variables of OST (≥95%), and baseline covariates included . All analyses were executed using SAS version 9.3 (SAS Institute Inc., Cary, North Carolina, USA).
We considered a range of sensitivity analyses to protect against model misspecification. Correct model specification applies to exposure, censoring and structural (weighted) models. In the context of the structural model, we considered subgroup analyses on those having no OST experience prior to HAART eligibility and those having OST experience. We also considered beginning follow-up at HAART initiation, as well as alternate definitions of OST receipt, including OST adherence at the minimum effect dose (60 mg/day) [28,29], OST adherence during periods of maintenance treatment (excluding periods of dose titration and tapering)  and alternate adherence thresholds for both HAART and OST adherence.
In the context of the exposure/censoring models used to construct the IPTW, first, we compared stabilized and unstabilized IPTW. We also considered truncated (stabilized) weights, assessing the stability of the IPTW (i.e. evaluating mean, standard deviation of IPTW) for appropriateness. We included variables statistically associated with HAART medication refill adherence for exposure/censoring models. For the assumption of positivity or correct model specification to hold, the mean of IPTWs should be approximately equal to one.
Characteristics of study cohort
Study sample selection is outlined in Fig. 2. Among 12 349 HIV-positive individuals in the STOP HIV/AIDS cohort, a total of 2928 (23.7%) HIV-positive individuals were identified as ever having received OST or having a history of injection drug use at HIV diagnosis. After excluding 584 individuals whose HAART eligibility status could not be verified and another 492 due to insufficient follow-up and unknown injection drug use status at HAART eligibility (i.e. identified as PWID through OST receipt after HAART eligibility only), the selected study cohort consisted of 1852 HIV-positive PWID.
The study sample was 34% female, median-aged 35 years and 27% of known Aboriginal ethnicity (Table 1). At the time of HAART eligibility or HAART initiation, 20% were on OST, and 39% had ever accessed OST in the past; individuals had a median of 0 (interquartile range: 0–1) prior OST episodes. The majority (46%) of the cohort became HAART eligible pre-2000. By the end of follow-up, 50% of the cohort accessed OST.
Longitudinal exposure to opioid substitution treatment and HAART
The cohort was followed up for a median duration of 5.5 years (2.2–10.1), and 25% of the cohort died during follow-up (Table 2). Individuals spent a median 44% of follow-up on HAART, 0% of follow-up on OST and 0% on both forms of treatment, jointly. Most importantly, individuals spent a median of 56% of the time they were in OST receiving HAART; when out of OST, individuals were receiving HAART only 38% of the time.
Inverse probability of treatment weight estimation
Variables that were associated with the outcome were selected for models for weight estimation. Excluded variables were baseline measures of aboriginal ethnicity, health authority of residence, CCI, as well as time-varying CCI. Stabilized weights were substantially more stable, with means at 1.00 and smaller standard deviations (IPTW model results and weights are presented in a supplementary appendix, https://links.lww.com/QAD/A670). Truncating weights at the first and 99th or fifth and 95th percentile resulted in smaller standard deviations with means close to 1. Stabilized, truncated weights (at the first and 99th percentiles) were thus used in baseline models.
Causal effect of opioid substitution treatment on HAART
Results of the structural model are presented in Table 3. The unadjusted odds of HAART medication refill adherence given OST was 1.54 [95% confidence interval (95% CI) 1.36–1.75]. Adjusting for baseline covariates in a GEE model (with logit link and binomial distribution) yielded an adjusted odds ratio (AOR) of 1.96 (1.72–2.24), with similar results after including time-varying covariates (AOR 1.91, 95% CI 1.68–2.19). Finally, adjusting for the potential effect of prior HAART medication refill adherence on OST through the use of IPTW, the marginal structural model yielded an AOR of 1.68 (95% CI 1.48–1.92). The time-varying effect of prior HAART adherence thus inflated the estimated odds ratio by a magnitude of 0.23 vs. the fully adjusted GEE model. The intermediate results of models estimated to construct the IPTW were presented in the supplementary appendix, https://links.lww.com/QAD/A670.
Sensitivity analysis on the effect of opioid substitution treatment on HAART
Baseline results indicating a positive and statistically significant effect of OST on HAART were supported by the results of our sensitivity analyses. We considered two changes in the classification of OST exposure: requiring treatment at the minimum effective dose, and requiring that individuals be at a stable maintenance dose to be classified as exposed. In each case, AORs in the structural model decreased compared with the baseline model formulation with truncated IPTW, but remained statistically significantly positive. Further, models with alternate adherence thresholds for both HAART and OST adherence had similar results compared with the baseline model. We also considered subgroups of individuals beginning follow-up at HAART initiation, and individuals with and without prior OST exposure at HAART eligibility. Odds ratios were higher than in the baseline model formulation for those with no prior exposure to OST.
Using a linked population-level database capturing HIV-positive PWID in a Canadian province, this study has demonstrated a 68% increase in the odds of HAART medication refill adherence due to OST exposure. The results were robust to a number of sensitivity analyses focusing on the assumptions of correct model specification inherent in the selected modelling approach.
These results largely confirm prior studies showing a positive association between OST and HAART adherence (measured within the same time period) in settings in which office-based OST is available; however the role of continuing injection drug use (in particular heroin use) was unclear [10,11,13]. Specifically, Palepu et al. controlled for time-varying heroin use and found that OST was independently positively associated with adherence, with an AOR of 1.52 (95% CI 1.16–2.00). In contrast, Roux et al. found that the odds of nonadherence for those engaged in OST and not injecting were no different from those currently opioid-abstinent, but patients reporting drug injection during OST had a two-fold risk of nonadherence compared with those currently abstinent.
These analyses were executed through prospective cohort studies that were able to ascertain self-reported drug use alongside time-varying measures of OST exposure and HAART medication refill adherence. However, in both studies, their measurement of OST exposure was coarse, capturing ‘ever accessing OST in the past 6 months’. Although we did not observe time-varying measures of the use of opioids, given the specificity with which we measured OST exposure, the high threshold for OST adherence and the blockade effect of methadone and buprenorphine against illicit opioid use, periods of OST exposure likely featured limited, to no illicit opioid use. Inclusion of this covariate would therefore constitute a violation of the assumption of positivity in this context .
Adjusting for the effect of prior HAART on OST resulted in a difference of 0.23 in the magnitude of the odds ratio, or 14% compared with unweighted (GEE) models. We believe this approximates the effect of healthcare system linkage through HAART on OST receipt. Although the remaining positive effect of OST on HAART cannot be separated into its hypothesized explanations (‘stabilized lifestyle’ and ‘linkage to the healthcare system’), we expect that healthcare system linkage through OST had a larger effect than the contrary. As opposed to antiretrovirals, certification is required to prescribe methadone or buprenorphine, making linkage to HIV care through OST prescribers perhaps more plausible than the contrary, as HAART prescribers with no OST exemption would have to refer clients elsewhere. The complex case-mix of this patient population, featuring a high prevalence of comorbid HIV, hepatitis C virus (HCV), mental health conditions and pain [31,32], is just one of the factors that challenge access and physician participation in the treatment of this patient population . Comprehensive and integrated healthcare, incorporating addiction treatment in primary care , is clearly critical for the successful treatment of HIV-positive individuals with opioid dependence. In this regard, integrated ‘directly observed therapy’ programmes, combining OST and HAART medication pickups, have shown mixed results, but continue to be assessed [35,36].
Our findings also have important implications for healthcare resource allocation decisions. The benefits in HAART medication refill adherence caused by exposure to OST lead directly to improved disease progression profiles , as well as reductions in infectivity through viral suppression [4–7]; these secondary benefits of OST should be incorporated directly into simulation models developed to determine the long-term cost-effectiveness of competing OST modalities  and population-level OST uptake [39,40].
We note that although maintenance-oriented treatment is recommended both federally and provincially in British Columbia, Canada, in some cases, patterns of medication dispensation resembling detoxification treatment have been observed . In the baseline model, we did not distinguish between these two modalities because there is no explicit indication within our health administrative databases to differentiate them. Nonetheless, we considered an alternate definition of reaching a minimum effective dose of 60 mg/day as an indirect measure of maintenance treatment receipt in sensitivity analysis. Furthermore, buprenorphine is not commonly used in BC, with the medication only gaining approval as a second-line treatment in 2010, initially only under special authority . As such, buprenorphine was only observed in only one patient month of follow-up within the study sample.
Our study offers the advantages of a complete population of HIV-positive PWID accessing OST and specific measurement of exposure and outcome, allowing for narrow repeated observation periods and a clear temporal ordering in cause and effect. Nonetheless, several limitations require careful consideration. First, although the study was population-based and featured complete follow-up data with precise measurement of exposure and outcome, the underlying policies and clinical practice in BC – specifically, maintenance-oriented OST with available office-based prescription and community pharmacy-based dispensation, and universal and comprehensive coverage for HIV medications and care – may differ substantially from other settings. Second, those identified as PWID at HIV diagnosis may have injected substances other than opioids exclusively, thus precluding indication for OST. Although we suspect the proportion of such individuals to be small, their inclusion may have attenuated the estimated effect of OST on HAART. Otherwise, although we judge the measurement of exposure and outcome to be at a high level, HAART adherence was measured according to pharmacy refill compliance records, which may overestimate true adherence, albeit equally during periods of OST exposure and nonexposure (nondifferential outcome misclassification), thus resulting in perhaps attenuated odds ratios. Finally, although we have tested the potential effects of model misspecification extensively, there remains potential unmeasured confounding due to depressive symptoms, alcohol use and potentially other drug use, previously identified as factors independently negatively associated with HAART adherence . The population under study likely included individuals with these behaviours; therefore, our results can be interpreted as a weighted average effect; samples with high levels of mental health comorbidity and/or polydrug use may have an attenuated effect, whereas nonexposed samples may have a larger effect. Otherwise, the causal interpretation of the results of marginal structural model relate to the affirmation of temporal ordering and control for key time-varying confounders; other interpretations of causal inferences consider a number of additional criteria, including repeated observation of the effect in different populations and study settings [41,42]. Cautious interpretations of the results are therefore recommended.
This study demonstrates the substantial benefits of OST in linking HIV-positive opioid-dependent individuals into HAART in a universal healthcare setting with freely available HAART. There is a priority to expand access to OST particularly within HIV-positive populations, to optimize HIV treatment uptake and adherence and its subsequent individual and population health and economic benefits.
We gratefully acknowledge the critical revisions suggested by Miguel Hernan on a previous draft of this article. We acknowledge BCMoH and Vancouver Coastal Health Decision Support Staff involved in data access and procurement, including Monika Lindegger, Clinical Prevention Services, BC Centre for Disease Control; Elsie Wong, Public Health Agency of Canada; Al Cassidy, BC Ministry of Health Registries and Joleen Wright and Karen Luers, Vancouver Coastal Health decision support. This study was funded by the National Institutes of Health/National Institute on Drug Abuse grant no. R01-DA032551 and the BC Ministry of Health-funded ‘Seek and treat for optimal prevention of HIV & AIDS’ pilot project.
B.N. and V.L. are Michael Smith Foundation for Health Research Scholars. J.S.G.M. has received grants from Abbott, Biolytical, Boehringer Ingelheim, Bristol-Myers Squibb, Gilead Sciences, Janssen, Merck and ViiV Healthcare. He is also supported by the Ministry of Health Services and the Ministry of Healthy Living and Sport, from the Province of British Columbia; through a Knowledge Translation Award from the Canadian Institutes of Health Research (CIHR); and through an Avant-Garde Award (No. 1DP1DA026182) from the National Institute of Drug Abuse, at the US National Institutes of Health. He has also received support from the International AIDS Society, United Nations AIDS Program, WHO, National Institute on Drug Abuse, National Institutes of Health Research-Office of AIDS Research, National Institute of Allergy & Infectious Diseases, The United States President's Emergency Plan for AIDS Relief (PEPfAR), Bill & Melinda Gates Foundation, French National Agency for Research on AIDS & Viral Hepatitis (ANRS), Public Health Agency of Canada.
The STOP HIV/AIDS Study Group is composed of the following: Rolando Barrios, MD, FRCPC, Senior Medical Director, VCH; Adjunct Professor, School of Population and Public Health, UBC; Patty Daly, MD, and Reka Gustafson, MD, Vancouver Coastal Health Authority; Mark Gilbert, Clinical Prevention Services, BC Centre for Disease Control; School of Population & Public Health, University of British Columbia; Perry R.W. Kendall, OBC, MBBS, MSc, FRCPC. Provincial Health Officer, British Columbia Ministry of Health; Clinical Professor, Faculty of Medicine UBC; Ciro Panessa, Gina McGowan and Nancy South, British Columbia Ministry of Health; Kate Heath, Julio S.G. Montaner, Robert S. Hogg and Bohdan Nosyk, BC Centre for Excellence in HIV/AIDS.
B.N. designed the study, assisted with analysis and wrote the first draft of the article. J.M. led the analysis and assisted with study design and interpretation of results. G.C. and B.Y. assisted with the analysis and study design. V.D.L., M.J.M., E.W. and J.S.G.M. aided in study design, interpretation of results and provided critical revisions to the article. J.S.G.M. led the procurement of the database. All authors approved the final draft.
The lead author (B.N.) affirms that the manuscript is an honest, accurate and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned have been explained.
Conflicts of interest
There are no conflicts of interest.
1. Mathers BM, Degenhardt L, Phillips B, Wiessing L, Hickman M, Strathdee SA, et al. 2007 Reference Group to the UN on HIV and Injecting Drug Use. Global epidemiology of injecting drug use and HIV among people who inject drugs: a systematic review
2. De Cock KM, Jaffe HW, Curran JW. The evolving epidemiology of HIV/AIDS
3. Mathers BM, Degenhardt L, Ali H, Wiessing L, Hickman M, Mattick RP, et al. 2009 Reference Group to the UN on HIV and Injecting Drug Use. HIV prevention, treatment, and care services for people who inject drugs: a systematic review of global, regional, and national coverage
4. Cohen MS, Chen YQ, McCauley M, Gamble T, Hosseinipour MC, Kumarasamy N, et al. Prevention of HIV-1 infection with early antiretroviral therapy
. N Engl J Med
5. Montaner JSG, Lima VD, Barrios R, Yip B, Wood E, Kerr T, et al. Association of highly active antiretroviral therapy coverage, population viral load, and yearly new HIV diagnoses in British Columbia, Canada: a population-based study
6. Nosyk B, Audoin B, Beyrer C, Cahn P, Granich R, Havlir D, et al. Examining the evidence on the causal effect of highly active antiretroviral therapy on transmission of human immunodeficiency virus using the Bradford Hill criteria
7. Wood E, Kerr T, Marshall BDL, Li K, Zhang R, Hogg RS, et al. Longitudinal community plasma HIV-1-RNA concentrations and incidence of HIV-1 among injecting drug users: prospective cohort study
8. Celentano DD, Galai N, Sethi AK, Shah NG, Strathdee SA, Vlahov D, Gallant JE. Time to initiating highly active antiretroviral therapy among HIV-infected injection drug users
9. Wood E, Hogg RS, Kerr T, Palepu A, Zhang R, Montaner JS. Impact of accessing methadone on the time to initiating HIV treatment among antiretroviral-naive HIV-infected injection drug users
10. Moatti JP, Carrieri MP, Spire B, Gastaut JA, Cassuto JP, Moreau J. Adherence to HAART in French HIV-infected injecting drug users: the contribution of buprenorphine drug maintenance treatment. The Manif 2000 study group
11. Palepu A, Tyndall MW, Joy R, Kerr T, Wood E, Press N, et al. Antiretroviral adherence and HIV treatment outcomes among HIV/HCV co-infected injection drug users: the role of methadone maintenance therapy
. Drug Alcohol Depend
12. Kapadia F, Vlahov D, Wu Y, Cohen MH, Greenblatt RM, Howard AA, et al. Impact of drug abuse treatment modalities on adherence to ART/HAART among a cohort of HIV seropositive women
. Am J Drug Alcohol Abuse
13. Roux P, Carrieri MP, Villes V, Dellamonica P, Poizot-Martin I, Ravaux I, Spire B. MANIF2000 cohort study group. The impact of methadone or buprenorphine treatment and ongoing injection on highly active antiretroviral therapy (HAART) adherence: evidence from the MANIF2000 cohort study
14. Roux P, Carrieri MP, Cohen J, Ravaux I, Poizot-Martin I, Dellamonica P, Spire B. Retention in opioid substitution treatment: a major predictor of long-term virological success for HIV-infected injection drug users receiving antiretroviral treatment
. Clin Infect Dis
15. Wood E, Hogg RS, Lima VD, Kerr T, Yip B, Marshall BD, Montaner JS. Highly active antiretroviral therapy and survival in HIV-infected injection drug users
16. Nosyk B, Geller J, Guh DP, Oviedo-Joekes E, Brissette S, Marsh DC, et al. The effect of motivational status on treatment outcome in the North American Opiate Medication Initiative (NAOMI) study
. Drug Alcohol Depend
17. Palepu A, Horton NJ, Tibbetts N, Meli S, Samet JH. Uptake and adherence to highly active antiretroviral therapy among HIV-infected people with alcohol and other substance use problems: the impact of substance abuse treatment
18. Nosyk B, Colley G, Yip B, Chan K, Heath K, Lima VD, et al. STOP HIV/AIDS Study Group. Application and validation of case-finding algorithms for identifying individuals with human immunodeficiency virus from administrative data in British Columbia, Canada
. PLoS One
19. Nosyk B, MacNab YC, Sun H, Fischer B, Marsh DC, Schechter MT, Anis AH. Proportional hazards frailty models for recurrent methadone maintenance treatment
. Am J Epidemiol
20. Carpenter CC, Fischl MA, Hammer SM, Hirsch MS, Jacobsen DM, Katzenstein DA, et al. Antiretroviral therapy for HIV infection in 1996. Recommendations of an international panel. International AIDS Society-USA
21. Thompson MA, Aberg JA, Hoy JF, Telenti A, Benson C, Cahn P, et al. Antiretroviral treatment of adult HIV infection: 2012 recommendations of the International Antiviral Society-USA panel
22. Hogg RS, Heath K, Lima VD, Nosyk B, Kanters S, Wood E, et al. Regional disparities in the burden of HIV/AIDS in Canada
. PLoS One
23. Schneeweiss S, Maclure M. Use of comorbidity scores for control of confounding in studies using administrative databases
. Int J Epidemiol
24. Hernan MA, Brumback BA, Robins JM. Estimating the causal effect of zidovudine on CD4 count with a marginal structural model for repeated measures
. Stat Med
25. Cole SR, Hernan MA. Constructing inverse probability weights for marginal structural models
. Am J Epidemiol
26. Robins JM, Hernan MA, Brumback B. Marginal structural models and causal inference in epidemiology
27. Harrell FE. SAS macros and data step programs useful in survival analysis and logistic regression. http://biostat.mc.vanderbilt.edu/wiki/pub/Main/SasMacros/survrisk.txt
[Accessed 15 December 2013].
28. Health Canada. Best practices, methadone maintenance treatment, Prepared by Jamieson, Beals, Lalonde and Associates Inc for the Office of Canada's Drug Strategy, Ottawa: Health Canada. 2002. http://guides.is.uwa.edu.au/content.php?pid=279406&sid=2301553
[Accessed 15 March 2014].
29. National Consensus Development on Effective Medical Treatment of Opiate Addiction. Effective medical treatment of opiate addiction
30. Nosyk B, Marsh DC, Sun H, Schechter MT, Anis AH. Trends in methadone maintenance treatment participation, retention, and compliance to dosing guidelines in British Columbia, Canada: 1996-2006
. J Subst Abuse Treat
31. Degenhardt L, Whiteford HA, Ferrari AJ, Baxter AJ, Charlson FJ, Hall WD, et al. Global burden of disease attributable to illicit drug use and dependence: findings from the Global Burden of Disease Study 2010
32. Dhingra L, Masson C, Perlman DC, Seewald RM, Katz J, McKnight C, et al. Epidemiology of pain among outpatients in methadone maintenance treatment programs
. Drug Alcohol Depend
33. Nosyk B, Anglin MD, Brissette S, Kerr T, Marsh D, Schackman BR, et al. A call for evidence-based medicine for the treatment of opioid dependence in North America
. Health Aff (Millwood)
34. Finkelstein R, Netherland J, Sylla L, Gourevitch MN, Cajina A, Cheever L. BHIVES Collaborative. Policy implications of integrating buprenorphine/naloxone treatment and HIV care
. J Acquir Immune Defic Syndr
2011; 56 (Suppl 1):S98–S104.
35. Hart JE, Jeon CY, Ivers LC, Behforouz HL, Caldas A, Drobac PC, Shin SS. Effect of directly observed therapy for highly active antiretroviral therapy on virologic, immunologic, and adherence outcomes: a meta-analysis and systematic review
. J Acquir Immune Defic Syndr
36. Lucas GM, Mullen BA, Galai N, Moore RD, Cook K, McCaul ME, et al. Directly administered antiretroviral therapy for HIV-infected individuals in opioid treatment programs: results from a randomized clinical trial
. PLoS One
37. Nosyk B, Min JE, Lima VD, Yip B, Hogg RS, Montaner JSG. on behalf of the STOP HIV/AIDS Study Group. Estimating HIV-1 disease progression to inform mathematical modeling and economic evaluation: an application using population-level data in British Columbia: 1996-2011
. J Acquir Immune Defic Syndr
38. Nosyk B, Guh DP, Bansback NJ, Oviedo-Joekes E, Brissette S, Marsh DC, et al. Cost-effectiveness of diacetylmorphine versus methadone for chronic opioid dependence refractory to treatment
39. Zarkin GA, Dunlap LJ, Hicks KA, Mamo D. Benefits and costs of methadone treatment: results from a lifetime simulation model
. Health Econ
40. Zaric GS, Barnett PG, Brandeau ML. HIV transmission and the cost-effectiveness of methadone maintenance
. Am J Public Health
41. Bradford-Hill A. The environment and disease: association or causation?
. Proc R Soc Med
42. Ward AC. The role of causal criteria in causal inferences: Bradford Hill's
. Epidemiol Perspect Innovat