Translator Disclaimer
1 May 2010 Precision of coherence analysis to detect cerebral autoregulation by near-infrared spectroscopy in preterm infants
Author Affiliations +
Coherence between spontaneous fluctuations in arterial blood pressure (ABP) and the cerebral near-infrared spectroscopy signal can detect cerebral autoregulation. Because reliable measurement depends on signals with high signal-to-noise ratio, we hypothesized that coherence is more precisely determined when fluctuations in ABP are large rather than small. Therefore, we investigated whether adjusting for variability in ABP (variabilityABP) improves precision. We examined the impact of variabilityABP within the power spectrum in each measurement and between repeated measurements in preterm infants. We also examined total monitoring time required to discriminate among infants with a simulation study. We studied 22 preterm infants (GA<30) yielding 215 10-min measurements. Surprisingly, adjusting for variabilityABP within the power spectrum did not improve the precision. However, adjusting for the variabilityABP among repeated measurements (i.e., weighting measurements with high variabilityABP in favor of those with low) improved the precision. The evidence of drift in individual infants was weak. Minimum monitoring time needed to discriminate among infants was 1.3-3.7 h. Coherence analysis in low frequencies (0.04-0.1 Hz) had higher precision and statistically more power than in very low frequencies (0.003-0.04 Hz). In conclusion, a reliable detection of cerebral autoregulation takes hours and the precision is improved by adjusting for variabilityABP between repeated measurements.



Near-infrared spectroscopy (NIRS) has been used in clinical research to study cerebral intravascular oxygenation in newborns for more than 20years .1 Recently, NIRS has also been introduced to detect cerebral autoregulation. Cerebral autoregulation is a physiological mechanism that protects the brain against a fluctuating arterial blood pressure (ABP) by maintaining an approximately constant level of cerebral blood flow.2 This is achieved through dilatation or constriction of cerebral blood vessels. The assessment of cerebral autoregulation is based on analysis of coherence between changes in cerebral intravascular oxygenation and changes in ABP. Using spontaneous fluctuations in ABP, this method offers a potential clinical tool to monitor cerebral autoregulation noninvasively and continuously in premature infants. Potentially, this can guide therapeutic interventions and therefore prevent brain damage.3

Coherence is a frequency domain estimate of linear correlation. It estimates the degree to which the variance in cerebral intravascular oxygenation at a certain frequency can be explained by the variance of the ABP at the same frequency. Coherence ranges from 0 to 1, with 1 indicating perfect correlation and 0 complete lack of correlation. High coherence between ABP and cerebral intravascular oxygenation indicates imperfect cerebral autoregulation, while low coherence indicates effective autoregulation (Fig. 1 ).4 Coherence provides a qualitative but not a quantitative measure of cerebral autoregulation. Quantitative assessment of autoregulation (i.e., the magnitude of pressure passivity) requires estimation of the operational gain or amplification factor between the two signals. However, quantification of operational gain is only reasonable if coherence is significant. Thus, an imprecise estimate of coherence makes analysis of gain inappropriate. Consequently, our study concentrates exclusively on optimizing the precision of coherence analysis to detect cerebral autoregulation.

Fig. 1

Example of simultaneous changes in the detrended signals of MAP and OI and their coherence-spectrum illustrating (a) high and (b) low coherence. The dotted line represents the shift from VLF to LF range.


A reliable detection of cerebral autoregulation requires a certain challenge to the autoregulatory system (i.e., a certain amount of ABP changes).4, 5, 6 In vulnerable patients such as premature infants, however, use of spontaneous ABP changes is preferable to induced changes. ABP with small changes has low values in its power spectral density (PSD). Detecting cerebral autoregulation on the basis of such a signal could be misleading and falsely indicate a low coherence and, therefore, effective cerebral autoregulation, whereas the variability of cerebral oxygenation is in fact dominated by “noise,” e.g., (i) physiological noise (such as variability of arterial content of carbon dioxide, arterial saturation or cerebral oxygen uptake) or (ii) instrumental noise (due to limitations of the optics, electronics or movements).

Previous reports using coherence to detect cerebral autoregulation have not adjusted for variability in ABP.7, 8, 9 This analytical approach has shown a lack of stability, with cerebral autoregulation that fluctuates (i.e., that comes and goes over short periods of time).8 Thus, we set out to investigate whether adjusting for the magnitude of ABP changes would lead to a more precise detection of cerebral autoregulation.

The hypothesis of our study was that coherence would be more precisely determined when ABP changes were large rather than small. In general, a reliable test of a functional relation requires that the stimulus itself exceeds the noise level of the measurements (i.e., a high signal-to-noise ratio). Hence, in the present study we examined the effect of ABP changes at specific frequencies [i.e., the effect of the shape of the power spectrum in each measurement (intra-measurement adjustment)]. We also examined the effect of ABP sometimes varying greatly and sometimes varying a little, by performing a weighted analysis where, out of several measurements in each infant, we weighted those with large ABP changes higher than those with small (intermeasurement adjustment).


Materials and Methods



