- Research article
- Open Access
- Open Peer Review
Adherence to oxygenation and ventilation targets in mechanically ventilated premature and sick newborns: a retrospective study
BMC Pediatricsvolume 13, Article number: 126 (2013)
Ventilator treatment exposes newborns to both hyperoxemia and hyperventilation. It is not known how common hyperoxemia and hyperventilation are in neonatal intensive care units in Norway. The purpose of this study was to assess the quality of current care by studying deviations from the target range of charted oxygenation and ventilation parameters in newborns receiving mechanical ventilation.
Single centre, retrospective chart review that focused on oxygen and ventilator treatment practices.
The bedside intensive care charts of 138 newborns reflected 4978 hours of ventilator time. Arterial blood gases were charted in 1170 samples. In oxygen-supplemented newborns, high arterial pressure of oxygen (PaO2) values were observed in 87/609 (14%) samples. In extremely premature newborns only 5% of the recorded PaO2 values were high. Low arterial pressure of CO2 (PaCO2) values were recorded in 187/1170 (16%) samples, and 64 (34%) of these were < 4 kPa. Half of all low values were measured in extremely premature newborns. Tidal volumes above the target range were noted in 22% of premature and 20% of full-term newborns.
There was a low prevalence of high PaO2 values in premature newborns, which increased significantly with gestational age (GA). The prevalence of low PaCO2 values was highest among extremely premature newborns and decreased with increasing GA. Further studies are needed to identify whether adherence to oxygenation and ventilation targets can be improved by clearer communication and allocation of responsibilities between nurses and physicians.
Clinical practice with respect to ventilator management, administration of oxygen, and assessment of oxygenation differs greatly among neonatal intensive care units (NICUs) [1, 2]. Ventilator treatment exposes newborns to both hyperoxemia and hyperventilation. The goal of ventilator treatment is to balance gas exchange while minimising trauma to the lung tissue . Adjusting oxygenation and ventilator therapy is challenging, and improved strategies are needed to minimise hyperoxemia  and hyperventilation with hypocarbia [5, 6] in preterm and full-term newborns. Appropriate oxygenation is achieved by titrating the fraction of inspired oxygen (FiO2) and the mean airway pressure (MAP). The aim of appropriate ventilation is to maintain an arterial pressure of carbon dioxide (PaCO2) of ~5.3 kPa (40 mm Hg) .
Nurses need a target range for oxygen saturation (SpO2) in order to titrate FiO2 appropriately. SpO2 in the 85–95% range excludes hyperoxia , but there is no consensus for the optimal saturation target ranges in premature newborns . A collaborative prospective meta-analysis of five ongoing trials in the USA, Australia, United Kingdom, New Zealand and Canada (NeOProM), aimed to establish the optimal SpO2 target ranges for extremely preterm newborns [9–11]. Interim results recommend SpO2 levels > 90% to avoid mortality . Previous recommendations suggested SpO2 levels < 85% [12, 13]. There is no evidence or consensus to guide the administration of oxygen in full-term newborns .
The Neovent study group found that time-cycled, pressure-limited ventilation was the most common mode currently used for neonatal ventilation. The tidal volume (TV) was usually targeted to 4–7 ml/kg . The same study group found that hypocarbia was relatively uncommon during neonatal ventilation, and they speculated that hypercarbia was more common because of the practice of permissive hypercarbia .
Few studies relate arterial oxygen tension values to SpO2 targets . It is not known how often the problem of hyperoxemia and hyperventilation occurs in NICUs in Norway. The present retrospective study is the first part of a larger study that aims to discover areas for quality improvement regarding oxygen and ventilator treatment of preterm and sick newborns. The purpose of this study was to investigate the documentation of oxygenation and ventilation among newborns receiving mechanical ventilation in a Norwegian NICU, to report on the following: (1) use of oxygen during ventilator treatment; (2) extent of charted deviations from oxygenation and ventilation targets; and (3) data associated with variations in MAP.
Patients and study design
The setting of this study was a level 4 NICU at Rikshospitalet, Oslo University Hospital, Norway. We retrospectively studied the documentation of oxygen and ventilator treatment practices between July 2010 and November 2011using intensive care charts. Patients who satisfied the inclusion criteria were identified through the NICU proprietary quality control database. Infants were eligible for inclusion if they had been mechanically ventilated for a minimum of 3 hours, and we chose to limit data collection to a maximum of 48 hours for each patient. The sample was grouped by gestational age (GA) into extremely premature (23–28 weeks GA), moderately premature (29–37 weeks GA), full-term (38–41 weeks GA), and newborns > 41 weeks post-conceptual age (defined as GA plus chronological age ). The principal diagnoses were categorised using ICD 10 (KITH–Health Affairs) and were: immature lungs, other lung/respiratory problems, circulatory problems, and infection.
Relevant variables for oxygenation and ventilator treatment were defined according to the literature, clinical practice, and discussion with experts in the field. Variables collected and reported in this study were GA, sex, birthweight, diagnoses (infection, lung problems, immature lungs and circulatory problems), peak inspiratory pressure (PIP), positive end-expiratory pressure (PEEP), inspiratory time (TI), FiO2, TV, MAP, highest preductal SpO2 value, arterial pressure of oxygen (PaO2), and PaCO2. Data were collected on expiratory TV measured from the ventilator, because this measures the effectiveness of mechanical ventilation . Blood gas analysis was carried out on arterial samples only, because capillary blood gases were deemed to have insufficient reliability for our purposes. This reliability problem relates to newborn infants crying as a reaction to vasopuncture, which frequently presents as rapid changes in PaO2 and PaCO2. Capillary blood gases were therefore analysed solely for comparison with arterial gases.
The limit of acceptable PaCO2 was set at 4.7–5.9 kPa . Normal limits for PaO2 were 6–10 kPa for premature newborns  and 8–10.7 kPa for full-term infants . Appropriate limits for SpO2 for newborns receiving supplemental oxygen were set at 88–93% in premature infants and not above 95% in full-term newborns, according to existing practice guidelines in our unit. The normal limits of TV were considered to be 4–6 ml/kg for premature and 5–8 ml/kg for full-term newborns .
Statistical analyses were performed using SPSS version 19.0 (SPSS Inc., Chicago, IL, USA). A research assistant checked data for accuracy. A power analysis resulted in the inclusion of charts from 138 premature and full term newborns. This was based on the assumption that the proportion of newborns with SpO2 outside the recommended limit, was approximately 10%, estimated with an accuracy of + / – 5% and calculated with a 95% confidence interval (CI).
Descriptive statistics and frequencies were calculated for GA, weight, number of newborns with measurements of arterial blood gases, 1-hour periods with and without arterial and capillary blood gas measurements (based on 24 1-hour periods per day), PaO2, PaCO2, and TV. The recorded values for TV were summarised with respect to the two limits that were applicable to premature and full-term newborns in the NICU. Preductal SpO2 values were analysed for the trend in mean values over time. The cutoff point for defining hypocarbia was PaCO2 < 4.7 kPa, and extreme hypocarbia was defined as PaCO2 < 4.0 kPa. When analysing the distribution of arterial blood gases, and comparing the prevalence of high versus low PaO2 and PaCO2 between the GA groups, we controlled for possible interdependence due to repeated measurements in each individual by using a generalised linear model and Wald’s analysis . We used partial correlation to analyse the correlation between PaCO2 and TV . A mixed linear model with repeated measures was used to analyse variations in MAP. Any p value < 0.05 was considered significant.
Approval by the Regional Committee for Medical Research Ethics in Norway was not required for this study because data collection was anonymous. Permission for the study was obtained from the Data Protection Officer at Oslo University Hospital and from the director of the NICU.
The documentation included 4978 hours of ventilator time for 138 newborns. Ventilator support consisted of 4702 hours of conventional mechanical ventilation and 276 hours of oscillation. The minimum duration of ventilator treatment was 5 hours and our predetermined maximum period of study was 48 hours. There were 85 male (62%) and 53 female (38%) infants. GA ranged from 23 to 52 weeks and weight at the time of study entry was 426–5345 g. There were 42 (30%) extremely premature, 42 (30%) moderately premature, 34 (25%) full-term and 20 (15%) newborns aged > 41 weeks. The leading diagnoses were lung immaturity [57 (41%)], other lung problems [43 (31%)], circulatory problems (including congenital heart disease) [32 (23%)] and infection [6 (5%)].
Use of oxygen during ventilator treatment
Figure 1 summarises the distribution of oxygen concentrations received by the newborns. In total, newborns received 2020 hours of ventilator support with FiO2 = 21%. Oxygen was given at a median FiO2 of 24% in premature newborns and a median of 21% in full-term newborns. The mean values were respectively 30% (CI 95% = 29.4 - 30.6) and 32% (CI 95% = 31.1 – 32.9).
Charting of deviations from oxygenation and ventilation targets
Blood samples were taken for arterial blood gas analysis from once an hour, to every 3 to 4 h. Table 1 shows the number of newborns with measurements of arterial blood gases.
Monitoring and assessment of oxygenation
Table 2 shows the arterial blood gas measurements in 1-hour periods obtained during ventilator treatment.
Capillary blood gas measurements were obtained in 47 patients. There were 226/4978 (5%) 1-hour periods with recorded capillary blood gases: with 39 (17%) in extremely premature, 94 (42%) in moderately premature, 56 (25%) in full-term, and 37 (16%) in newborns > 41 weeks old. Analysis comparing the frequency of capillary and arterial blood gas sampling showed that arterial blood gas sampling was most common. Wald’s analysis indicated significant differences between the GA groups in the frequency of using arterial blood gases (p = 0.02). Further analysis showed that samples for arterial blood gas analysis were withdrawn significantly more often from the extremely premature versus moderately premature infants (p = 0.004).
Table 3 shows the distribution of low, normal and high PaO2 values according to GA in newborns who received oxygen supplementation.
Although extremely premature infants had relatively few high PaO2 values, the percentage of high value samples increased with GA. There was a significant difference between GA groups with regard to the risk of low versus normal PaO2 values (p < 0.001) and normal versus high values (p < 0.001). Further calculations showed that the odds ratio (OR) of normal PaO2 (compared to high PaO2) in newborns of GA 23–28 weeks was 16 (95% CI = 6.5–40) times greater than that in newborns of GA > 41 weeks.
The mean SpO2 trends over time are shown in Figure 2.
Monitoring and assessment of ventilation
Low PaCO2 values were recorded in 187/1170 (16%) samples (70 patients). When the cutoff for low PaCO2 was set at 4 – 4.7 kPa, there were 42 (34%) low PaCO2 values in extremely premature, 30 (24%) in moderately premature, and 35 (28%) in full-term newborns, and 16 (13%) in newborns aged > 41 weeks. The analysis showed no significant difference in the occurrence of low PaCO2 between GA groups (p = 0.639). However, when the cutoff value was reduced to < 4 kPa (n = 64), there was a significant difference (p = 0.015) between the groups. The OR of PaCO2 < 4 kPa in newborns of GA 23–28 weeks was 5 (95% CI = 1.6–15) times greater than that in newborns of GA > 41 weeks.
For premature infants of GA 23–37 weeks, there were 2374 recorded TV measurements in 1-hour periods, with a median TV of 4.7 ml/kg. Less than half of the measurements (n = 1058) were within the normal reference range of 4–6 ml/kg, and 520 (22%) were > 6 ml/kg. There was no significant correlation between PaCO2 and TV (r = 0.007, p = 0.87). Full-term newborns had a total of 1535 TV measurements in 1-hour periods, with a median TV of 6.1 ml/kg. More than half of all measurements (n = 883) were in the normal range, while 301 (20%) were > 8 ml/kg. There was a significant weak negative correlation between PaCO2 and TV (r = −0.12, p = 0.03).
Data associated with variations in MAP
A linear mixed effects model was used to estimate variations in MAP using the independent variables PEEP, TI and PIP (Table 4).
Table 4 shows that there was a significant correlation between MAP variability, TI and PIP, which decreased with the passage of time. There was no significant correlation between MAP variability and PEEP. This suggests that changing PEEP was not used as a strategy to adjust MAP and thereby affect oxygenation. The two groups of premature infants received PEEP with a mean pressure of 4.7 cm H2O (CI 95% = 4.67 – 4.73), and full-term newborns received PEEP with a mean pressure of 4.9 cm H2O (CI 95% = 4.85 – 4.95).
The purpose of oxygen administration in the NICU is to prevent free radical damage . In the present study, the median FiO2 in premature newborns was lower (24%) than in the study by van Kaam et al., who found a median FiO2 value of 28% during conventional mechanical ventilation . Furthermore, we found that newborns achieved acceptable mean SpO2 during the observation period (Figure 2). Targeting saturation to 88–93% in premature newborns and not more than 95% in full-term newborns has been a goal for several years [25, 26]. However, these targets have not yet been fully agreed upon, because the results of ongoing randomised trials testing high versus low SpO2 targets are still pending [9, 10, 27]. Our calculation of mean SpO2 may have masked episodes of hypoxemia or hyperoxia when nurses titrated FiO2. Nevertheless, it suggests that the goal for saturation targets was met to a large extent. Using SpO2 alone to guide decisions concerning oxygen administration is not evidence-based practice. Thus, nurses also have to assess skin colour, heart rate, and values from blood gases as well transcutaneous O2 measurements .
Although it has been suggested that entrenched clinical practices and cultures make it difficult to change the use of oxygen , our results showed that extremely premature newborns did not receive excessive amounts of oxygen (Tables 1 and 3). The incidence of hyperoxemia increased with GA. Because of the haemoglobin oxygen dissociation curve , sick newborns with pulmonary hypertension are at higher risk of developing high levels of hyperoxemia when they are treated with an oxygen saturation of ~95%. It is common practice in NICUs that nurses wait until an infant is stable before withdrawing arterial blood for gas analysis, and this practice will affect the results. In addition nurses’ workload is an important factor in the achievement of SpO2 goals and appropriate oxygen management in the NICU. There is evidence that compliance with saturation targets is improved with higher nurse: patient ratios . The good results for oxygen management in the NICU in the present study may in part have been due to the practice of having a 1:1 nurse:patient ratio for all infants on mechanical ventilation.
Our study revealed that low PaCO2 values occurred most commonly in extremely premature infants. Hypocarbia may cause cerebral vasoconstriction, resulting in decreased oxygen delivery to the brain [7, 10]. Moreover, in extremely preterm infants, PaCO2 in the normal range seems to yield the best electroencephalography activity . PaCO2 can be regulated by controlling the minute ventilation with TV, or the ventilator rate . It is therefore noteworthy that we found high TVs in 22% of preterm infants and 23% of newborns aged >37 weeks. Nevertheless, the median TV of 4.7 ml was lower than the 5.3 ml recorded in the study of van Kaam et al. . Mechanical ventilation using high TVs is known to cause lung damage . We did not find any significant correlation between PaCO2 and TV in the premature infants, and only a weak negative correlation in full-term infants. This suggests that the relationship between lung physiology and what happens during respirator treatment may not be simple.
Our results showed that variations in PEEP had no significant effect on MAP variability (Table 4). Changing PEEP is often an effective way to adjust oxygenation  and to regulate MAP for experienced clinicians who can accurately assess changes in measured data and calculate the impact of any adjustments . The premature newborns in our study had a PEEP mean value of 4.7 cm H2O which was similar to that observed in the study of van Kaam et al., although they suggested that PEEP values >7 cm H2O might be protective for the lungs .
It is suggested that maintaining appropriate oxygenation is hindered by insufficient communication of unit policies as well as personal bias about the best practice . Control of oxygenation and ventilation is crucial during mechanical ventilation. Therefore, further studies should identify how communication and allocation of responsibilities between nurses and physicians can reduce the incidence of hypocarbia and hyperoxemia.
One limitation of this study was that we had no record of fluctuations in oxygenation levels, nor observations of how soon after blood gas analyses adjustments were made. Regarding the analysis of the mean SaO2, it would have been helpful to include an analysis that indicated the uncertainty in the results. However, because of the volume of repeated measurements for each individual, a standard box plot could not be used. Another weakness of the study was related to the high TV values, which were presented as a general occurrence regardless of how many high values there were from each individual patient.
We showed that, in general, premature newborns were treated within the desired limits of SpO2, and few high PaO2 values were noted. The occurrence of high PaO2 values increased significantly with GA. Many recorded TVs were too high and hypocarbia during ventilation was more common in the extremely premature infants.
MTS, RN, Intensive Care Nurse Specialist, Master of Nursing Science, PhD student, Faculty of Medicine, University of Oslo, Norway.
ITB, RN, Master of Nursing Science, PhD. Professor Institute of Nursing and Health Sciences, Faculty of Medicine, University of Oslo, Norway. ITB is a supervisor of MTS.
TWRH, MD, PhD, Professor, Department of Neonatology, Women’s & Children’s Division, Oslo University Hospital – Rikshospitalet, Norway and Institute for Clinical Medicine, Faculty of Medicine, University of Oslo, Norway. TWRH has been the senior supervisor of MTS.
Fraction of inspired oxygen
Mean airway pressure
Neonatal intensive care unit
Positive end-expiratory pressure
Peak inspiratory pressure
Sola A, Saldeno YP, Favareto V: Clinical practices in neonatal oxygenation: where have we failed? What can we do?. J Perinatol. 2008, 28 (Suppl 1): S28-S34.
Solberg MT, Hansen TW, Bjork IT: Nursing assessment during oxygen administration in ventilated preterm infants. Acta Paediatr. 2011, 100 (2): 193-197. 10.1111/j.1651-2227.2010.02094.x.
Muhlethaler V, Malcolm G: Mechanical ventilation in the newborn; a simplified approach.Part 1: Intermittent positive pressure ventilation. J Paediatr Child Health. 2012, 48 (8): 649-652. 10.1111/j.1440-1754.2010.01872.x.
Deuber C, Terhaar M: Hyperoxia in very preterm infants: a systematic review of the literature. J Perinat Neonat Nur. 2011, 25 (3): 268-274. 10.1097/JPN.0b013e318226ee2c.
Groenendaal F, De Vooght KMK, van Bel F: Blood gas values during hypothermia in asphyxiated term neonates. Pediatrics. 2009, 123 (1): 170-172. 10.1542/peds.2008-1955.
Curley G, Laffey JG, Kavanagh BP: Bench-to-bedside review: carbon dioxide. Crit Care. 2010, 14 (2): 220-10.1186/cc8926.
Thome UH, Ambalavanan N: Permissive hypercapnia to decrease lung injury in ventilated preterm neonates. Semin Fetal Neonat M. 2009, 14 (1): 21-27. 10.1016/j.siny.2008.08.005.
Quine D, Stenson BJ: Arterial oxygen tension (Pao2) values in infants< 29 weeks of gestation at currently targeted saturations. Arch Dis Child-Fetal. 2009, 94 (1): F51-F53.
Askie LM: CURRENT OPINION optimal oxygen saturations in preterm infants: a moving target. Curr Opin Pediatr. 2013, 25: 000-000.
Askie LM, Brocklehurst P, Darlow BA, Finer N, Schmidt B, Tarnow-Mordi W, Ne OCG: NeOProM: neonatal oxygenation prospective meta-analysis collaboration study protocol. BMC Pediatr. 2011, 11: 6-
Carlo W, Finer N, Walsh M, Rich W, Gantz M, Laptook A, Yoder B, Faix R, Das A, Poole W: Target ranges of oxygen saturation in extremely preterm infants. New Engl J Med. 2010, 362 (21): 1959-
Sola A, Rogido MR, Deulofeut R: Oxygen as a neonatal health hazard: call for detente in clinical practice. Acta Paediatr. 2007, 96 (6): 801-812. 10.1111/j.1651-2227.2007.00287.x.
Sweet DG, Carnielli V, Greisen G, Hallman M, Ozek E, Plavka R, Saugstad OD, Simeoni U, Speer CP, Halliday HL: European consensus guidelines on the management of neonatal respiratory distress syndrome in preterm infants–2010 update. Neonatology. 2010, 97 (4): 402-417. 10.1159/000297773.
Stenson BJ, Orme JA: The twists and turns of neonatal oxygen therapy. Early Hum Dev. 2012, 88 (12): 961-963. 10.1016/j.earlhumdev.2012.09.003.
van Kaam AH, De Jaegere AP, Rimensberger PC, Debeer A, Chemin A, Norbert K, Autret F, Andreou A, Kroon A, Minić A: Incidence of hypo-and hyper-capnia in a cross-sectional European cohort of ventilated newborn infants. Arch Dis Child-Fetal. 2013, 98: F323-F326. 10.1136/archdischild-2012-302649.
Goldstein MR: Ventilator management: what does it all mean?. Newborn Infant Nurs Rev. 2006, 6 (2): 79-86. 10.1053/j.nainr.2006.03.003.
Kramme R, Hoffmann K-P, Pozos RS: Springer Handbook of Medical Technology. 2012, Berlin, Heidelberg: Springer Berlin Heidelberg
Gardner SL, Merenstein GB: Handbook of neonatal intensive care. 2006, St. Louis: Mosby
Deulofeut R, Critz A, Adams-Chapman I, Sola A: Avoiding hyperoxia in infants < or = 1250 g is associated with improved short- and long-term outcomes. J Perinatol. 2006, 26 (11): 700-705. 10.1038/sj.jp.7211608.
Bhutani VK, Sivieri EM: Clinical use of pulmonary mechanics and waveform graphics. Clin Perinatol. 2001, 28 (3): 487-503. 10.1016/S0095-5108(05)70102-0. v
Fitzmaurice GM, Laird NM, Ware JH: Applied longitudinal analysis. 2011, Chichester: Wiley
Cohen J: Statistical power analysis for the behavioral sciences. 1988, Hillsdale, N. J.: Laurence Erlbaum
Raghuveer TS, Bloom BT: A paradigm shift in the prevention of retinopathy of prematurity. Neonatology. 2011, 100 (2): 116-129. 10.1159/000322848.
van Kaam AH, Rimensberger PC, Borensztajn D, De Jaegere AP: Ventilation practices in the neonatal intensive care unit: a cross-sectional study. J Pediatr. 2010, 157 (5): 767-771. 10.1016/j.jpeds.2010.05.043. e763
VanderVeen DK, Mansfield TA, Eichenwald EC: Lower oxygen saturation alarm limits decrease the severity of retinopathy of prematurity. J American Assoc Pediat Ophthal Strabismus. 2006, 10 (5): 445-448. 10.1016/j.jaapos.2006.04.010.
Chow LC, Wright KW, Sola A: Can changes in clinical practice decrease the incidence of severe retinopathy of prematurity in very low birth weight infants?. Pediatrics. 2003, 111 (2): 339-345. 10.1542/peds.111.2.339.
Saugstad OD, Aune D: In search of the optimal oxygen saturation for extremely low birth weight infants: a systematic review and meta-analysis. Neonatology. 2011, 100 (1): 1-8. 10.1159/000322001.
Dawson JA, Vento M, Finer NN, Rich W, Saugstad OD, Morley CJ, Davis PG: Managing oxygen therapy during delivery room stabilization of preterm infants. J Pediatr. 2012, 160 (1): 158-161. 10.1016/j.jpeds.2011.07.045.
Goldsmith J gJ: Neonatal intensive care unit oxygen management: a team effort. Pediatrics. 2007, 119 (119): 1195-
Sink DW, Hope SAE, Hagadorn JI: Nurse: patient ratio and achievement of oxygen saturation goals in premature infants. Arch Dis Child-Fetal. 2011, 96 (2): F93-F98. 10.1136/adc.2009.178616.
Wikstrom S, Lundin F, Ley D, Pupp IH, Fellman V, Rosen I, Hellstrom-Westas L: Carbon dioxide and glucose affect electrocortical background in extremely preterm infants. Pediatrics. 2011, 127 (4): e1028-e1034. 10.1542/peds.2010-2755.
Habre W: Neonatal ventilation. Best Pract Res Clin Anaesthesiol. 2010, 24 (3): 353-364.
Morley C: Treatment of respiratory failure: Mechanical ventilation. Neonatology A practiacal approach to neonatal diseases. edn. Edited by: Buonocore G, Bracci R, Weindling M. 2012, Milano: Springer, 497-508.
Bhutani VK, Benitz WE: Pulmonary Function and Graphics. Assisted ventilation of the neonate. Edited by: Goldsmith JP, Karotkin EH. 2011, St. Louis, Mo: Elsevier/Saunders, 1 online resource (xvii, 631 s.): ill, 5
Bancalari E, Claure N: Control of oxygenation during mechanical ventilation in the premature infant. Clin Perinatol. 2012, 39 (3): 563-10.1016/j.clp.2012.06.013.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2431/13/126/prepub
We are grateful to Professor Magne Thoresen for his statistical advice. This study was supported by the Norwegian Premature Association and the Norwegian Extra Foundation for Health and Rehabilitation.
The authors declare that they have no competing interests.
MTS formulated the hypotheses and designed and conducted the research, analyzed the data, and wrote the paper. ITB participated in the formulation of the hypotheses and study design, co-authored the paper. TWRH participated in the formulation of the hypotheses and study design, offered advice on data collection, and co-authored the paper. All authors read and approved the final manuscript.