Our prospective observational pilot study was conducted with approval from the Hospital Clinic of Barcelona IRB (Ethics Committee's approval/ID: HCB 2020/0666). Written informed consent was obtained from participating patients. The study was performed in accordance with relevant guidelines and regulations.
Subject selection
We enrolled adults of any race and ethnicity who were admitted to a routine nursing ward with a diagnosis of SARS-CoV-2 infection, confirmed by a reverse-transcription polymerase chain reaction test. We excluded patients younger than 18 years old and pregnant women. Because we targeted patients at risk of clinical deterioration who were eligible for non-invasive and invasive mechanical ventilation support, we also excluded patients in whom escalation to life-sustaining respiratory support was not anticipated because treatment was unlikely to prove helpful.
Protocol
We included patients who had COVID-19 and were hypoxemic, defined as arterial pressure of oxygen (PaO2) < 80 mmHg or pulse oximeter saturation (SpO2) < 90% while breathing ambient air. All were hospitalized for treatment and required supplemental oxygen, to ensure a SpO2 over 91%.
Patients were considered to be in respiratory failure and candidates for escalation to non-invasive ventilation when they developed moderate-to-severe shortness-of-breath and tachypnoea of ≥ 30 bpm; when the PaO2/Inspired fraction of oxygen (FiO2) ratio was ≤ 200; if a FiO2 > 0.4 was needed to maintain a SpO2 > 92%; or when pH was < 7.35 accompanied by a PaCO2 exceeding 45 mmHg.
High-flow nasal cannula oxygen (HFNCO) was generally the initial treatment, followed by Non-Invasive Mechanical Ventilation (NIMV). Patients were escalated to invasive mechanical ventilation if they did not improve and had a ROX index [(SpO2/FiO2) / respiratory rate] < 3 within 2 hours, < 3.5 within 6 hours, or < 4 within 12 hours after HFNCO and NIVM initiation.
Measurements
Demographic, morphometric, and clinical characteristics, as well as nurse-recorded vital signs including blood pressure, heart rate, respiratory rate, oxygen saturation, ventilation mode, supplemental oxygen concentration and delivery system were recorded.
Pulmonary management was recorded, including use of HFNCO, NIMV, invasive mechanical ventilation, and clinically detected deterioration. We also recorded episodes when nurses requested non-routine physician or ICU team evaluation, transfer to a higher-acuity unit, and discharge disposition including vital status. Arterial blood gas analysis was not mandated per protocol, but available values were extracted from medical records.
The Circadia Contactless Breathing Monitor (model C100, Circadia Technologies, Ltd., London, United Kingdom) is an FDA-cleared non-contact device that uses radar to measure patient respiratory rate. Accuracy of the Circadia System has been validated by direct comparison against manually scored reference end-tidal CO2 capnography and polygraphy ventilatory effort, both during spot measurements and continuous monitoring in awake and sleeping patients.[18]
We began Circadia monitoring as soon as practical after patients consented to participate in our study. Circadia Monitors were positioned about 1.5 meters from patients. Monitoring continued for 15 days, unless aborted for hospital discharge, initiation of invasive mechanical ventilation, or death. Respiratory rate data from the Circadia Monitor were transmitted wirelessly to cloud-based storage for offline analysis and was not available to clinicians who managed participating patients. Respiratory rate from the Circadia Monitor was recorded at 3-second intervals while patients were within the detection range of the Circadia Monitor (typically while in bed).
Exposure and outcomes
Our primary outcome was pulmonary decompensation, defined as an escalation in the pulmonary care beyond Venturi-mask support. The following levels of respiratory support were defined, from less to more severe: 1) ambient air; 2) nasal cannula; 3) Venti-mask; 4) HFNCO; 5) NIMV; 6) invasive mechanical ventilation; and 7) death. Our primary aim was to assess whether the information extracted from respiratory rate continuous monitoring is associated with respiratory care escalation 24 hours in advance in hospitalized COVID-19 patients.
The exposure was short-term and long-term respiratory patterns associated with respiratory care escalation. Respiratory patterns (variability and fluctuations over time) were captured through features, which included respiratory rate mean, respiratory rate standard deviation, kurtosis of the respiratory rate, skewness of the respiratory rate, and trend (extracted from the slope of a linear regression model), all computed over three rolling periods, 30 min, 24 hours and 72 hours in duration, occurring 24 hours prior to respiratory care escalation (Supplemental Table 1).
Data analysis
Analysis was restricted to subjects who had successful Circadia monitoring for at least half of their hospitalization. We considered patients to have escalated respiratory care when they required an increment in their ventilatory support category, with others falling into the non-escalating category. Patients could experience multiple escalation episodes within a day. However, we considered only the initial escalation within each 48-hour period to avoid situations where a prior de-escalation of treatment was quickly unsuccessful, rather than representing a true new worsening of the underlying pulmonary function. We also excluded patients who de-escalated ventilatory support after joining the study. We considered a patient to be de-escalating if they required HFNCO or NIMV at study inclusion and thereafter progressively improved (Supplemental Fig. 1).
Patients who required ventilatory support escalation were matched to three non-escalating patients on baseline variables (Table 1) using minimum Euclidean distance. Escalation episodes in each patient were matched to a similar point in the three matched non-escalation patients based on time since admission. When the exact corresponding time was unavailable because the non-escalation patient was already discharged at this time point, the last available time window was used. A non-escalated patient could be matched with more than 1 escalated patient (Supplemental Table 2).
Table 1
Patient demographic, morphometric, and clinical characteristics
| Before Matching | After Matching |
Escalation | | Escalation | |
| YES (N = 13) | NO (N = 112) | P-value | YES (N = 13) | NO (N = 26) | P-value |
Age, years | 67 ± 13 | 61 ± 13 | 0.06 | 67 ± 17 | 63 ± 12 | 0.322 |
Female, n (%) | 8 (32) | 50 (36) | 0.049 | 3 (23) | 6 (23) | 1 |
BMI, kg/m2 | 31 ± 6 | 28 ± 5 | < 0.01 | 31 ± 5 | 29 ± 4 | 0.111 |
Medical History | | | | | | |
Arterial hypertension, n (%) | 16 (10) | 68 (48) | 0.058 | 8 (61) | 19 (73) | 0.713 |
Diabetes mellitus, n (%) | 10 (6) | 22 (16) | 0.021 | 3 (23) | 6 (23) | 1 |
Ischemic cardiac disease, n (%) | 2 (1) | 9 (6) | 0.487 | 1 (8) | 2 (8) | 1 |
Chronic heart failure, n (%) | 2 (1) | 2 (1) | 0.318 | 2 (15) | 1 (4) | 0.399 |
COPD, n (%) | 2 (1) | 11 (8) | 0.275 | 3(2) | 11(23) | 1 |
Asthma, n (%) | 1 (1) | 11 (8) | 0.206 | 1 (8) | 1(4) | 1 |
Neurological disease, n (%) | 1 (1) | 11 (8) | 0.588 | 1 (8) | 1(4) | 1 |
Variables were summarized as means ± SDs and N (%). BMI, body mass index; COPD, Chronic Obstructive Pulmonary Disease. T-tests were used for comparisons. P < 0.05 was considered significant |
Respiratory rate was sampled at 3-second intervals. Using a 30-minute rolling analysis window, 30-minute respiratory features were obtained through various aggregation methods, including the mean. The mean respiratory rate, extracted over 30 minutes, was further aggregated into windows of 24 and 72 hours, at 30-minute increments, to obtain the 24-hour and 72-hour features. Extracted respiratory features included mean, standard deviation, minimum, maximum, skewness, kurtosis, and trend.
Since some patients had multiple escalation events, and to adjust for confounders, we estimated the standard error and confidence intervals using bootstrap resampling of the logistic regression model (10,000 runs, sampling episode ID with replacement and using the percentile method). We also used the variance inflation factor (VIF) to estimate collinearity between respiratory features. Higher values of VIF (above 10) indicate the existence of multicollinearity, meaning that such variables are highly correlated and should not simultaneously be included in the model.[19] The model was fitted using the glm package in R (version 3.6.0).
We used a multivariable logistic regression model to assess the degree to which various respiratory rate features were associated with respiratory care escalation occurring 24 h later. The criterion for a predictor to be significant was a p-value < 0.05 in the multivariable logistic regression model, uncorrected for multiple comparisons. Performance of the model was measured using traditional measures of diagnostic accuracy, specifically c-statistic from the logistic regression model (area under the receiver operating characteristic curve (ROC-AUC)), sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV). We assumed a decision threshold > 0.5 from the logistic regression model predicted escalation.
Sample size considerations
There was no a priori sample size estimate for our pilot study. Instead, we enrolled all qualifying and consenting patients over the study period. Given the observed confidence intervals, our actual sample size appears sufficient to evaluate early prediction of respiratory escalation.