Article Text

Original research
Diagnostic accuracy of subjective dyspnoea in detecting hypoxaemia among outpatients with COVID-19: a retrospective cohort study
  1. Linor Berezin1,
  2. Alice Zhabokritsky1,
  3. Nisha Andany1,2,
  4. Adrienne K Chan1,2,
  5. Jose Estrada-Codecido2,
  6. Andrea Gershon1,3,
  7. Philip W Lam1,2,
  8. Jerome A Leis1,2,
  9. Scott MacPhee4,
  10. Samira Mubareka2,5,
  11. Andrew E Simor1,2,
  12. Nick Daneman1,2
  1. 1Faculty of Medicine, University of Toronto, Toronto, Ontario, Canada
  2. 2Division of Infectious Diseases, Sunnybrook Health Sciences Centre, University of Toronto, Toronto, Ontario, Canada
  3. 3Division of Respirology, Sunnybrook Health Sciences Centre, University of Toronto, Toronto, Ontario, Canada
  4. 4Department of Nursing, Sunnybrook Health Sciences Centre, University of Toronto, Toronto, Ontario, Canada
  5. 5Department of Laboratory Medicine and Pathology, University of Toronto, Toronto, Ontario, Canada
  1. Correspondence to Dr Nick Daneman; Nick.Daneman{at}sunnybrook.ca

Abstract

Objectives The majority of patients with mild-to-moderate COVID-19 can be managed using virtual care. Dyspnoea is challenging to assess remotely, and the accuracy of subjective dyspnoea measures in capturing hypoxaemia have not been formally evaluated for COVID-19. We explored the accuracy of subjective dyspnoea in diagnosing hypoxaemia in COVID-19 patients.

Methods This is a retrospective cohort study of consecutive outpatients with COVID-19 who met criteria for home oxygen saturation monitoring at a university-affiliated acute care hospital in Toronto, Canada from 3 April 2020 to 13 September 2020. Dyspnoea measures were treated as diagnostic tests, and we determined their sensitivity (SN), specificity (SP), negative/positive predictive value (NPV/PPV) and positive/negative likelihood ratios (+LR/−LR) for detecting hypoxaemia. In the primary analysis, hypoxaemia was defined by oxygen saturation <95%; the diagnostic accuracy of subjective dyspnoea was also assessed across a range of oxygen saturation cutoffs from 92% to 97%.

Results During the study period, 89/501 (17.8%) of patients met criteria for home oxygen saturation monitoring, and of these 17/89 (19.1%) were diagnosed with hypoxaemia. The presence/absence of dyspnoea had limited accuracy for diagnosing hypoxaemia, with SN 47% (95% CI 24% to 72%), SP 80% (95% CI 68% to 88%), NPV 86% (95% CI 75% to 93%), PPV 36% (95% CI 18% to 59%), +LR 2.4 (95% CI 1.2 to 4.7) and −LR 0.7 (95% CI 0.4 to 1.1). The SN of dyspnoea was 50% (95% CI 19% to 81%) when a cut-off of <92% was used to define hypoxaemia. A modified Medical Research Council dyspnoea score >1 (SP 98%, 95% CI 88% to 100%), Roth maximal count <12 (SP 100%, 95% CI 75% to 100%) and Roth counting time <8 s (SP 93%, 95% CI 66% to 100%) had high SP that could be used to rule in hypoxaemia, but displayed low SN (≤50%).

Conclusions Subjective dyspnoea measures have inadequate accuracy for ruling out hypoxaemia in high-risk patients with COVID-19. Safe home management of patients with COVID-19 should incorporate home oxygenation saturation monitoring.

  • COVID-19
  • infectious diseases
  • public health
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • This is the first study to evaluate the diagnostic accuracy of subjective dyspnoea in detecting hypoxaemia in the setting of COVID-19.

  • The diagnostic accuracy of patient-reported presence of dyspnoea in capturing hypoxaemia was evaluated across a range of oxygen saturation (SpO2) cutoffs from 92% to 97% and stratified based on age, presence of lung disease and date of symptom onset.

  • Objective dyspnea scales, including the modified Medical Research Council dyspnoea scale and Roth score were also evaluated.

  • Methodological limitations of the study include the retrospective study design and small sample size.

  • This study was limited to patients who were considered high risk for severe COVID-19, and the data collected for this study were from single patient assessments and did not assess whether changes in dyspnoea correlate with changes in SpO2 over time.

