Data
The present study used data from multiple Danish national registers, all of which are continuously updated.
The Civil Registration System was established in 1968 and includes unique individual identification number, name, gender, date of birth, place of birth and residence, citizenship, identity of parents and spouses (23).
The Registry of Drug Abusers in Treatment has been recording information on people seeking treatment for DUD in publically funded treatment centers under the Danish social services since 1996 (24).
The National Patient Register was established in 1977 and contains personal and admission data for secondary care. The personal data include the unique identification number, municipality, and region of residence. The admission data include hospital and department codes, admission type, patient contact type (inpatient, outpatient, or emergency department), referral information, contact reason, and dates of admission and discharge (25).
The Psychiatric Central Research Register has recorded episodes of psychiatric care since 1970, and contains information for all outpatient, inpatient, and emergency contacts at psychiatric hospitals, including dates of beginning and end of treatment, diagnoses, type of referral, place of treatment, place of residence, and mode of admission (26).
The Danish Registry for Causes of Death contains information on dates and causes of death based up on the death certificate. Since 1875, the Danish National Board of Health has maintained the registers covering deaths among all Danish residents dying in Denmark, and since 1970 such records have been computerized (27).
The Central Criminal Register contains information about offenses and offenders in criminal cases for use in criminal procedures since November 1978. Information is transferred from the central crime register to Statistics Denmark (28).
All registers were linked, using the unique identification number assigned to each individual up on birth or first entry to Denmark as an immigrant.
The data for this study are stored on secure servers at Statistics Denmark, and all procedures were approved by the Danish Data Protection Agency. Since the data used for this study were collected and stored for monitoring and quality assurance, no ethics evaluation was needed under Danish law.
Inclusion criteria
Patients were included in the study if they had been enrolled in a publicly funded outpatient treatment facility for DUD in Denmark between 2000 and 2010, and were in the age range of 18 to 75 years at time of admission. Patients were excluded if their date of death was invalid. Less than five cases were omitted from the analyses, as their date of death was recorded as January 1st 1960, although they had been in treatment after January 1st 2000 (the exact number cannot be given due to data protection rules prohibiting the download of microdata). We followed the patients over the entire observation period, beginning from first registered treatment enrolment to completed suicide or December 31, 2010, whichever occurred first.
Measures
Outcome variable. The outcome in the study was defined as time from the first registered enrollment at a treatment center for DUD to completed suicide. Dates and causes of death were identified using the Danish Register of Cause of Death. We used ICD-10 codes to identify all completed suicides ascribed to intentional self-harm (X60-X84) or the sequelae of intentional self-harm (Y87.0) (29).
Predictors. The information on all substances used by the patients in the twelve months prior to enrolment in treatment were extracted from The Registry of Drug Abusers Undergoing Treatment. These variables were dummy-coded for the following types of substances: any opioids; central stimulants; cannabis; any recorded problem drinking in the database; use of benzodiazepines; methylenedioxymethamphetamine; and intravenous drug use. Further, we used a categorical predictor indicating previous drug treatment versus no previous drug treatment, or missing information on previous treatment based on the admission form. Using the code of reason for contact from the National Patient Register, we constructed a dummy variable representing any record of admission to a hospital in Denmark due to self-harm within the past 12 months leading up to the first registered admission to drug use disorder treatment. Using the Central Criminal Register, a dummy variable was created to indicate whether a person had been charged with a crime within the past 12 months, leading up to the first registered treatment admission. The sociodemographic variables include gender, age, civil status (living without a partner or not), not being in education, employment or training, and immigrant status (born in Denmark or not).
Analyses
Descriptive statistics are reported as percentages for dichotomous variables and means with standard deviations for all other variables. Comparison between groups was done using Nelson-Aalen curves of the cumulative hazard and estimated cumulative incidence functions. Time-to-event analysis for completed suicide was conducted using Fine and Gray’s competing risks analysis (30), in which the cumulative incidence function (CIF); i.e., Ce (t) gives the proportion of patients at time t who have experienced event e, while accounting for the fact that patients can experience another event that prevents event e from happening, labeled the competing event or competing risk (e.g., death not caused to suicide will rule out later death caused by suicide).
In the analyses, subjects were considered to have experienced the event if they died due to suicide, to have experienced the competing event, if they died from any other cause, and to be censored if no event took place by December 31st, 2010.
All p-values were 2-tailed, and level of significance was assessed as a Type I error with rate of alpha 0.05. All statistical analyses were performed using Stata 15 (31).
In order to compare the rate of suicide among people treated for DUD with the rate in the total population, we drew matched controls from a representative sample of the national population without a history of public funded treatment for DUD or alcohol use disorders. The sample size of the control group is five randomly drawn individuals from the total population (n=139,710), proportionate to each individual enrolled at a DUD treatment center within the time period of this study, using similar gender and age group at enrolment date at the treatment centers. The age groups were 18-27, 28-37, 38-47, and 48-75 years at the time of admission to treatment.
A random enrolment date was generated for the control group. The only restriction made in the matching was that the individuals in the control group did not die before the enrolment date, the year that the follow-up started, and did not have a record of treatment for a drug or alcohol use disorder in the databases. In the control group, the time-to-event outcome variable was created after generating a random enrolment date with a beta distribution.
We analyzed the standardized mortality ratios (SMRs) to estimate the suicide gap between individuals with DUD, recent psychiatric history, or both, compared with a general population sample with neither. Comparisons were made using three categories: individuals in the control group who had past year psychiatric care history, patients treated for DUD with past year psychiatric care, and patients treated for DUD without past year psychiatric care. We used individuals in the control group without past year psychiatric care as the reference group.