We studied 22 premature infants with a gestational age of <32weeks with an indwelling arterial catheter, which had been inserted for clinical reasons at birth. The study was carried out at the Neonatal Care Unit, Copenhagen University Hospital, Rigshospitalet, Copenhagen. The study was conducted in parallel with an ongoing study on the link between fetal inflammation and cerebral autoregulation, and hence, only children born at Rigshospitalet were eligible. Infants with major malformations or diagnosis of brain injury at recruitment were excluded. We aimed to examine infants within the first 24h after birth. Median age at examination was 17.4 (4.043.5)h . Clinical data were obtained from the infants’ medical records. The latest blood gas value taken before the examination was also recorded [median time interval between blood sample and examination: 2.0 (06.4)h ]: Demographic and clinical data of the infants are summarized in Table 1 . The Danish Local Ethical Committee approved the study (Journal No. H-A-2007-0109), and written informed parental consent was obtained in all cases.

Table 1

Demographic data.

No. of infants22
Sex [no]10 males/12 females
Gestational age [weeks]27.5 [24.1–29.4]
Birth weight [g]940 [460–1266]
Apgar at 5min 10 [6–10]
Umbilical cord pH7.32 [7.03–7.45]
Circulatory support in the first day of life
Inotropics [no (%)]2 (9%)
Volume treatment [no (%)]6 (27%)
Intropics or volume treatment [no (%)]6 (27%)
Clinical state at measurement
Postnatal age [h]17.4 [4.0–43.5]
Inotropics [no (%)]2 (9%)
Nasal-CPAP [no (%)]20 (91%)
Mechanical ventilation [no (%)]2 (9%)
Inspiratory oxygen concentration [%]25 [21–40]
Arterial oxygen saturation [%]95 [85–99]
PaCO2 [kPa]5.8 [4.4–8.2]
PaO2 [kPa]8.0 [4.0–14.0]
Haemoglobin [mM]9.2 [5.5–11.9]
MAP [mm Hg]35 [26–42]
Later outcome
Intraventricular haemorrhage (all) [no (%)]5 (23%)
Intraventricular haemorrhage (Gr. III/IV) [no (%)]2 (9%)
Cystic perventricular leucomalacia [no (%)]0
Neonatal mortality [no (%)]2 (9%)
Median and [range] unless otherwise stated.


Continuous NIRS and Arterial Blood Pressure Recording

NIRS data were recorded using a NIRO-300 spectrophotometer (Hamamatsu Photonics, Hamamatsu City, Japan). The probes were fixed in a nontransparent, soft probe holder and secured to the frontotemporal or frontoparietal region of the head with a flexible bandage. The interoptode distance was 4cm . Because the skin of premature infants is very thin and vulnerable, no adhesive tape was used. Changes in the relative concentrations of oxygenated and deoxygenated haemoglobin were recorded and used to calculate the oxygenation index (OI). OI is the difference between oxygenated and deoxygenated haemoglobin divided by a factor of 2. Changes in OI have been validated to represent changes in cerebral blood flow in several animal studies.10, 11, 12

ABP waveforms were measured with a pressure transducer connected to the indwelling arterial catheter. Before the mean ABP (MAP) was calculated, the ABP data were filtered by a first-order low-pass filter with a cutoff frequency of 0.3Hz .

NIRS and MAP data were sampled simultaneously and synchronized with arterial saturation (assessed by pulse oximetry) at 2Hz by an analog-to-digital converter, and stored on a laptop for further analysis. Using in-house written software (Labview) recording was automatically stopped if changes in arterial saturation exceeded 5%. The infants were clinically stable during measurements. One investigator (G.H.H.) observed the infants and the NIRS instrument continuously during the measurements and interrupted the recording if movement artefacts caused a baseline shift (e.g., radiographs and manipulation of the infant due to suctioning, etc.) or if there was loss of the ABP signal due to blood-gas sampling. To eliminate episodes of artifacts, the data were divided into 10-min measurements of uninterrupted recordings.


Frequency Analysis

Transfer function analysis of coherence between continuous measurements of OI and MAP were performed in two frequency ranges. The very low-frequency (VLF) range ( 0.0030.04Hz , corresponding to a cycle length of 25300s ) and the low-frequency (LF) range ( 0.040.1Hz , corresponding to a cycle length of 1025s ). Each 10-min measurement was subdivided into three segments of 5min with 50% overlap and a Hanning window was applied to minimize spectral leakage. Both time signals were then detrended and transformed into PSDs by Fourier transform (Matlab, MathWorks, Natick, Massachusetts).


Adjustment for Variability in ABP

We adjusted for variability in ABP in two steps: (i) within each measurement (intrameasurement adjustment) and (ii) between repeated measurements (intermeasurement adjustment) (Fig. 2 ). For intrameasurement adjustment, we calculated coherence in three different ways, as follows:

  • 1. Mean value of the coherence spectrum (Cohmean) , i.e., no adjustment for variability in ABP

    Eq. 1

    where N is the number of PSD values within the predefined frequency range and fi is the i ’th frequency within the frequency range.

  • 2. Intrameasurement ABP weighted mean value according to the distribution of PSDMAP (CohMAPwmean) , thereby attaching most weight to the frequencies with highest variability in ABP,

    Eq. 2

    where wi=PSDMAP(fi)i=1NPSDMAP(fi) and i=1Nwi=1 .

  • 3. Point measurement at the frequency with maximum PSDMAP (CohMAPmax) ,

    Eq. 3

    where fmaxMAP corresponds to the frequency at which PSDMAP has the maximum value within the predefined frequency range.

