Study setting
The pandemic significantly disrupted access to health services across various settings in Australia, notably affecting elective surgeries, routine care, and mental health services (16, 17). To address these challenges, the Australian government rapidly expanded telehealth services, which included funding for general practitioners, specialists, and allied health professionals to provide care through phone or video consultations. This rapid expansion towards counteracting the reduction in In-person consultations and ensuring continued access to health services, particularly mental health care. To understand the events that occurred around March 2020, we may state that two interventions occurred, one natural (COVID-19) and another one (the expansion of telehealth MBS items) meant to counter the effect of the natural one. This study aimed at evaluating the impact of the March 2020 COVID-19 pandemic and the subsequent telehealth policy change on general practitioners’ (GP) consultations for mental health conditions in Australia. The analysis focuses on three distinct periods: i) The immediate impact period (March 2020 – February 2021); ii) The pandemic recovery period (March 2021 – February 2022); and iii) The post-pandemic period (March 2022 – February 2023).
This study hypothesizes that: i) The pandemic had no significant effect on In-person mental health consultations by GPs; ii) There was no significant change in the overall Total GP Consultations by GPs for mental health conditions in Australia before and after the implementation of telehealth services expansion. For clear understanding of the effects of telehealth policy on overall mental health consultations by GPs, we first tested whether the COVID-19 significantly caused a decline in In-person mental health consultations. Then, an evaluation of the telehealth changes was performed to determine whether it sustained the overall mental health consultations by GPs at its pre-pandemic level.
Data Preparation
We analyzed the publicly available data on GP mental health consultations in Australia from March 2017 to February 2023, indicating 36 months period before and after the introduction of telehealth services which occurred in March 2020, accessed on the website of Services Australia(1). The datasets were obtained from the Services Australia website and included monthly counts of three types of consultations: In-person, Phone, and Video consultations, represented by MBS item numbers – 2713, 92115, and 92127 respectively. These are defined as follows: i) 2713- This refers to In-person professional consultation by a GP for a mental condition (not less than 20 minutes) including taking a relevant history, identifying the presenting issue (if not previously documented), providing treatment and advice, and making referrals for additional services or treatments if necessary, as well as recording the consultation's outcomes (Jones & Smith, 2020); ii) 92115- This indicates telehealth consultation by a GP for a mental health condition (not less than 20 minutes) involving taking a relevant history, identifying the presenting issue (if not previously documented), providing treatment and advice, and making referrals for additional services or treatments if necessary, along with documenting the consultation's outcomes(Doe & Roe, 2019); and iii) 92127 – This represents phone consultation by a GP for a mental disorder lasting at least 20 minutes, including taking a relevant history, identifying the presenting issue (if not previously documented), providing treatment and advice, and making referrals for additional services or treatments if necessary, and recording the outcomes of the consultation (Brown & Green, 2018).
Guidelines and details for the highlighted MBS items (2713, 92115, and 92127) can are provided on the Services Australia website on Medicare Benefits Schedule (Department of Health, 2023). These descriptions above depict a clear outline of the types of consultations available for mental health disorders provided by GPs to address and document patient needs (Department of Human Services, 2023).
Insights into the dataset were gained through a detailed descriptive analysis of the dataset. The dataset was analyzed for the separate consultations including: i) In-person consultations only; ii) Phone consultations only; iii) Video consultations only; and iv) Total GP mental health consultations (consisting of In-person, Phone, and Video consultations). Summary statistics, and overall distribution of the consultation data were described.
Statistical analysis
'Total GP Consultations', derived from the sum of In-person, Phone, and Video consultations for each month was created to provide an understanding of the overall trend in GP consultations. We performed ANOVA to compare the mean number of the modes of the three modes of consultations and the Total GP Consultations across the three specified periods. The ANOVA tested the null hypothesis that there were no differences in the mean number of consultations across the periods against the alternative hypothesis that at least one period differed. The ANOVA model involved fitting an ordinary least squares (OLS) model and then performing the ANOVA using the fitted model.. If the ANOVA results indicated significant differences, a post-hoc analysis using Tukey's Honestly Significant Difference (HSD) test was conducted to identify which specific periods differed from each other. The post-hoc analysis provided pairwise comparisons between periods, highlighting where significant differences were present.
Furthermore, we evaluated the telehealth policy changes of the Australian government through the use Interrupted Time Series (ITS) analysis via Ordinary Least Squares (OLS) regression. The ITS Ordinary Least Squares method is simple and accessible, with easy-to-interpret coefficients indicating intervention impacts [10–11]. To test this hypothesis, an Interrupted Time Series (ITS) model was employed to analyze the trend of mental health consultations pre- and post- the introduction of telehealth services. The intervention point for the analysis was March 2020 when several telehealth MBS items (including for mental health services) were expanded and funded by the Government of Australia.
The ITS model specification is as follows:
Yt represents the number of mental health consultations at time t.
Tt denotes the time since the start of the observation period.
Dt is a dummy variable representing the intervention period (0 before intervention, 1 after intervention).
Pt represents the post-intervention trend (interaction term between time and intervention period),
β0 is the intercept.
β1 is the coefficient for the underlying pre-intervention trend.
β2 measures the immediate effect of the intervention.
β3 represents the change in trend post-intervention; and
ϵt is the error term.
The key hypotheses to be tested are: i) Immediate effect - The coefficient β2 is expected to be significantly positive, indicating an immediate increase in mental health consultations following the introduction of telehealth services; ii) Change in trend: The coefficient β3 is anticipated to be significantly positive, suggesting an improved or sustained trend in mental health consultations post-intervention. These outcomes will support the alternate hypothesis that telehealth mitigated the reduction in In-person consultations, ensuring continued access to mental health services during the pandemic.
The dataset was segmented into pre- and post-intervention periods using a binary indicator, with March 2020 as the period when the intervention was implemented. The date column was converted to a datetime format and set as the index to facilitate time series analysis. The following dependent variables were examined: i) In-person consultations – This indicates the services offered through physical contact with patients in GP offices; and ii) Total GP Consultations - This was calculated by summing the service counts of In-person, Phone, and Video consultations for each month. The primary variables used in the study were: i) Time Trend (T) - A continuous variable representing time, starting from March 2017 in the dataset to February 2023; ii) Intervention Dummy (D): This is the binary variable denoting the intervention period, coded as 0 for months before March 2020 and 1 for March 2020 and onwards; iii) Post-intervention Trend (P) – This connotes an interaction term between the time trend (T) and the intervention dummy (D), representing the change in trend post-intervention.
The OLS regression model was fitted, incorporating time trend, intervention dummy, post-.the binary intervention indicator, and their interaction to capture immediate level changes and trend (slope) changes post-intervention. The model was fitted to the data, and actual versus fitted values were plotted with the intervention point clearly marked. A detailed model summary was generated to provide statistical insights into the coefficients and overall model fit, facilitating a comprehensive evaluation of the intervention's effectiveness (Penfold & Zhang, 2013). Counterfactual predictions were also generated from the OLS models. Diagnostic checks, including the Durbin-Watson statistic and plots of the Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF), were performed to ensure no significant autocorrelation in the residuals. We did not model the Phone and Video Consultations as they had zero as their pre intervention values. This is because developing a model including zero as baseline data in an interrupted time series (ITS) analysis may be misleading as it violates the assumptions of statistical models, leading to inflated variance estimates and biased coefficient estimates (Bernal et al., 2017; Linden, 2015). The robust framework used for this evaluation was essential to ensure a comprehensive understating of the impact of the pandemic and policy changes on GP consultations for mental health conditions in Australia.
Furthermore, an evaluation of the impact of the telehealth policy changes was conducted, considering the periods defined above using an Analysis of Variance (ANOVA). The ANOVA analysis aimed to determine whether there were statistically significant differences in the mean number of consultations across specified periods affected by the pandemic and subsequent policy changes. To this analysis, the dataset was filtered to include only the periods from March 2020 to February 2021, March 2021 to February 2022, and March 2022 to February 2023. These periods were chosen to capture the immediate impact of the pandemic and the policy changes, as well as the longer-term effects. The filtering process ensured that only relevant data points were included in the analysis, providing a focused view of the changes over time.
All analyses were conducted using Python, leveraging the statsmodels library for conducting ANOVA and time series modeling and analysis (Seabold & Perktold, 2010); matplotlib library for data visualization and trend plotting (Hunter, 2007); pandas library for data preparation and manipulation (McKinney, 2010); and Seaborn for plotting the actual consultations and counterfactual values (Waskom et al., 2017).