Next Article in Journal
A New Electro-Optical Switch Modulator Based on the Surface Plasmon Polaritons of Graphene in Mid-Infrared Band
Next Article in Special Issue
Wearable Sensor-Based Exercise Biofeedback for Orthopaedic Rehabilitation: A Mixed Methods User Evaluation of a Prototype System
Previous Article in Journal
Imaging for Small UAV-Borne FMCW SAR
Previous Article in Special Issue
Data Analytics of a Wearable Device for Heat Stroke Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Breathing Parameters Using an Inertial Measurement Unit (IMU)-Based System

1
Dipartimento di Elettronica, Informazione e Bioingegneria, Politecnico di Milano, 20133 Milan, Italy
2
Scientific Institute, IRCCS E. Medea, Bioengineering Lab, 23842 Bosisio Parini, Lecco, Italy
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(1), 88; https://doi.org/10.3390/s19010088
Submission received: 24 October 2018 / Revised: 10 December 2018 / Accepted: 20 December 2018 / Published: 27 December 2018
(This article belongs to the Special Issue Data Analytics and Applications of the Wearable Sensors in Healthcare)

Abstract

:
Breathing frequency (fB) is an important vital sign that—if appropriately monitored—may help to predict clinical adverse events. Inertial sensors open the door to the development of low-cost, wearable, and easy-to-use breathing-monitoring systems. The present paper proposes a new posture-independent processing algorithm for breath-by-breath extraction of breathing temporal parameters from chest-wall inclination change signals measured using inertial measurement units. An important step of the processing algorithm is dimension reduction (DR) that allows the extraction of a single respiratory signal starting from 4-component quaternion data. Three different DR methods are proposed and compared in terms of accuracy of breathing temporal parameter estimation, in a group of healthy subjects, considering different breathing patterns and different postures; optoelectronic plethysmography was used as reference system. In this study, we found that the method based on PCA-fusion of the four quaternion components provided the best fB estimation performance in terms of mean absolute errors (<2 breaths/min), correlation (r > 0.963) and Bland–Altman Analysis, outperforming the other two methods, based on the selection of a single quaternion component, identified on the basis of spectral analysis; particularly, in supine position, results provided by PCA-based method were even better than those obtained with the ideal quaternion component, determined a posteriori as the one providing the minimum estimation error. The proposed algorithm and system were able to successfully reconstruct the respiration-induced movement, and to accurately determine the respiratory rate in an automatic, position-independent manner.

1. Introduction

Continuous monitoring of respiratory parameters such as breathing frequency (fB), inspiratory time (TI) and expiratory time (TE) could foster early diagnosis of a wide range of respiratory disorders and help to track a patient’s condition, discriminating between stable and at-risk patients [1,2]. Conditions of interest could be sleep breathing disorders, sudden infant death syndrome, chronic obstructive pulmonary disease (COPD) and neuromuscular disorders. The current gold standard for measuring fB is to count the number of breaths in one minute, through auscultation or observation [3,4]. Other methods for breathing function assessment currently used in clinical practice are spirometry or pneumotachograph based on airflow measurement by using mouthpiece or facemask. In overnight polysomnography, breathing activity is assessed both by measuring respiratory flow, through pressure transducer or thermistors near the nostrils, and respiratory efforts (breathing-derived chest-wall movements), by strain-gauge belts. Also, exhaled carbon dioxide sensors, transthoracic inductance and impedance plethysmography and ECG—or PPG—derived fB have been used to measure breathing signal. Despite their accuracy, these methods are uncomfortable and intrusive, and are not suitable for continuous monitoring in the clinical environment and at home. An emerging area of interest is to use motion sensors to detect the small breathing-derived movements/orientation changes of the chest wall. This method is particularly suitable for long-term breathing monitoring because it is unobtrusive, tolerable, and low-cost. The principle was first presented with a single-axis accelerometer in animal model (dog) using a pressure transducer in the trachea as reference [5]. Starting from this point, a variety of studies demonstrated the feasibility of using one accelerometer placed on the chest wall to derive respiratory signal and/or breathing frequency in different positions [6,7,8,9,10,11,12,13,14]. Morillo et al. [8] combined a piezoelectric single-axis accelerometer and a polarized capacitive microphone placed on the suprasternal notch to collect information of the cardiac, respiratory, and snoring activities for the screening of patients affected by Sleep Apnea-Hypopnea Syndrome. Measurements were limited to the supine position, that was selected to increase the sensitivity of the single-axis accelerometer, limiting the generality of the findings. The analysis method was based on the estimation of breathing frequency through the identification of the peak of the spectrum or autocorrelation; the main limitation of this approach is that, when the breathing is irregular, a main peak may not exist, and individual breaths must be identified and counted. Hung et al. [7] moved from single-axis to biaxial accelerometers. The aim of their study was to evaluate the reliability of the device in terms of detection of the onsets of expiration and expiration, and to assess the feasibility of differentiating between different breathing patterns (normal breathing, apnea, deep breathing). The signals from both axes (anteroposterior and longitudinal) of the accelerometer were summed, limiting the analysis to the sagittal plane, in sitting and lying positions. An adaptive band-pass filter was applied with a variable passband centered at the detected dominant breathing frequency.
As emerged by these studies, single or dual-axis accelerometers can be used to derive breathing signal when appropriately aligned with the major axis of rotation, which changes when the subject move from a posture to another. Contrarily, the use of a tri-axial accelerometer allows measuring inclination changes due to breathing regardless of orientation. In this case, the problem lies in the identification of the accelerometer axis to consider when posture changes. Bates et al. [13] proposed a method to track the major axis of rotation as it changes, to continuously monitor angular motion due to breathing also when subject change position/orientation. An alternative possibility to the best axis selection is fusing the axes. Jin et al. [12] proposed a posture-independent signal processing method based on three possible algorithms for accelerometer axes fusion. They demonstrated that methods based on Principal Component Analysis (PCA) obtained the highest performance in terms of Signal-to-Noise Ratio (SNR), but no results were provided about breathing rate estimation or validation against a reference method.
With the entry of tri-axial accelerometers new opportunities opened for the monitoring of breathing frequency using inertial sensors, but their use was still confined to static conditions since, when the subject is moving, the degree of the movement-related signal would exceed that due to breathing. One possible approach is to identify non-breathing motion, as proposed by Bates et al [13]. In a successive study, Mann et al. [11] furtherly developed the method proposed by Bates et al. [13], by adding activity tracking, and allowing identification of asymmetric breaths, that was not possible in the original method. An attempt to remove motion artifacts by using signal processing was made by Liu et al. [14]. They proposed an elegant method based on PCA-fusion of the three axes of an accelerometer and on filtering of the first principal component by using an adaptive filter that varied according to the energy expenditure derived by the same accelerometer. To overcome the problems of using a single accelerometer in dynamic conditions, a possibility is to fuse data from accelerometers and from other sensors, such as gyroscopes. Yoon et al. [15] investigated the feasibility of measuring breathing-related motions also during dynamic activities of the subject, by fusing data from a tri-axial accelerometer and a gyroscope and applying Kalman filter. They found that, during dynamic exercises, fusion of accelerometer and gyroscope data provided benefits in terms of reduction of estimation error. Gollee et al. used a more complex system, an inertial measurement unit (IMU) fusing accelerometer, gyroscope and magnetometer but considered only static conditions [16]. Another approach to overcome the problems related to motion artefacts is modularity. Lapi et al. [17] tried to overcome limitations deriving from the use of a single accelerometer by proposing a system based on a couple of 3-axis accelerometers placed bilaterally on the chest. Using two accelerometers permitted to detect respiration-related chest-wall movements regardless of sensor positioning with respect to the gravity vector; secondly, the breathing frequency can be obtained even when one of the two sensors is silenced by postural constraints. Recently, Gaidhani at al. [18] proposed a method that uses two IMUs composed by a 3-axis accelerometer, a 3-axis gyroscope and a 3-axis magnetometer, placed on the anterior and posterior side of the chest to decompose the motions experienced by the two IMUs into trunk movements and breathing actions. This paper presents an automatic processing algorithm to derive breathing frequency and other breathing temporal parameters from quaternion-based orientation signals recorded simultaneously at thoracic and abdominal level by using a modular, wireless, IMU-based device [19]. An important step of the processing algorithm is dimension reduction (DR) that allows the extraction of a single respiratory signal starting from 4-component quaternion data. Three different methods of DR are proposed and compared; two of them are based on the selection of one quaternion component, the third one is based on PCA-fusion of the 4 quaternion components. Results obtained using the IMU-based device, with the three different methods, are validated against optoelectronic plethysmography, an already established method to evaluate ventilation through an external measurement of the chest-wall surface motion [20,21,22,23,24,25].

2. Materials and Methods

2.1. Device Architecture and Hardware Description

The system used in this study is composed by three IMU-sensor units that communicate via Bluetooth with a smartphone; here data are pre-processed and saved. Two of the three sensor units (peripheral units) are dedicated to the recording of chest-wall respiratory-related movements and are placed on the thorax and on the abdomen to record respiratory information about both the compartments; the third sensor unit (reference central unit) is placed on a body area that is integral with the chest wall, but not involved in respiratory movements (e.g., coccyx or anterior superior iliac crest). The measurement of chest-wall movements, related to both abdominal and thoracic compartments, allows the consideration of the two-degree-of-freedom (DoF) model of chest-wall breathing movements [26], that considers abdomen and rib cage (thorax) as acting independently. Moreover, the compartmental contribution to total chest-wall volume changes according to posture and the breathing strategy adopted by each subject. Thus, the recording of chest-wall movements at different levels provides on the one hand, a more accurate estimation of the breathing signal, and on the other hand allows investigation of asynchronies between compartments, typical of different pathological conditions. The reference central unit, in addition to performing a central role within the Bluetooth piconet, can be used to discriminate between static and dynamic conditions and to map the activity state of the subject. Moreover, it could be used to reduce movement information not linked to breathing by means of frequency domain analysis or by referring orientation change experienced by the peripheral units to the coordinate frame of the reference unit. Each unit is composed by a printed circuit board equipped with a low-power microcontroller, a Bluetooth Low Energy (BLE) module, a 9-DoF IMU (3-axis accelerometer, a 3-axis gyroscope, and a 3-axis magnetometer) and lithium polymer rechargeable battery. A voltage regulator circuit, and Li–Po battery recharge circuit with mini USB port are also included in the design. Differently from the peripheral units, the reference central unit is equipped with a different BLE module, able to support simultaneous central/peripheral role and also brings a Micro Secure Digital (SD) Memory Card Connector for data logging. The dimensions of each peripheral unit, comprehensive of the 3D-printed housing, are 41 mm × 33 mm× 19 mm (LWH), and the weight is 25 g, including the battery, while the reference central unit measures 45 mm × 45 mm × 15 mm (LWH), and weighs 35 g. A prototypal version of this device has been described in [19].