Fig. 2

Example of ten 10-min measurements in one of the infants (infant number 20) in the VLF range. (b) Illustrates the effect of adjusting for ABP variability among repeated measurements (intermeasurement adjustment). In this example, the measurements with high ABP variability (i.e., mean power of MAP, illustrated with bars) have the highest CoherenceST (cross), and hence, the adjusted overall estimate of CoherenceST is higher than the unadjusted. In this case, the estimate even changes from being statistically insignificant to statistically significant (CoherenceST0) . (a) Illustrates the effect of adjusting for ABP variability within measurement number 4 (intrameasurement adjustment). Mean coherence (Cohmean) [solid line], intrameasurement ABP weighted mean coherence (CohMAPwmean) [dashed-doted line], and coherence at the frequency with maximum PSDMAP (CohMAPmax) [dashed line] are shown.


For intermeasurement adjustment, we adjusted each infant’s overall estimate of coherence for the varying degree of variability in ABP between repeated measurements. This was performed by a weighted analysis. We calculated a weighting factor for each measurement as the proportion of mean PSDMAP in the measurement concerned compared to the overall sum of mean PSDMAP in each infant. This factor, reaching the sum of 1 in each infant, was used to weight each repeated measurement. The weighted model is referred to as “intermeasurement ABP weighting.”


Threshold of Significant Coherence

We estimated the 95% confidence limit for each coherence measurement mentioned above using Monte Carlo simulations. This confidence limit was then used as a threshold (ThresholdCOH) to determine whether a measurement is coherent or not (i.e., coherent if the coherence measurement is above ThresholdCOH ). In the Monte Carlo simulations, OI and MAP signals were both reshuffled randomly to obtain uncorrelated signals with signal variances that were the same as the original signal. The same procedure was repeated 10,000 times, and 10,000 values resulted for each coherence measurement. These 10,000 values were then sorted in ascending order, and the 95% percentile was taken as the 95% confidence limit of ThresholdCOH . We standardized each measurement of coherence by subtracting ThresholdCOH from it yielding standardized coherence (CoherenceST) . CoherenceST0 indicates imperfect cerebral autoregulation.


Measure of Precision

We applied different analytical approaches to the same data and compared their precision in terms of repeatability between repeated measurements. We hypothesized that cerebral autoregulation stayed stable between two successive 10-min measurements, and hence, that a high repeatability reflects a high signal-to-noise ratio (i.e., a low amount of measurement error). Statistically, we looked for the analytical approach with the best discrimination among infants based on their level of coherence. This means that we looked for the analytical approach with the highest variability in level of coherence among infants (i.e., a large range of estimates among infants) that at the same time had a small range of variation between repeated measurements in the same infant [i.e., a small standard error (SE)]. These conditions are combined in the P value. Hence, we chose to favor the analytical approach with the most significant difference in level of coherence among our infants.


Statistical Methods

Averages of coherence, CoherenceST , ThresholdCOH, and variability in MAP were calculated for each infant, and from these averages, we calculated a grand mean. Variability in MAP in each measurement was calculated as the square root of mean PSDMAP in each frequency range. To estimate the variability in mean PSDMAP between measurements, we calculated each infant’s coefficient of variation.

One-way analysis of variance for repeated measurements with spatial power as covariance structure was used to describe each infants overall estimate of CoherenceST (“proc mixed” in SAS 9.1, SAS Institute). This covariance structure takes the nonequidistant nature of the measurements due to artifacts, blood-gas sampling, or changes in arterial saturation exceeding 5% into account. We used two models. Model 1 used infant identification as the only covariate, allowing each infant to have its own level of CoherenceST with no effect of time. Model 2 added time and interaction between time and infant identification, thereby allowing each infant to have an individual level and trend of CoherenceST . We used model 2 to identify a possible drifting nature of cerebral autoregulation. We used model 1 to calculate each infant’s overall estimate of coherence and to divide the infants into the two following groups: effective autoregulation (i.e., CoherenceST<0 ) and imperfect autoregulation (i.e., CoherenceST0 ). An independent-samples T-test was conducted to compare gestational age (GA) and postnatal age for infants in these two groups. A nonparametric Fisher’s exact test was used to test if clinical state (use of inotropics or mechanical ventilation), incidence of intraventricular haemorrhage (IVH) grade 3–4, or death in the first month of life differed between these groups (SPSS statistics 17.0). P<0.05 was considered statistically significant.

In order to estimate the minimum monitoring time needed to reach a robust estimate of cerebral autoregulation, we examined statistical power of Cohmean , CohMAPwmean , and CohMAPmax by means of a simulation study. Using estimated values of CoherenceST , the within infant variation, the autocorrelation, and the residual variation from model 1 with an autoregressive structure, we generated data sets with 22 infants for each of the three measurements of coherence. The number of measurements varied between 5 and 35. Statistical power according to number of measurements was then estimated by the empirical rejection rate in the 1000 simulated data sets.




