Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Pediatrics. Author manuscript; available in PMC 2011 May 1.
Published in final edited form as:
PMCID: PMC2884182

Comparative Safety of Antidepressant Agents for Children and Adolescents Regarding Suicidal Acts

Sebastian Schneeweiss, MD, ScD,a Amanda R. Patrick, MS,a Daniel H. Solomon, MD, MPH,a Colin R. Dormuth, MA, MS, ScD,b Matt Miller, MD, ScD,c Jyotsna Mehta, MS,a Jennifer C. Lee, BS,a and Philip S. Wang, MD, DrPHa,d



The objective of this study was to assess the risk of suicide attempts and suicides after initiation of antidepressant medication use by children and adolescents, for individual agents.


We conducted a 9-year cohort study by using population-wide data from British Columbia. We identified new users of antidepressants who were 10 to 18 years of age with a recorded diagnosis of depression. Study outcomes were hospitalization attributable to intentional self-harm and suicide death.


Of 20 906 children who initiated antidepressant therapy, 16 774 (80%) had no previous antidepressant use. During the first year of use, we observed 266 attempted and 3 completed suicides, which yielded an event rate of 27.04 suicidal acts per 1000 person-years (95% confidence interval [CI]: 23.9–30.5 suicidal acts per 1000 person-years). There were no meaningful differences in the rate ratios (RRs) comparing fluoxetine with citalopram (RR: 0.97 [95% CI: 0.54–1.76]), fluvoxamine (RR: 1.05 [95% CI: 0.46–2.43]), paroxetine (RR: 0.80 [95% CI: 0.47–1.37]), and sertraline (RR: 1.02 [95% CI: 0.56–1.84]). Tricyclic agents showed risks similar to those of selective serotonin reuptake inhibitors (RR: 0.92 [95% CI: 0.43–2.00]).


Our finding of equal event rates among antidepressant agents supports the decision of the Food and Drug Administration to include all antidepressants in the black box warning regarding potentially increased suicidality risk for children and adolescents beginning use of antidepressants.

Keywords: antidepressant medications, fluoxetine, suicide, intentional self-harm, cohort study, epidemiology

WHAT'S KNOWN ON THIS SUBJECT: Antidepressants may be associated with increased risk of suicidality for children, but variations in risk among individual antidepressant agents are unknown.

WHAT THIS STUDY ADDS: This study found no significant variations in rates of suicidal acts among antidepressant agents, which supports the decision of the Food and Drug Administration to require a black box warning regarding potentially increased suicidality risk for all antidepressants.

In October 2004, the US Food and Drug Administration issued an advisory that antidepressants may be associated with increased risks of suicidal thoughts and behaviors among children and adolescents.1 These warnings were prompted by a meta-analysis of all available randomized trials of antidepressants in this age group, in which patients assigned randomly to receive antidepressants had nearly twice the rate of suicidal ideation or behavior, compared with those who received placebo.2 However, the interpretation of the findings of that meta-analysis is limited by a number of factors, including the brief duration of trials, few suicide attempts and almost no completed suicides, varying definitions of suicidal ideation, non-comparable doses, and heterogeneous patient populations.

Epidemiological studies of the safety of antidepressants in pediatric populations are very limited. Studies comparing antidepressant users with nonusers have generally suggested an increased risk of suicidal acts associated with antidepressant use,36 but these results could be caused by the channeling of antidepressant agents to children and adolescents who are more likely to attempt suicide, leading to intractable confounding by indication.7 Only 2 studies compared suicide risk across antidepressant agents,8,9 and questions remain regarding whether those studies had adequate statistical power or adequately controlled for prescribing biases caused by preferential avoidance of the use of tricyclic agents for patients at high risk for suicidality.

Although they raised important concerns, Food and Drug Administration meta-analyses and advisories have not provided patients, clinicians, or policymakers with adequate guidance for the treatment decisions they face. In the current study, we sought to address whether the putative link between antidepressants and suicidality applies equally to all antidepressant classes, agents, and durations of use or whether there are particular regimens with safety advantages that should be prescribed preferentially in pediatric and adolescent populations.


Patients and Data Source

We conducted a cohort study of all British Columbia residents 10 to 18 years of age who initiated use of an antidepressant medication between January 1, 1997, and December 31, 2005. Children <10 years of age were excluded on the basis of the assumption that they would not attempt suicide.10,11 Initiation was defined as filling an antidepressant prescription without having filled one in the preceding year. We analyzed only the first treatment episode for each patient during the study period. We required evidence of depression, as indicated by a diagnosis of depression recorded during 2 office visits or as a hospital discharge diagnosis during the 6 months before through 2 months after the initiation date. We allowed subjects with diagnoses after the initiation date to be included in the cohort because some practitioners may prescribe an antidepressant as part of diagnosing the condition, in which case the diagnosis may follow the prescription. To ensure complete ascertainment of previous drug use, we required that subjects be residents of British Columbia, as evidenced by enrollment in the provincial Medical Services Plan during the year before initiation. Bupropion was not considered because of its potential use for smoking cessation, and escitalopram was not considered because it was not marketed until December 2004, the end of the study period.