Introduction

As of 19 January 2021, there have been more than 93 million laboratory-confirmed novel COVID-19 cases and 2 million deaths documented globally.1 The spectrum of disease of COVID-19 ranges from asymptomatic or mild symptoms, to severe respiratory failure and death.2 Approximately 20% of patients with COVID-19 experience dyspnoea, which is more commonly associated with severe disease.3 Fatal cases of COVID-19 have higher rates of dyspnoea, lower blood oxygen saturation (SpO2) and greater rates of complications such as acute respiratory distress syndrome.4 5

In an effort to reduce avoidable hospitalisations, healthcare contacts and transmission, most patients with COVID-19 can be managed in the community using virtual healthcare platforms, and transferred to hospital only if they develop progressive respiratory disease.6 Subjective dyspnoea can be assessed remotely using patient interview, and augmented by surrogate measures such as the Roth score7 and modified Medical Research Council (mMRC) dyspnoea Scale.8 However, the accuracy of these measures has not been formally evaluated in the context of COVID-19.6 7 Of great concern is the risk of false reassurance if patients develop hypoxaemia without the subjective sensation of dyspnoea. ‘Silent hypoxaemia’, or low SpO2 in the absence of dyspnoea, has been reported in the setting of COVID-19 and clinicians have speculated that it may be associated with increased out-of-hospital mortality9; case reports have described patients presenting to hospital with rapid deterioration and respiratory failure without signs of respiratory distress.10–12

Previous studies of the utility of dyspnoea measurement in diagnosing hypoxaemia in other respiratory conditions such as chronic obstructive pulmonary disease (COPD), congestive heart failure and lung cancer have yielded conflicting results.13–16 The association has not been studied during the COVID-19 pandemic despite the highly publicised concept of ‘silent hypoxaemia’. Therefore, we sought to determine the diagnostic accuracy of subjective dyspnoea measures in diagnosing hypoxaemia among a cohort of outpatients with COVID-19.

Methods

Study participants

All consecutive patients with laboratory-confirmed COVID-19 followed as outpatients by the Sunnybrook Health Sciences Centre COVID-19 Expansion to Outpatients (COVIDEO) virtual care service from 3 April 2020 to 13 September 2020 were included in this retrospective cohort study. The patients were diagnosed based on a positive mid-turbinate or nasopharyngeal swab for COVID-19 RNA detected by real-time PCR. COVIDEO is a virtual care model for monitoring of outpatients with COVID-19 at Sunnybrook Health Sciences Centre, and is the basis of similar programmes at other hospitals.17 Patients were contacted by an infectious diseases physician for assessment and monitoring either by telephone or through the Ontario Telemedicine Network virtual platform.

A portable pulse oximeter was delivered to the homes of high-risk patients as defined by age >60 years, pregnancy, extensive comorbidities or presence of cardiorespiratory symptoms, such as chest pain or dyspnoea. The requirement for informed consent was waived.

Patient and public involvement

No patient involved.

Data collection

The demographic characteristics, clinical data, measures of subjective dyspnoea (presence of shortness of breath, mMRC dyspnoea scale score, Roth score), physical examination findings and SpO2 readings for study participants were collected from electronic medical records by one investigator (either AZ or SM). For the analysis, values were obtained from the patient’s initial virtual care assessment with a pulse oximeter, and subjective dyspnoea measures were taken at the same time as the objective measure of hypoxia.

Predictor variables

The primary predictor of interest was patient-reported presence of dyspnoea. Secondary predictor variables were patient-reported breathing faster at rest, breathing harder than normal, feeling more breathless today than yesterday, as well as dyspnoea as measured by the mMRC dyspnoea scale and the Roth score.

The mMRC dyspnoea scale has been studied extensively in a variety of respiratory conditions.8 It is composed of five categories that describe the degree of activity limitation due to worsening breathlessness. The participant assigns themselves a score ranging from 0 to 4 based on their perception of which activities result in dyspnoea, with higher scores indicating a greater impairment in their ability to perform daily activities.