General Features of Blood Pressure, Coherence, and Threshold of Coherence

We monitored 22 infants continuously with a median duration of 2.1 (1.23.7)h yielding 215 measurements of 10min [median: 10, (4–17)]. Grand mean variability in MAP was 4.8mm Hg (SD: 1.5) in VLF and 1.8mm Hg (SD: 0.7) in LF. Mean coefficient of variation in mean PSDMAP between measurements was 62% (SD: 21%) and 50% (SD: 21%) for VLF and LF, respectively ( P : 0.06).

CohMAPmax had the highest ThresholdCOH in both frequency bands [0.80 (SD: 0.002) in both VLF and LF] followed by CohMAPwmean [VLF: 0.50 (SD: 0.001) and LF: 0.47 (SD: 0.007)] and Cohmean [VLF: 0.47 (SD: 0.006) and LF: 0.45 (SD: 0.005)].

The variability between infants in single measurements of Cohmean was 0.06 (Table 2 ), while the variability in Cohmean within infants was 0.1. This means that the imprecision of a measurement of coherence from a single 10-min measurement exceeded the variability in grand mean coherence between infants.

Table 2

Grand mean of each infant’s average coherence, threshold of coherence, and standardized coherence (CoherenceST) .

BandCoherence measurement N CoherenceThreshold of coherenceStandardized coherence
Cohmean 220.450.0620.47 0.03 0.062
VLF CohMAPwmean 220.490.0650.50 0.01 0.065
CohMAPmax 220.510.0900.80 0.29 0.090
cohmean 220.440.0560.45 0.01 0.056
LF CohMAPwmean 220.460.0670.47 0.00 0.067
CohMAPmax 220.540.1200.80 0.26 0.121


Intrameasurement Adjustment for Blood Pressure Variability

Overall, Cohmean had the lowest SE of estimated CoherenceST . CohMAPwmean had a slightly higher SE of estimated CoherenceST , while CohMAPmax had the highest. The range of estimated CoherenceST showed the same, however, less pronounced, tendency (Table 3 ). Figure 3 illustrates this tendency on infant level. These observations are combined in the P values of a significant difference in CoherenceST among infants (Table 3), where Cohmean turned out with the most significant difference. CohMAPwmean was only able to distinguish among infants significantly in the LF-band, whereas no significant differences among infants were seen for CohMAPmax .

Fig. 3

Estimated CoherenceST and number of measurements (in square brackets) in each infant illustrating the effect of intrameasurement adjustment for ABP variability in VLF and LF. Mean coherence (Cohmean) [△], intrameasurement ABP weighted mean coherence (CohMAPwmean) [◻], and coherence at the frequency with maximum PSDMAP (CohMAPmax) [○] are shown. Error bars represent 1.96× SE of the estimate. Infants treated with Dopamine during the measurement are marked with (D). Infants who developed IVH grade 3–4 are marked with (H). Infants who died in the first month of life are marked with () .


Table 3

Effect of adjusting for ABP variability within each measurement.

BandCoherence measurementEstimated standardized coherenceSE of estimatep-value
Cohmean 0.214 (0.1100.104) 0.040.00131
VLF CohMAPwmean 0.220 (0.1320.088) 0.050.0865
CohMAPmax 0.376 (0.4400.064) 0.090.5050
Cohmean 0.221 (0.0810.140) 0.040.00121
LF CohMAPwmean 0.264 (0.0900.175) 0.050.01221
CohMAPmax 0.495 (0.4780.017) 0.120.5512

a P<0.05 .


Intermeasurement Adjustment for Blood Pressure Variability

We found a positive but insignificant correlation between each infant’s mean variability in MAP and estimates of CoherenceST based on Cohmean (VLF: r2=0.07 , P=0.2 ; LF: r2=0.001 , P=0.9 ). However, residual plots of estimated CoherenceST (model 1) against mean PSDMAP showed a pattern of decreasing residuals with increasing mean PSDMAP (data not shown). We therefore weighted the estimated Cohmean with mean PSDMAP in each measurement (intermeasurement ABP weighting). This weighting resulted in a wider range of estimated CoherenceST among infants (VLF: 0.28, LF: 0.26) accompanied by almost no change in the SE (VLF: 0.03, LF: 0.04). Consequently, the P values were diminished (i.e., discrimination among infants increased). In VLF the P value changed from 0.0013 to 0.0007 and in LF from 0.0012 to 0.0001).


Level of CoherenceST

Because Cohmean and intermeasurement adjustment for blood pressure variability performed best in terms of precision, we used these conditions in model 1 to estimate each infant’s overall level of CoherenceST in both frequency bands. Indication of imperfect cerebral autoregulation [i.e., CoherenceST0 , was found in 41% (n=9) and 50% (n=11) of the infants for VLF and LF, respectively]. There was no significant difference in GA, postnatal age, or PaCO2 between infants with an imperfect (CoherenceST0) and infants with an effective autoregulation (CoherenceST<0) . Two infants received treatment with Dopamine, two infants had mechanical ventilation, two infants developed IVH grade 3–4, and two infants died in the first month of life (Table 1). None of these parameters differed significantly between the two groups.


