All patients with an invasive cancer diagnosis were included. Patients diagnosed only through autopsy or death certificate (<1.5% of patients) were excluded. Data were extracted for patients with more than 100,000 person-years or more of survival time; thus, certain uncommon and aggressive cancers were excluded, including Kaposi’s sarcoma, multiple myeloma, hepatobiliary cancer, male breast cancers, and mesotheliomas. Since certain cancers represent a heterogeneous group of disease (e.g. leukemia, lymphoma), these cancers were grouped for certain parts of the analysis, so they could be reported accurately.
Mortality codes in SEER are assigned from death certificates, completed by the doctor caring for the patient at the time of demise. Patients were considered to have committed suicide if the cause of death was coded as: suicide and self-inflicted injury (50220). Patients with other causes of death, including accident, homicide, and legal intervention were excluded from the present analysis. Survival time in SEER is measured in months, and the smallest nonzero value is 1 month, which was the minimum time to any event.
Notably, SEER does not code comorbidities or diagnoses associated with suicide, including suicidal ideation, previous suicide attempts, or use of anti-depressive medications. The observed associations between cancer and suicide may be confounded by psychiatric disorders and the use of medications, but we are unable to control for these factors in the current work. These are limitations of the analysis and limit the interpretability of the results.
For objective 1, we calculated standardized mortality ratios (SMRs), which provide the relative risk of death for patients with cancer as compared to all US residents, stratified by cancer subgroup3 (link),22 ,27 . Data were characterized with SMRs adjusted by age, race, and sex to the US population over the same time. Five-year age categories were used for standardization using SEER*Stat 8.2.1 and Microsoft Excel 15.0.4 (Microsoft, Redmond, WA)27 –29 (link).
For objective 2, odds ratios (ORs) with 95% CIs were calculated based on the number of observed events per patient subgroup. Further, the absolute and relative number of suicides per patient age group (at time of diagnosis) were calculated. We also performed a survival analysis using a Cox proportional hazards model to calculate hazard ratios (HRs), with the survival time being from diagnosis until suicide, and non-suicide deaths plus living patients being censored.