2.2. Quaternion-Based Orientation Estimation and Fusion Algorithm

The final goal is to derive breathing signal by measuring orientation changes during the respiratory movements, both at thoracic and abdominal level. The IMUs provide 3D-acceleration, 3D-magnetic field, and 3D-angular rate. These measures are combined to provide accurate 3D orientation data aboard each unit. The orientation is represented with quaternions [27,28], that even though may suffer from problems of interpretation in terms of meaningfully physical angles, are interesting mathematical entities (four-dimensional complex number (q = [q0 q1 q2 q3]), since they require less computing time and avoid the singularity problems (i.e., “gimbal lock”) typical of other orientation descriptors, e.g., Euler angles. The fusion of the data collected from the sensors is done by using the sensor fusion algorithm proposed by Madgwick et al. [29], based on an analytically derived and optimized gradient descent algorithm enabling levels of accuracy exceeding that of the Kalman-based algorithm, with low computational (277 scalar arithmetic operations each filter update) load and low sampling rates (e.g., 10 Hz); this orientation filter also provides an online magnetic distortion compensation algorithm and gyroscope bias drift compensation. The sensors data were collected at 40 Hz and the fusion algorithm was updated with the same rate, but due to limited buffer of the BLE module and to the stricter timings used for the Bluetooth communication, just one quaternion out of 4 computed is considered (10 Hz); nevertheless, the final sampling rate was considered appropriate given the relative low frequency of the respiratory signal [0.1 ÷ 1 Hz]. Thus, the microprocessor of each unit, receives data from accelerometer, gyroscope and magnetometer that are on board and implements Madgwick fusion filter [29] to compute a quaternion representing the change of orientation of each unit relative to the earth frame ( q ^ E a r t h T h ,   q ^ E a r t h A b ,   q ^ E a r t h R e f ), or more correctly the change of orientation of the earth relative to each unit frame [29]. In fact, in quaternion form, an arbitrary orientation of a coordinate frame B relative to coordinate frame A, achieved through a rotation of angle θ around an axis Ar (rx, ry, rz) defined in frame A, is univocally represented through the normalized quaternion q ^ B A defined by Equation (1):
q ^ B A = [ q 0   q 1   q 2   q 3 ] = [ cos θ 2 r x   sin θ 2 r y   sin θ 2   r z   sin θ 2 ]

2.3. Quaternion-Derived Breathing Frequency

All the elaborations and computations needed to extract breathing parameters from data collected by the device were performed offline using MATLAB, the processing took on average 1.027 ± 0.129 seconds for the analysis of signals of 1071 ± 270 samples. A signal processing procedure was designed to extract the breathing frequency starting from quaternions representing the change of orientation of each unit relative to the earth frame ( q ^ E a r t h T h ,   q ^ E a r t h A b ,   q ^ E a r t h R e f ). The block diagram of the signal processing part is presented in Figure 1. The algorithm is divided into 4 main blocks: (i) pre-processing, (ii) DR, (iii) spectrum analysis, and (iv) processing.
Pre-processing block includes the preliminary steps that leads to chest-wall respiratory-related orientation change signals. The orientations changes of thoracic and abdominal units were referred to the reference unit frame (that in turn represents orientation changes of trunk) applying Equations (2) and (3) respectively:
q ^ R e f T h = q ^ E a r t h T h     q ^ * E a r t h R e f   = q ^ E a r t h T h     q ^ R e f E a r t h ,
q ^ R e f A b = q ^ E a r t h A b     q ^ * E a r t h R e f   = q ^ E a r t h A b     q ^ R e f E a r t h ,
These two quaternions represent the outputs of the pre-processing block and the input of the DR block.
Dimension-reduction block takes the quaternions obtained from Equations (2) and (3), that are composed by 4 components each, and provides as output 2 single-component signals (1 for the abdomen and 1 for the thorax) representing chest-wall respiratory-related orientation change signals. These two signals represent the input of the power spectrum block and of the processing block. To reduce dimension from 4 components to 1, two possibilities were investigated as shown in Figure 2:
(i).
Best quaternion component selection
(ii).
PCA-based fusion of the quaternion components
To select the best component among the 4 components representing the orientation quaternion, two different methods were proposed, both based on spectrum analysis. The idea was to choose the component with the highest breathing information, computing the power spectral density estimate (PSD) between 0.5–2 Hz for each component and selecting the component with: (1) maximum PSD peak (“Peak” method) or (2) maximum area under the PSD (“Area” method). To assess the goodness of these two methods in predicting the best quaternion component, the ideal component (“Ideal”) was determined a posteriori, case by case, based on minimum breathing frequency estimation error (see Section 2.5).
Since more than one quaternion component is supposed to convey breathing information, the possibility to maximize this information fusing the 4 components of the quaternion by means of PCA was investigated. PCA is a mathematical procedure that transforms an original set of correlated variables into a (smaller) number of uncorrelated variables by determining a set of orthogonal vectors called principal components, which are defined by a linear combination of the original variables [30,31]. To do this, the directions in the data with the most variation, i.e., the eigenvectors corresponding to the largest eigenvalues of the covariance matrix, are computed and the data are projected onto these directions. To compute the eigenvectors, data were arranged into a two-dimensional matrix X(m × n), where m was the number of observations of the time series and n the number of variables (quaternion components). Then, the univariate means were subtracted from the n columns, to center the data. Singular Value Decomposition (SVD) was used to compute the eigenvectors (V = [v1, v2, v3, v4]) and corresponding eigenvalues (λ1, λ2, λ3, λ4). Original data were finally projected in the new coordinate system (Y = XV) and the first principal component, accounting for the largest possible variance, was selected and passed to other blocks [30,31].
Spectrum Analysis block include a set of steps needed to optimize the subsequent processing phase. The two signals representing chest-wall (abdominal and thoracic) respiratory-related orientation obtained downstream of the dimension-reduction block underwent the following steps (Figure 1):
(i).
A low-frequency threshold (fLOW) was determined based on a first estimate of the breathing frequency (fB). The rough estimate of fB was done by identifying maxima points of the signal and computing the fB, breath by breath, as reciprocal of the temporal distance between consecutive maxima points. Then, the mean (fB_Rough) and the standard deviation (fB_Rough_SD) of the fB over the entire trial were computed. To facilitate maxima points identification, signals were at first band-bass filtered using a first-order infinite impulse response (IIR) Butterworth filter [0.05 Hz–2 Hz] and smoothed with a third-order Savitzky–Golay [32] finite impulse response (FIR) filter (fixed window length = 31 samples). Low thresholds fLOWAb and fLOWTh were determined for the abdominal and thoracic signals respectively as difference fB_Rough − fB_Rough_SD. Then the minimum value between fLOWAb and fLOWTh was chosen as final low-frequency threshold, named fLOW, and it was used in the next step.
(ii).
PSD estimate (Welch’s method, Hamming window size: 300 samples, overlapping: 50 samples) was computed and the spectrum frequency corresponding to the breathing rate was identified, both for the thorax (fpeak_T) and the abdomen (fpeak_A), by looking for the local peak of the PSD within the window [fLOW ÷ 2 Hz]. The use of a low threshold, based on a rough estimate of the breathing frequency, supports the selection of the PSD peak linked to breathing rate and avoid selecting wrong peaks, often related to low-frequency oscillation artifacts.
(iii).
The breathing frequency derived by the spectrum was used to set an adaptive band-pass filter, as proposed in a previous study [7], centered on fpeak frequency. For the abdomen, upper (fU) and lower (fL) cut-off frequency points for the band-pass filter were defined, by applying Equations (4) and (5) respectively [7]:
fU_A = fpeak_A + 0.04,
fL_A = max (0.05, (fpeak_A − 0.04)),
For the thorax, Equations (6) and (7) were applied:
fU_T = fpeak_T + 0.04,
fL_T = max (0.05, (fpeak_T − 0.04)),
Moreover, based on fpeak, a set of parameters was selected to optimize subsequent smoothing and minima/maxima detection phases of the processing block.
Processing block includes all the steps needed to extract breathing frequency and temporal parameters from the signals obtained downstream of the dimension-reduction block. Chest-wall respiratory-related orientation change signals (abdominal and thoracic) underwent the following steps:
(i).
Adaptive band-pass filter. The signals were band-pass filtered (first-order IIR Butterworth filter), with fU and fL cut-off frequency points determined within the spectrum analysis block.
(ii).
Smoothing. Filtered signals were furtherly smoothed (third-order Savitzky–Golay FIR filter) to simplify subsequent identification of maxima and minima points. The level of smoothing (window length) was automatically selected based on fpeak, i.e., increasing window length for decreasing fpeak. Relation between optimal window length values and fpeak values has been determined empirically.
(iii).
Minima and maxima points detection. A set of optimized parameters (i.e., minimum peak distance (MPD) and minimum prominence threshold (MPT)) was automatically selected based on fpeak to optimize recognition of minima and maxima points of the smoothed signals. Optimal MPD and MPT values depending on fpeak were experimentally determined.
(iv).
Breathing frequency extraction. Breath by breath, inspiratory time (TI) was computed as the temporal distance between a minimum point (mi) and the consecutive maximum point (Mi); Expiratory time (TE) was computed as the temporal distance between the maximum point (Mi) and the consecutive minimum point (mi + 1); total time (TTOT) was computed as T TOT = T I + T E [s], duty cycle (DC) was computed as T I T TOT × 100 [%] and breathing frequency was computed as 60 ( T TOT ) [breaths/minute]. A mean value for each of the above-mentioned parameter was computed for each trial (~3 min).

2.4. Experimental Setup

To evaluate the capability of the device and of the proposed methods to correctly estimate breathing frequency (and temporal parameters) 8 healthy volunteers (4 males, 4 females) were enrolled. All subjects gave their informed consent for inclusion before they participated in the study. The study was conducted in accordance with the Declaration of Helsinki, and the protocol (Project identification code n° 534) was approved by the Ethics Committee of Scientific Institute IRCCS Medea (date of approval: 25 January 2018). Chest-wall movements during breathing, in seated and supine position, were measured using the proposed device and optoelectronic plethysmography (OEP) simultaneously. OEP [21] is a technique based on a similar functioning principle of the proposed device; in fact, it allows assessment of ventilatory and breathing pattern by measuring chest-wall movements related to breathing, by using motion capture principles. The system is composed of eight infrared video cameras working at a sampling rate of 60 Hz. It can compute the 3D coordinates of retro-reflective markers positioned on the chest wall in specific anatomic points. From the three-dimensional coordinates of the markers, it is possible to obtain the volume enclosed by the chest-wall surface, by applying the Gauss’s theorem. The chest wall is modelled by a bicompartmental model, composed of rib cage and abdomen, and thus it is possible to investigate the contribution of both the compartments. This is an advantage for the validation of the proposed device, in fact, using OEP as reference method it is possible to compare the data recorded with the thoracic and abdominal units of the device with those obtained by using OEP for the thoracic and abdominal compartments, respectively. OEP has been widely validated against spirometer, in healthy subjects, in different conditions and positions, also during submaximal and maximal exercise on cycle ergometer, obtaining discrepancies in tidal volume measurements always <5% [22,23,24,33,34].
The subjects were prepared, placing the reflective markers according to the 89-marker protocol (previously described in [24,35]) used for seated position and the 52-marker protocol (previously described in [36,37]) used for supine position or, more generally, when a back support is present. Then peripheral IMU-units were placed on the thorax and on the abdomen, while reference IMU-unit was placed on the coccyx in seated position, and on the bed in supine position (Figure 3).
Subjects were then asked to seat or lie on a bed and were invited to perform a slow vital capacity maneuver (SVC) and then to start breathing with the following patterns: (I) quiet breathing (QB), (II) increasing fB but same tidal volume of QB (↑fB, VT=), (III) increasing fB and reducing tidal volume (↑fB, VT↓), (IV) decreasing fB with the same tidal volume of QB (↓fB, VT=), (V) decreasing fB increasing tidal volume of QB (↓fB, VT↑). QB trial was repeated two times, thus, each subject performed 6 trials of the duration of 3 min each. The SVC maneuver was used to align OEP signal and device signals during data analysis, since it is generally recognizable with respect to QB. In fact, SVC requires a maximal inspiration followed by a complete expiration without forced or rapid effort.
The subjects were asked to maintain the same breathing pattern (namely, QB, ↑fB, VT=, ↑fB, VT↓, ↓fB, VT=, ↓fB, VT↑) until the end of the trial; in case of fatigue they were asked to perform a second SVC before returning to QB. This procedure was repeated in seated position and in supine position.

2.5. Statistical Analysis

For each trial, mean values of fB, TI, TE and DC were extracted from the best quaternion components identified online by using “Area” and “Peak” methods, and from the signal obtained with the PCA-based fusion method, both for the thoracic and abdominal tracings. Moreover, to evaluate the performance of the selection methods (“Area” and “Peak”) and their ability to select the best component, the same parameters were obtained for all the quaternion components (q0, q1, q2, q3) and compared with those obtained by OEP, on the abdominal and thoracic compartment, respectively. The “Ideal” quaternion component was identified a posteriori, trial by trial, as the one providing the minimum estimation error of the breathing frequency. Obviously, the “Ideal” component cannot be identified during online analysis, or when a reference method is not present. Thus, for each trial, 5 sets of parameters were available:
  • fB_OEP, TI_OEP, TE_OEP and DC_OEP
  • fB_Peak, TI_Peak, TE_Peak and DC_Peak
  • fB_Area, TI_Area, TE_Area and DC_Area
  • fB_PCA, TI_PCA, TE_PCA and DC_PCA
  • fB_Ideal, TI_Ideal, TE_Ideal and DC_Ideal
Among the entire set of trials, those with fB_OEP < 6 breaths/minute or fB_OEP > 60 breaths/minute were discarded. Then, the absolute (Equation (8)) and relative (Equation (9)) errors of estimation in static conditions (supine and seated position) were computed for each parameter:
A b s o l u t e   E r r o r   ( E ) = | D e v i c e O E P | ,
R e l a t i v e   E r r o r   ( E % ) = | D e v i c e O E P | O E P × 100
For all the dimension-reduction methods (“Area”, “Peak”, “PCA”), mean and standard deviation (SD) were computed for E and E% considering all the subjects and all the trials, for the supine and seated position and compared with those obtained considering the “Ideal” component. The error obtained with the “Ideal” component, identified a posteriori, is thus the minimum error obtainable using a single quaternion component, and represents the performance that the other methods (“Area”, “Peak”, and PCA-fusion) should achieve or beat. For E% obtained in fB and DC estimation, non-parametric alternative to the one-way Analysis of variance (ANOVA) with repeated measures (Friedman test) was performed to assess if significant differences between methods (“Area”, “Peak”, “PCA”) and “Ideal” component occurred, “Ideal”); post-hoc analysis was done performing Wilcoxon signed-rank tests on the different combinations of related methods, applying the correction for multiple comparisons using false discovery rate (FDR) method [38,39].
For fB, TI, TE, linear regression analysis and correlation analysis (Pearson’s product-moment correlation rP, or Spearman’s rank-order correlation rS, if data were not normally distributed) were performed between measurements obtained with the device and measurements obtained with the OEP, for the supine and seated position, respectively.
To assess the agreement between measurements obtained with the device and with the OEP, Bland–Altman analysis was performed plotting the difference of the two paired measurements (device–OEP) against the mean of the two measurements [40,41,42]. Mean of the differences (d) and limits of agreement (LOA: from d − (1.9 × SD) to d + (1.9 × SD)) were calculated. The presence of heteroscedasticity was always examined to assess the presence of proportional biases and/or the correlation between differences and mean values. As proposed by Brehm et al. [43], to determine if data were heteroscedastic a visual inspection of Bland–Altman plots was performed at first. If the errors (y-axes: absolute differences) increased with increasing measured values (x-axes: mean), the data were suspected of being heteroscedastic. Then Kendall’s tau (τ) correlation between the absolute differences and the corresponding means was computed to assess the degree of heteroscedasticity. Data were denoted heteroscedastic when a positive, significant correlation (τ > 0.1 and p-value < 0.05) was found, for other cases data were considered homoscedastic [43].
When heteroscedasticity was present the “classical” 95% confidence and tolerance limits cannot be constructed, thus the approach based on the construction of V-shaped limits was applied: the regression line (ordinary least squares (OLS) best fit) was constructed for differences on mean values and the V-shaped confidence limits (upper confidence limit: UCL, lower confidence limit: LCL) were constructed modelling the variability in the SD of the differences directly as a function of the level of the measurement, using a method based on absolute residuals from a fitted regression line [44,45].

