Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The role of geography and distance on physician follow-up after a first hospitalization with a diagnosis of a schizophrenia spectrum disorder: A retrospective population-based cohort study in Ontario, Canada

  • Martin Rotenberg ,

    Roles Conceptualization, Formal analysis, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing

    martin.rotenberg@camh.ca

    Affiliations Centre for Addiction and Mental Health, Toronto, Ontario, Canada, Department of Psychiatry, University of Toronto, Toronto, Ontario, Canada, ICES, Toronto, Ontario, Canada

  • Peter Gozdyra,

    Roles Methodology, Resources, Software, Supervision, Validation, Writing – review & editing

    Affiliation ICES, Toronto, Ontario, Canada

  • Kelly K. Anderson,

    Roles Methodology, Supervision, Writing – review & editing

    Affiliations ICES, Toronto, Ontario, Canada, Department of Epidemiology and Biostatistics, Western University, London, Ontario, Canada, Department of Psychiatry, Western University, London, Ontario, Canada

  • Paul Kurdyak

    Roles Conceptualization, Methodology, Supervision, Writing – review & editing

    Affiliations Centre for Addiction and Mental Health, Toronto, Ontario, Canada, Department of Psychiatry, University of Toronto, Toronto, Ontario, Canada, ICES, Toronto, Ontario, Canada

Abstract

Background

Timely follow-up after hospitalization for a schizophrenia spectrum disorder (SSD) is an important quality indicator. We examined the proportion of individuals who received physician follow-up within 7 and 30 days post-discharge by health region and estimated the effect of distance between a person’s residence and discharging hospital on follow-up.

Methods

We created a retrospective population-based cohort of incident hospitalizations with a discharge diagnosis of a SSD between 01/01/2012 and 30/03/2019. The proportion of follow-up with a psychiatrist and family physician within 7 and 30 days were calculated for each region. The effect of distance between a person’s residence and discharging hospital on follow-up was estimated using adjusted multilevel logistic regression models.

Results

We identified 6,382 incident hospitalizations for a SSD. Only 14.2% and 49.2% of people received follow-up care with a psychiatrist within 7 and 30 days of discharge, respectively, and these proportions varied between regions. Although distance from hospital was not associated with follow-up within 7 days of discharge, increasing distance was associated with lower odds of follow-up with a psychiatrist within 30 days.

Conclusion

Post-discharge follow-up is poor across the province. Geospatial factors may impact post-discharge care and should be considered in further evaluation of quality of care.

Introduction

People in the early phase of a schizophrenia spectrum disorder (SSD) can benefit from timely interventions and follow-up to reduce the risk of poor outcomes, including hospitalization, poor treatment adherence, substance misuse and mortality [1, 2]. Those who are hospitalized at the time of first diagnosis may be at even greater risk of poor outcomes, as recent population-based data has described these people are over 13 times more likely to die within 5 years of diagnosis compared with age-matched, population-based controls due to suicide and unnatural deaths [3]. Quality standards in both the USA [4] and Canada [5] have focused on post-discharge quality of care to address the problem of lengthy times to post-discharge follow up after a SSD-related hospitalization. In Ontario (Canada) timely outpatient follow-up within 7-days of an inpatient hospitalization discharge for a SSD is a quality indicator [5]. Timely follow-up with a physician is important to monitor response and tolerability to treatment, assess safety, and support recovery as a person transitions back to the community following hospitalization.

At the present time we do not know what proportion of people who are discharged from hospital after a first diagnosis of SSD are seen within 7 days. Previous work from Ontario, Canada has estimated that only 45% of people were seen by a psychiatrist within 30 days of discharge [6]. There is limited data on how other factors, specifically geographic factors may impact timely follow-up and quality of care for people with a SSD.

A previous meta-analysis of health service use and distance has shown that geographical factors impact health service use, and people who live closer to healthcare facilities have better health related outcomes [7]. Although only three of the 108 studies included heterogenous samples of a variety of mental disorder including SSDs, the three studies included [810] all found a distance decay effect, whereby health service use outcomes in heterogenous samples with mental disorders were inversely associated with how far a person resides from services. Additional studies that have recently been published focusing on people with a SSD in Nigeria [11], Taiwan [12] and the USA [13] which have found distance to hospital impact the likelihood of missing a follow-up clinic appointment post-discharge [11], impacted length of stay [12], and be associated with rehospitalization [13]. None of the previous studies focused on first episode psychosis. There is a dearth of research examining how geospatial factors like distance to services may impact service use and quality of care during an early phase of a SSD, particularly during the potentially high-risk period following hospitalization discharge.

The objectives of this study were to 1) estimate the proportion of people discharged with SSD in Ontario who received timely-follow up while accounting for geographical factors and 2) examine the association between distance from a person’s place of residence to hospital and timely follow-up while adjusting for potential confounding factors.

For the purpose of this study, we examined distance from hospital to a person’s residence (as opposed to a local outpatient service or provider) as based on previous work more than half of people in previous population-based cohorts in the province were not seen by outpatient services after discharge [6]. Moreover, nearly all the specialized early psychosis intervention programs in the province are either hospital-based or located very close to a local hospital [14]. Hospitals would also be most familiar with relevant services that are within local health regions and have more established referral pathways to services that are closer.