The PharmaNet database (British Columbia Ministry of Health) used to identify antidepressant initiators includes the name, dose, and dispensed quantity for all prescription drugs dispensed in British Columbia pharmacies. This information is entered by pharmacists through a province-wide network that ensures minimal under-reporting and misclassification, and it is recorded for all dispensings, independent of the payor. The British Columbia Ministry of Health also maintains linkable data on all physician services and hospitalizations for all persons in its publicly funded health care system. Up to 25 diagnoses for hospital discharges and 1 diagnosis for each medical service are recorded, with good specificity and completeness.12 These administrative data were further linked to vital statistics data that included information on cause of death, including suicide.

Antidepressant Medication Exposure

Antidepressant medications were grouped into the following categories: selective serotonin reuptake inhibitors (SSRIs), including citalopram, fluoxetine, fluvoxamine, paroxetine, and sertraline; serotonin-norepinephrine reuptake inhibitors (SNRIs), including venlafaxine; tricyclic antidepressants, including amitriptyline, amoxapine, clomipramine, desipramine, doxepin, imipramine, maprotiline, nortriptyline, protriptyline, and trimipramine; other newer and atypical agents, including mirtazapine, nefazodone, and trazodone; and monoamine oxidase inhibitors (MAOIs), including moclobemide, phenelzine, and tranylcypromine. Duloxetine, a SNRI, was not marketed in Canada during the study period.

Exposure status was assigned on the basis of the initiated medication. We determined the availability of antidepressant supply for each patient day by combining consecutive antidepressant dispensings for each patient on the basis of dispensing dates and reported “days supply.”13 When a dispensing occurred before the previous dispensing should have run out, use of the new dispensing was assumed to have begun the day after the end of the old dispensing. Each patient's exposure risk window ended when the subject had been without antidepressant supply for 14 days. Patients underwent censoring at the end of their exposure risk window, the date of switching to or augmenting treatment with another antidepressant medication, emigration from the province, occurrence of a study outcome, death, the end of 1 year, or the end of the study period, whichever came first.

Study End Points

Study outcomes were attempted suicide, completed suicide, and the composite outcome of suicidal acts. Suicide attempt was defined as a hospitalization with an International Classification of Diseases, Ninth Revision (ICD-9), external-cause-of-injury code (E-code) for deliberate self-harm (E950.x–E 958.x). Completed suicides were deaths with an ICD-9 E-code of E950.x to E958.x or an International Classification of Diseases, 10th Revision (ICD-10), diagnosis code of X60 to X84 listed as the cause of death. Violent suicide attempts (ICD-9 E-codes) and completed suicides (ICD-10 codes) were defined as those involving hanging (E953.0; X70), gunshot/explosion (E955; X72–X75), jumping/lying in front of a moving object (E957, E958.0; X80, X81), stabbing/blunt trauma (E956; X78, X79), vehicle collision (E958.5, E958.6; X82), electrocution (E958.4; X83), or self-immolation (E958.1; X76). Assignment and recording of an E-code are mandatory in Canada whenever a diagnosis of poisoning or injury is coded.14 A study of the validity of deliberate self-harm codes for identifying suicide attempts reported a positive predictive value of 86% for deliberate self-harm E-codes.15

Patient Characteristics

Patient characteristics were assessed at treatment initiation, on the basis of medical claims during the year preceding cohort entry. These characteristics included age, gender, calendar year of cohort entry, adjusted family income status (below C$16 000, $16 000–28 000, or above C$28 000, as defined by premium subsidy levels),16 enrollment in an income assistance plan, number of acute hospitalizations, Charlson comorbidity score,17 and number of distinct generic entities prescribed. Psychiatric disorders were defined as the presence of 2 recorded outpatient diagnoses or 1 inpatient diagnosis. Diagnoses included anxiety and obsessive-compulsive disorders, sleep disorders, mania, attention-deficit/hyperactivity disorder (ADHD), substance abuse, psychotic disorders, delirium, personality disorders, and other mental disorders (hysteria, acute reactions to stress, sexual deviations, and disturbances of emotions specific to childhood or adolescence). Additional indicators of psychiatric severity were the number of psychiatric hospitalizations, number of psychiatric visits with a depression diagnosis, suicide attempt in the past year, concurrent stimulant use, and number of psychiatric drug classes (antipsychotic agents, benzodiazepines, and other psychotropic drugs) taken simultaneously. Concurrent stimulant use was defined as having stimulant supply available on the antidepressant initiation date. In addition to psychiatric comorbidities, we measured a number of general medical comorbidities, including malignancy, pain requiring opiates, poisoning or drug toxicity, injury other than poisoning, and seizure disorders. A complete list of covariates and their definitions is provided in Appendix 1.