The Roth score is a tool for quantifying the severity of dyspnoea, in which the patient is asked to count audibly to 30 in their native language, and the maximal count and counting time are recorded. A prior validation study demonstrated a strong positive correlation between pulse oximetry measurement and both counting time (r=0.59; p<0.001) and maximal count (r=0.67; p<0.001) achieved in one breath.7

Outcomes

The reference measure was SpO2 as measured by a ChoiceMMed pulse oximeter (model MD300C20). In the primary outcome definition, hypoxaemia was considered to be present if SpO2 was <95% in order to provide sufficient power to estimate the diagnostic test characteristics. Secondary outcomes included hypoxaemia cut-offs varying from 92% to 97%. Patients received instructions on correct pulse oximeter use and were told to wait 5–10 s for readings to calibrate prior to recording the SpO2 measurements.

Statistical analysis

In the primary analysis, the subjective dyspnoea measures were treated as diagnostic tests, and the specificity (SP), sensitivity (SN), positive predictive value (PPV), negative predictive value (NPV), positive likelihood ratio (+LR) and negative LR (−LR) were determined in order to evaluate the predictive value in detecting hypoxaemia. For the continuous predictors, the test characteristics were provided across a range of different score thresholds.

Diagnostic test characteristics of the primary dyspnoea measure were also determined in subgroups stratified based on patient characteristics, including (1) age <60 vs >60 years, (2) presence versus absence of underlying lung disease, and date from symptom onset (<7 vs >7 days). The Wilson method with continuity correction was used to calculate 95% CIs in order to avoid a negative lower limit.

A secondary analysis examined the strength of association between the presence of dyspnoea and hypoxaemia with a χ2 test or Fisher’s exact test (for sample sizes <5) with dyspnoea treated as a binary variable (present or absent). This relationship is represented by a violin plot. In additional analyses, a correlation coefficient was calculated to assess whether there was an association between the participants’ Roth scores and their SpO2 measurements. These associations were displayed graphically with a scatter plot (online supplemental figures 1 and 2). The relationship between the mMRC dyspnoea scale and hypoxaemia was analysed with a χ2 test and represented by violin plot.

All analyses were conducted in SAS Statistical Software V.9.3. For all statistical analyses, p<0.05 was considered significant.

Sample size calculation

The primary test characteristic of interest was the SN of dyspnoea as a test for hypoxaemia. The sample size was estimated based on a test of single proportion, namely SN. It was estimated that 62 patients would be required in order to estimate SN with a 10% margin of error and 95% CI if the true SN was 80%.

Results

Demographic and clinical characteristics of outpatients with COVID-19

From 3 April to 13 September 2020, a total of 501 patients with COVID-19 were followed by COVIDEO. Of these patients, 89 (17.8%) met criteria for home SpO2 monitoring (age >60 years, pregnancy, extensive comorbidities or presence of cardiorespiratory symptoms). One patient was lost to follow-up after provision of the oxygen monitoring device.

Overall, the median age of patients was 52 years (IQR 38–64 years) and 57 patients (64%) were female (table 1). The median number of days from symptom onset to clinical assessment was 6 (IQR 3–8). Among these patients, the most common comorbidities were hypertension (36%), obesity (20%), diabetes (17%), asthma (16%) and malignancy (16%). Twenty-nine patients (33%) had no comorbidities. The most common symptoms reported on intake assessment were fatigue/malaise (66%), cough (63%) and myalgias (45%). While the patients were being followed by COVIDEO, 11 (12%) required hospitalisation, with a median duration of hospitalisation of 3 days (IQR 2.5–7). Five (6%) patients were admitted to the intensive care unit (ICU), and the median length of ICU stay was 6 days (IQR 2–11). One patient was intubated, and no patients died within 30 days of their COVID-19 diagnosis.

Table 1

Demographics and clinical characteristics among outpatients with COVID-19 monitored with home oxygen saturation devices

Association of Dyspnoea measurements with detection of hypoxaemia