Methods

Study design

We used administrative health data from ICES to construct a retrospective cohort of Ontario residents ages 16 to 40 years who were discharged from a hospitalization with an incident diagnosis of SSD between January 1, 2012 and March 31, 2019. People were followed up for 30 days after hospital discharge. These ages were used as they 1) overlap with many of the early psychosis intervention program age criteria in Ontario, and 2) provide coverage beyond age 30 which is closer to the actual the mean age of diagnosis of a SSD in Ontario (32.5 years) [15]. The follow-up period was prior to the onset of the COVID-19 pandemic during which there was a rapid update of virtual care [16].

Setting

Ontario is a Canadian province with a population of approximately 14 million people and has a single-payer universal healthcare system. The province is approximately 1 million km2 and has a large geographic spread. Historically there have been 14 Local Health Integration Network (LHIN) regions (S1 Appendix) that range in size and contain 1 to 2 million people and have been used for the purposes of integration, funding, and planning health care. The LHINs are further subdivided into 97 primary sub-LHINs. For this study, we used LHIN regions and sub-LHIN regions for geospatial analyses and multi-level modeling.

Data

The datasets used were linked using unique encoded identifiers and analyzed at ICES. ICES is an independent, non-profit research institute whose legal status under Ontario’s health information privacy law allows it to collect and analyze health care and demographic data, without consent, for health system evaluation and improvement.

The data sources included the Registered Persons Database; the Ontario Mental Health Reporting System (OMHRS); the Canadian Institute for Health Information Discharge Abstract Database; the National Ambulatory Care Reporting System; the Ontario Health Insurance Plan (OHIP) claims database; the ICES Physician Database; information about Ontario health care institutions funded by the Ministry of Health and Long-Term Care; Immigration, Refugees, and Citizenship Canada Permanent Residents database; the Ontario Marginalization Index; the Postal Code Conversion File; Ontario Road Network; and the Canadian Census of Population. Further details on the data sources are presented in S2 Appendix.

Case definition

We ascertained incident or “first-episode” cases of SSD during a hospitalization by ensuring that people included in the cohort had no previous SSD diagnosis. We included all first hospitalizations with a discharge diagnosis of Schizophrenia, Schizoaffective Disorder and Psychotic Disorder NOS based on DSM-IV / ICD-9 codes (295,298) and ICD-10-CA codes (F20, F25, F29) [17]. All people included in the cohort were residents of Ontario and eligible for OHIP coverage, had at least a 5-year lookback period to ascertain if cases were incident, had a hospitalization length of stay greater than 3 days to ensure sufficient time for diagnosis, and had a length of stay no longer than 90 days as long stays are not common in this clinical population. We also excluded forensic hospitalizations as outpatient care pathways differ in the forensic system. Finally, we restricted our regression models to a subset of the cohort who were hospitalized in designated mental health beds that reported detailed clinical assessment to OMHRS, because this contains information on clinical factors important for understanding physician follow-up [11, 1822].

Exposure: Distance

The exposure of interest was the straight line (Euclidian) distance between the discharging hospital and a person’s place of residence. This was calculated based on location of residence at the time of hospitalization discharge and the location of the hospital.

Outcomes: Physician follow-up

Our primary outcome was defined as follow-up vs. no follow-up with a psychiatrist within 7 days post-discharge. We focused on the 7-day time window to reflect Ontario’s quality standard for follow-up for people with SSD post-discharge.

Our secondary outcomes included: follow-up vs. no follow-up with a psychiatrist within 30 days post-discharge, follow-up vs. no follow-up with any physician (psychiatrist or FP) within 7 days post-discharge, and follow-up vs. no follow-up with any physician within 30 days post-discharge.

FP follow-up was included as part of our secondary outcomes as FPs play an important role in providing mental health services to this population in Ontario [23]. Considering FPs may provide mental health services as part of a general medical visit [24] specific mental health billing codes were not used, and any visit to a FP during the follow-up period was included.

Covariates

We included covariates that were known to impact physician follow-up in this population [6, 11, 1822, 25] and were available in the health administrative data. Further details on the rational for covariate inclusion, covariate categorization and the conceptual framework are included in S3 Appendix.

Sociodemographic factors included age, sex, immigrant status, housing stability, and living arrangements.

Socioenvironmental factors included rural residence based on the Rurality Index of Ontario [26] and area-level marginalization based on the Ontario Marginalization Index (ON-Marg) which includes material deprivation, residential instability, dependency and ethnic concentration dimensions [27]. We also assigned all people to one of the 14 LHINs and 97 primary sub-LHINs based on residence at the time of discharge.

Clinical factors included diagnosis (dichotomized as either schizophrenia/schizoaffective disorder or psychosis NOS), length of stay in hospital, as well as history of problematic substance use, Positive Symptom Scale (PSS), involuntary admission, and insight (full, limited or no insight) which were based on RAI-MH assessment data.

At the hospital level we differentiated between community and academic hospitals. From a service use perspective, we identified any previous contact in the year prior to hospitalization with a psychiatrist, with a FP, any previous ED visits for any mental health reason, and any mental health hospitalization for a reason other than a SSD. Temporally, we also captured the year of discharge to account for potential changes in the structure and funding of mental health services during the study period.