Statistical Analyses

We used Cox proportional-hazards regression analysis to estimate the effect of antidepressant use on the composite outcome of suicidal acts (attempted or completed suicide). Fluoxetine was used as the reference group for comparisons among the SSRIs, because this drug is widely used and has demonstrated efficacy in children.14,15 For between-class comparisons, SSRIs served as the reference group. To reduce the potential for confounding by indication, a bias that arises when prognostic factors, such as a patient's perceived suicide risk, influence treatment decisions, we adjusted for a number of patient demographic and clinical characteristics (Appendix 1).7 Methods of adjustment included traditional, multivariate, adjusted outcome models, models adjusted for a propensity score including the prespecified covariates in Appendix 1, and models adjusted for a high-dimensional propensity score including empirically identified covariates. These methods are described in Appendix 2.

We plotted Kaplan-Meier curves for suicide attempt–free survival as a function of the duration of continuous antidepressant use. We also created multivariate, adjusted, Kaplan-Meier plots by weighting data for each subject by the inverse of his or her probability for treatment, as estimated in the propensity score analysis described above.18


From 1997 to 2005, a total of 20 906 children and adolescents initiated antidepressant therapy. As shown in Table 1, girls accounted for 63% of the population. ADHD (5.1%), manic-depressive illness (2.1%), substance abuse (1.7%), and other mental disorders (8.1%) were the most commonly recorded psychiatric comorbidities; 7.1% of the population had been hospitalized with a psychiatric condition, 2.0% had attempted suicide, and 2.6% were using stimulants at the time of initiation of antidepressant use. SSRIs were the most common class of antidepressants prescribed, accounting for 78% of antidepressant use. Paroxetine, citalopram, and sertraline, the most commonly prescribed drugs, accounted for 25%, 17%, and 17% of overall use, respectively. A total of 16 774 subjects (80%) had no antidepressant use in the past 3 years and contributed to the treatment-naive main analysis. Children and adolescents taking fluvoxamine (5.1%), MAOIs (0.2%), and tricyclic drugs (7.2%) were more likely to be male, to have a diagnosis of ADHD, and to be taking stimulants concurrently, compared with children and adolescents given other agents. There were no other meaningful differences between treatment groups.

Distribution of Patient Baseline Characteristics

During the first year of use, we identified a total of 268 children and adolescents who attempted (n = 266) or completed (n = 3) suicide, which yielded a rate of 27.0 suicidal acts (attempted and completed suicides) per 1000 person-years (95% confidence interval [CI]: 23.9–30.5 suicidal acts per 1000 person-years) (Table 2). The suicidal act rates for subgroups of treatment-naive subjects and subjects without a previous suicide attempt were similar to the overall rate (Table 2). Although the hazards were proportional throughout the 1-year follow-up period, most of the events occurred in the first 6 months after initiation (Figs 1 and and2).2). We identified a total of 43 violent suicide attempts or suicides during the follow-up period, which resulted in a rate of 4.3 cases per 1000 person-years (95% CI: 3.1–5.8 cases per 1000 person-years).

Suicide and suicide attempt-free probability according to SSRI agent and time since treatment initiation. Because of the relatively low incidence of suicidal acts, the y-axis scale was inflated to allow for comparisons between agents.
Suicide and suicide attempt-free probability according to antidepressant class and time since treatment initiation. AD indicates antidepressant. Because of the relatively low incidence of suicidal acts, the y-axis scale was inflated to allow for comparisons ...
Suicide and Attempted Suicide Counts and Rates During 1-Year Follow-up Period

In our primary analysis of treatment-naive children and adolescents, there were no clinically meaningful differences between fluoxetine and the other SSRIs (Table 3). In unadjusted analyses, tricyclic agents were associated with a reduced risk of attempted or completed suicide, compared with SSRIs (rate ratio [RR]: 0.59 [95% CI: 0.28–1.27]). This effect was attenuated to 0.71 (95% CI: 0.33–1.52) in a propensity score-adjusted analysis and was further reduced to 0.92 (95% CI: 0.43–2.00) after adjustment for the high-dimensional propensity score. MAOIs were not investigated because of the small numbers of exposed subjects. Truncation of follow-up time at 6 months after initiation of antidepressant use yielded quantitatively similar results (Appendix 3). Citalopram seemed to be associated with a lower risk of violent suicides and attempts (RR: 0.23 [95% CI: 0.04–1.29]), but the estimate was unstable because of few violent events.