A total of 17 (19.1%) patients were diagnosed with hypoxaemia. hypoxaemia was significantly associated with the presence of dyspnoea (p=0.046), mMRC dyspnoea scale score over 0 (p=0.014), over 1 (p=0.001) and over 2 (p=0.001) (table 2). Weak associations were identified between patients’ Roth scores and their SpO2 measurements for maximum count (r=0.29; p=0.23) and counting time (r=0.12; p=0.617), respectively. The distribution of SpO2 (%) values in COVID-19 outpatients who reported dyspnoea and with various mMRC dyspnoea scale scores is shown in figure 1.

Figure 1

Comparison of SpO2 and measures of subjective dyspnoea. (A) Violin plots showing the distribution of SpO2 (%) values in COVID-19 outpatients who reported dyspnoea and those who did not. (B) Violin plots showing the distribution of SpO2 (%) values in COVID-19 outpatients with various mMRC dyspnoea scale scores. The width of each plot is proportional to the number of patients with the respective SpO2 (represented by black dots). The median SpO2 is indicated by the central horizontal black line and the dotted lines correspond to the IQR. mMRC, modified Medical Research Council; SpO2, oxygen saturation.

Table 2

Association of dyspnoea measurements with detection of hypoxaemia

Diagnostic accuracy of dyspnoea measurements in the detection of hypoxaemia

The presence or absence of subjective dyspnoea had an SN 47% (95% CI 24% to 72%), SP 80% (95% CI 68% to 88%), NPV 86% (95% CI 75% to 93%), PPV 36% (95% CI 18% to 59%), +LR 2.4 (95% CI 1.2 to 4.7), −LR 0.7 (95% CI 0.4 to 1.1) for diagnosing hypoxaemia (table 3). The presence of subjective dyspnoea had lower SN (25% (95% CI 16% to 37%]), 27% (95% CI 17% to 41%], 40% (95% CI 23% to 59%) for detecting hypoxaemia as defined by thresholds of <97%,<96% and <95%, respectively. At a lower SpO2 threshold of <92%, the SN increased only slightly to 50% (95% CI 19% to 81%). The other binary measures of subjective dyspnoea, including breathing faster at rest, breathing harder than normal, and feeling more breathless than the day before had lower SN (0% (95% CI 0% to 24%], 0% (95% CI 0% to 24%) and 6.2% (95% CI 0.3% to 32%), respectively), and higher SP (96% (95% CI 87% to 99%], 97% (95% CI 89% to 100%) and 96% (95% CI 87% to 99%), respectively) (table 3).

Table 3

Diagnostic accuracy of dyspnoea measurements in the detection of hypoxaemia: estimation of test characteristics including sensitivity (SN), specificity (SP), negative (NPV) and positive (PPV) predictive value and negative (−LR) and positive (+LR)

mMRC dyspnoea scale scores were recorded and available for 63 patients (70.8%). An mMRC dyspnoea scale score of greater than 0 was determined to have an SN of 54% (95% CI 26% to 80%), SP 82% (95% CI 68% to 91%), NPV 87% (95% CI 74% to 95%), PPV 44% (95% CI 21% to 69%), −LR 0.6 (95% CI 0.3 to 1.0) and +LR 3.0 (95% CI 1.4 to 6.5) for the detection of hypoxaemia. At higher cut-off values, the SN of the mMRC dyspnoea scale was reduced to 39% (95% CI 15% to 68%) for scores greater than 1 and 2% and 8% (95% CI 0.4% to 38%) for scores greater than 3. The SP for mMRC dyspnoea scale scores greater than 1, 2 and 3 at capturing hypoxaemia was 98% (95% CI 88% to 100%).

Roth scores were available for 19 patients (29.7%). The Roth score had a higher SN at higher cut-off values for counting time. A maximum count of less than 12 had an SN of 25% (95% CI 1.3% to 78%), SP 100% (95% CI 75% to 100%), NPV 83% (95% CI 58% to 96%), PPV 100% (95% CI 6% to 100%) and a −LR of 0.75 (95% CI 0.4 to 1.3).