Analysis

We first calculated the crude proportion of people in the entire cohort and in each LHIN for each outcome. Considering that crude proportions do not account for differences in population size across the LHINs, we used the 2016 population estimate in each LHIN based on the 2016 Canadian Census to calculate standardized event ratios (SER). We visualized SER for each outcome of interest using choropleth maps at the LHIN-level for the entire province. A SER less than 1 indicates there was a lower proportion of people with physician follow-up in a region than expected, and a SER greater than 1 indicates there was a larger proportion of people with follow-up than expected.

We described the characteristics of the analytic cohort stratified by our primary outcome using standardized differences [28].

We used a log base 2 transformation of the distance variable due to the wide range and right-skewed distribution. The odds ratio (OR) of this transformation can be interpreted as the change in the odds for every 2-fold increase in distance. Considering the transformed variable was linearly associated with the logit of the outcome, did not change the statistical significance of the exposure estimates, and did not change the estimates of other variables in the models, we reported the log base 2 transformation for all models.

We first fit multilevel logistic regression models for all outcomes of interest with random intercepts at the sub-LHIN region level and no predictor variables. Random intercepts were fit at the sub-LHIN region level to account for the nested structure (people receiving care within health administrative regions) and clustering within these smaller geographic regions. All further multilevel models were fit with random intercepts at the sub-LHIN level.

We used the null models with only the random intercept and no predictor variables to estimate general contextual effects, which provide information on between-cluster variances at the sub-LHIN region level for all outcomes of interest. We calculated intraclass correlation coefficients (ICC) to quantify the variation in the outcomes attributable to the sub-LHIN level [29]. Median odds ratio (MOR) were also calculated and can be interpreted as the median increase in the odds ratio as one moves from an area with a low proportion of follow-up to a high proportion area [30, 31].

We then proceeded to fit unadjusted multilevel logistic regression models and then fit adjusted models for all outcomes of interest with our exposure of interest and covariates. The null models and unadjusted models fit with the primary exposure were used to assess the fit of the fully adjusted models. All effect estimates were presented as ORs with 95% confidence intervals. Estimates were considered significant at the p<0.05 level when the confidence intervals did not overlap with unity.

We conducted sensitivity analyses to assess the impact of potential exposure misclassification, which included comparing model estimates with alternative ways of measuring the primary exposure of interest, including drive time and driving distance from hospital to a person’s residence. Additional sensitivity analyses restricted the sample based on housing stability and extreme distances.

We used R 3.1.2 [32] to conduct descriptive analyses and to visualize data, Stata 15.1 [33] for regression analyses, and ArcGIS 10.2 [34] for network analyses for the sensitivity analyses.

Ethics approval

The study was approved by the Centre for Addiction and Mental Health Research Ethics Board, Toronto, Canada. Individual informed consent was waived. The study data was de-identified prior to analysis.

Results

We identified 25,637 total hospitalizations for SSD between January 1, 2012 and March 31, 2019. Of these, 6,382 hospitalizations were for an incident diagnosis of SSD and met the inclusion criteria and were included in the descriptive analysis. Of these hospitalizations 5,882 were discharges from a psychiatric inpatient unit and were included in the analytic sample that used detailed clinical data. The characteristics of the full cohort are presented in Table 1, and the characteristics of the analytic sample stratified by the primary outcome are presented in S4 Appendix.

thumbnail
Table 1. Characteristics of all people with an incident diagnosis of a schizophrenia spectrum disorder (SSD) diagnosed during a first hospitalization in Ontario between Jan 1, 2012 and March 30, 2019 (n = 6,382).

https://doi.org/10.1371/journal.pone.0287334.t001

Only 14% of the analytic cohort received psychiatric follow-up within 7 days post-discharge. The median distance between a person’s residence and the discharging hospital was 8.0 km (IQR: 3.8 km, 18.4 km) in the group that received follow-up with a psychiatrist within 7 days post-discharge, compared to 9.71 km (IQR: 3.9 km, 24.09 km) in the group that did not receive follow-up.

A greater median length of stay was observed for people who had follow-up with a psychiatrist post-discharge. We observed a greater proportion of people who were admitted to hospital involuntarily in the group that received follow-up by a psychiatrist. Furthermore, a greater proportion of people in the group who were seen by a psychiatrist within 7 days resided in urban areas with greater levels of ethnic concentration. We also observed a greater proportion of people admitted to academic hospitals who were seen by a psychiatrist within 7 days post-discharge. With respect to previous service use, a greater proportion of people who received follow-up with a psychiatrist had previous contact with a psychiatrist and FPs within 1 year prior to hospitalization.

Follow-up by Local Health Integration Network

Only 14.2% of the entire cohort had a visit with a psychiatrist within 7 days post-discharge. The crude proportion of the cohort receiving follow-up differed between LHINs and ranged between 4.8% to 27.6% of the cohort receiving follow-up with a psychiatrist within 7 days post-discharge. Almost half of the entire cohort (49.6%) were seen by a psychiatrist within 30 days post-discharge, and this ranged between 17.7% to 63.2% based on region.