Event RRs for Suicidal Acts and Violent Suicidal Acts During 1-Year Follow-up Period

We did not find evidence of multiplicative effect-measure modification by concurrent stimulant use. There was some evidence of sertraline effect-measure modification by ADHD status. The presence of ADHD among sertraline users increased the rate of suicidal acts substantially, relative to the rate observed for fluoxetine users. Sertraline was associated with a 22% reduced risk (RR: 0.68 [95% CI: 0.39–1.18]), compared with fluoxetine, among subjects without ADHD, whereas it was associated with a large increase in risk (RR: 4.70 [95% CI: 0.56–39.25]) among subjects with ADHD, although both estimates lacked precision. Similarly, we found some indication that the initiation of tricyclic antidepressant use increased the risk of suicide and suicide attempt more among patients with previous psychiatric hospitalization (RR: 2.80 [95% CI: 0.84–9.27]) than among those with no such hospitalization (RR: 0.45 [95% CI: 0.17–1.22]).


In a population of 20 906 children and adolescents 10 to 18 years of age who were beginning antidepressant therapy in British Columbia between 1997 and 2005, we observed no clinically meaningful variation in the risk of suicidal acts according to antidepressant agent within the class of SSRIs or between antidepressant classes. We observed some numeric variations in the effect sizes between individual agents, but the 95% CIs were substantially overlapping. Because of the small number of patients receiving MAOIs, we did not attempt to compare the risk of suicidal acts between MAOIs and SSRIs. Sertraline might have demonstrated an increased risk of suicidality among children with a recorded diagnosis of ADHD, although this finding was very unstable because of the few events. The majority of events occurred in the first half-year after treatment initiation. Our results are consistent with the findings of several other observational studies that reported small or no differences in rates of suicides and suicide attempts between antidepressant classes,8,9 as well as meta-analyses of randomized, controlled trial data that found no heterogeneity among SSRIs.1921

Studies that do not find statistically significant differences merit special attention to their statistical power to detect a clinically meaningful difference.22 Because of sample size, we had the greatest power to detect a difference in suicidal act rates between venlafaxine and the SSRI drug class. We could exclude a hazard reduction of >11% and an increase of >108% with a confidence level of 95% for venlafaxine users versus SSRI users. Our estimates of the comparative safety of fluvoxamine versus fluoxetine were the least precise. For the fluvoxamine versus fluoxetine comparison, we could exclude a hazard reduction of >54% and an increase of >143%.

During the 9-year study period, the suicide risk among all British Columbia residents 13 to 17 years of age, as reported by the British Columbia Coroner's Office and British Columbia Statistics, averaged ~0.052 suicide deaths per 1000 people.2325 The rate we observed after initiation of antidepressant use was 5 times higher, which likely reflects the current depressed state and greater degree of psychiatric comorbidity in our population.

A major strength of our study is that the large, stable, study population allowed us to examine a variety of medications and important subgroups. We had an adequate sample size to restrict our study to new initiators of antidepressants with documented diagnoses of depression. An incident user design reduces the likelihood of missing early adverse events, allows for an evaluation of risks over time, ensures that the assessment of patient baseline characteristics is not influenced by effects of antidepressant treatment, and reduces the likelihood that current treatment assignment is influenced by past drug-related experiences, such as adverse affects and refractory symptoms. Our decision to censor data for subjects at treatment discontinuation and to use a proportional-hazards analysis inherently controls for differences in treatment persistence. To obtain clearly identified exposure groups, we compared monotherapies and censored patient follow-up data as soon as the patient switched drugs or augmented therapy. This analytic strategy makes treatment groups more comparable with respect to initial health state and avoids analytic difficulties associated with comparisons of patients who escalate or change treatment in response to treatment failure or adverse effects with patients who do not. This reduces the generalizability of our findings to patients receiving monotherapy; however, this reduction in generalizability is outweighed by the improvement in validity. The extent of generalizability would depend on the extent to which prescribing by British Columbia physicians and underlying suicide risk and comorbidity patterns among British Columbia children and adolescents are similar to those in other regions and health care settings.

Nonrandomized studies using health care utilization data are particularly scrutinized for their control of confounding and the potential for misclassifying diagnoses.26 Confounding would occur if certain antidepressants were more likely to be given to patients with a greater background risk of suicide. Therefore, we controlled for sociodemographic, clinical, and health care utilization factors likely to be independent predictors of suicidality, by using traditional multivariate and high-dimensional propensity score techniques. However, our ability to adjust fully for mental health status was limited by the measurement and reporting of mental health conditions as ICD-9 codes on insurance claims to the provincial government. Random misclassification of confounders in health care utilization databases leads to incomplete adjustment of confounding bias.27 For example, our finding that expanding the population to include non–treatment-naive subjects increased the apparent protective effect of tricyclic antidepressants, compared with SSRIs, suggests that there may be residual confounding by the presence of treatment-resistant depression or other suicide risk factors in this analysis.