3. Results

3.1. Breathing Patterns

Table 1 presents the mean and SD of breathing rate for each breathing pattern (QB1 e QB2, ↑fB, VT=, ↑fB, VT↓, ↓fB, VT=, ↓fB, VT↑) estimated with OEP and device, using “PCA”, “Area”, and “Peak” methods and the “Ideal” component, for all subjects, in supine and seated position. Sample size (n) of each condition is reported in Table 1 for the breathing pattern ↑fB, VT↓, just one thoracic tracing was available for seated position (n = 1). It can be noticed that each subject demonstrated a different breathing frequency for each breathing pattern and SD in the forced breathing patterns is higher than those obtained for QB, meaning that subjects interpreted the required speed differently.

3.2. Accuracy Errors

Relative errors of estimation in supine and seated position computed for best component-selection methods (“Peak”, “Area”), for PCA-fusion method and for the “Ideal” quaternion component for fB and DC are presented in Figure 4. For what concerns fB estimation in supine position, relative errors obtained using PCA were similar or even better than those provided by the “Ideal” component; on the contrary, both component-selection methods, namely “Area” and “Peak”, provided errors higher than 10%, both for the abdominal and the thoracic compartments. Errors obtained with PCA resulted significantly lower than those obtained with the “Area” method, both for the abdominal (Wilcoxon post-hoc test FDR-adjusted, p = 0.038) and thoracic compartment (Wilcoxon post-hoc test FDR-adjusted, p = 0.015); also, PCA was significantly better than “Peak” method considering abdominal compartment (Wilcoxon post-hoc test FDR-adjusted, p = 0.038). Errors obtained with “Ideal” component resulted significantly lower than those obtained with the component-selection methods both for the abdominal (Wilcoxon post-hoc test FDR-adjusted, Ideal vs. Area p = 0.038, Ideal vs. Peak p = 0.038) and thoracic (Wilcoxon post-hoc test FDR-adjusted, Ideal vs. Area p = 0.000, Ideal vs. Peak p = 0.020) compartments.
In seated position, fB estimation errors obtained with component-selection methods were lower on average than those obtained in supine position, while PCA performances declined. This led to a sort of equalization effect, confirmed also by the statistical analysis: significant differences remained only for comparisons “Ideal” vs. “Area” method (Wilcoxon post-hoc test FDR-adjusted, AB: p = 0.102, TH: p = 0.015) and “Ideal” vs. “Peak” method (Wilcoxon post-hoc test FDR-adjusted, AB: p = 0.006, TH: p = 0.015).
Regarding duty cycle, relative errors of estimation obtained with the different methods are comparable, with exception of those provided by “Ideal” component that are on average lower, both in supine (Wilcoxon post-hoc test FDR-adjusted, AB: Ideal vs. Peak p = 0.042, Ideal vs. Area p = 0.042; TH: Ideal vs. Area p = 0.006) and seated position (Wilcoxon post-hoc test FDR-adjusted, AB: Ideal vs. Peak p = 0.015, Ideal vs. PCA p = 0.006; TH: Ideal vs. Area p = 0.015, Ideal vs. Peak p = 0.015, Ideal vs. PCA p = 0.05).
Absolute estimation errors of fB, TI and TE obtained with the device using different methods (Area, Peak, Ideal, PCA) relative to OEP are reported in Table 2.