Visualization of the SERs for all outcomes by LHIN are presented in Fig 1. Overall, there was a similar pattern across all outcomes, with a greater than expected level of follow-up in the LHINs with the largest urban centers (Toronto Central LHIN, followed by the Champlain and the Hamilton-Niagara-Haldimand-Brant LHINs).

thumbnail
Fig 1. Maps of standardized event ratios of follow-up with a psychiatrist or any physician within 7 and 30 days of hospitalization discharge by Local Health Integration Network region.

Notes: LHINs with the top 10% SERs are highlighted in red, SER = standardized event ratios, PSYCH7 = Follow-up with a psychiatrist within 7 days of hospitalization discharge, PYCH30 = Follow-up with a psychiatrist within 30 days of hospitalization discharge, ANY7 = Follow-up with any physician within 7 days of hospitalization discharge, ANY30 = Follow-up with any physician within 30 days of hospitalization discharge. *This figure contains information licensed under the Open Government Licence–Ontario, https://www.ontario.ca/page/open-government-licence-ontario.

https://doi.org/10.1371/journal.pone.0287334.g001

Contextual effects at the regional level

An ICC of 0.10 (95% CI = 0.05 to 0.16) indicated that 10% of the variance occurred at the sub-LHIN region level for the primary outcome, follow-up with a psychiatrist within 7 days. A MOR of 1.80 indicated that the sub-LHIN regions with the greatest odds of follow-up had an 80% increased odds of follow-up when compared to the regions with the lowest odds of follow-up. ICCs for secondary outcomes ranged between 0.15 (95% CI = 0.11 to 0.21) for follow-up with a psychiatrist within 30 days and 0.03 (95% CI = 0.01 to 0.06) for FP follow-up within 7 days. MORs ranged between 2.07 for follow-up with a psychiatrist within 30 days and 1.33 for FP follow-up within 7 days.

The association between distance and follow-up

The unadjusted and fully adjusted estimates of the primary effect of distance from residence to hospital on the primary and secondary outcomes are presented in Table 2. Full model estimates, including secondary effects, are presented in Tables 36.

thumbnail
Table 2. Unadjusted and adjusted odds ratios of log base 2 transformation of distance from all models (n = 5,590).

https://doi.org/10.1371/journal.pone.0287334.t002

thumbnail
Table 3. Unadjusted and adjusted mixed logistic regression models of factors associated with psychiatrist follow-up within 7 days of hospitalization discharge (n = 5,590).

https://doi.org/10.1371/journal.pone.0287334.t003

thumbnail
Table 4. Unadjusted and adjusted mixed logistic regression models of factors associated with psychiatrist follow-up within 30 days of hospitalization discharge (n = 5,590).

https://doi.org/10.1371/journal.pone.0287334.t004

thumbnail
Table 5. Unadjusted and adjusted mixed logistic regression models of factors associated with any physician follow-up within 7 days of hospitalization discharge (n = 5,590).

https://doi.org/10.1371/journal.pone.0287334.t005

thumbnail
Table 6. Unadjusted and adjusted mixed logistic regression models of factors associated with any physician follow-up within 30 days of hospitalization discharge (n = 5,590).

https://doi.org/10.1371/journal.pone.0287334.t006

Primary effects

Distance was not associated with the primary outcome of follow-up with a psychiatrist within 7 days post-discharge in either the unadjusted (OR = 0.99, 95% CI = 0.93 to 1.02) or adjusted models (aOR = 1.00, 95% CI = 0.95 to 1.05). However, distance was independently associated with the secondary outcome of follow-up with a psychiatrist within 30 days post-discharge. In the fully adjusted model, there was a 4% reduction in the odds of follow-up for every 2-fold increase in distance between a person’s residence and location of hospitalization (aOR 0.96, 95% CI = 0.92 to 0.99). Distance was not associated with follow-up with any physician within 7 days post-discharge (aOR = 1.01, 95% CI = 0.97 to 1.05) or 30 days post-discharge (aOR = 0.98, 95% CI = 0.95 to 1.02).

Secondary effects

The primary effect estimates (between distance and follow-up) were obtained from models constructed to obtain unbiased effect estimates of this relationship. Unadjusted and adjusted estimates for all covariates are also presented for completeness (Tables 36). These secondary effects must be interpreted with caution, as these estimates may be biased due to the models not being constructed to obtain unconfounded estimates of the relationship between the secondary exposures and outcomes for people with a history of substance use [35, 36].

In the adjusted model estimating the odds of psychiatric follow-up within 7 days of discharge (Table 3) there was a statistically significant increase in the odds of follow-up in areas with moderate levels of dependency (Q3), areas with higher levels of ethnic concentration (Q3,Q4), for people who were admitted voluntarily, admitted to a teaching hospital, had a longer length of stay and had previous psychiatric contact. There was a statistically significant reduction in the odds of follow-up with a psychiatrist within 7 days of discharge if a person had a history of substance use and with each increasing year in the cohort during which an admission occurred.