In terms of the validity of our suicide attempt definition, intentional self-harm E-codes do not distinguish between suicidality and self-harm without suicidal intent, despite their widespread use in research and surveillance to identify suicide attempts. However, a study of the validity of deliberate self-harm E-codes to identify suicide attempts reported a positive predictive value of 86% for these codes, relative to the standard method of medical chart review.15 To the extent that suicide attempts are assigned the correct E-codes, these admissions should be captured in the British Columbia hospital discharge database, because this database had an E-coding completeness rate of >95% for the years evaluated. Although it is possible that some suicide attempts were recorded as unintentional injuries or injuries of unknown intent, a Canadian study found that E-codes and medical chart review yielded consistent estimates of the proportions of poisoning-related hospitalizations among subjects 18 to 24 years of age that were attributable to intentional self-harm, and a US study reported 95% agreement between E-codes and physician opinions regarding the intent of injuries that resulted in hospitalization. We did not have data on suicide attempts that were treated in an emergency department without inpatient admission. However, the majority of pediatric and adolescent patients with nonfatal self-harm seen in the emergency department are admitted to the hospital,11 and our inclusion of only more-severe attempts that required hospitalization is likely to increase the specificity of our definition, at some cost in sensitivity. High specificity is desirable, in that relative risk estimates are unbiased if outcomes are assessed with 100% specificity, even if sensitivity is far lower.28


Our analysis supports the decision of the Food and Drug Administration to include all antidepressants in the black box warning regarding increased suicidality risk for children and adolescents initiating use of antidepressants. Once a decision to initiate pharmacotherapy is made, treatment decisions should be made on the basis of efficacy and less so on the basis of safety. Clinicians should be vigilant in monitoring children and adolescents for whom use of any antidepressant agents is initiated.


This project was funded by a grant from the National Institute of Mental Health (grant RO1-MH078708).

Funded by the National Institutes of Health (NIH).


attention-deficit/hyperactivity disorder
confidence interval
external-cause-of-injury code
International Classification of Diseases, Ninth Revision
International Classification of Diseases, 10th Revision
monoamine oxidase inhibitor
rate ratio
selective serotonin reuptake inhibitor
serotonin-norepinephrine reuptake inhibitor


Baseline Covariates and Definitions

Nonpsychiatric DisordersaDefinition
Chronic lung disease491.xx
Diabetes mellitus250.xx
Plus ≥1 prescription for any of the following: insulin, acarbose, miglitol, metformin, acetahexamide, chlorpropamide, glimepiride, glipizide, glyburide, repaglinide, tolazamide, tolbutamide, pioglitazone, rosiglitazone, troglitazone, or nateglinide
Gastrointestinal hemorrhage578.x
531.0, 531.2, 531.4, 531.6
532.0, 532.2, 532.4, 532.6
533.0, 533.2, 533.4, 533.6
534.0, 534.2, 534.4, 534.6
Ventricular arrhythmia427.4x, 427.41, 427.42, 427.5x, 427.1x, 427.69
Other arrhythmia427.xx or 426.xx except for 427.4x, 427.41, 427.42, 427.5x, 427.1x, 427.69
Congenital heart diseases745.xx, 746.xx
Other congenital anomalies of cardiovascular system747.xx
Rheumatic heart diseases393.xx–398.xx except for 398.91
Hypertensive heart and kidney diseases402.xx–404.xx except for 402.01, 402.11, 402.91, 404.01, 404.11, 404.91, 404.03, 404.13, 404.93
Nonrheumatic valve424.xx
Endocarditis, pericarditis, and myocarditis420.xx, 421.xx, 422.xx
Rheumatic fever390.xx, 391.xx, 392.xx
Essential hypertension401.xx or use of thiazide
Secondary hypertension405.xx
Hypothyroidism243, 244.x
Injury other than poisoning800.xx–959.xx
Colorectal cancer153.xx, 230.3x, 154.xx(except 154.2, 154.3 and 154.4), 230.4
Lung cancer162.xx, 231.2
Breast cancer174.xx, 233.0
Prostate cancer185.xx, 233.4
Other malignancy140.xx–208.xx (hospital diagnosis), excluding lung cancer, breast cancer, prostate cancer, and colorectal cancer codes
Other heart disease411.xx–414.xx
Pain requiring high-potency opiatesUse of diacetylmorphine HCl, fentanyl, hydromorphone, hydroxyephedrine/normethadone, meperidine, methadone, or morphine
Pain requiring middle-potency opiatesUse of codeine, hydrocodone, oxycodone, pseudoephedrine, drocodone, or propoxyphene
Peptic ulcer disease and related conditions (without hemorrhage)531.1, 531.3, 531.5, 531.7, 531.9
532.1, 532.3, 532.5, 532.7, 532.9
533.1, 533.3, 533.5, 533.7, 533.9
534.1, 534.3, 534.5, 534.7, 534.9
Pneumonia480.x, 481, 482.xx, 483.x, 484.x, 485, 486
Poisoning or drug toxicity960.xx–979.xx
Rheumatoid arthritis714.x plus disease-modifying antirheumatic drug or steroid use
Seizure disorders345, 345.xx
Psychiatric disordersb
    Depressive illnesses293.83