The diagnostic test with the highest SN for diagnosing hypoxaemia was a Roth score maximum counting time of less than 25 s, which still had an SN of only 75% (95% CI 22% to 99%), and inadequate SP 13% (95% CI 2.3% to 42%), NPV 67% (95% CI 13% to 98%), PPV 19% (95% CI 5.0% to 46%), −LR 1.88 (95% CI 0.2 to 16) and +LR 0.87 (95% CI 0.48 to 1.6). When all subjective dyspnoea predictors are combined in a single variable, the SN is 59% (95% CI 34% to 81%), SP 67% (95% CI 55% to 78%), NPV 87% (95% CI 75% to 94%), PPV 30% (95% CI 16% to 49%), −LR 0.6 (95% CI 0.3 to 1.1) and +LR 1.8 (95% CI 1.0 to 3.0).

The diagnostic accuracy of dyspnoea presence in the detection of hypoxaemia was most impacted when stratified by the presence of underlying lung disease. In the patients with underlying lung disease, the SN and SP of the presence of dyspnoea in detecting hypoxaemia was 100% (95% CI 20%– to 100%) and 80% (95% CI 51%– to 95%), respectively. A lower SN (22% (95% CI 3.9% to 60%)) and high SP (96% (95% CI 79% to 100%)) was observed for patients over 60 years when results were stratified based on age. Stratifying based on days from symptom onset did not impact the diagnostic accuracy of dyspnoea in detecting hypoxaemia (table 4).

Table 4

Diagnostic accuracy of the presence of dyspnoea in the detection of hypoxaemia stratified based on patient characteristics, including age, presence versus absence of underlying lung disease and date from symptom onset

Discussion

To our knowledge, this is the first study to evaluate the diagnostic accuracy of subjective dyspnoea in detecting hypoxaemia in the setting of COVID-19. Self-reported shortness of breath has very limited utility for detecting hypoxaemia, with an SN of only 47% and SP of only 80% for detecting SpO2 levels below 95%. Using a lower SpO2 threshold of less than 93% did not meaningfully improve the SN of subjective dyspnoea in diagnosing hypoxaemia (SN 50%). Other binary measures of subjective dyspnoea, including breathing faster at rest, breathing harder than normal and feeling more breathless than yesterday offered high SP. Similarly, an mMRC dyspnoea score exceeding 1, a Roth maximal count less than 12 and Roth counting time under 8 s offered high SP and +LR to rule in hypoxaemia. Identifying patients with these features may be helpful in the remote assessment of COVID-19 outpatients. However, none of these measures offered sufficient SN or –LR to help rule out hypoxaemia—which is the more clinically important consideration for these patients. Even when all variables were combined into a single maximally sensitive predictor, the SN was just 59%.

Previous studies examining the correlation between subjective dyspnoea and hypoxaemia in other respiratory conditions have yielded inconsistent findings. The strongest confirmation of the potential diagnostic utility of dyspnoea emerged from a study of 76 patients admitted to the emergency department with acute exacerbations of COPD, in which dyspnoea scores exceeding 3 or 4 on a five-stage scoring system were found to have a sensitivity of 93.5% for detecting hypoxaemia.14 Additionally, the mMRC dyspnoea scale has been found to be significantly correlated with SpO2 in measurements obtained during exercise among patients with idiopathic pulmonary fibrosis.18 Conversely, several other studies have shown no correlation between perceived dyspnoea and hypoxaemia in conditions such as advanced lung cancer, COPD and palliative care patients.15 16 19 While previous studies show variable relationships between dyspnoea and hypoxaemia in various respiratory pathologies, our study shows that neither binary measures of subjective dyspnoea, the mMRC dyspnoea scale, or the Roth score can be used to rule out hypoxaemia in the setting of COVID-19.