Statistical Power

We estimated statistical power by the percentage among 1000 simulated data sets, where the null hypothesis (i.e., no significant difference in CoherenceST among infants) was rejected and plotted power against number of measurements for the three coherence measurements in VLF and LF (Fig. 4 ). The increase of statistical power was highest for Cohmean with the line approaching a power of 80% by 22 (3.7h) and 8 (1.3h) measurements, for VLF and LF, respectively. The line for CohMAPwmean had a much smaller slope and only reached a power of approximately 35% and 60% when simulating 35 measurements (5.8h) , for VLF and LF, respectively. For CohMAPmax the simulation study showed no change in power with an increasing number of measurements.

Fig. 4

Analysis of statistical power of the three coherence measurements in the (a) very low-frequency and (b) low-frequency range. Power is calculated as percentage of simulations out of 1000, where 22 infants’ standardized coherence differed significantly. Mean coherence (Cohmean) (continuous line), intrameasurement ABP weighted mean coherence (CohMAPwmean) (dotted line), and coherence at the frequency with maximum PSDMAP (CohMAPmax) (dashed line) are shown.



Time Trend

A time trend ( P value0.05 in model 2) was only demonstrated for Cohmean in the VLF-range (data not shown). No significant drift was found in the LF range, where statistical power was almost three times bigger. Consequently, our results do not provide compelling evidence of a drift in coherence over time.


Frequency Dependency

We found no evident concordance between estimated CoherenceST in VLF and LF (estimates are adjusted for inter-measurement variability in ABP) (Fig. 5 ). Only 12 infants had concordant results when classified according to CoherenceST<0 or CoherenceST0 . Even though the scatter plot is dominated by noise (i.e., wide confidence intervals of estimated CoherenceST in individual infants), it seems as if noise alone does not explain the lack of concordance between the two frequency ranges. However, the difference between estimated CoherenceST in VLF and LF was not significant ( P=0.3 in paired sampled T-test).

Fig. 5

Scatter plot of estimated CoherenceST (inter-measurement adjusted for variability in ABP) in VLF and LF illustrating the amount of concordance. Error bars represent 1.96× SE of each infant’s estimated CoherenceST [●]. The shaded areas indicate concordance between estimated CoherenceST in VLF and LF.




The principal finding of this study is that the precision of coherence analysis to detect cerebral autoregulation by NIRS is improved when the magnitude of variability in ABP is taken into account. Second, it apparently requires several hours of recordings to obtain a reliable estimate of cerebral autoregulation to characterize an individual infant.


Repeatability as a Measure of Precision

We set out to find a method that could provide a more precise detection of cerebral autoregulation based on coherence analysis of spontaneous variations in ABP and cerebral NIRS and used repeatability as a proxy for precision. We based this assumption on the fact that we only included clinically stable infants. Furthermore, as cerebral autoregulation reflects the overall activity of millions of vascular smooth muscle cells, it may be assumed that changes in autoregulatory capacity occur gradually. In an experimental study, autoregulation was abolished by 20min of severe hypoxia and recovered over 57h .13 Hence, because the infants we studied were clinically stable, we assumed that cerebral autoregulation stayed reasonably stable between two successive 10-min measurements. This was in accordance with the lack of clear evidence of drift in coherence over time in individual infants.


Adjustment for Variability in Blood Pressure

Contrary to our expectations, we found that the mean coherence (Cohmean) was best to discriminate among infants (i.e., had the highest precision). Thus, in other words, adjusting for variability in ABP within each measurement did not increase the signal-to-noise ratio. The same tendency was seen in our simulation study where fewer measurements (less monitoring time) were needed to discriminate among infants based on the mean coherence (Cohmean) . We ascribe the lower precision of the CohMAPwmean and CohMAPmax compared to Cohmean to the effect of the signal-to-noise ratio in MAP. The noise is transferred to the power spectrum by the Fourier transform and gives rise to uncertainty in determining the weighting and the frequency with maximum power, respectively.

However, in line with our hypothesis, precision improved when we adjusted for the varying degree of variability in ABP between measurements in each infant (intermeasurement ABP weighting). Thus, in other words, weighting measurements with large fluctuations in ABP in favor of those with small actually increased the signal-to-noise ratio. Unfortunately, a power study using simulated data sets could not be conducted to explore how this weighting would affect the monitoring time because it would require simulating an ABP signal.

To our knowledge, no previous studies have investigated the influence of variability in ABP on estimated coherence between spontaneous fluctuations in ABP and cerebral NIRS measurements. Our results, however, are in line with previous studies of cerebral autoregulation based on transfer function analysis between Doppler measurements of cerebral blood flow velocity and spontaneous variations in ABP.14, 15 It has been proposed to set a minimum variability in ABP for a more precise estimation of cerebral autoregulation from spontaneous fluctuations in ABP.5, 14 However, this approach will result in data loss, because in practice, the variability in ABP can be quite low over a considerable length of time. It is important to use data as effectively as possible when research involves vulnerable patients, such as premature infants. Hence, our approach with weighting coherence with variability in ABP in each measurement seems relevant for research studies comparing cerebral autoregulation among vulnerable patients.