Anxiety disorders300.0x
Sleep disorders307.4x
347, 347.xx
Manic-depressive illness296.0x
ADHD312.xx, 314.xx
Alcohol abuse291.0x
303, 303.xx
Drug abuse or dependence292.0x
304, 304.xx
Psychotic disorders290.8x
295, 295.xx
297, 297.xx
298, 298.1x–298.9x
299, 299.xx
Other psychiatric disorders
Personality disorders301, 301.xx
Other mental disorders300.1x
302, 302.xx
306, 306.xx
307.0x–307.3x, 307.5x–307.9x
308, 308.xx
310, 310.xx
313, 313.xx
315, 315.xx
316, 316.xx
aDefined as having ≥1 hospitalization with the condition as any diagnosis or ≥1 physician claim with the condition recorded unless otherwise noted.
bDefined as having ≥1 hospitalization with the condition as any diagnosis or ≥2 physician claims with the condition recorded.


To reduce the potential for confounding by indication, a bias that arises when prognostic factors, such as a patient's perceived suicide risk, influence treatment decisions, we adjusted for a number of patient demographic and clinical characteristics (Appendix 1).7 Because of the anticipated small number of outcomes, we reduced these potential confounders into an exposure propensity score.29 A propensity score is a patient's probability of receiving treatment (in this case, another SSRI versus fluoxetine or a non-SSRI versus a SSRI), adjusted for his or her measured characteristics, and typically is estimated by using a logistic regression model.30 Treatment effects can be estimated within propensity score levels, which define groups of patients who are similar in their probabilities of having received treatment. We estimated propensity scores by using logistic regression models including the patient covariates in Appendix 1; no further variable selection was performed. After plotting and comparing the distribution of propensity scores for the exposed and unexposed patients, we truncated our study population to the area of overlap. Because relative risk estimates varied within deciles of propensity scores, we used exposure propensity score decile as a categorical variable in the Cox proportional-hazards models.

In addition to constructing a propensity score that included the covariates described above, in a secondary analysis we constructed a high-dimensional propensity score that included empirically identified covariates. This technique examined all dispensed drugs, recorded diagnoses, and performed procedures reported for a patient before the initiation of antidepressant therapy and identified codes that might potentially bias the exposure and outcome in question. All identified codes were sorted according to their potential for confounding and were used to create a propensity score that was based on the top 500 identified variables in addition to the variables identified by the investigators (Appendix 1). This new technique has been shown to improve adjustment for confounding.31

We also conducted a multivariate adjusted analysis. This model was constructed by using a backward selection algorithm with the P value for covariate exclusion criteria set to .2, as recommended by simulation studies on variable selection in epidemiology.32 To decrease further the likelihood of confounding by indication,7 we performed the primary analysis among treatment-naive subjects initiating antidepressant use, defined as children and adolescents with no antidepressant use in the past 3 years. This design reduces the likelihood that we analyzed data for children and adolescents initiating use of a medication for recurrent depression. We conducted a subgroup analysis with children and adolescents without a previous suicide attempt. In addition, we tested for effect-measure modification by various covariates, such as concurrent use of stimulants, previous ADHD diagnosis, previous psychiatric hospitalization, and previous suicide attempt.


Event RRs for Suicidal Acts and Violent Suicidal Acts During 6-Month Follow-up Period

