Paul S. Addison, Rui Wang, Alberto A. Uribe, Sergio D. Bergese, "On Better Estimating and Normalizing the Relationship between Clinical Parameters: Comparing Respiratory Modulations in the Photoplethysmogram and Blood Pressure Signal (DPOP versus PPV)", Computational and Mathematical Methods in Medicine, vol. 2015, Article ID 576340, 6 pages, 2015. https://doi.org/10.1155/2015/576340
On Better Estimating and Normalizing the Relationship between Clinical Parameters: Comparing Respiratory Modulations in the Photoplethysmogram and Blood Pressure Signal (DPOP versus PPV)
DPOP (ΔPOP or Delta-POP) is a noninvasive parameter which measures the strength of respiratory modulations present in the pulse oximeter waveform. It has been proposed as a noninvasive alternative to pulse pressure variation (PPV) used in the prediction of the response to volume expansion in hypovolemic patients. We considered a number of simple techniques for better determining the underlying relationship between the two parameters. It was shown numerically that baseline-induced signal errors were asymmetric in nature, which corresponded to observation, and we proposed a method which combines a least-median-of-squares estimator with the requirement that the relationship passes through the origin (the LMSO method). We further developed a method of normalization of the parameters through rescaling DPOP using the inverse gradient of the linear fitted relationship. We propose that this normalization method (LMSO-N) is applicable to the matching of a wide range of clinical parameters. It is also generally applicable to the self-normalizing of parameters whose behaviour may change slightly due to algorithmic improvements.
An example of the least-squares fit for the DPOP-PPV data pairs is shown in Figure 2. However, we suspect that this is not the best estimate of the true relationship between DPOP and PPV as examples of tight linear DPOP-PPV relationships may be observed in the literature [3, 8, 10]. These other reported studies often involve manually selected data, where great lengths have been taken to ensure the best quality data. Figure 2, however, represents all data acquired during the clinical trial. This is more representative of the data that will be fed into a pulse oximeter in practice and hence contains significant noise. It is this kind of data we have a specific interest in within our development of robust parameters for use in medical devices. We also know that when there are no respiratory modulations present in the blood pressure (BP) waveform, there will be no respiratory modulations in the pulse oximeter waveform. Hence, assuming no other signal components, the relationship should pass through the origin. The regression line plotted in Figure 2 clearly does not do this, and the standard least-squares fitted line appears to be dragged upwards in the region of the origin due to the asymmetry in signal noise components.
In an attempt to better understand the asymmetry in the distribution of noise in the DPOP measure, we examined more closely segments of signal where such errors occurred. Often they coincided with significant up slopes or down slopes in the baseline due to vasomotion. An example of this baseline slope is shown in Figure 3 where long-term (i.e., longer than the respiratory cycle) increases or decreases in the baseline are observed. We examined the effect of persistent gradients in the baseline on the calculated DPOP. A simple numerical model was developed whereby a uniformly occurring distribution of amplitude modulations ranging from 0 to and a range of baseline gradient induced amplitude measurement errors (as measured from the trough to pulse peak shown in Figure 1) from to were input into (1). The plot of the resulting distribution of errors in DPOP is shown in Figure 4(a). The asymmetry is obvious in the plot. This may be compared to the actual distribution of DPOP errors from the least-median-of-squares line shown in Figure 4(b). There is a resemblance between the two figures with both possessing a distinctly positive skew. Although a very simple gradient induced noise model was used, the parametric bounds of the model can be varied quite widely (in terms of amplitude modulations and gradients) with the same conclusion: an asymmetric distribution with a positive skew results when gradient noise is introduced into the DPOP measure. Note that other noise will be present in the experimental data (motion, drug effects, etc.) and hence we do not expect an exact match with the synthetic data. In practice other noise sources could include those of a physiological nature (e.g., venous blood movement [16–18], autonomic modulation , motion artefact , and arrhythmia ) and device specific noise (e.g., amplification and filtering effects  or arterial line flushes ).
2.3. Regression Line Approximations of the DPOP-PPV Relationship
Given that the data set shown in Figure 2 is relatively noisy, exhibiting the asymmetry described above, the question is, how may we numerically approximate this relationship line better to account for the effects of the asymmetric spread of data points and thus improve the representation of the underlying association? We propose a number of steps in this regard.
First, we tackled the intent to have the fitted line pass through the origin. For least-squares linear regression this was actually a simple task. The mathematics of traditional least-squares regression was reworked with the extra criterion that the line intercept equalled zero . This resulted in the modified regression line shown in Figure 5(a). We see that such an approximation appears to fit the original data better in the vicinity of the origin. However, as noted above, the noise manifests itself predominantly in terms of a positive error in DPOP which tends to drag the fitted line upwards from the suspected relationship. Many of these positive errors in DPOP are, in fact, quite large and thus dominate the least-squares technique. In an attempt to counteract the disproportional effect such outliers have on the traditional least-squares method we applied a least-median-of-squares (LMS) regression  to the data. The corresponding regression line for the data set is shown in Figure 5(b). We can see that this line appears to better approximate the suspected true relationship. Finally, we merged these two techniques, combining the desire to both account for outliers and pass through the origin. Figure 5(c) plots an LMS regression line which has been forced to pass through the origin (the “LMSO” line).
2.4. Normalizing the Fluid Responsiveness Parameters
The steps described in the previous section were aimed at improving our estimate of the underlying relationship between the two parameters. However, in practice, either one or the other parameter is used to determine the fluid responsiveness of the patient (they are not normally used at the same time). The determination of fluid responsiveness is performed by comparing the parameter value against a predefined threshold, for example, 13% for PPV [3, 10, 19, 23, 24]. If the patient’s PPV value is less than the threshold, then the patient may be deemed unresponsive to volume expansion and a decision was made not to administer fluid. Alternatively, if the patient’s PPV value is greater than the threshold then this may indicate that the patient will respond to volume expansion. However, it is well known that this threshold is different for PPV and DPOP. This is also evidenced by the fact that these two parameters do not have 1 : 1 correspondence in the current data. This lack of parity is to be expected as the signals are quite different, with additional biomechanics at play between the BP generated pulse in the fluid (blood) column at the arm, where PPV is measured, and the pressure-mechanical coupling through the vessel walls and body tissue matrix “downstream” at the peripheral finger site where the pleth is measured [25, 26].
However, now that we have determined the linear relationship that passes through the origin (the LMSO method of Figure 5(c)), it is a trivial task to normalize the two measures. And, as we require DPOP to be the proxy for PPV, we rescale the DPOP values to match their corresponding PPV values. This is performed by multiplying each DPOP value by the inverse of the gradient of our best fit line. This rescaling, or normalizing, of the data using the LMSO line of Figure 5(c) is shown in Figure 6 (the “LMSO-N” method).
We have described a number of simple steps to improve the determination of the underlying relationship between two clinical parameters. These involved a least-median-of-squares estimator that is more robust to outlying data points and incorporating the additional criterion that the relationship must pass through the origin. The noise appeared asymmetric in nature and although we had no reference truth in this study, we know from the literature that tight linear DPOP-PPV relationships might be expected where great lengths have been taken to ensure the best quality data. If the DPOP parameter is being used to determine whether to administer fluids, then an error could possibly cause the wrong action to be taken (all medical monitoring devices exhibit some error and hence may be susceptible to this kind of issue to some degree). Whether a higher or lower DPOP value occurs at a threshold value of PPV (e.g., 13% ), indicative of the threshold for administering fluids, depends on the rotation and translation of the best fit line according to the regression and normalization used, which itself is data dependent. Further work on larger data sets is required to fully understand this effect. This area is ripe for more sophisticated analysis of this kind described here and the proposed method is perhaps a first step in this regard, although the authors do feel that it strikes a balance between improvement and oversophistication.
In addition, we have detailed a technique for normalizing the parameters using the inverse of the gradient of the best fit line. This method of normalizing DPOP in terms of PPV is advantageous in that the clinician may change between the two parameters seamlessly as there is no need for separate parameter-dependent thresholds. The LMSO-N method is also useful in self-normalizing any parameter against itself. For example, the method could be used to normalize the algorithm output of the parameter to the previous version of the same parameter, thus maintaining consistency across algorithm revisions.
Although we have discussed our methods in terms of the relationship between DPOP and PPV, other relationships both in terms of fluid responsiveness parameters (e.g., stroke volume variation and systolic pressure variation) and other medical parameters may benefit from this technique.
- M. Feissel, J.-L. Teboul, P. Merlani, J. Badie, J.-P. Faller, and K. Bendjelid, “Plethysmographic dynamic indices predict fluid responsiveness in septic ventilated patients,” Intensive Care Medicine, vol. 33, no. 6, pp. 993–999, 2007.
- G. A. Westphal, E. Silva, A. R. Gonçalves, M. Caldeira Filho, and L. F. Poli-de-Figueiredo, “Pulse oximetry wave variation as a noninvasive tool to assess volume status in cardiac surgery,” Clinics, vol. 64, no. 4, pp. 337–343, 2009.
- L. Ø. Høiseth, I. E. Hoff, Ø. Skare, K. A. Kirkebøen, and S. A. Landsverk, “Photoplethysmographic and pulse pressure variations during abdominal surgery,” Acta Anaesthesiologica Scandinavica, vol. 55, no. 10, pp. 1221–1230, 2011.
- J. R. Chandler, E. Cooke, C. Petersen et al., “Pulse oximeter plethysmograph variation and its relationship to the arterial waveform in mechanically ventilated children,” Journal of Clinical Monitoring and Computing, vol. 26, no. 3, pp. 145–151, 2012.
- M. Bland, An Introduction to Medical Statistics, Oxford University Press, Oxford, UK, 1997.
- P. S. Addison, R. Wang, A. A. Uribe, and S. D. Bergese, “Increasing signal processing sophistication in the calculation of the respiratory modulation of the photoplethysmogram (DPOP),” Journal of Clinical Monitoring and Computing, 2014.
- P. S. Addison, R. Wang, S. J. McGonigle, A. A. Uribe, and S. D. Bergese, “Calculation of the respiratory modulation of the photoplethysmogram (DPOP) incorporating a correction for low perfusion,” Anesthesiology Research and Practice, vol. 2014, Article ID 980149, 9 pages, 2014.
- J. P. Phillips, A. Belhaj, K. Shafqat, R. M. Langford, K. H. Shelley, and P. A. Kyriacou, “Modulation of finger photoplethysmographic traces during forced respiration: venous blood in motion?” in Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '12), pp. 3644–3647, San Diego, Calif, USA, 2012.
- K. H. Shelley, D. Tamai, D. Jablonka, M. Gesquiere, R. G. Stout, and D. G. Silverman, “The effect of venous pulsation on the forehead pulse oximeter wave form as a possible source of error in SpO2 calculation,” Anesthesia and Analgesia, vol. 100, no. 3, pp. 743–747, 2005.
- K. H. Shelley, M. Dickstein, and S. M. Shulman, “The detection of peripheral venous pulsation using the pulse oximeter as a plethysmograph,” Journal of Clinical Monitoring, vol. 9, no. 4, pp. 283–287, 1993.
- S. A. Landsverk, L. O. Hoiseth, P. Kvandal, J. Hisdal, O. Skare, and K. A. Kirkeboen, “Poor agreement between respiratory variations in pulse oximetry photoplethysmographic waveform amplitude and pulse pressure in intensive care unit patients,” Anesthesiology, vol. 109, no. 5, pp. 849–855, 2008.
- C. Sandroni, F. Cavallaro, C. Marano, C. Falcone, P. de Santis, and M. Antonelli, “Accuracy of plethysmographic indices as predictors of fluid responsiveness in mechanically ventilated adults: a systematic review and meta-analysis,” Intensive Care Medicine, vol. 38, no. 9, pp. 1429–1437, 2012.
- J. G. Eisenhauer, “Regression through the origin,” Teaching Statistics, vol. 25, no. 3, pp. 76–80, 2003.
- D. L. Massart, L. Kaufman, P. J. Rousseeuw, and A. Leroy, “Least median of squares: a robust method for outlier and model error detection in regression and calibration,” Analytica Chimica Acta, vol. 187, pp. 171–179, 1986.
- F. L. Poli de Figueiredo, E. Silva, and M. Rocha, “Pulse oximetry wave respiratory variations for the assessment of volume status in patients under mechanical ventilation,” Journal of Critical Care, vol. 32, article A96, 2004.
- G. Natalini, A. Rosano, M. E. Franceschetti, P. Facchetti, and A. Bernardini, “Variations in arterial blood pressure and photoplethysmography during mechanical ventilation,” Anesthesia and Analgesia, vol. 103, no. 5, pp. 1182–1188, 2006.
- M. Briet, P. Boutouyrie, S. Laurent, and G. M. London, “Arterial stiffness and pulse pressure in CKD and ESRD,” Kidney International, vol. 82, no. 4, pp. 388–400, 2012.
- T. Khamdaeng, J. Luo, J. Vappou, P. Terdtoon, and E. E. Konofagou, “Arterial stiffness identification of the human carotid artery using the stress-strain relationship in vivo,” Ultrasonics, vol. 52, no. 3, pp. 402–411, 2012.