In the adjusted model for psychiatric follow-up within 3 days of discharge (Table 4) there was a statistically significant increase in the odds of follow-up in areas with higher levels of ethnic concentration (Q3,Q4), for people with a greater length of stay, and had previous contact with a psychiatrist prior to the hospitalization. As reported in the primary effects section, there was statistically significant and independent association between distance and follow-up with a psychiatrist within 30 days of follow-up, whereby there was a reduction in the in the odds of follow-up for every 2-fold increase in distance between a person’s residence and location of hospitalization. Furthermore, in the adjusted model there was also a statistically significant reduction in the odds of follow-up if a person lived alone, lived in a rural area, used substances, had no insight and had a previous mental health hospitalization not for a SSD, and with each increasing calendar year during the cohort period.

In the adjusted model estimating the odds of any physician follow-up within 7 days of discharge (Table 5) there was a statistically significant increase in the odds of follow-up in areas with moderate to high levels of ethnic concentration (Q3,Q4,Q5), moderate residential instability (Q3), in people with a longer length of stay in hospital, admitted to a teaching hospital and had previous contact with a family physician in the year prior to hospitalization. There was a statistically significant reduction in the odds of follow-up if a person lived alone and with each increasing calendar year during the cohort period.

In the adjusted model estimating any physician follow-up within 30 days of discharge there was a statistically significant increase in the odds of follow-up in areas with higher levels of ethnic concentration (Q4,Q5), and if a person was admitted as a voluntary patient, spent a longer time in hospital, and had contact with a psychiatrist or family physician for mental health reasons in the year prior to admission. There was a statistically significant reduction in the odds of follow-up if a person lived alone, in areas with higher levels of deprivation (Q5), if a person uses substances, has limited or no insight, and had a mental health hospitalization for a disorder other than a SSD in the year prior to the admission.

Sensitivity analyses

We found no substantive change in estimates or significance when we compared models using Euclidian distance to models with alternative definitions of the exposure specifically, driving distance and driving time. Moreover, model estimates from sensitivity analyses that restricted the sample based on housing stability and extreme distances did not suggest misclassification of the primary exposure. Further details on the sensitivity analyses are presented in S5 Appendix.

Discussion

Less than 3 out of every 20 people who are discharged from hospital in Ontario after a new diagnosis of SSD received follow-up with a psychiatrist within 7 days. Outpatient follow-up within 7 days post-discharge is the quality standard in Ontario for all people with schizophrenia during any phase of the illness [5]. Within 30 days post-discharge, slightly less than half of all people in the cohort were seen by a psychiatrist. In a similar clinical sample from the same jurisdiction but from an earlier time period (1999 and 2009) 44.8% of all people with a new diagnosis of SSD received outpatient follow-up with a psychiatrist within 30 days post-discharge [6]. Although our current study found a slight improvement in the proportion of people who received follow-up with a psychiatrist within 30 days post-discharge, it remains that less than half of all people who are hospitalized at first diagnosis with SSD receive follow-up care with a psychiatrist within 30 days post-discharge, and this trend has persisted over two decades.

Our findings also highlight disparities in post-discharge quality of care based on where a person resides in the province. The greatest proportion of people who received timely follow-up across all outcomes live in the health regions which include the largest metropolitan centers in the province (i.e., Toronto, Ottawa, and Hamilton).

The low proportion of follow-up we observed outside of the large metropolitan centers is in keeping with the unequal distribution of psychiatrists across the province [37]; however, other factors including the structure and delivery of mental health services must also be considered. From a policy perspective further consideration of the incentivization of quality care would be reasonable to examine despite initial negative studies [38]. Unlike previous work which has focused on broader populations, it would be warranted to study incentives targeted to specific clinical population based on service need and geography. It is unlikely that incentivization of care alone will improve quality, and alternative models to standard mental health care delivery that are both more efficient and responsive to the care needs of service users should also be considered.

We found no association between the distance from a person’s place of residence and the discharging hospital and the primary outcome of follow-up with a psychiatrist within 7 days post-discharge. However, we found a modest association between our secondary outcome of psychiatric follow-up within 30 days post-discharge with distance. These findings are suggestive of a distance decay effect when examined for specific types of physician follow-up within 30 days of hospitalization discharge. This finding is consistent with the literature on this topic from other areas of health care [7] and in similar clinical populations [810]. The fact that we did not find a significant association between distance and the outcomes that assessed care contacts within 7 days of discharge may be due to only a small proportion of the sample receiving follow-up within 7 days of hospitalization discharge. Moreover, it is possible that care received within 7 days of hospitalization discharge may be driven by factors other than distance, specifically factors like past service use, which may be more strongly associated with follow-up.

Further replication and examination is warranted to determine the degree of clinical significance of this distance decay finding, considering other factors (e.g., past service use and clinical factors) may be more strongly associated with follow-up. It may also be important to consider how distance may impact length of stay in hospital and the risk of rehospitalization, as previously examined in non-first episode cohorts of people with SSD diagnoses [12, 13]. Clinically, our results support the fact that greater care should be taken when discharging people from a first inpatient hospitalization who live further away from hospital and in specific regions where there may be poor access to services.

Several secondary exposures (ranging from area level marginalization, rurality, past psychiatric follow-up, family physician follow-up, substance use, level of insight, length of stay and year of admission) were associated with the odds of follow-up, it is important to note the adjusted models were built to obtain unbiased estimates of the primary exposure of distance on the outcomes of interest. These findings are not surprising and are in keeping with previous literature which has found sociodemographic, socioenvironmental, clinical, service use, hospital and temporal factors to all be associated with follow-up in this clinical population [6, 11, 1822, 25]. Although the secondary exposure estimates may be biased, these associations warrant further examination in future work.