Frequency-Dependent Functional Relationship

Simulations to estimate the statistical power of the various signal analytical approaches showed a remarkable difference in statistical power in the two frequency ranges. Statistical power was almost three times bigger in LF compared to VLF. From a practical point of view, this indicates that a three-times-longer recording time is required in VLF compared to LF. We speculate that this finding might be attributed to the fact that other regulatory mechanisms for cerebral hemodynamics have more influence in the VLF band (i.e., 0.0030.4Hz , 25300scycle length) compared to the LF band (i.e., 0.040.1Hz , 1025scycle length). This will result in a lower signal-to-noise ratio. Because 91% of our infants were breathing spontaneously, irregular fluctuations in PaCO2 might have contributed to the lower statistical power in the VLF band.

The lack of concordance between VLF and LF (Fig. 5) implies that they might represent different aspects or time scales of cerebral autoregulation. Because the autoregulatory ability is controlled by vascular smooth muscle cells, we speculate that an injurious stimulus initially reduces the response time and then eventually destroys the ability to react to changes in blood pressure. Consequently, an infant might actually have a bad autoregulation at higher frequencies and a good one at lower frequencies (upper left quadrant on Fig. 5). However, a higher amount of physiological noise in VLF compared to LF could also contribute to the lack of concordance by falsely indicating a good autoregulation at low frequencies. Still, further analyses are needed to clarify if the two frequency bands contribute with shared or different information.


Nature of Cerebral Autoregulation

Our results indicated a range of levels of coherence, even in these relatively healthy preterm infants. Furthermore, the drift in individual infants was minimal. In an earlier work by Soul,8 who used the same technique to detect cerebral autoregulation in preterm infants, data were analyzed and interpreted differently. The authors computed a pressure-passive index (i.e., percentage of all 10-min measurements that were classified as coherent) to characterize autoregulation in each infant. The computation of Soul’s study was based on an assumption of autoregulation being “on” or “off” according to the mean coherence being below or above the threshold of statistical significance. Our data show that single values of coherence from 10-min measurements have an imprecision that exceeds the variability among infants. Consequently, we believe that the precision of this measuring method is too low to get a meaningful estimate from a single measurement. Similarly, using the maximum value of coherence (i.e., picking a single measurement among several) to characterize each infant’s level of autoregulation is not an effective nor unbiased way of data analysis. In two previous studies, this was done and the maximum coherence was associated with mortality7 and occurrence of IVH in premature infants.9

An imperfect cerebral autoregulation (i.e., coherence exceeding the limit of statistical significance) was detected in approximately one-half of our stable infants, and no correlation was found to disease (i.e., IVH, neonatal death, or use of Dopamine). The fact that an imperfect cerebral autoregulation did not correlate with disease was unexpected. However, the present study was not powered to test for this correlation. The proportion of infants with an imperfect cerebral autoregulation is in line with earlier studies using the same technique to study “dynamic” cerebral autoregulation in premature infants.7, 9 Contrarily, studies using a “static” measuring approach have, in general, shown an intact cerebral autoregulation in stable premature infants and an impaired autoregulation in sick premature infants. The discrepancy between studies based on continuous cerebral NIRS and studies based on single measurements of cerebral blood flow by means of either xenon-133 clearance technique16, 17 or cerebral NIRS18, 19 may be apparent. The conclusion that autoregulation is present in stable premature infants in the latter studies came from a statistically insignificant relation between arterial blood flow and cerebral blood flow. In fact, the 95%-confidence limit was +2%mm Hg (i.e., compatible with an imperfect autoregulation that would give rise to a statistically significant coherence with the techniques discussed here).


Strength and Limitations


Appropriate threshold of significance

Strength of our technique is that we used a measurement-specific threshold of significant coherence based on Monte Carlo simulations. Previous reports in this field using mean analysis of coherence have either used a conventional threshold of significant coherence of 0.57, 9, 20 or a calculated threshold8 based on a published mathematical definition.21 However, a conventional threshold cannot be generalized because it does not account for differences in analytical conditions such as length of measurements and segments, number of harmonics in the frequency range, type of tapering window applied, or degree of overlap between segments. These differences affect the degrees of freedom of the estimated coherence and thus, the threshold of significant coherence.5, 22 Besides, the published mathematical definition is only valid for point measurements.21 Our ThresholdCOH for Cohmean and CohMAPwmean happened to be close to the conventional limit of 0.5. ThresholdCOH for our point measurement (CohMAPmax) was almost identical to a calculated threshold of 0.78 based on our analytical conditions [i.e., sampling rate (2Hz) , length of measurements ( 10min=1200 samples), and length of segments ( 5min=600 samples)]. The fact that our ThresholdCOH showed very low standard deviations indicates that it may not be necessary to run Monte Carlo simulations on all of the recorded data. A previous study23 has claimed an advantage of this method because it maintains the statistical distribution of the original data. However, we found that the threshold is only changed significantly by the way the coherence is calculated. Hence, we propose that a reliable threshold of significant coherence can be achieved by running Monte Carlo simulations on a single recording representing the analytical conditions.


Oxygenation versus blood flow