Discrepancies between respiratory rate and SpO2 in COVID-19 patients with acute respiratory failure have been highlighted previously, suggesting that a normal respiratory rate may belie profound hypoxaemia in this setting.20 High levels of anxiety may contribute to feelings of dyspnoea in patients who are non-hypoxaemic. There are also a growing number of case reports documenting silent hypoxaemia among COVID-19 patients, where patients present with hypoxaemia in the absence of respiratory symptoms.10 11 21 The underlying mechanism responsible for severe hypoxaemia in the absence of dyspnoea is not well elucidated. It has been postulated that this clinical picture may be consistent with a phenotype of COVID-19 pneumonia (L-phenotype) characterised by low elastance, low ventilation-perfusion ratio and near normal compliance.12 The relatively high compliance results in preserved gas volumes, while hypoxaemia may result due to a ventilation–perfusion mismatch caused by impaired lung perfusion regulation and loss of hypoxic vasoconstriction.22 23 Additionally, the absence of dyspnoea despite severe hypoxaemia may reflect pulmonary vaso-occlusive disease, whereby patients develop clinically silent microvascular thrombi in early stages of the disease, which if left untreated, results in worsening hypercoagulability and rapid clinical deterioration due to a thromboinflammatory cascade.24 25 While at this point the exact mechanism remains speculative, our data suggest that the discrepancy between dyspnoea and hypoxaemia makes it difficult to accurately assess patients remotely and emphasises the importance of SpO2 monitoring in order to avoid missing patients with developing respiratory failure.

This study has several limitations. The data collected for this study were from patients’ initial pulse oximeter assessment and did not assess whether changes in dyspnoea correlate with changes in SpO2 over time. This is a potentially important notion when monitoring patients who are (or are not) becoming increasingly dyspneic while self-isolating in their homes. While the number of patients included was sufficient for the primary analysis, they were insufficient for precise estimates of subgroups stratified by age, presence of lung disease, date of symptom onset and for calculation of the diagnostic test characteristics at lower SpO2 cutoffs. In our study, less than 20% of included patients were diagnosed with hypoxaemia. While this represents a small sample of patients with hypoxaemia, it is clear that in order to prevent missing any patients with hypoxaemia who require admission all high risk patients require SpO2 monitoring.

Additionally, our study was limited to patients who were considered at high risk of severe disease, and it is possible that the diagnostic test characteristics might differ in younger and healthier patients. Lastly, pulse oximeters may have variable accuracy as individuals become increasingly hypoxic and are further impacted by individual patient characteristics; however, a perfect reference standard of invasive blood oxygen measurement would be neither practical nor ethical in the outpatient setting.26

Conclusions

Our findings indicate that subjective dyspnoea does not accurately capture hypoxaemia in patients with COVID-19. Although some dyspnoea scores have high specificity and +LR for identifying hypoxaemia, none of these measures have sufficient sensitivity to rule out hypoxaemia. Therefore, relying on surrogate measures of dyspnoea alone is not sufficient to remotely monitor high-risk outpatients with COVID-19. Home SpO2 monitoring should be a mandatory component of remote management all high-risk outpatients with COVID-19.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors The authors all stand behind the conclusions of this manuscript, agree to be accountable for all aspects of the work, and support its publication. LB contributed to the planning, conception and study design, data analysis, interpretation of the data, and reporting of the work. AZ contributed to data analysis, interpretation of the data, and reporting of the work. NA contributed to the planning, conception, study design, study conduct, data acquisition and reporting of the work. AKC contributed to the planning, conception, study design, study conduct, data acquisition and reporting of the work. JE-C contributed to data analysis, interpretation of the data and reporting of the work. AG contributed to the interpretation of the data and reporting of the work. PWL contributed to the planning, conception, study design, study conduct, data acquisition and reporting of the work. JAL contributed to the planning, conception, study design, study conduct, data acquisition and reporting of the work. SMP contributed to the planning, conception, study design, study conduct, data acquisition and reporting of the work. SM contributed to the planning, conception, study design, study conduct, data acquisition and reporting of the work. AES contributed to the planning, conception, study design, study conduct, data acquisition and reporting of the work. ND contributed to the planning, conception, and study design, study conduct, data acquisition, data analysis, interpretation of the data and reporting of the work. All authors contributed to the manuscript preparation and have given approval for its submission.

  • Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

  • Patient consent for publication Not required.

  • Ethics approval This study was approved by the institutional review board at Sunnybrook Health Sciences as minimal risk research, using data collected for routine clinical practice.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data availability statement Data are available on reasonable request. All data relevant to the study are included in the article or uploaded as online supplemental information. Deidentified data are available on reasonable request.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.