RR (95% CI)
Suicidal Acts
Violent Suicidal Acts
UnadjustedAdjusted for Age, Gender, and Calendar YearAdjusted for Propensity Score DecileAdjusted for Propensity Score Decile
Children and adolescents with no antidepressant use in past 3 y
        Citalopram1.02 (0.62–1.70)1.00 (0.58–1.74)1.02 (0.58–1.80)0.23 (0.04–1.29)
        Fluvoxamine1.36 (0.65–2.84)1.29 (0.60–2.79)1.16 (0.54–2.51)0.00 (0.00–0.00)
        Paroxetine1.07 (0.66–1.76)0.88 (0.52–1.48)0.86 (0.51–1.45)1.41 (0.48–4.13)
        Sertraline1.14 (0.67–1.93)1.01 (0.58–1.74)0.95 (0.55–1.65)1.10 (0.34–3.53)
        SNRIs1.07 (0.71–1.62)1.12 (0.74–1.70)1.09 (0.72–1.67)0.29 (0.04–2.12)
        Other newer and atypical agents0.37 (0.12–1.17)0.39 (0.12–1.21)0.36 (0.11–1.14)1.37 (0.32–5.87)
        Tricyclic drugs0.63 (0.29–1.34)0.70 (0.33–1.51)0.75 (0.35–1.62)0.56 (0.07–4.17)
All children and adolescents
        Citalopram0.87 (0.56–1.37)0.90 (0.53–1.53)0.91 (0.53–1.58)0.23 (0.04–1.30)
        Fluvoxamine1.38 (0.79–2.43)1.34 (0.75–2.40)1.15 (0.64–2.07)1.12 (0.21–5.95)
        Paroxetine1.06 (0.71–1.58)0.97 (0.64–1.47)0.96 (0.63–1.47)2.14 (0.78–5.86)
        Sertraline1.02 (0.66–1.58)0.94 (0.60–1.48)0.88 (0.56–1.38)1.62 (0.56–4.66)
        SNRIs0.98 (0.66–1.45)1.07 (0.72–1.60)1.05 (0.70–1.57)0.50 (0.12–2.10)
        Other newer and atypical agents0.58 (0.26–1.30)0.61 (0.27–1.37)0.56 (0.25–1.27)1.61 (0.49–5.30)
        Tricyclic drugs0.40 (0.19–0.84)0.43 (0.20–0.92)0.48 (0.22–1.02)0.35 (0.05–2.60)
Children and adolescents with no previous suicide attempt
        Citalopram0.85 (0.54–1.34)0.88 (0.52–1.51)0.89 (0.51–1.55)0.22 (0.04–1.22)
        Fluvoxamine1.40 (0.80–2.46)1.37 (0.77–2.44)1.12 (0.62–2.02)0.97 (0.18–5.30)
        Paroxetine1.01 (0.68–1.52)0.92 (0.60–1.40)0.91 (0.60–1.40)1.99 (0.72–5.50)
        Sertraline1.00 (0.65–1.56)0.92 (0.59–1.45)0.86 (0.55–1.36)1.56 (0.54–4.49)
        SNRIs0.97 (0.65–1.45)1.06 (0.70–1.59)1.04 (0.69–1.57)0.53 (0.13–2.21)
        Other newer and atypical agents0.59 (0.26–1.33)0.62 (0.27–1.40)0.58 (0.26–1.30)1.57 (0.48–5.18)
        Tricyclic drugs0.40 (0.19–0.85)0.45 (0.21–0.95)0.49 (0.23–1.04)0.35 (0.05–2.57)

The reference group for the SSRI comparison was fluoxetine; the reference group for the antidepressant drug class comparisons was SSRIs.


FINANCIAL DISCLOSURE: The authors have indicated they have no financial relationships relevant to this article to disclose.


