This was a historical cohort study. A total of 1396 patients admitted between March 2011 and March 2012 in 16 MMT clinics (13 private and 3 state clinics) in Rafsanjan, entered the study and their retention rate was examined for one year.
The patients’ data abstracted from their medical records using checklists and collected by clinics’ staff. The checklists included demographic information, criminal history, pattern of drug use and the average dosage of daily methadone use during treatment.
Retention duration in MMT was the outcome variable. Retention was defined by the number of days a patient remained in MMT and calculated from the first MMT entry up to the date patients dropped out of treatment.
Predictor variables of retention were examined in two groups of individual and program related factors. Individual related factors included age, gender, marital status, education, employment, physical and mental illness history, addiction treatment history, polysubstance abuse, opiates use duration and prison history. Polysubstance abuse is defined as history of using two or more groups of addictive substances prior to treatment. Program related factors included type of clinic (private or state clinics), average dosage of methadone and experience of clinic providing methadone treatment. Clinic experience was defined as duration that clinic had provided methadone treatment.
3.1. Statistical Analysis
Results presented as mean ± SD (standard deviation) for numeric variables and summarized by absolute frequencies and percentages for categorical variables.
A multivariable Cox proportional hazards model with forward variable selection process for exploring variables associated with retention in methadone treatment was constructed. Variables incorporated into the multivariable model were age, gender, marital status, education, employment, physical and mental illness history, previous drug treatment, polysubstance abuse, duration of opioid dependence, prison history, type of clinic (private or state clinics), average dosage of methadone and experience of clinic providing methadone treatment.
The associations were presented as hazard ratios (HRs) with 95% confidence intervals (CIs). The overall fit of the final model was checked by Cox-Snell residuals and graphical inspections were also performed to assess any departures from proportionality. To further examine, we checked the proportional hazards assumption by adding a time-dependent covariate adjusted for other covariates (
32). No statistically significant violation of the assumption of proportionality was detected in multivariable analysis. Retention rate was also estimated using the Kaplan-Meier method, the standard estimator of the survival function, as univariable analyses.
For statistical analysis, statistical software SPSS version 15.0 for windows (SPSS Inc., Chicago, IL) and statistical package SAS version 9.1 for windows (SAS Institute Inc., Cary, NC, USA) were used. All p-values were 2-tailed, with statistical significance defined by P ≤ 0.05.