Abstract
Importance Pulse oximetry, a ubiquitous vital sign in modern medicine, has inequitable accuracy that disproportionately affects Black and Hispanic patients, with associated increases in mortality, organ dysfunction, and oxygen therapy. Although the root cause of these clinical performance discrepancies is believed to be skin tone, previous retrospective studies used self-reported race or ethnicity as a surrogate for skin tone.
Objective To determine the utility of objectively measured skin tone in explaining pulse oximetry discrepancies.
Design, Setting, and Participants Admitted hospital patients at Duke University Hospital were eligible for this prospective cohort study if they had pulse oximetry recorded up to 5 minutes prior to arterial blood gas (ABG) measurements. Skin tone was measured across sixteen body locations using administered visual scales (Fitzpatrick Skin Type, Monk Skin Tone, and Von Luschan), reflectance colorimetry (Delfin SkinColorCatch [L*, individual typology angle {ITA}, Melanin Index {MI}]), and reflectance spectrophotometry (Konica Minolta CM-700D [L*], Variable Spectro 1 [L*]).
Main Outcomes and Measures Mean directional bias, variability of bias, and accuracy root mean square (ARMS), comparing pulse oximetry and ABG measurements. Linear mixed-effects models were fitted to estimate mean directional bias while accounting for clinical confounders.
Results 128 patients (57 Black, 56 White) with 521 ABG–pulse oximetry pairs were recruited, none with hidden hypoxemia. Skin tone data was prospectively collected using 6 measurement methods, generating 8 measurements. The collected skin tone measurements were shown to yield differences among each other and overlap with self-reported racial groups, suggesting that skin tone could potentially provide information beyond self-reported race. Among the eight skin tone measurements in this study, and compared to self-reported race, the Monk Scale had the best relationship with differences in pulse oximetry bias (point estimate: −2.40%; 95% CI: −4.32%, - 0.48%; p=0.01) when comparing patients with lighter and dark skin tones.
Conclusions and relevance We found clinical performance differences in pulse oximetry, especially in darker skin tones. Additional studies are needed to determine the relative contributions of skin tone measures and other potential factors on pulse oximetry discrepancies.
Question Can skin tone capture information beyond race to help explain pulse oximetry discrepancies?
Findings Pulse oximetry bias across races seems to persist across skin tone when measured using administered visual scales, reflectance colorimetry, or reflectance spectrophotometry. Among the eight skin tone measurements in this study, and compared to self-reported race, the Monk Scale seemed to best correlate with pulse oximetry bias when comparing patients with lighter and dark skin tones.
Meaning Compared to self-reported race, skin tone is associated with some pulse oximetry discrepancies; we recommend using skin tone to assist the regulatory clearance of equitable pulse oximeters.
Background
Racial and ethnic bias in pulse oximetry stands out as a quintessential health inequity, whereby the same medical devices that guide clinical decision-making may fail to function equally well for all patients.1 The reliability of pulse oximetry has been a reason for concern for decades,2–6 but it was not until the COVID-19 pandemic when Sjoding and colleagues’ seminal paper reported racial bias in pulse oximetry measurements that pulse oximetry became a health equity issue.7
Followed by other studies,8–15 oxygen saturation measured by pulse oximetry (SpO2) is widely reported to overestimate the “true” arterial oxygen (SaO2), measured by arterial blood gas (ABG), disproportionately affecting Black and Hispanic patients. A seemingly small discrepancy is associated with higher rates of “hidden hypoxemia” among these patients,2, 8, 12 with associated inequities in oxygen therapies 9, 16 and increases in mortality and organ dysfunction.8
Pulse oximeters estimate arterial oxygen saturation by measuring light absorption of oxyhemoglobin and deoxyhemoglobin in capillary blood.17, 18 Previous studies have shown that skin tone can independently affect light absorption, causing discrepant readings, especially among darker-skinned individuals. 19–21As such, previous retrospective studies share a fundamental limitation: self-reported race or ethnicity is used as a surrogate for skin tone, although the root cause of these discrepancies is believed to be skin tone.22
In this cohort study, we prospectively collected skin tone data from critically ill patients in various body locations using different devices. We paired this data with pulse oximetry measurements, ABG, and other Electronic Health Records (EHR) data to investigate the utility of skin tone data in explaining pulse oximetry performance.
As a pilot study, our objectives were 2-fold: first, to provide a framework to conduct larger clinical studies that assess the association between different skin tone measurements and pulse oximetry discrepancies; and second, to provide evidence that can support recent discussions from the Food and Drug Administration (FDA)23–26 in pursuit of guidelines to evaluate pulse oximetry performance in a more inclusive spectrum of patients.
Pulse oximetry performance is commonly assessed by the FDA using the accuracy root mean square (ARMS), with a threshold of ARMS ≤ 3.0% for transmittance pulse oximeters between SaO2 and SpO2 measurements ranging from 70–100%.27 ARMS can be decomposed into the average difference in magnitude between SaO2 and SpO2, often referred to as “bias”, and the variability of those differences, which can be measured as the standard deviation of the difference between SaO2 and SpO2.
Methods
This study was approved by the Duke Health IRB under Pro00110842, following the American Medical Association’s recommendations on health equity language and adhering to the STROBE statement.28, 29
Cohort Selection
Patients admitted to an adult intensive care unit (ICU) at Duke University Hospital were screened. Standard-of-care ABG up to 5 minutes after a pulse oximetry measurement was required for eligibility, resulting in SaO2–SpO2 pairs.
Exclusion criteria included unremovable fingernail polish, admission for a vascular complication (e.g., grafting or stenting), amputation, and large areas of skin discoloration where the accuracy of skin tone measurements could be affected. Pairs containing either a SaO2 or a SpO2 measurement out of the 70–100% range were excluded.8 27
Data Collection and Processing
All data and patients’ consent were stored in Duke Health’s REDCap, with data processing performed in Python 3.10.
The mathematical definitions of mean directional bias, variability of bias, and ARMS are in Supplemental Formulas 2, 3, and 4.
Skin Tone
Three types of skin tone assessment were conducted using different devices: administered visual scales (Fitzpatrick, Monk30, and Von Luschan scales, visible in Supplemental Figures 6a, 6b and 6c); reflectance colorimetry (Delfin SkinColorCatch); and reflectance spectrophotometry (Variable Spectro 1 Pro Bridge Set and Konica Minolta CM-700D Spectrophotometer).
All the skin tone data are collected within 7 days of the SaO2–SpO2 pairs and with controlled lighting to ensure reproducibility. Using the L*, individual typology angle (ITA) and Melanin Index (Melanin Index) color spaces, eight different skin tone measures were collected in this study, as detailed in Supplemental Table 2. Further details are demonstrated in Supplemental Text.
Data merging
Pulse oximetry values and ABG panel data were merged into SaO2-SpO2 pairs and recorded in REDCap. Demographic data was merged from the EHR system. Three race groups were defined, including “Black”, “Other” and “White” patients. The group “Other” captures minority patients who self-identify as Asian, American Indian / Alaskan natives, more than two races, and unknown race–groups that separately represent a small proportion of patients. Vital signs captured within 4 hours prior to the SaO2–SpO2 pair were merged from the EHR system. Mean arterial pressure (MAP) from the arterial line was preferred when available, otherwise cuff values were used. Laboratory test values from the previous 24 hours, relative to the SaO2– SpO2 pair, were merged (listed in Supplemental Table 1).
Missingness
Missing data occurred occasionally in two skin tone measurements, Variable L* and Konica Minolta L*, due to technical issues or patient refusal. Missingness rates for skin tone measurements can be found in Table 1. In the merged EHR clinical data, missingness occurred in vital signs and laboratory test values when no value was found within the set windows. Missingness rates for these covariates are in Supplemental Table 1. Patients with missing data that would be necessary for modeling were dropped from the study, and sensitivity analyses were run to assess the robustness of this design choice.
Measurement variability
The standard deviation was computed across the different values of the same measure and location, and compared with the average standard deviations across all locations. Average standard deviations across palm and finger locations were also computed, as depicted in Supplemental Table 3.
Statistical Analysis
Exploratory data analysis was performed using Python 3.1031, as described in the Supplemental Text, and summarized using the tableone package,32 in Table 1. For each tertile of skin tone, mean directional bias, variability of bias, and ARMS were computed, as reported in Figure 2. Statistical analysis was conducted in R 4.3.1,33 using the packages nlme for the mixed-effects analysis.34, 35
Linear mixed-effects models
Linear mixed-effects models with patient identifiers as a random effect were fitted and adjusted for potential confounders (race and clinical features). Pairs with missing data were dropped for the analysis.
As a baseline, we built a model to assess the effect of self-reported race in pulse oximetry mean directional bias, adjusting for pH, SaO2, heart rate (HR), and MAP (Supplemental Formula 5). The following models, documented in Supplemental Formula 6, included these same covariates, as well as the skin tone variables, separate per model (eight models were built in total to assess the individual effect of each skin tone variable).
Lastly, to investigate the combined effect of all the skin tone measurements on pulse oximeter mean directional bias, we fitted two linear mixed-effects models (Supplemental Formula 7 and Supplemental Formula 8). The first model included six skin tone variables, excluding Konica Minolta L* and Variable L* due to missingness. The second model included all eight skin tone measurements, as a sensitivity analysis. These differences in design resulted in a lower sample size for the second model. Table 2 summarizes the built models. All the significance levels in linear mixed-effects models are calculated using likelihood ratio tests (LRT) using chi-squared statistics.
Results
Cohort Characteristics
From January 1, 2023 to June 30, 2023, a total of 1,167 admitted inpatients with qualifying SaO2-SpO2 pairs were screened at Duke University Hospital. Out of 301 patients who met our inclusion criteria and were approached, 134 patients consented to this study (see the flow diagram in Figure 1). After 6 exclusions due to withdrawal, missing location data, or incomplete data, 128 patients were considered for analysis (39.8% female, 43% Black, see Table 1).
After excluding readings that did not fall into the 70–100% range, a total of 521 SaO2–SpO2 pairs were obtained from this cohort. SpO2 values ranged from 82% to 100%, and SaO2 values ranged from 83.8% to 99.0%. The difference between SaO2 and SpO2 ranged from −9.0% to 8.8% – see Supplemental Table 1 for further pair-level characteristics.
Measurement variability on skin tone scales
Supplemental Table 3 shows that objective scales resulted in lower standard deviations when compared to subjective scales. Palm averages are found to be more stable, when compared to other locations, supporting the design choice of taking palm averages as the preferred measurement for subsequent analyses.
Skin tone and race data
Exploratory data analysis showed that the skin tone measurements yielded differences among each other and overlap with self-reported racial groups. Figure 2 depicts the unadjusted mean directional bias, variability of bias, and ARMS, per skin tone tertile. Further detailed text can be found in the Supplemental Text.
Linear mixed effect model on mean directional bias
To understand the effect of self-reported race and each of the skin tones on pulse oximetry discrepancy, the fitted linear mixed-effects model did not find race to have a significant effect (□2 = 0.90, p = 0.64), Nevertheless, the obtained coefficient is in the expected direction (- 0.23%, 95% CI: −0.76%, 0.30% for Black patients). Among eight skin tone measurements, only the Monk scale yielded a significant effect on the mean directional bias (point estimate: −2.40%; 95% CI: −4.32%, −0.48%; p-value: 0.01). – see Table 2 for a full report of the models.
To examine the combined effect of the six skin tone variables, excluding Konica Minolta L* and Variable L* due to missing data (the goal being to maximize sample size), the LRT (Supplemental Formula 7) rejected the null hypothesis (□2 = 14.98, p = 0.02). This suggests that at least one of these six skin tone measurements affects the mean directional bias. However, in a sensitivity analysis that included all eight skin tone variables, the LRT (Supplemental Formula 8) showed borderline insignificance (□2 = 14.86, p = 0.06)
Discussion
The objective of this study was to investigate the relationship between skin tone measurements and pulse oximetry discrepancies. We prospectively collected skin tone data from 128 critically-ill patients, comprising a total of 521 pairs of pulse oximetry-ABG data, and leveraging six different tools across sixteen body sites. We addressed the fundamental limitation of previous studies on pulse oximetry racial discrepancies,7–10, 12–14, 16 which solely relied on racial or ethnic groups as a proxy for skin tone. As a prospective cohort study where we enrolled patients in a comprehensive screening process, we minimized assumptions associated with secondary data analysis36 and obtained a cohort with equal representation of Black and White patients.
The collected skin tone measurements were shown to yield differences among scales and when compared to self-reported race, suggesting that skin tone data carries information beyond self-reported race. When assessing the relation of skin tone with pulse oximetry bias, racial and ethnic disparities7 do seem to persist, whereby darker patients show a higher degree of mean directional bias (see Figure 2 and Table 2). These findings are aligned with a recent similar report.37 As opposed to a model that solely relies on self-reported race (and clinical confounders) to account for pulse oximetry bias – where the effect of self-reported race is found to be nonsignificant – models that accounted for skin tone found at least one of the measures to be significant (Supplemental Formula 7, results in Table 2). This finding suggests that skin tone is related, beyond self-reported race, to pulse oximetry bias.
In this study, we tested different devices for skin tone assessment that ranged from a negligible cost (color-printed scales) to thousands of dollars (Konica Minolta’s spectrophotometer). The measurement variability was non-negligible, but lower for objective scales, as expected. However, in exploring pulse oximetry bias, the Monk scale – designed to be easy-to-use for evaluation of technology, while representing a broad range of skin tones 30, 38 – was found to yield the strongest association with pulse oximetry bias (Table 2). As this study does not show clear evidence that more sophisticated and expensive devices (colorimetry or spectrophotometry) add value in this application of skin tone characterization, we believe that further investigation is necessary, including a broader range of skin tone measurement devices and a larger sample size.
In the exploratory analysis, both across self-reported race and skin tone measurements, darker-pigmented patients observed a lower SaO2–SpO2 variability (see Supplemental Table 4 and Figure 2). To assess precision while adjusting for confounding, we identified the within-stratum variation by modeling heterogeneous variance in the linear mixed-effects models (Supplemental Formula 9) using the tertiles defined in Figure 2. Supplemental Table 6 lists the built models, where the within-stratum variance is shown to be consistently lower in the darkest tertile of all eight skin tone measurements, and higher in the lightest tertile, except for Delfin Melanin Index. Consequently, the darker-pigmented patients of our cohort seem to have more consistently wrong pulse oximetry readings, despite having a lower ARMS. Although this finding is not aligned with a recent report,37 concerns about the reliability of ARMS among different racial groups in pulse oximetry performance evaluation and potential inequitable treatment delay had already been raised before.39 Considering the far-reaching implications of these findings in pulse oximetry regulation, further investigation is necessary.
Considering the prohibitive cost of replacing existing pulse oximeters,40 our work stands as a fundamental milestone to any interim solution that may tackle pulse oximetry inaccuracies leveraging existing technologies. If the finding that pulse oximeters are “more consistently wrong” among darker-skinned patients stands in follow-up studies, one could argue that clinical algorithms that perform a holistic correction – and not simple race corrections41 – could be more attainable, due to the observed lower variance among darker-skinned patients. Given this finding, which suggests that skin tone may provide additional information beyond race, we propose that incorporating skin tone measurements could help mitigate residual confounding in algorithms solely reliant on self-reported race.
Implications in regulation and pulse oximetry clearance
In response to FDA’s recent discussions on pulse oximetry performance discrepancies,23–25 we believe that this pilot study provides initial evidence to support the suggested need of thoughtfully collecting and assessing skin tone data in pulse oximetry clearances. Besides being an important factor in pulse oximetry miscalibration, and a more objective measure than self-reported race, skin tone data seems to yield utility in pulse oximetry discrepancies. Consequently, besides requiring racial and ethnic diversity for pulse oximetry clearance,26 we recommend the FDA to require the quantification and representation of a full spectrum of skin tones, while not disregarding the potential impact of other unmeasured confounders. Recognizing Beer-Lambert’s law’s impact on light transmission, we underline the importance of assessing other potential confounding variables such as perfusion, skin thickness, systemic vascular resistance, or local vascular resistance, for which further investigation is necessary as these are not commonly measured in medicine.
Moreover, our findings on increased variability of bias and better ARMS among darker-skinned patients, despite worsened mean directional bias, raise questions about FDA’s conventional reliance on ARMS for pulse oximetry clearance. These considerations require larger follow-up studies and are not necessarily aligned with other reports.37 Considering that bias in the direction of overestimation of SaO2 may carry more downstream clinical harm than bias in the opposite direction, we would like to build upon previous concerns39 and bring to debate the question: “What is an equitable performance assessment metric for pulse oximetry clearance?”.
Limitations and Future Work
As our skin tone data presented non-negligible measurement variability across sites and examiners, we considered the average of the left and right, dorsal and ventral palm readings in 125 patients, out of 128, to obtain more stable skin tone measurements. Due to our limited sample size, most of the findings are not significant despite a reasonable effect size. In the future, we would like to run a larger study with more patients. Although this study’s cohort had over 40% Black patients, the darkest skin tones are rarely observed, which might be due to the population skin tone distribution of the community where our study was based. Moreover, we would like to enroll more patients with hypoxemia (i.e., SaO2 < 88%), to potentially investigate the impact of skin tone in hidden hypoxemia phenomena. Additionally, we would like to examine other potential covariates that may contribute to pulse oximetry disparities. Finally, as our prospective study suggests that, despite its effect, skin tone is unlikely to be the sole contributor to pulse oximetry discrepancies, we advocate for the need of further investigation on other unmeasured confounders.
Conclusion
This pilot study analyzed skin tone measurements with pulse oximetry performance discrepancies among critically ill patients. We prospectively collected skin tone assessments via administered visual scales, reflectance colorimetric, and spectrophotometric devices. Pulse oximetry varied across skin pigmentation and, similarly to previous reports, darker-skin-toned patients yielded a greater bias, independently of clinical confounders. However, with a large variation in pulse oximetry data, skin tone is unlikely to be the sole contributor to performance discrepancies in pulse oximetry. While these findings necessitate a larger sample size to be further validated and select the best method(s) for skin tone measurement, we hope this paper provides a framework for future similar studies, as well as initial evidence to support FDA’s discussions on regulation changes towards more equitable pulse oximeters.
Data Availability
An IRB revision is underway to evaluate this possibility but legal counsel has suggested that re-consent may be necessary for submission to a public access-controlled archive.
Footnotes
↵* co-first authors
Conflicts of interest AIW holds equity and management roles in Ataia Medical.
Funding AIW is supported by the Duke CTSI by the National Center for Advancing Translational Sciences (NCATS) of the National Institutes of Health under UL1TR002553 and REACH Equity under the National Institute on Minority Health and Health Disparities (NIMHD) of the National Institutes of Health under U54MD012530.
JWG is a 2022 Robert Wood Johnson Foundation Harold Amos Medical Faculty Development Program and declares support from RSNA Health Disparities grant (#EIHD2204), Lacuna Fund (#67), Gordon and Betty Moore Foundation, and NIH (NIBIB) MIDRC grant under contracts 75N92020C00008 and 75N92020C00021