1. Food and Drug Administration Summary minutes of the September 13–14, 2004 Center for Drug Evaluation and Research Pharmachopharmacologic Drugs Advisory Committee and the FDA Pediatric Advisory Committee. [February 1, 2005]. Available at:
2. Hammad TA, Laughren T, Racoosin J. Suicidality in pediatric patients treated with anti-depressant drugs. Arch Gen Psychiatry. 2006;63(3):332–339. [PubMed]
3. Tiihonen J, Lönnqvist J, Wahlbeck K, Klaukka T, Tanskanen A, Haukka J. Antidepressants and the risk of suicide, attempted suicide, and overall mortality in a nationwide cohort. Arch Gen Psychiatry. 2006;63(12):1358–1367. [PubMed]
4. Didham RC, McConnell DW, Blair HJ, Reith DM. Suicide and self-harm following prescription of SSRIs and other antidepressants: confounding by indication. Br J Clin Pharmacol. 2005;60(5):519–525. [PubMed]
5. Olfson M, Marcus SC, Shaffer D. Antidepressant drug therapy and suicide in severely depressed children and adults: a case-control study. Arch Gen Psychiatry. 2006;63(8):865–872. [PubMed]
6. Valuck RJ, Libby AM, Sills MR, Giese AA, Allen RR. Antidepressant treatment and risk of suicide attempt by adolescents with major depressive disorder: a propensity-adjusted retrospective cohort study. CNS Drugs. 2004;18(15):1119–1132. [PubMed]
7. Walker AM. Confounding by indication. Epidemiology. 1996;7(4):335–336. [PubMed]
8. Jick H, Kaye J, Jick S. Antidepressants and the risk of suicidal behaviors. JAMA. 2004;292(3):338–343. [PubMed]
9. Martinez C, Rietbrock S, Wise L, et al. Antidepressant treatment and the risk of fatal and non-fatal self-harm in first episode of depression: nested case-control study. BMJ. 2005;330(7488):389. [PMC free article] [PubMed]
10. British Columbia Injury Research and Prevention Unit BCIRPU online data tool. [March 18, 2009]. Available at:
11. Centers for Disease Control and Prevention Web-based Injury Statistics Query and Reporting System (WISQARS) [March 18, 2009]. Available at:
12. Williams JI, Young W. Inventory of Studies on the Accuracy of Canadian Health Administrative Databases. Institute for Clinical Evaluative Sciences; Toronto, Canada: 1996.
13. Dormuth CR, Glynn RJ, Neumann P, Maclure M, Brookhart AM, Schneeweiss S. Impact of two sequential drug cost-sharing policies on the use of inhaled medications in older patients with chronic obstructive pulmonary disease or asthma. Clin Ther. 2006;28(6):964–978. [PubMed]
14. Canadian Institute for Health Information Abstracting Manual. Canadian Institute for Health Information; Ottawa, Canada: 1999.
15. Iribarren C, Sidney S, Jacobs DR, Jr, Weisner C. Hospitalization for suicide attempt and completed suicide: epidemiological features in a managed care population. Soc Psychiatry Psychiatr Epidemiol. 2000;35(7):288–296. [PubMed]
16. Warburton RN. Takeup of income-tested health-care premium subsidies: evidence and remedies for British Columbia. Can Tax J. 2005;53(1):1–28.
17. Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chron Dis. 1987;40(5):373–383. [PubMed]
18. Robins JM. Marginal structural models. In: 1997 Proceedings of the American Statistical Association, Section on Bayesian Statistical Science. American Statistical Association; Alexandria, VA: 1998. [May 5, 2009]. pp. 1–10. Available at:
19. Mosholder AD, Willy M. Suicidal adverse events in pediatric randomized, controlled clinical trials of antidepressant drugs are associated with active drug treatment: a meta-analysis. J Child Adolesc Psychopharmacol. 2006;16(1–2):25–32. [PubMed]
20. Dubicka B, Hadley S, Roberts C. Suicidal behaviour in youths with depression treated with new-generation antidepressants: meta-analysis. Br J Psychiatry. 2006;189(5):393–398. [PubMed]
21. Hetrick S, Merry S, McKenzie J, Sindahl P, Proctor M. Selective serotonin reuptake inhibitors (SSRIs) for depressive disorders in children and adolescents. Cochrane Database Syst Rev. 2007;(3):CD004851. [PubMed]
22. Freiman JA, Chalmers TC, Smith H, Jr, Kuebler RR. The importance of β, the type II error and sample size in the design and interpretation of the randomized control trial: survey of 71 “negative” trials. N Engl J Med. 1978;299(13):690–694. [PubMed]
23. British Columbia Coroner's Office Teen deaths: 1997 to 2004. [November 18, 2008]. Available at: www.
24. British Columbia Coroner's Office Suicide statistics: 1997 to 2004. [November 18, 2008]. Available at:
25. BC Stats British Columbia population by selected age groups. [November 18, 2008]. Available at:
26. Schneeweiss S, Avorn J. A review of uses of health care utilization databases for epidemiologic research on therapeutics. J Clin Epidemiol. 2005;58(4):323–337. [PubMed]
27. Greenland S, Robins J. Confoundingandmisclassification. Am J Epidemiol. 1985;122(3):495–506. [PubMed]
28. Kelsey JL, Whittemore AS, Evans AS, Thompson WD. Methods in Observational Epidemiology. 2nd ed. Oxford University Press; New York, NY: 1996.
29. Braitman LE, Rosenbaum PR. Rare outcomes, common treatments: analytic strategies using propensity scores. Ann Intern Med. 2002;137(8):693–695. [PubMed]
30. Rosenbaum P, Rubin D. The central role of the propensity score in observational studies for causal effects. Biometrika. 1983;70(1):41–55.
31. Schneeweiss S, Rassen JA, Glynn RJ, Avorn J, Mogun H, Brookhart MA. High-dimensional propensity score adjustment in studies of treatment effects using health care claims data. Epidemiology. 2009;20(4):512–522. [PMC free article] [PubMed]
32. Maldonado G, Greenland S. Simulation study of confounder-selection strategies. Am J Epidemiol. 1993;138(11):923–936. [PubMed]