It is also important to note that the delivery of outpatient services has rapidly changed since March 2020 in response to the SARS-CoV-2 pandemic. The pandemic has generated a swift reconfiguration of service delivery to virtual care, including telepsychiatry [39] which has historically been underutilized [40]. Prior to March 2020, telepsychiatry had been identified as an avenue to improve poor rates of post-discharge care to people with mental health needs [40] and as a means to improve disparities in physician supply that exist between regions. Yet, challenges and barriers to the use of virtual care [41] may also be determined in part by geography (e.g., access to highspeed internet to successfully use a virtual platform, having a suitable location to privately meet with a provider virtually). Further examination of the use of virtual care in populations with SSD in relation to post-discharge follow-up should be undertaken.

Strengths and limitations

To our knowledge, this is the first study that examines the impact of geospatial factors and distance on post-discharge quality of care among people with SSD. The focus on people diagnosed after first hospitalization is also a specific population that may be at a greater risk of negative outcomes. We also conducted sensitivity analyses to ensure our exposure of interest was not misclassified.

Although we excluded people with prior diagnoses of SSD, we were not able to account for people with affective psychoses (e.g., bipolar I disorder), who are at risk of diagnostic transition [42] and may be diagnosed with SSD during a hospitalization. This is important as people with affective psychoses are eligible for specialized mental health services like early psychosis intervention programs [43] that also provide services to people with SSD, which may impact service use post-discharge. We also restricted our regression analyses to people who had been discharged from a psychiatric bed (92.2% of all hospitalizations), as non-psychiatric beds do not report data to OMHRS which includes important clinical data that may confound the exposure-outcome relationship we were interested in.

We focused on physician services and could not account for contacts with non-physicians including nurses, case managers, social workers, and psychologists who all provide important aspects of care to this clinical population. Outpatient contacts with non-physician providers post-discharge would also meet the Ontario quality standard. Moreover, we were not able to ascertain whether a follow-up visit was arranged (which is the standard of care) but the person failed to attend, or if an appointment was never booked using the health administrative data available. Future studies may consider linking institutional and clinic-level data to examine these factors further. We were also unable to account for follow-up care which may have occurred outside of Ontario, which may be more likely to occur in regions bordering other provinces.

This study focused on one distance exposure, however other distance exposures may be important to consider in the context of post-discharge quality of care (e.g., the distance between a person’s residence to their outpatient provider). As noted, nearly all early psychosis intervention programs in the province are either hospital based or located near a hospital [14] and many people in this clinical population do not receive any outpatient care [6]. The exposure used for this study is a reasonable one to focus on considering the study’s emphasis on post-hospitalization discharge quality of care as well as the novel nature of this work in this clinical population and jurisdiction. This study highlights that it may be feasible to further examine other geospatial exposures in future work.

Conclusions

Very few people who are diagnosed with SSD at the time of a first hospitalization received timely follow-up with a psychiatrist post-discharge in Ontario, Canada. Whether a person received follow-up care differs based on where they reside in the province. We did not find distance to be associated with psychiatric follow-up within 7 days post-discharge, however, increasing distance was independently associated with a reduction in the odds of psychiatric follow-up within 30 days post-discharge.

Our findings support the need for further work to evaluate the impact of geospatial factors on post-discharge quality of care and downstream outcomes in this clinical population. We need to consider developing, implementing, and evaluating interventions to improve post-discharge follow-up for people who have been diagnosed with SSD during a hospitalization.

Supporting information

S1 Appendix. Location of the 14 Local Health Integration Network (LHIN) regions in Ontario.

https://doi.org/10.1371/journal.pone.0287334.s001

(DOCX)

S3 Appendix. Covariate description, categorization and conceptual framework.

https://doi.org/10.1371/journal.pone.0287334.s003

(DOCX)

S4 Appendix. Characteristics of the analytic cohort who received follow-up with psychiatrist within 7 days of hospitalization discharge vs. those who did not (n = 5,882).

https://doi.org/10.1371/journal.pone.0287334.s004

(DOCX)

Acknowledgments

Parts of this material are based on data and information compiled and provided by: the Ontario Ministry of Health (MOH), Canadian Institute for Health Information (CIHI). The analyses, conclusions, opinions and statements expressed herein are solely those of the authors and do not reflect those of the funding or data sources; no endorsement is intended or should be inferred. Parts of this material are based on data and/or information compiled and provided by CIHI. However, the analyses, conclusions, opinions and statements expressed in the material are those of the author(s), and not necessarily those of CIHI. Parts or whole of this material are based on data and/or information compiled and provided by Immigration, Refugees and Citizenship Canada (IRCC). However, the analyses, conclusions, opinions and statements expressed in the material are those of the author(s), and not necessarily those of IRCC.