We used the oxygenation index as a surrogate measure of cerebral blood flow. This assumption is only reasonable if arterial oxygenation is stable during measurements. We kept the impact of variations in arterial oxygenation to a minimum by rejecting all data with changes in arterial saturation exceeding 5%. The cerebrovascular transit time (i.e., the time delay in reaching equilibrium of oxygenation as response to changes in cerebral blood flow) represents another potential limitation when oxygenation index is used as a surrogate measure of cerebral blood flow.24 However, because cerebrovascular transit time in neonates is 10s ,25 this limitation would not affect analysis in our frequency range. The fact that cerebral oxygenation is “volume weighted” and not “flow weighted”26 represents another problem. This is relevant if the volume of tissue is heterogeneous, as is the case with the brain. Change in mean blood flow is dominated by changes to high flow regions, whereas change in mean (vascular) oxygenation is dominated by changes in high blood volume regions. Consequently, changes in oxygenation will only truly reflect changes in blood flow if the following conditions are met: (i) changes in blood flow to all parts of the volume of tissue are proportional and (ii), changes in blood flow and blood volume are proportional. These are reasonable assumptions in normal conditions and during small changes, but it is difficult to know.


Other factors influencing cerebral blood flow

Another potential limitation of our approach lies in the very concept of using coherence to estimate cerebral autoregulation. Coherence assumes linearity between ABP and cerebral blood flow (i.e., it assumes ABP to drive cerebral blood flow), thereby ignoring other potential input signals. Thus, it does not account for a possible effect of other physiological variables on cerebral blood flow (e.g., PaCO2 ). We tried to keep this problem to a minimum by only investigating infants who were clinically stable.


Statistical versus physiological significance

Strictly speaking a coherence exceeding the threshold of statistical significance only indicates that autoregulation is not working perfectly (i.e., that the slope of cerebral blood flow versus ABP is not precisely zero). The quantitative assessment of this slope (i.e., the actual value of autoregulation) requires estimation of the operational gain (i.e., magnitude of change in cerebral oxygenation per change in ABP). Because autoregulation is always imperfect to some degree, a statistical significant coherence could correspond to an imperfect but not impaired cerebral autoregulation.



A reliable detection of cerebral autoregulation by means of analysis of coherence between cerebral NIRS and spontaneous fluctuations in ABP requires several hours of monitoring. The precision can be improved by adjusting for the varying degree of blood pressure variability between repeated measurements. These findings are relevant for research studies comparing cerebral autoregulation among vulnerable patients.


We thank the Ludvig and Sara Elsass Foundation and the University of Copenhagen, Denmark, for funding the study. We also thank the participating infants and their parents for contributing to this study. Finally, we thank the staff of the Neonatal Department and the Medical Engineering Department at Rigshospitalet for practical and technical assistance, respectively.



J. E. Brazy, D. V. Lewis, M. H. Mitnick, and V. Jobsis, “Noninvasive monitoring of cerebral oxygenation in preterm infants: preliminary observations,” Pediatrics, 75 (2), 217 –225 (1985). 0031-4005 Google Scholar


G. Greisen, “Autoregulation of cerebral blood flow,” NeoReviews, 8 (1), 22 –31 (2007). Google Scholar


J. J. Volpe, “Neurobiology of periventricular leukomalacia in the premature infant,” Pediatr. Res., 50 (5), 553 –562 (2001). 0031-3998 Google Scholar


C. A. Giller and D. G. Iacopino, “Use of middle cerebral velocity and blood pressure for the analysis of cerebral autoregulation at various frequencies: the coherence index,” Neurol. Res., 19 (6), 634 –640 (1997). 0160-6412 Google Scholar


R. B. Panerai, “Cerebral autoregulation: from models to clinical applications,” Cardiovasc. Eng., 8 (1), 42 –59 (2008). 1567-8822 Google Scholar


E. G. Ramos, D. M. Simpson, R. B. Panerai, J. Nadal, J. M. Lopes, and D. H. Evans, “Objective selection of signals for assessment of cerebral blood flow autoregulation in neonates,” Physiol. Meas, 27 (1), 35 –49 (2006). 0967-3334 Google Scholar


F. Y. Wong, T. S. Leung, T. Austin, M. Wilkinson, J. H. Meek, J. S. Wyatt, and A. M. Walker, “Impaired autoregulation in preterm infants identified by using spatially resolved spectroscopy,” Pediatrics, 121 (3), e604 –e611 (2008). 0031-4005 Google Scholar


J. S. Soul, P. E. Hammer, M. Tsuji, J. P. Saul, H. Bassan, C. Limperopoulos, D. N. Disalvo, M. Moore, P. Akins, S. Ringer, J. J. Volpe, F. Trachtenberg, and A. J. du Plessis, “Fluctuating pressure-passivity is common in the cerebral circulation of sick premature infants,” Pediatr. Res., 61 (4), 467 –473 (2007). 0031-3998 Google Scholar


M. Tsuji, J. P. Saul, P. A. du, E. Eichenwald, J. Sobh, R. Crocker, and J. J. Volpe, “Cerebral intravascular oxygenation correlates with mean arterial pressure in critically ill premature infants,” Pediatrics, 106 (4), 625 –632 (2000). 0031-4005 Google Scholar