3.3. Linear Regression and Correlation Analysis

Scatter plots showing the relationship between measurements obtained with the OEP and with the device, using “Area”, “Peak”, “Ideal” components, and PCA-fusion respectively are presented for fB (Figure 5), TI (Figure 6) and TE (Figure 7). For each scatter plot the regression line is computed, both for thorax and abdomen, and the relative equations are reported.
Correlation coefficients for the comparisons Device vs. OEP are reported in Table 3. Regarding the main parameter, fB, results obtained with correlation analysis confirmed what emerged from estimation error analysis: in supine position, PCA exhibited the best performances in terms of correlation with OEP measurements both in terms of regression line and correlation coefficient. In seated position, “Ideal” component was the one with the highest correlation with OEP measurements, followed by PCA.
With reference to TI estimation in supine position, “Ideal” component provided the best performances, followed by PCA method; “Peak” and “Area” methods provided comparable, poor performances. In seated position, measurements of TI provided by component-selection methods were on average more correlated with OEP measurements than measurements obtained using PCA-fusion method. The “Ideal” component presented the best results, followed by “Area” and “Peak” methods; PCA provided the worst performance considering the abdominal compartment, while correlation between measurements obtained with the thoracic unit and OEP measurements was good.
Estimation of TE was on average more problematic. In terms of regression lines in fact, slope values were far from the unity for all the considered methods, highlighting a proportional error leading to an overestimation for low values of expiratory time and an underestimation at high expiratory times, as shown in Figure 7. For what concerns supine position, correlation coefficients were good both for “Ideal” component and PCA-fusion method; on the contrary, correlation coefficients were low both for “Area” and “Peak” methods. Also, in seated position correlation coefficients provided by “Ideal” and PCA-fusion method were higher than those provided by “Area” and “Peak” methods, especially with respect to the thoracic compartment.

3.4. Bland–Altman Analysis

Bland–Altman plots representing the agreement between measurements obtained with the OEP and with the device, using “Area”, “Peak”, “Ideal” components, and PCA-fusion respectively are presented for fB (Figure 8), TI (Figure 9) and TE (Figure 10). In Bland–Altman plots, the difference of the two paired measurements (device–OEP) is plotted against the mean of the two measurements (device+OEP)⁄2. Results of agreement analysis, including evaluation of heteroscedasticity (Kendall’s τ correlation and relative p-value) are reported in Table 4. As shown there, for homoscedastic data, the mean of the differences representing the fixed bias, and LOAs were computed. On the other hand, for heteroscedastic data, OLS line of best fit representing the proportional bias and upper and lower 95% V-shape confidence limits (UCL and LCL) are reported.
With respect to the main parameter (fB), agreement between OEP and the device is very strong when the “Ideal” component or the PCA-fusion are used, both in supine and seated position. In relation to time estimation, the agreement decreases for all the methods considered. In particular, for what concerns inspiratory times, a significant relationship between errors and mean value emerged, with a general increase of the difference (device–OEP) at higher time values (overestimation of the device), both in supine and seated position. Also, for expiratory times absolute errors increased with increasing time values, but in this case the device underestimated at high time values (negative slope of the OLS).

3.5. Quaternion Component Selection

Considering the quaternion components selected by the “Area” and “Peak” methods as best component or identified as “Ideal” component, a clear rule did not emerge. In fact, there was not a quaternion component that was selected as best component with a considerable frequency. Relative frequencies of quaternion component selection with the different methods are presented in Figure 11. It is interesting to notice that quaternion component q0 was never selected by “Area” and “Peak” methods, while the “Ideal” component was q0 in 14.86 % of cases (n = 74) in supine position and 6.76% of cases (n = 74) in seated position. In seated position, the component q1 was selected more frequently as best component both by using “Area” (44.59) and “Peak” (39.19%) methods. In contrast, in supine position, the components q2 (“Area” 51.35%, “Peak”: 41.89) and q3 (“Area”: 39.19% and “Peak”: 40.54%) were selected more frequently.
Excluding q0 component, that was clearly less selected, the other quaternion components were almost equally selected as “Ideal” component considering all the trials, both in supine position (q0: 14.86%, q1: 22.97%, q2: 32.43%, q3: 29.73%), and seated position (q0: 6.76%, q1: 33.78, q2: 22.97%, q3: 36.49%).
In regard to the ability of the two component-selection methods (“Area” and ”Peak”) to identify the “Ideal” component, i.e., the component providing the minimum fB estimation error, in supine position, the “Area” method was able to identify the “Ideal” component in 45.94% of the cases (relative frequency for the event “the component selected by “Area” method and the “Ideal” component corresponded”), the “Peak” method identified the “Ideal” component in 52.70% of the cases, while for the 45.94% of the cases neither the “Area” method nor the “Peak” method were able to identify the “Ideal” component. In 44.59% of cases, “Area” method and “Peak” method selected the same quaternion component, that was also identified as “Ideal” component.

4. Discussion

In this study, we presented an automatic, position-independent processing algorithm to derive breathing signal, and subsequently breathing temporal parameters, from chest-wall orientation changes acquired using an IMU-based device previously developed by our group [19], composed of three sensor units. Even if the modular configuration of the device was designed to reduce non-breathing movements, the aim of this work was neither to demonstrate the effectiveness of this approach nor to support the presence of the reference unit. On the contrary the focus is on the analysis algorithm, which uses quaternion form to represent orientation, avoiding singularity problem that affects Euler representation; thus, thoracic and abdominal orientation change signals are 4-dimensional entities. The proposed algorithm includes therefore a dimension-reduction block to obtain a 1-dimension signal representing chest-wall orientation changes due to breathing activity.
Another aim of this study was, therefore, to compare three different dimension-reduction methods. The first two methods (“Peak” and “Area”) are based on the selection of the quaternion component with the highest breathing information. The third method is based on the fusion of the 4 components of the quaternion using PCA.
The PCA-fusion method performed better than the best component-selection methods (“Peak” and “Area”) as regards breathing frequency estimation, both in supine and seated position. In supine position, it provided better results than “Ideal” component, while in seated position it provided closer performances to those obtained by using the “Ideal” component with respect to “Area” and “Peak” methods.
About estimation of other temporal parameters (TI, TE and duty cycle), “Ideal” component provided the best results, while PCA-fusion method gave results comparable to the best component-selection methods. Thus, a quaternion component providing the best performance exists, the problem lies in its a priori identification. In fact, both “Area” and “Peak” methods failed to identify it on the basis of spectral analysis (in 45.94% of the cases neither the “Area” method nor the “Peak” method were able to identify the “Ideal” component), and no quaternion component emerged as the most selected as “Ideal” component (supine q0:14.86%, q1: 22.97%, q2: 32.43%, q3: 29.73%; seated q0: 6.76%, q1: 33.78, q2: 22.97%, q3: 36.49%).
Geometrical or morphological considerations to determine which quaternion component is more involved in breathing movement are problematic when considering quaternions, and are position- and IMU-placement dependent, thus not suitable in dynamic conditions. On the contrary, PCA-fusion method represents an interesting solution to this problem because it fuses the information of the four quaternion components regardless the position of the subject or the IMUs placement, avoiding the necessity to select a best component/axes, as reported in previous studies [7,11,13].
In this study, we found that PCA-fusion method provided the best fB estimation performance in terms of mean absolute errors (<2 breaths/minute), correlation (r > 0.963) and agreement (see Table 4) with the reference method. Comparing our results in terms of accuracy errors with those obtained by previous studies is difficult because in most cases only relative errors were reported, but these errors depend on the breathing frequency adopted. Liu et al. [14] reported a mean absolute error of 15.45 breaths/minute (thus about 7 times higher than the error obtained in this study) during quiet sitting. Bates et al. [13] obtained an RMS error of 0.38 breaths/minute and a peak error of 3 breaths/minute in a postoperative patient during sleep. Considering comparable conditions (abdominal compartment in supine position) we obtained an RMS error of 1.51 breaths/minute, using PCA-fusion method, but in our study different, forced, breathing patterns were included, leading to higher mean error.
Regarding correlation between fB measurements obtained with the proposed method and OEP, our results are comparable to those obtained by Bates et al. [13] that reported a correlation coefficient equal to 0.985 between measurements of fB obtained with the accelerometer and nasal cannula. Mann et al. [11] obtained a correlation coefficient of 0.97 between measurements of fB obtained with a tri-axial accelerometer and with a system based on oxygen consumption measurement (Oxycon Mobile). In both cases [11,13], breath-by-breath analysis was not possible. Liu et al. [14] reported correlation coefficients lower than 0.6 between fB computed with a 3-axis accelerometer and with the reference (Airflow CO2 analysis).
There are few studies in the literature that performs Bland–Altman analysis to assess the agreement between breathing frequency measurements obtained by using inertial sensors and other validated methods. In the study from Morillo et al. [8] agreement analysis using Bland Altman plots was done against PSG thermistor. The authors reported a mean difference (fixed bias) of 0.02 (SD = 1.09) breaths/minute and LOAs from −3.05 to +2.11 breaths/minute in the range ~12 ÷ 35 breaths/minute. In that case, the use of a single-axis accelerometer, prevent the use of that method during postural changes. Dehkrodi et al. [9] used a tri-axial accelerometer placed on the suprasternal notch extending the validation presented in [8] to different sleeping positions and breathing conditions (Deep: 13.5 ± 4.3, Normal: 16.5 ± 5.2 and Shallow: 39.7 ± 30.3 breath/minute). They reported a mean difference (fixed bias) of 0.042 breaths/minute and LOAs from –0.65 to 0.74 breaths/minute. Lapi et al. [17] performed agreement analysis with Bland Altman plots for measurements of breathing frequency (range 12 ÷ 26 breaths/minute) obtained with the accelerometer and with the standard method (counting breaths by visual inspection), in supine position. They reported a mean difference (fixed bias) of 0.33 breaths/minute and LOA from −1.92 to 2.60 with 3.2% of data outside that range. For all the above-mentioned studies heteroscedasticity of data was not considered or reported making it difficult compare them directly with our results. In fact, taking into account heteroscedasticity of data, for fB in supine position, we built Bland–Altman plot with proportional bias (OLS: y = 0.008x + 0.130, thus going from 0.18 (at x = 6 breaths/minute) to 0.61 (at x = 60 breaths/minute) breaths/minute) and V-shaped limits (LCL: y = −0.038x − 2.039 thus the lower limit goes from −2.26 (at x = 6 breaths/minute) to −4.32 (at x = 60 breaths/minute) breath/minute; UCL: y = 0.054x + 2.299; thus the upper limit goes form 2.62 (at x = 6 breaths/minute) to 5.539 (at x = 60 breaths/minute) breaths/minute) ). Thus, considering comparable breathing frequency ranges our results are closer to those obtained by Morillo et al. [8] and Lapi et al. [17]. On the contrary, Dehkrodi et al. [9] obtained better results; unfortunately, the steps to obtain the acceleration derived respiratory (ADR) signal are not described in detail.
For the best of our knowledge, this is the first study that provides a detailed analysis of respiratory timing measurements obtained by using inertial sensor systems, validating them against an established method.