References

  1. 1. Anderson KK, Norman R, MacDougall A, Edwards J, Palaniyappan L, Lau C, et al. Effectiveness of Early Psychosis Intervention: Comparison of Service Users and Nonusers in Population-Based Health Administrative Data. Am J Psychiatry. 2018;175: 443–452. pmid:29495897
  2. 2. Kane JM, Robinson DG, Schooler NR, Mueser KT, Penn DL, Rosenheck RA, et al. Comprehensive Versus Usual Community Care for First-Episode Psychosis: 2-Year Outcomes From the NIMH RAISE Early Treatment Program. Am J Psychiatry. 2015;173: 362–372. pmid:26481174
  3. 3. Kurdyak P, Mallia E, de Oliveira C, Carvalho AF, Kozloff N, Zaheer J, et al. Mortality After the First Diagnosis of Schizophrenia-Spectrum Disorders: A Population-based Retrospective Cohort Study. Schizophr Bull. 2021;47: 864–874. pmid:33459778
  4. 4. National Committee for Quality Assurance. Follow-Up After Hospitalization for Schizophrenia (7- and 30-day). 2012. Available: https://www.qualityforum.org/QPS/
  5. 5. Health Quality Ontario. Quality Standards Schizophrenia: Care for Adults in Hospitals. 2016.
  6. 6. Anderson KK, Kurdyak P. Factors Associated with Timely Physician Follow-up after a First Diagnosis of Psychotic Disorder. Can J Psychiatry Rev Can Psychiatr. 2017;62: 268–277. pmid:27738261
  7. 7. Kelly C, Hulme C, Farragher T, Clarke G. Are differences in travel time or distance to healthcare for adults in global north countries associated with an impact on health outcomes? A systematic review. BMJ Open. 2016;6: e013059. pmid:27884848
  8. 8. Joseph AE, Boeckh JL. Locational variation in mental health care utilization dependent upon diagnosis: a Canadian example. Soc Sci Med [D]. 1981;15: 395–340. pmid:7323872
  9. 9. Skarsvåg K, Wynn R. Travel Time and the Use of Psychiatric Outpatient Clinic Services in Coastal Northern Norway. Can J Psychiatry. 2004;49: 153–154. pmid:15065754
  10. 10. McCarthy JF, Blow FC, Valenstein M, Fischer EP, Owen RR, Barry KL, et al. Veterans Affairs Health System and Mental Health Treatment Retention among Patients with Serious Mental Illness: Evaluating Accessibility and Availability Barriers. Health Serv Res. 2007;42: 1042–1060. pmid:17489903
  11. 11. Nwefoh E, Aguocha CM, Achor J, Uwakwe R, Onyeama G, Bakare MO, et al. Missed post-hospitalisation clinic appointment in a psychiatric hospital in Southeast Nigeria. Ment Health Relig Cult. 2018;21: 564–577.
  12. 12. Tseng K-C, Hemenway D, Kawachi I, Subramanian SV, Chen WJ. Travel Distance and the Use of Inpatient Care among Patients with Schizophrenia. Adm Policy Ment Health. 2008;35: 346–356. pmid:18512144
  13. 13. Hudson CG. Five-year rehospitalization experience of a state-wide cohort of persons with schizophrenia. Soc Psychiatry Psychiatr Epidemiol. 2019;54: 861–870. pmid:30603806
  14. 14. Early Psychosis Intervention Ontario Network. 2020 [cited 24 Jun 2020]. Available: https://help4psychosis.ca/looking-for-help/where-to-get-help/
  15. 15. Rotenberg M, Tuck A, Anderson KK, McKenzie K. The Incidence of Psychotic Disorders and Area-level Marginalization in Ontario, Canada: A Population-based Retrospective Cohort Study. Can J Psychiatry. 2021; 07067437211011852. pmid:33896210
  16. 16. Saunders NR, Kurdyak P, Stukel TA, Strauss R, Fu L, Guan J, et al. Utilization of Physician-Based Mental Health Care Services Among Children and Adolescents Before and During the COVID-19 Pandemic in Ontario, Canada. JAMA Pediatr. 2022;176: e216298. pmid:35129604
  17. 17. Kurdyak P, Lin E, Green D, Vigod S. Validation of a Population-Based Algorithm to Detect Chronic Psychotic Illness. Can J Psychiatry Rev Can Psychiatr. 2015;60: 362–368. pmid:26454558
  18. 18. Marcus SC, Chuang C-C, Ng-Mak DS, Olfson M. Outpatient Follow-Up Care and Risk of Hospital Readmission in Schizophrenia and Bipolar Disorder. Psychiatr Serv Wash DC. 2017;68: 1239–1246. pmid:28669289
  19. 19. Olfson M, Marcus SC, Doshi JA. Continuity of Care After Inpatient Discharge of Patients With Schizophrenia in the Medicaid Program: A Retrospective Longitudinal Cohort Analysis. J Clin Psychiatry. 2010;71: 0–0. pmid:20441730
  20. 20. Marino L, Wissow LS, Davis M, Abrams MT, Dixon LB, Slade EP. Predictors of outpatient mental health clinic follow-up after hospitalization among Medicaid-enrolled young adults. Early Interv Psychiatry. 2016;10: 468–475. pmid:25639939
  21. 21. Adeponle AB, Baduku AS, Adelekan ML, Suleiman GT, Adeyemi SO. Prospective Study of Psychiatric Follow-up Default and Medication Compliance after Discharge at a Psychiatric Hospital in Nigeria. Community Ment Health J. 2009;45: 19–25. pmid:18633705
  22. 22. Compton MT. Barriers to Initial Outpatient Treatment Engagement Following First Hospitalization for a First Episode of Nonaffective Psychosis: A Descriptive Case Series. J Psychiatr Pract. 2005;11: 62–69. pmid:15650625
  23. 23. Anderson KK, Archie S, Booth RG, Cheng C, Lizotte D, MacDougall AG, et al. Understanding the role of the family physician in early psychosis intervention. BJPsych Open. 2018;4: 447–453. pmid:30450223
  24. 24. Craven M, Cohen M, Campbell D, Williams J, Kates N. Mental Health Practices of Ontario Family Physicians: A Study Using Qualitative Methodology. Can J Psychiatry. 1997;42: 943–949. pmid:9429064
  25. 25. Kurdyak P, Vigod SN, Newman A, Giannakeas V, Mulsant BH, Stukel T. Impact of Physician Follow-Up Care on Psychiatric Readmission Rates in a Population-Based Sample of Patients With Schizophrenia. Psychiatr Serv. 2018;69: 61–68. pmid:28859584
  26. 26. Glazier R, Zagorski B, Rayner J. Comparison of primary care models in Ontario by demographics, case mix and emergency department use, 2008/09 to 2009/10. Toronto, ON: Institute for Clinical Evaluative Sciences; 2012. Available: https://www.ices.on.ca/flip-publication/comparison-of-primary-care-models-in-ontario-by-demographics/files/assets/basic-html/page15.html
  27. 27. Matheson F, Dunn JR, Smith KLW, Moineddin R, Glazier RH. Development of the Canadian Marginalization Index: A New Tool for the Study of Inequality. Can J Public Health Rev Can Santee Publique. 2012;103: S12–S16.
  28. 28. Austin PC. Using the Standardized Difference to Compare the Prevalence of a Binary Variable Between Two Groups in Observational Research. Commun Stat—Simul Comput. 2009;38: 1228–1234.
  29. 29. Austin PC, Merlo J. Intermediate and advanced topics in multilevel logistic regression analysis. Stat Med. 2017;36: 3257–3277. pmid:28543517
  30. 30. Austin PC, Wagner P, Merlo J. The median hazard ratio: a useful measure of variance and general contextual effects in multilevel survival analysis. Stat Med. 2017;36: 928–938. pmid:27885709
  31. 31. Austin PC, Stryhn H, Leckie G, Merlo J. Measures of clustering and heterogeneity in multilevel Poisson regression analyses of rates/count data. Stat Med. 2018;37: 572–589. pmid:29114926
  32. 32. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2020. Available: https://www.R-project.org/
  33. 33. Stata Statistical Software: Release 15.1. College Station, TX: StataCorp LLC; 2017.
  34. 34. ArcGIS version 10.2. Redlands, CA: Environmental Systems Research Institute; 2018.
  35. 35. Bandoli G, Palmsten K, Chambers CD, Jelliffe-Pawlowski LL, Baer RJ, Thompson CA. Revisiting the Table 2 Fallacy: A Motivating Example Examining Preeclampsia and Preterm Birth. Paediatr Perinat Epidemiol. 2018;32: 390–397. pmid:29782045
  36. 36. Westreich D, Greenland S. The Table 2 Fallacy: Presenting and Interpreting Confounder and Modifier Coefficients. Am J Epidemiol. 2013;177: 292–298. pmid:23371353
  37. 37. Kurdyak P, Stukel TA, Goldbloom D, Kopp A, Zagorski BM, Mulsant BH. Universal coverage without universal access: a study of psychiatrist supply and practice patterns in Ontario. Open Med. 2014;8: e87–e99. pmid:25426177
  38. 38. Rudoler D, Oliveira C de, Cheng J, Kurdyak P. Payment incentives for community-based psychiatric care in Ontario, Canada. CMAJ. 2017;189: E1509–E1516. pmid:29229712
  39. 39. Gratzer D, Torous J, Lam RW, Patten SB, Kutcher S, Chan S, et al. Our Digital Moment: Innovations and Opportunities in Digital Mental Health Care. Can J Psychiatry. 2021;66: 5–8. pmid:32603188
  40. 40. Serhal E, Crawford A, Cheng J, Kurdyak P. Implementation and Utilisation of Telepsychiatry in Ontario: A Population-Based Study. Can J Psychiatry. 2017;62: 716–725. pmid:28541753
  41. 41. Crawford A, Serhal E. Digital Health Equity and COVID-19: The Innovation Curve Cannot Reinforce the Social Gradient of Health. J Med Internet Res. 2020;22. pmid:32452816
  42. 42. Consoli A, Brunelle J, Bodeau N, Louët E, Deniau E, Perisse D, et al. Diagnostic transition towards schizophrenia in adolescents with severe bipolar disorder type I: An 8-year follow-up study. Schizophr Res. 2014;159: 284–291. pmid:25217364
  43. 43. Durbin J, Selick A, Langill G, Cheng C, Archie S, Butt S, et al. Using Fidelity Measurement to Assess Quality of Early Psychosis Intervention Services in Ontario. Psychiatr Serv. 2019;70: 840–844. pmid:31159664