M. Tsuji, A. duPlessis, G. Taylor, R. Crocker, and J. J. Volpe, “Near infrared spectroscopy detects cerebral ischemia during hypotension in piglets,” Pediatr. Res., 44 (4), 591 –595 (1998). 0031-3998 Google Scholar


J. S. Soul, G. A. Taylor, D. Wypij, A. J. Duplessis, and J. J. Volpe, “Noninvasive detection of changes in cerebral blood flow by near-infrared spectroscopy in a piglet model of hydrocephalus,” Pediatr. Res., 48 (4), 445 –449 (2000). 0031-3998 Google Scholar


A. Pryds, J. Tonnesen, O. Pryds, G. M. Knudsen, and G. Greisen, “Cerebral pressure autoregulation and vasoreactivity in the newborn rat,” Pediatr. Res., 57 (2), 294 –298 (2005). 0031-3998 Google Scholar


A. Tweed, J. Cote, H. Lou, G. Gregory, and J. Wade, “Impairment of cerebral blood flow autoregulation in the newborn lamb by hypoxia,” Pediatr. Res., 20 (6), 516 –519 (1986). 0031-3998 Google Scholar


J. Liu, D. M. Simpson, and R. Allen, “High spontaneous fluctuation in arterial blood pressure improves the assessment of cerebral autoregulation,” Physiol. Meas, 26 (5), 725 –741 (2005). 0967-3334 Google Scholar


D. M. Simpson, R. B. Panerai, E. G. Ramos, J. M. Lopes, M. N. Marinatto, J. Nadal, and D. H. Evans, “Assessing blood flow control through a bootstrap method,” IEEE Trans. Biomed. Eng., 51 (7), 1284 –1286 (2004). 0018-9294 Google Scholar


D. Jayasinghe, A. B. Gill, and M. I. Levene, “CBF reactivity in hypotensive and normotensive preterm infants,” Pediatr. Res., 54 (6), 848 –853 (2003). 0031-3998 Google Scholar


O. Pryds, G. Greisen, H. Lou, and B. Friis-Hansen, “Heterogeneity of cerebral vasoreactivity in preterm infants supported by mechanical ventilation,” J. Pediatr., 115 (4), 638 –645 (1989). Google Scholar


L. Tyszczuk, J. Meek, C. Elwell, and J. S. Wyatt, “Cerebral blood flow is independent of mean arterial blood pressure in preterm infants undergoing intensive care,” Pediatrics, 102 (2), 337 –341 (1998). 0031-4005 Google Scholar


M. J. Munro, A. M. Walker, and C. P. Barfield, “Hypotensive extremely low birth weight infants have reduced cerebral blood flow,” Pediatrics, 114 (6), 1591 –1596 (2004). 0031-4005 Google Scholar


H. Bassan, K. Gauvreau, J. W. Newburger, M. Tsuji, C. Limperopoulos, J. S. Soul, G. Walter, P. C. Laussen, R. A. Jonas, and A. J. du Plessis, “Identification of pressure passive cerebral perfusion and its mediators after infant cardiac surgery,” Pediatr. Res., 57 (1), 35 –41 (2005). 0031-3998 Google Scholar


J. A. Taylor, D. L. Carr, C. W. Myers, and D. L. Eckberg, “Mechanisms underlying very-low-frequency RR-interval oscillations in humans,” Circulation, 98 (6), 547 –555 (1998). 0009-7322 Google Scholar


R. Bortel and P. Sovka, “Approximation of statistical distribution of magnitude squared coherence estimated with segment overlapping,” Signal Process., 87 (5), 1100 –1117 (2007). 0165-1684 Google Scholar


R. B. Panerai, P. J. Eames, and J. F. Potter, “Multiple coherence of cerebral blood flow velocity in humans,” Am. J. Physiol. Heart Circ. Physiol., 291 (1), H251 –H259 (2006). 0363-6135 Google Scholar


F. Y. Wong, M. Nakamura, T. Alexiou, V. Brodecky, and A. M. Walker, “Tissue oxygenation index measured using spatially resolved spectroscopy correlates with changes in cerebral blood flow in newborn lambs,” Intensive Care Med., 35 (8), 1464 –1470 (2009). 0342-4642 Google Scholar


C. E. Elwell, M. Cope, and D. T. Delpy, “An analytical method for determining cerebrovascular transit time using near infrared spectroscopy,” Adv. Exp. Med. Biol., 428 561 –568 (1997). 0065-2598 Google Scholar


M. Smith and C. Elwell, “Near-infrared spectroscopy: shedding light on the injured brain,” Anesth. Analg., 108 (4), 1055 –1057 (2009). Google Scholar
©(2010) Society of Photo-Optical Instrumentation Engineers (SPIE)
Gitte Holst Hahn, Karl Bang Christensen, Terence S. Leung, and Gorm Greisen "Precision of coherence analysis to detect cerebral autoregulation by near-infrared spectroscopy in preterm infants," Journal of Biomedical Optics 15(3), 037002 (1 May 2010).
Published: 1 May 2010

Back to Top