5. Conclusions

PCA-fusion method provided overall best performances with respect to selecting the best quaternion component identified based on spectrum analysis. In supine position results obtained fusing the 4 quaternion components were even better than those obtained with the “Ideal” component, identified a posteriori considering the minimum breathing frequency estimation error. Performance in seated position were worse than those obtained in supine position, probably because subjects were seated without the back support and some oscillations of the trunk were more likely to occur. This could particularly affect PCA-based method where the first principal component selected for further analysis is the one with the largest variance, and thus more subject to larger body motions. This must be taken into account in dynamic conditions; signal baseline removal prior to PCA-fusion should be considered in this case.
The analysis algorithm proposed in this work, applying PCA-fusion as dimension-reduction method, can be used to analyze further data. In fact, an extended validation of the proposed device and method is needed also in dynamic conditions, during daily activities, considering not only healthy subjects but also patients that could particularly take advantage of this system (e.g., COPD, neuromuscular patients, sleep apnea, etc.). This would also allow study of asynchronies of thoraco-abdominal compartments taking advantage of the modular configuration of the device. Another key step will be the adaptation of our method, currently implemented as an offline analysis, to online monitoring, moving the computation process aboard the smartphone. This enhancement could allow immediate computation of an average breathing frequency over a certain period (e.g., 60 s) directly aboard the smartphone, fostering the use of the device in other applications such as sport and fitness, exercise testing, breathing training to use different respiratory muscles, rehabilitation protocols and treatment evaluation where respiratory assessment could be of great interest.

6. Patents

The present work is partially described in the International Patent application n° PCT/IB2018/054956, priority date 11 July 2017, title “A wearable device for the continuous monitoring of the respiratory rate”. Inventors: Ambra Cesareo, Andrea Aliverti, Assignee: Politecnico di Milano.

Author Contributions

A.C. and A.A. contributed IMU circuit design and fabrication; A.C. conceived and designed the experiments; A.C. and Y.P. performed the experiments and analyzed the data; A.C. and E.B. prepared the original draft; E.B. and A.A. reviewed the paper; E.B. and A.A. supervised the project and contributed to funding acquisition.

Funding

This research received no external funding.

Acknowledgments

This work was supported by “FoRST -Fondazione per la Ricerca Scientifica Termale”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Subbe, C.; Davies, R.; Williams, E.; Rutherford, P.; Gemmell, L. Effect of Introducing the Modified Early Warning Score on Clinical Outcomes, cardio-pulmonary Arrests and Intensive Care Utilisation in Acute Medical Admissions. Anaesthesia 2003, 58, 797–802. [Google Scholar] [CrossRef] [PubMed]
  2. Cretikos, M.A.; Bellomo, R.; Hillman, K.; Chen, J.; Finfer, S.; Flabouris, A. Respiratory Rate: The Neglected Vital Sign. Med. J. Aust. 2008, 188, 657. [Google Scholar]
  3. World Health Organization. Acute Respiratory Infections in Children: Case Management in Small Hospitals in Developing Countries, a Manual for Doctors and Other Senior Health Workers; World Health Organization: Geneva, Switzerland, 1990. [Google Scholar]
  4. Karlen, W.; Gan, H.; Chiu, M.; Dunsmuir, D.; Zhou, G.; Dumont, G.A.; Ansermino, J.M. Improving the Accuracy and Efficiency of Respiratory Rate Measurements in Children using Mobile Devices. PLoS ONE 2014, 9, e99266. [Google Scholar] [CrossRef] [PubMed]
  5. Torres, A.; Fiz, J.; Galdiz, B.; Gea, J.; Morera, J.; Jané, R. Assessment of Respiratory Muscle Effort Studying Diaphragm Movement Registered with Surface Sensors. Animal Model (Dogs). In Proceedings of the 26th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (IEMBS’04), San Francisco, CA, USA, 1–4 September 2004; pp. 3917–3920. [Google Scholar]
  6. Reinvuo, T.; Hannula, M.; Sorvoja, H.; Alasaarela, E.; Myllyla, R. Measurement of Respiratory Rate with High-Resolution Accelerometer and EMFit Pressure Sensor. In Proceedings of the 2006 IEEE Sensors Applications Symposium, Houston, TX, USA, 7–9 February 2006; pp. 192–195. [Google Scholar]
  7. Hung, P.; Bonnet, S.; Guillemaud, R.; Castelli, E.; Yen, P.T.N. Estimation of Respiratory Waveform using an Accelerometer. In Proceedings of the 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro (ISBI 2008), Paris, France, 14–17 May 2008; pp. 1493–1496. [Google Scholar]
  8. Morillo, D.S.; Ojeda, J.L.R.; Foix, L.F.C.; Jiménez, A.L. An Accelerometer-Based Device for Sleep Apnea Screening. IEEE Trans. Inf. Technol. Biomed. 2010, 14, 491–499. [Google Scholar] [CrossRef] [PubMed]
  9. Dehkordi, P.K.; Marzencki, M.; Tavakolian, K.; Kaminska, M.; Kaminska, B. Validation of Respiratory Signal Derived from Suprasternal Notch Acceleration for Sleep Apnea Detection. In Proceedings of the 2011 Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Boston, MA, USA, 30 August–3 September 2011; pp. 3824–3827. [Google Scholar]
  10. Fekr, A.R.; Janidarmian, M.; Radecka, K.; Zilic, Z. A Medical Cloud-Based Platform for Respiration Rate Measurement and Hierarchical Classification of Breath Disorders. Sensors 2014, 14, 11204–11224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Mann, J.; Rabinovich, R.; Bates, A.; Giavedoni, S.; MacNee, W.; Arvind, D. Simultaneous Activity and Respiratory Monitoring using an Accelerometer. In Proceedings of the 2011 International Conference Body Sensor Networks (BSN), Dallas, TX, USA, 23–25 May 2011; pp. 139–143. [Google Scholar]
  12. Jin, A.; Yin, B.; Morren, G.; Duric, H.; Aarts, R.M. Performance Evaluation of a Tri-Axial Accelerometry-Based Respiration Monitoring for Ambient Assisted Living. In Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC 2009), Minneapolis, MN, USA, 2–6 September 2009; pp. 5677–5680. [Google Scholar]
  13. Bates, A.; Ling, M.J.; Mann, J.; Arvind, D. Respiratory Rate and Flow Waveform Estimation from Tri-Axial Accelerometer Data. In Proceedings of the 2010 International Conference on Body Sensor Networks (BSN), Singapore, 7–9 June 2010; pp. 144–150. [Google Scholar]
  14. Liu, G.; Guo, Y.; Zhu, Q.; Huang, B.; Wang, L. Estimation of Respiration Rate from Three-Dimensional Acceleration Data Based on Body Sensor Network. Telemed. e-Health 2011, 17, 705–711. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Yoon, J.; Noh, Y.; Kwon, Y.; Kim, W.; Yoon, H. Improvement of Dynamic Respiration Monitoring through Sensor Fusion of Accelerometer and Gyro-Sensor. J. Electr. Eng. Technol. 2014, 9, 334–343. [Google Scholar] [CrossRef]
  16. Gollee, H.; Chen, W. Real-Time Detection of Respiratory Activity using an Inertial Measurement Unit. In Proceedings of the 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS 2007), Lyon, France, 23–26 August 2007; pp. 2230–2233. [Google Scholar]
  17. Lapi, S.; Lavorini, F.; Borgioli, G.; Calzolai, M.; Masotti, L.; Pistolesi, M.; Fontana, G.A. Respiratory Rate Assessments using a Dual-Accelerometer Device. Respir. Physiol. Neurobiol. 2014, 191, 60–66. [Google Scholar] [CrossRef]
  18. Gaidhani, A.; Moon, K.S.; Ozturk, Y.; Lee, S.Q.; Youm, W. Extraction and Analysis of Respiratory Motion using Wearable Inertial Sensor System during Trunk Motion. Sensors 2017, 17, 2932. [Google Scholar] [CrossRef]
  19. Cesareo, A.; Gandolfi, S.; Pini, I.; Biffi, E.; Reni, G.; Aliverti, A. A Novel, Low Cost, Wearable Contact-Based Device for Breathing Frequency Monitoring. In Proceedings of the 2017 39th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Seogwipo, Korea, 11–15 July 2017; pp. 2402–2405. [Google Scholar]
  20. Romei, M.; Mauro, A.L.; D’angelo, M.; Turconi, A.; Bresolin, N.; Pedotti, A.; Aliverti, A. Effects of Gender and Posture on Thoraco-Abdominal Kinematics during Quiet Breathing in Healthy Adults. Respir. Physiol. Neurobiol. 2010, 172, 184–191. [Google Scholar] [CrossRef]
  21. Aliverti, A.; Pedotti, A. Opto-electronic plethysmography. In Mechanics of Breathing; Springer: Milano, Italy, 2002; pp. 47–59. [Google Scholar]
  22. Aliverti, A.; Dellacà, R.; Pelosi, P.; Chiumello, D.; Gattinoni, L.; Pedotti, A. Compartmental Analysis of Breathing in the Supine and Prone Positions by Optoelectronic Plethysmography. Ann. Biomed. Eng. 2001, 29, 60–70. [Google Scholar] [CrossRef] [PubMed]
  23. Aliverti, A.; Dellaca, R.; Pelosi, P.; Chiumello, D.; Pedotti, A.; Gattinoni, L. Optoelectronic Plethysmography in Intensive Care Patients. Am. J. Respir. Crit. Care Med. 2000, 161, 1546–1552. [Google Scholar] [CrossRef] [PubMed]
  24. Cala, S.; Kenyon, C.; Ferrigno, G.; Carnevali, P.; Aliverti, A.; Pedotti, A.; Macklem, P.; Rochester, D. Chest Wall and Lung Volume Estimation by Optical Reflectance Motion Analysis. J. Appl. Physiol. 1996, 81, 2680–2689. [Google Scholar] [CrossRef] [PubMed]
  25. Ferrigno, G.; Carnevali, P.; Aliverti, A.; Molteni, F.; Beulcke, G.; Pedotti, A. Three-Dimensional Optical Analysis of Chest Wall Motion. J. Appl. Physiol. 1994, 77, 1224–1231. [Google Scholar] [CrossRef] [PubMed]
  26. Konno, K.; Mead, J. Measurement of the separate volume changes of rib cage and abdomen during breathing. J. Appl. Physiol. 1967, 22, 407–422. [Google Scholar] [CrossRef] [PubMed]
  27. Hamilton, W.R. XI. On Quaternions; Or on a New System of Imaginaries in Algebra. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1848, 33, 58–60. [Google Scholar] [CrossRef]
  28. Rosenfeld, B. The History of Non-Euclidean Geometry: Evolution of the Concept of a Geometrical Space; Springer-Verlag: Berlin, Germany, 1988. [Google Scholar]
  29. Madgwick, S.O.; Harrison, A.J.; Vaidyanathan, R. Estimation of IMU and MARG Orientation using a Gradient Descent Algorithm. In Proceedings of the 2011 IEEE International Conference on Rehabilitation Robotics (ICORR), Zurich, Switzerland, 29 June–1 July 2011; pp. 1–7. [Google Scholar]
  30. Pearson, K. LIII. On Lines and Planes of Closest Fit to Systems of Points in Space. Lond. Edinb. Dublin Philos. Mag. J. Sci. 1901, 2, 559–572. [Google Scholar] [CrossRef]
  31. Hotelling, H. Analysis of a Complex of Statistical Variables into Principal Components. J. Educ. Psychol. 1933, 24, 417. [Google Scholar] [CrossRef]
  32. Savitzky, A.; Golay, M.J. Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Anal. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  33. Kenyon, C.M.; Cala, S.J.; Yan, S.; Aliverti, A.; Scano, G.; Duranti, R.; Pedotti, A.; Macklem, P.T. Rib cage mechanics during quiet breathing and exercise in humans. J. Appl. Physiol. 1997, 83, 1242–1255. [Google Scholar] [CrossRef]
  34. Layton, A.M.; Moran, S.L.; Garber, C.E.; Armstrong, H.F.; Basner, R.C.; Thomashow, B.M.; Bartels, M.N. Optoelectronic plethysmography compared to spirometry during maximal exercise. Respir. Physiol. Neurobiol. 2013, 185, 362–368. [Google Scholar] [CrossRef] [PubMed]
  35. LoMauro, A.; Cesareo, A.; Agosti, F.; Tringali, G.; Salvadego, D.; Grassi, B.; Sartorio, A.; Aliverti, A. Effects of a Multidisciplinary Body Weight Reduction Program on Static and Dynamic Thoraco-Abdominal Volumes in Obese Adolescents. Appl. Physiol. Nutr. Metab. 2016, 41, 649–658. [Google Scholar] [CrossRef] [PubMed]
  36. Lo Mauro, A.; D’Angelo, M.G.; Romei, M.; Motta, F.; Colombo, D.; Comi, G.P.; Pedotti, A.; Marchi, E.; Turconi, A.C.; Bresolin, N.; et al. Abdominal Volume Contribution to Tidal Volume as an Early Indicator of Respiratory Impairment in Duchenne Muscular Dystrophy. Eur. Respir. J. 2010, 35, 1118–1125. [Google Scholar] [CrossRef] [PubMed]
  37. Cesareo, A.; LoMauro, A.; Santi, M.; Biffi, E.; D’Angelo, M.G.; Aliverti, A. Acute Effects of Mechanical Insufflation-Exsufflation on the Breathing Pattern in Stable Subjects with Duchenne Muscular Dystrophy. Respir. Care 2018. [Google Scholar] [CrossRef] [PubMed]
  38. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  39. Benjamini, Y.; Yekutieli, D. The Control of the False Discovery Rate in Multiple Testing Under Dependency. Ann. Stat. 2001, 1165–1188. [Google Scholar]
  40. Altman, D.G.; Bland, J.M. Measurement in Medicine: The Analysis of Method Comparison Studies. Statistician 1983, 307–317. [Google Scholar] [CrossRef]
  41. Bland, J.M.; Altman, D. Statistical Methods for Assessing Agreement between Two Methods of Clinical Measurement. Lancet 1986, 327, 307–310. [Google Scholar] [CrossRef]
  42. Bland, J.M.; Altman, D.G. Measuring Agreement in Method Comparison Studies. Stat. Methods Med. Res. 1999, 8, 135–160. [Google Scholar] [CrossRef]
  43. BREHM, M.; Scholtes, V.A.; Dallmeijer, A.J.; Twisk, J.W.; Harlaar, J. The Importance of Addressing Heteroscedasticity in the Reliability Analysis of ratio-scaled Variables: An Example Based on Walking energy-cost Measurements. Dev. Med. Child Neurol. 2012, 54, 267–273. [Google Scholar] [CrossRef]
  44. Bland, J. How Do I Estimate Limits of Agreement When the Mean or SD of Differences Is Not Constant? Available online: https://www-users.york.ac.uk/~mb55/meas/glucose.htm (accessed on 4 October 2018).
  45. Ludbrook, J. Confidence in Altman–Bland Plots: A Critical Review of the Method of Differences. Clin. Exp. Pharmacol. Physiol. 2010, 37, 143–149. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Block Diagram of the Analysis algorithm that allows derivation of breathing temporal parameters (fB. TI, TE) from quaternion-based orientation change signals recorded on Thorax, Abdomen and Reference point.
Figure 1. Block Diagram of the Analysis algorithm that allows derivation of breathing temporal parameters (fB. TI, TE) from quaternion-based orientation change signals recorded on Thorax, Abdomen and Reference point.
Sensors 19 00088 g001
Figure 2. Dimension-reduction block in detail. Starting from the 4 components [q0, q1, q2, q3] of each quaternion (Abdominal: Ab and Thoracic: Th), three methods are applied to obtain a single-component signal: two methods based on best quaternion component selection (“Area” and “Peak”) and one method based on the fusion of the 4 components through Principal Component Analysis (PCA). “Area” method selects the quaternion component with the larger area under the Power Spectral Density (PSD) estimate, while “Peak” method selects the quaternion component with the highest PSD’s peak. PCA-fusion method selects the first principal component (PC_1) that accounts for the largest variance in the data.
Figure 2. Dimension-reduction block in detail. Starting from the 4 components [q0, q1, q2, q3] of each quaternion (Abdominal: Ab and Thoracic: Th), three methods are applied to obtain a single-component signal: two methods based on best quaternion component selection (“Area” and “Peak”) and one method based on the fusion of the 4 components through Principal Component Analysis (PCA). “Area” method selects the quaternion component with the larger area under the Power Spectral Density (PSD) estimate, while “Peak” method selects the quaternion component with the highest PSD’s peak. PCA-fusion method selects the first principal component (PC_1) that accounts for the largest variance in the data.
Sensors 19 00088 g002
Figure 3. Experimental setup. Retroreflective-marker configuration for optoelectronic plethysmography (OEP) and IMU-unit (Ab: Abdomen, Th: Thorax, Ref: Reference) placement in supine (A and B panels) and seated (C and D panels) positions. Panel E shows the experimental setup and the OEP Lab; Infrared cameras of the motion capture system are also noticeable.
Figure 3. Experimental setup. Retroreflective-marker configuration for optoelectronic plethysmography (OEP) and IMU-unit (Ab: Abdomen, Th: Thorax, Ref: Reference) placement in supine (A and B panels) and seated (C and D panels) positions. Panel E shows the experimental setup and the OEP Lab; Infrared cameras of the motion capture system are also noticeable.
Sensors 19 00088 g003
Figure 4. Relative errors (E%) of estimation of breathing frequency (a,b) and Duty Cycle (c,d) in supine (a,c) and seated (b,d) positions, computed for each method (Peak, Area and PCA) and for the “Ideal” component with respect to the reference (OEP). Errors are computed both for the Thoracic and abdominal compartments. Horizontal blue lines indicate statistical significance of difference (post-hoc analysis, Wilcoxon test FDR corrected).
Figure 4. Relative errors (E%) of estimation of breathing frequency (a,b) and Duty Cycle (c,d) in supine (a,c) and seated (b,d) positions, computed for each method (Peak, Area and PCA) and for the “Ideal” component with respect to the reference (OEP). Errors are computed both for the Thoracic and abdominal compartments. Horizontal blue lines indicate statistical significance of difference (post-hoc analysis, Wilcoxon test FDR corrected).
Sensors 19 00088 g004
Figure 5. Comparisons of breathing frequency (fB expressed as breaths/minuteute) measurements by using the proposed device and by using Optoelectronic plethysmography (OEP) presented as regression analysis, in supine (top panels) and seated (bottom panels) positions. For what concerns fB measurements obtained with the IMU-device, three dimension-reduction methods were considered: Area, Peak and PCA-fusion. The performance obtained by using these three methods is benchmarked against that obtained with the Ideal quaternion component determined a posteriori based on the minimum estimation error. The regression line between measurements done by OEP and the proposed device is plotted, and the relative equation presented, both for the thorax and the abdomen.
Figure 5. Comparisons of breathing frequency (fB expressed as breaths/minuteute) measurements by using the proposed device and by using Optoelectronic plethysmography (OEP) presented as regression analysis, in supine (top panels) and seated (bottom panels) positions. For what concerns fB measurements obtained with the IMU-device, three dimension-reduction methods were considered: Area, Peak and PCA-fusion. The performance obtained by using these three methods is benchmarked against that obtained with the Ideal quaternion component determined a posteriori based on the minimum estimation error. The regression line between measurements done by OEP and the proposed device is plotted, and the relative equation presented, both for the thorax and the abdomen.
Sensors 19 00088 g005
Figure 6. Comparisons of inspiratory time (TI expressed as seconds) measurements by using the proposed device and by using Optoelectronic plethysmography (OEP) presented as regression analysis, in supine (top panels) and seated (bottom panels) positions. For what concerns TI measurements obtained with the IMU-device, three dimension-reduction methods were considered: Area, Peak and PCA-fusion. The performance obtained by using these three methods is benchmarked against that obtained with the Ideal quaternion component determined a posteriori based on the minimum estimation error. The regression line between measurements done by OEP and the proposed device is plotted, and the relative equation presented, both for the thorax and the abdomen.
Figure 6. Comparisons of inspiratory time (TI expressed as seconds) measurements by using the proposed device and by using Optoelectronic plethysmography (OEP) presented as regression analysis, in supine (top panels) and seated (bottom panels) positions. For what concerns TI measurements obtained with the IMU-device, three dimension-reduction methods were considered: Area, Peak and PCA-fusion. The performance obtained by using these three methods is benchmarked against that obtained with the Ideal quaternion component determined a posteriori based on the minimum estimation error. The regression line between measurements done by OEP and the proposed device is plotted, and the relative equation presented, both for the thorax and the abdomen.
Sensors 19 00088 g006
Figure 7. Comparisons of expiratory time (TE, expressed as seconds) measurements by using the proposed device and by using Optoelectronic plethysmography (OEP) presented as regression analysis, in supine (top panels) and seated (bottom panels) positions. Regarding TE measurements obtained with the IMU-device, three dimension-reduction methods were considered. Area, Peak and PCA-fusion. The performance obtained by using these three methods is benchmarked against that obtained with the Ideal quaternion component determined a posteriori based on the minimum estimation error. The regression line between measurements done by OEP and the proposed device is plotted, and the relative equation presented, both for the thorax and the abdomen.
Figure 7. Comparisons of expiratory time (TE, expressed as seconds) measurements by using the proposed device and by using Optoelectronic plethysmography (OEP) presented as regression analysis, in supine (top panels) and seated (bottom panels) positions. Regarding TE measurements obtained with the IMU-device, three dimension-reduction methods were considered. Area, Peak and PCA-fusion. The performance obtained by using these three methods is benchmarked against that obtained with the Ideal quaternion component determined a posteriori based on the minimum estimation error. The regression line between measurements done by OEP and the proposed device is plotted, and the relative equation presented, both for the thorax and the abdomen.
Sensors 19 00088 g007
Figure 8. Agreement analysis between OEP and the IMU-based device for breathing frequency (fB, expressed as breaths/minuteute) measurements, in supine (top panels) and seated (bottom panels) position. In each Bland–Altman plot the differences between measurements of fB obtained by using the IMU-based device and by using OEP are plotted against the mean of the two measurements. For homoscedastic data, the mean of the differences (bias: —) and limits of agreement (black dotted line) from mean − 1.96 s to mean + 1.96 s are represented by lines parallel to the X axis. For heteroscedastic data, the proportional bias (—) is represented by the ordinary least squares (OLS) line of best fit for the difference on mean values; V-shaped upper and lower 95% confidence limits (- - -) are calculated according to Bland [44].
Figure 8. Agreement analysis between OEP and the IMU-based device for breathing frequency (fB, expressed as breaths/minuteute) measurements, in supine (top panels) and seated (bottom panels) position. In each Bland–Altman plot the differences between measurements of fB obtained by using the IMU-based device and by using OEP are plotted against the mean of the two measurements. For homoscedastic data, the mean of the differences (bias: —) and limits of agreement (black dotted line) from mean − 1.96 s to mean + 1.96 s are represented by lines parallel to the X axis. For heteroscedastic data, the proportional bias (—) is represented by the ordinary least squares (OLS) line of best fit for the difference on mean values; V-shaped upper and lower 95% confidence limits (- - -) are calculated according to Bland [44].
Sensors 19 00088 g008
Figure 9. Agreement analysis between OEP and the IMU-based device for inspiratory time (TI, [s]) measurements, in supine (top panels) and seated (bottom panels) position. In each Bland–Altman plot the differences between measurements of TI obtained by using the IMU-based device and by using OEP are plotted against the mean of the two measurements. For homoscedastic data, the mean of the differences (bias: —) and limits of agreement (- - -) from mean − 1.96 s to mean + 1.96 s are represented by lines parallel to the X axis. For heteroscedastic data, the proportional bias (—) is represented by the OLS line of best fit for differences on mean values; V-shaped upper and lower 95% confidence limits (- - -) are calculated according to Bland [44].
Figure 9. Agreement analysis between OEP and the IMU-based device for inspiratory time (TI, [s]) measurements, in supine (top panels) and seated (bottom panels) position. In each Bland–Altman plot the differences between measurements of TI obtained by using the IMU-based device and by using OEP are plotted against the mean of the two measurements. For homoscedastic data, the mean of the differences (bias: —) and limits of agreement (- - -) from mean − 1.96 s to mean + 1.96 s are represented by lines parallel to the X axis. For heteroscedastic data, the proportional bias (—) is represented by the OLS line of best fit for differences on mean values; V-shaped upper and lower 95% confidence limits (- - -) are calculated according to Bland [44].
Sensors 19 00088 g009
Figure 10. Agreement analysis between OEP and the IMU-based device for expiratory time (TE, [s]) measurements, in supine (top panels) and seated (bottom panels) position. In each Bland–Altman plot the differences between measurements of TE obtained by using the IMU-based device and by using OEP are plotted against the mean of the two measurements. For homoscedastic data, the mean of the differences (bias: —) and limits of agreement (- - -) from mean − 1.96 s to mean + 1.96 s are represented by lines parallel to the X axis. For heteroscedastic data, the proportional bias (—) is represented by the OLS line of best fit for differences on mean values; V-shaped upper and lower 95% confidence limits (- - -) are calculated according to Bland [44].
Figure 10. Agreement analysis between OEP and the IMU-based device for expiratory time (TE, [s]) measurements, in supine (top panels) and seated (bottom panels) position. In each Bland–Altman plot the differences between measurements of TE obtained by using the IMU-based device and by using OEP are plotted against the mean of the two measurements. For homoscedastic data, the mean of the differences (bias: —) and limits of agreement (- - -) from mean − 1.96 s to mean + 1.96 s are represented by lines parallel to the X axis. For heteroscedastic data, the proportional bias (—) is represented by the OLS line of best fit for differences on mean values; V-shaped upper and lower 95% confidence limits (- - -) are calculated according to Bland [44].
Sensors 19 00088 g010
Figure 11. Relative frequencies of quaternion component (q0, q1, q2, q3) selection using Area and Peak methods and of quaternion component selection as Ideal component, in supine and seated position. Each portion of the rings represents the ratio between the number of times that each quaternion component has been selected (by Area and Peak methods or as Ideal component respectively) and the total number of trials (n = 74).
Figure 11. Relative frequencies of quaternion component (q0, q1, q2, q3) selection using Area and Peak methods and of quaternion component selection as Ideal component, in supine and seated position. Each portion of the rings represents the ratio between the number of times that each quaternion component has been selected (by Area and Peak methods or as Ideal component respectively) and the total number of trials (n = 74).
Sensors 19 00088 g011
Table 1. Breathing frequencies for different breathing patterns.
Table 1. Breathing frequencies for different breathing patterns.
SupineQB 1↑fB, VT=↓fB, VT=↑fB, VT↓fB, VTQB 2
AB (n = 8)TH (n = 7)AB (n = 8)TH (n = 8)AB (n = 6)TH (n = 6)AB (n = 3)TH (n = 2)AB (n = 5)TH (n = 6)AB (n = 8)TH (n = 8)
OEP17.13 ± 2.2317.62 ± 2.2139.49 ± 10.2639.48 ± 10.2511.17 ± 2.6411.14 ± 2.6348.11 ± 13.2947.66 ± 13.958.38 ± 2.408.61 ± 2.2415.29 ± 4.3415.24 ± 4.48
DeviceArea17.20 ± 2.1518.53 ± 3.9639.19 ± 10.8438.16 ± 14.6715.52 ± 8.1419.29 ± 9.4234.47 ± 34.5830.01 ± 25.5210.11 ± 3.9110.62 ± 4.4421.28 ± 8.1416.55 ± 3.62
Peak17.20 ± 2.1516.63 ±2.2938.92 ± 11.1238.96 ± 13.0115.52 ± 8.1416.97 ± 9.5934.47 ± 34.5853.39 ± 7.5510.11 ± 3.919.32 ± 2.1221.25 ± 8.1516.55 ± 3.62
PCA17.23 ± 2.0917.03 ± 1.9840.06 ± 9.2240.63 ± 8.4511.51 ± 2.5911.56 ± 3.0448.61 ± 12.6449.93 ± 11.948.61 ± 2.198.95 ± 2.0615.11 ± 4.8215.23 ± 3.79
Ideal17.57 ± 2.5216.34 ± 2.0039.81 ± 10.5040.94 ± 9.9712.11 ± 1.9913.60 ± 2.6149.53 ± 12.3148.12 ± 13.899.34 ± 3.338.83 ± 2.0617.74 ± 4.7015.37 ± 3.82
SeatedQB 1↑fB, VT=↓fB, VT=↑fB, VT↓fB, VTQB 2
AB (n = 8)TH (n = 8)AB (n = 8)TH (n = 7)AB (n = 8)TH (n = 8)AB (n = 2)TH (n = 1)AB (n = 6)TH (n = 5)AB (n = 7)TH (n = 6)
OEP16.99 ± 2.6516.94 ± 2.7742.74 ± 10.9746.49 ± 7.5814.57 ± 13.4014.49 ± 13.4133.20 ± 16.0122.0110.04 ± 3.5610.05 ± 3.8118.47 ± 4.3218.47 ± 4.49
DeviceArea17.15 ± 2.8615.78 ± 3.6243.06 ± 11.2240.71 ± 13.5915.71 ± 13.5217.06 ± 13.0532.63 ± 12.8727.6210.83 ± 2.9312.58 ± 4.0819.80 ± 5.1517.04 ± 2.09
Peak16.48 ± 3.5816.27 ± 3.7643.06 ± 11.2240.71 ± 13.5916.53 ± 13.5017.15 ± 12.9932.63 ± 12.8727.6210.83 ± 2.9312.58 ± 4.0817.23 ± 2.1916.54 ± 1.95
PCA16.54 ± 3.0116.58 ± 2.1042.26 ± 11.3845.31 ± 8.5616.07 ± 13.0714.91 ± 13.3933.79 ± 14.5027.2710.76 ± 2.9011.35 ± 3.7316.59 ± 2.1715.47 ± 1.69
Ideal16.68 ± 2.5916.61 ± 2.3442.64 ± 11.2145.42 ± 7.9815.29 ± 13.3311.07 ± 2.3633.38 ± 13.9224.5310.40 ± 3.5310.29 ± 3.6619.06 ± 4.8918.73 ± 5.04
Across subject mean ± SD of the breathing frequency (fB, [breaths/minute]) measurements with OEP and the device, using best component-selection methods (“Area” and “Peak”), PCA-fusion method and “Ideal” component, for the requested patterns. Data are reported for the supine and seated positions, subdivided in abdominal and thoracic contributions.
Table 2. Absolute errors of breathing frequency (E_fB), Inspiratory time (E_TI) and expiratory time (E_TE) obtained for the device with respect to OEP, using best component-selection methods (“Area” and “Peak”), PCA-fusion method and “Ideal” component.
Table 2. Absolute errors of breathing frequency (E_fB), Inspiratory time (E_TI) and expiratory time (E_TE) obtained for the device with respect to OEP, using best component-selection methods (“Area” and “Peak”), PCA-fusion method and “Ideal” component.
AreaPeakPCAIdeal
E_fB [breaths/minute]supineAB3.64 ± 7.463.64 ± 7.461.00 ± 1.241.39 ± 2.76
TH5.46 ± 8.893.17 ± 4.921.55 ± 1.511.56 ± 1.96
seatedAB2.19 ± 2.492.12 ± 2.741.71 ± 2.251.04 ± 1.24
TH3.35 ± 5.683.31 ± 5.691.79 ± 2.040.96 ± 0.22
E_TI [s]supineAB0.48 ± 0.730.48 ± 0.730.33 ± 0.510.20 ± 0.38
TH0.43 ± 0.520.41 ± 0.490.47 ± 0.670.17 ± 0.25
seatedAB0.33 ± 0.580.36 ± 0.560.46 ± 0.710.16 ± 0.27
TH0.43 ± 0.480.44 ± 0.490.42 ± 0.350.17 ± 0.25
E_TE [s]supineAB0.58 ± 0.820.58 ± 0.820.43 ± 0.580.29 ± 0.52
TH0.79 ± 0.940.67 ± 0.920.46 ± 0.630.36 ± 0.71
seatedAB0.43 ± 0.560.43 ± 0.560.43 ± 0.550.22 ± 0.31
TH0.56 ± 0.660.56 ± 0.660.39 ± 0.410.24 ± 0.36
Data are reported as mean ± SD, in supine and seated position for thoracic (TH) and abdominal (AB) compartments.
Table 3. Correlation outcomes across subjects and breathing patterns. Coefficient of correlation (r) between measurements obtained using Device vs. OEP are reported for fB. TI. and TE using best component-selection methods (“Area” and “Peak”), PCA-fusion method and “Ideal” component, in supine (Thorax: n = 37. Abdomen: n = 37) and seated (Thorax: n = 35. Abdomen: n = 39) position.
Table 3. Correlation outcomes across subjects and breathing patterns. Coefficient of correlation (r) between measurements obtained using Device vs. OEP are reported for fB. TI. and TE using best component-selection methods (“Area” and “Peak”), PCA-fusion method and “Ideal” component, in supine (Thorax: n = 37. Abdomen: n = 37) and seated (Thorax: n = 35. Abdomen: n = 39) position.
SupineSeated
ThoraxAbdomenThoraxAbdomen
fBArea0.580 $0.706 $0.748 $0.915 $
Peak0.833 $0.706 $0.759 $0.861 $
PCA0.963$0.985$0.953 $0.924 $
Ideal0.935 $0.931 $0.974$0.977$
TIArea0.727 #0.665 $0.812 #0.812 #
Peak0.785 #0.659 $0.809 #0.824 #
PCA0.783 #0.874#0.926 #0.731 #
Ideal0.943#0.862 $0.951#0.948#
TEArea0.600 $0.713 #0.682 #0.818 #
Peak0.687 $0.712 #0.723 #0.835 #
PCA0.891#0.864 #0.888 #0.824 #
Ideal0.874 #0.966$0.938#0.951#
Correlations are all significant (p-value < 0.001). $ Spearman correlation coefficient; # Pearson correlation coefficient; Bold: best correlation result.
Table 4. Agreement analysis outcomes across subjects and different breathing patterns. Bland and Altman plot statistics for measurements of fB. TI. and TE using best component-selection methods (“Area” and “Peak”), PCA-fusion method and “Ideal” component and, in supine (n = 74) and seated (n = 74) position.
Table 4. Agreement analysis outcomes across subjects and different breathing patterns. Bland and Altman plot statistics for measurements of fB. TI. and TE using best component-selection methods (“Area” and “Peak”), PCA-fusion method and “Ideal” component and, in supine (n = 74) and seated (n = 74) position.
τp-ValueHeteroscedastic?Fixed Bias a/OLSLOA c/V-Shape Limits d
fB supineArea0.1590.045Yesy = −0.054x + 2.316 bUCL: y = 0.085x + 10.907 d
LCL: y = −0.192x − 6.275
Peak0.1420.074No1.380 aFrom −11.95 to 14.72 c
PCA0.2110.008Yesy = 0.008x + 0.130 bUCL: y = 0.054x + 2.299 d
LCL: y = −0.038x − 2.039
Ideal0.0380.631No0.884 aFrom −4.171 to 5.940 c
fB seatedArea0.1420.074No0.084 aFrom −9.635 to 9.803 c
Peak0.1320.096No−0.121 aFrom −9.931 to 9.688 c
PCA0.1080.174No−0.23 aFrom −5.474 to 5.010 c
Ideal0.1960.014Yesy = −0.021x + 0.597 bUCL: y = 0.028x + 2.057 d
LCL: y = −0.071x − 0.864
TI supineArea0.3020.000Yesy = 0.084x − 0.019 bUCL: y = 0.618x + 0.095 d
LCL: y = −0.450x − 0.132
Peak0.3340.000Yesy = 0.104x − 0.021 bUCL: y = 0.638x + 0.093 d
LCL: y = −0.430x − 0.135
PCA0.3750.001Yesy = 0.283x − 0.175 bUCL: y = 0.926x − 0.840 d
LCL: y = −0.390x + 0.354
Ideal0.2920.000Yesy = −0.090x + 0.121 bUCL: y = 0.158x + 0.163 d
LCL: y = −0.338x + 0.078
TI seatedArea0.4300.000Yesy = 0.1022x − 0.0141 bUCL: y = 0.834x + 0.112 d
LCL: y = −0.618x − 0.409
Peak0.4220.000Yesy = 0.220x − 0.075 bUCL: y = 0.642x − 0.173 d
LCL: y = −0.438x + 0.197
PCA0.4890.000Yesy = 0.171x − 0.0332 bUCL: y = 0.7182x − 0.226 d
LCL: y = −0.377x + 0.160
Ideal0.3130.000Yesy = −0.059x + 0.129 bUCL: y = 0.211x + 0.069 d
LCL: y = −0.329x + 0.189
TE supineArea0.4210.000Yesy = −0.170x + 0.166 bUCL: y = 0.508x + 0.358 d
LCL: y = −0.847x − 0.026
Peak0.4050.000Yesy = −0.138x + 0.144 bUCL: y = 0.496x + 0.306 d
LCL: y = −0.771x − 0.017
PCA0.5220.000Yesy = −0.209x + 0.328 bUCL: y = 0.667x + 0.037 d
LCL: y = −1.084x + 0.620
Ideal0.4840.000Yesy = −0.153x + 0.148 bUCL: y = 0.3987x − 0.185 d
LCL: y = −0.705x + 0.481
TE seatedArea0.3840.000Yesy = −0.216x + 0.364 bUCL: y = 0.303x + 0.532 d
LCL: y = −0.735x + 0.197
Peak0.3960.000Yesy = −0.231x + 0.413 bUCL: y = 0.226x + 0.666 d
LCL: y = −0.657x + 0.101
PCA0.4220.000Yesy = −0.127x + 0.320 bUCL: y = 0.284x + 0.498 d
LCL: y = −0.538x + 0.142
Ideal0.3160.000Yesy = −0.058x + 0.054 bUCL: y = 0.383x + 0.337 d
LCL: y = −0.500x − 0.228
τ: Kendall’s τ correlation coefficient and relative p-value (heteroscedasticity test). a Fixed Bias. obtained as the mean of differences (device – OEP). for homoscedastic data. b OLS: ordinary least squares line of best fit (proportional bias) for heteroscedastic data. c LOA: limits of agreement. computed as mean difference ± 1.96SD (for homoscedastic data). d V-shape limits: UCL and LCL 95% confidence limits. calculated according to Bland [44].

Share and Cite

MDPI and ACS Style

Cesareo, A.; Previtali, Y.; Biffi, E.; Aliverti, A. Assessment of Breathing Parameters Using an Inertial Measurement Unit (IMU)-Based System. Sensors 2019, 19, 88. https://doi.org/10.3390/s19010088

AMA Style

Cesareo A, Previtali Y, Biffi E, Aliverti A. Assessment of Breathing Parameters Using an Inertial Measurement Unit (IMU)-Based System. Sensors. 2019; 19(1):88. https://doi.org/10.3390/s19010088

Chicago/Turabian Style

Cesareo, Ambra, Ylenia Previtali, Emilia Biffi, and Andrea Aliverti. 2019. "Assessment of Breathing Parameters Using an Inertial Measurement Unit (IMU)-Based System" Sensors 19, no. 1: 88. https://doi.org/10.3390/s19010088

APA Style

Cesareo, A., Previtali, Y., Biffi, E., & Aliverti, A. (2019). Assessment of Breathing Parameters Using an Inertial Measurement Unit (IMU)-Based System. Sensors, 19(1), 88. https://doi.org/10.3390/s19010088

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop