- Research
- Open access
- Published:
A multivariate time-frequency method to characterize the influence of respiration over heart period and arterial pressure
EURASIP Journal on Advances in Signal Processing volume 2012, Article number: 214 (2012)
Abstract
Respiratory activity introduces oscillations both in arterial pressure and heart period, through mechanical and autonomic mechanisms. Respiration, arterial pressure, and heart period are, generally, non-stationary processes and the interactions between them are dynamic. In this study we present a methodology to robustly estimate the time course of cross spectral indices to characterize dynamic interactions between respiratory oscillations of heart period and blood pressure, as well as their interactions with respiratory activity. Time-frequency distributions belonging to Cohen’s class are used to estimate time-frequency (TF) representations of coherence, partial coherence and phase difference. The characterization is based on the estimation of the time course of cross spectral indices estimated in specific TF regions around the respiratory frequency. We used this methodology to describe the interactions between respiration, heart period variability (HPV) and systolic arterial pressure variability (SAPV) during tilt table test with both spontaneous and controlled respiratory patterns. The effect of selective autonomic blockade was also studied. Results suggest the presence of common underling mechanisms of regulation between cardiovascular signals, whose interactions are time-varying. SAPV changes followed respiratory flow both in supine and standing positions and even after selective autonomic blockade. During head-up tilt, phase differences between respiration and SAPV increased. Phase differences between respiration and HPV were comparable to those between respiration and SAPV during supine position, and significantly increased during standing. As a result, respiratory oscillations in SAPV preceded respiratory oscillations in HPV during standing. Partial coherence was the most sensitive index to orthostatic stress. Phase difference estimates were consistent among spontaneous and controlled breathing patterns, whereas coherence was higher in spontaneous breathing. Parasympathetic blockade did not affect interactions between respiration and SAPV, reduced the coherence between SAPV and HPV and between respiration and HPV. Our results support the hypothesis that non-autonomic, possibly mechanically mediated, mechanisms also contributes to the respiratory oscillations in HPV. A small contribution of sympathetic activity on HPV-SAPV interactions around the respiratory frequency was also observed.
Introduction
Non-stationary signal processing techniques applied to electrophysiological signals have been successful in the assessment of important features of cardiovascular control physiology.
Although it is well known that respiratory activity affects cardiovascular regulation[1–4], the mechanisms responsible for the coordination between the autonomic control of circulation and respiration are still unclear and currently matter of debate[5–7]. Respiration affects both heart period length and arterial pressure, which in turn are mutually related. Interconnections between these three physiological parameters can be broadly reduced to four interactions: (i) Blood pressure decreases during inspiration and increases during expiration, following the changes in intrathoracic pressure. (ii) Heart period decreases during inspiration and increased during expiration, a phenomenon called respiratory sinus arrhythmia (RSA)[1, 2], that may be due to direct central nervous modulation, reflex phenomena and mechanical influence of respiration[5–7]. (iii) A decrease (increase) in blood pressure provokes a decrease (increase) in heart period, via the baroreflex, a negative feedback mechanism that buffers short-term fluctuations in arterial pressure. (iv) A change in heart period causes a change in blood pressure through direct mechanical effects. As a result, the respiratory induced oscillations in the heart period and arterial pressure variability are mutually related.
To assess the functioning of the cardiovascular system in both normal and impaired conditions, it is important to reliably characterize the dynamic interactions between physiological signals involved in its regulation, i.e., heart period, arterial pressure and respiration. Cross spectral analysis has been widely used to study the cardiovascular system[3, 8–10]. Cross spectral indices include spectral coherence, partial spectral coherence, phase difference, and time delay. Spectral coherence measures the strength of the coupling between the spectral components of two signals. Partial coherence measures the coherence between two signals after removing the influence that a third signal has over them[11, 12]. Phase difference and time delay measure the latencies between two correlated oscillations. A limitation of spectral analysis is that it requires signals to be stationary and does not provide a dynamic characterization of the interactions among them. The purpose of this study is to present a methodology to robustly estimate the time course of cross spectral indices to characterize dynamic interactions between cardiovascular signals in non-stationary conditions. In particular, an assessment of the coupling between respiratory induced oscillations of heart period and blood pressure, as well as their interactions with respiratory activity is proposed. The characterization is based on the estimation of the time course of coherence, partial coherence, phase difference and time delay. The purpose of this study is to present a methodology to robustly estimate the time course of cross spectral indices to characterize dynamic interactions between respiration, heart period and arterial pressure, simultaneously. The characterization of the cardiovascular function is based on the estimation of the time course of coherence, partial coherence, phase difference, and time delay. Recently, distributions belonging to the Cohen’s class were used to estimate time-frequency (TF) coherence between cardiovascular signals[13], and time-varying baroreflex sensitivity[14]. In this article, we present a comprehensive framework which also includes TF phase difference and TF partial coherence to simultaneously characterize the dynamic interactions between these signals during tilt table test in spontaneous and controlled breathing, and after autonomic blockade.
Study populations
In this study, we used data from two experiments. In the first experiment, subjects performed a tilt table test under spontaneous respiration[15]. In the second experiment, subjects performed both a tilt table test and autonomic blockade under controlled respiration[3]. Tilt table test is one of the most well-established test for the assessment of cardiovascular autonomic function[9], while autonomic blockade data offers the unique possibility of studying the influence of sympathetic and parasympathetic activity separately[3]. As explained in the following, in this study we divided these data in three groups.
Tilt table test during spontaneous respiration (TTSR)
Fourteen subjects (9 males, 24–34 yr, median 28 yr) without any previous cardiovascular history underwent a head-up tilt table test according to the following protocol: 4 min in early supine position (TES), 5 min tilted head-up to an angle of 70° (THT) and 4 min back to later supine position (TLS)[15].
12-lead ECG were recorded with Biopac MP150 system, with a sampling frequency of 1000 Hz. Blood pressure was measured at the finger by means of Finometer device with sampling frequency of 250 Hz. During the procedure, the Finometer was recalibrated at the beginning of THT and TLS. The recalibration took few seconds and introduced artefacts which were detected and corrected by interpolation. Arterial pressure from the finger was not corrected for the hydrostatic gradient change during tilt. Respiratory signal was recorded with a sampling frequency of 125 Hz, by using TSD201 transducer which measure thoracic expansion while breathing. This signal gives a measure correlated with lung volume changes.
Tilt table test during controlled respiration (TTCR)
Data set TTCR[3, 4], composed of 14 subjects undergoing a tilt table test, is used as control for the results obtained in TTSR. Fourteen subjects (aged 19–38 year, median 21 years) participated in the study according to the following protocol: 13 min in supine position and 13 min head-up tilted. During the entire experiment, subjects were instructed to breath at a controlled rate, by following auditory cues. The time between the onset of two consecutive inspirations was randomly chosen from a Poisson distribution with mean 0.2 Hz. These data were originally used in[3], with the purpose of assessing the frequency response of the cardiovascular system. The rationale for randomizing the respiratory cycle was to whiten the power spectral density of the respiratory signal, which is seen as the input of the cardiovascular system.
Lung volume changes were measured with a two belt chest-abdomen inductance plethysmograph (Respitrace Systems). Arterial pressure was measured with a 22-gauge Teflon catheter in the radial artery of the nondominant hand. For each participant, intervals of about 7 min, TES and THT, recorded during supine and head-up tilt conditions, respectively, were selected for the analysis.
Autonomic blockade during controlled respiration (ABCR)
After having been tilted, subjects of TTCR data set were returned to the supine position and given either atropine or propranolol. Atropine (0.03 mg/kg) was used for parasympathetic blockade in seven subjects, while propranolol (0.2 mg/kg) was used to sympathetic blockade in the other seven subjects. Physiological data were recorded in supine position following the same controlled respiratory protocol used in TTCR data set. Intervals of about 7 min, where stationarity was assumed were selected for the analysis of blockade condition. Data from TTCR recorded in supine position (TES) are used as control for autonomic blockade conditions.
Methods
Signal processing
QRS complexes in the ECG were detected by using a methodology described in[16]. Instantaneous heart rate was derived by integral pulse frequency modulation (IPFM) model, which also accounts for the presence of ectopic beats[17], and then evenly resampled at 4 Hz, using spline interpolation. Instantaneous heart period was then obtained as the reciprocal of instantaneous heart rate. The heart period variability (HPV) signal, xH(t), was obtained by high pass filtering the heart period signal with a cut-off frequency of 0.03 Hz. The systolic arterial pressure series was obtained by taking the maximum of the pressure signal within a short interval following a QRS detection. The time series were subsequently interpolated at the time of occurrence of the systolic peak by splines with a sampling frequency of 4 Hz, and the SAPV signal, xS(t), was obtained by high-pass filtering with a cut-off frequency of 0.03 Hz. The signal from the thoracic belt was decimated to 4 Hz to get the respiratory signal, xR(t), which gives a measure correlated with instantaneous lung volume.
Cross time-frequency analysis
In the following, signals {xi(t),xk(t)} ∈ {xH(t),xS(t), − xR(t)}, indicate the complex analytical signal representation of HPV, SAPV and inverse respiratory signal, respectively. Owing to the inversion of the respiratory signal, and under physiological conditions, all the members of the triplet {xH(t),xS(t),−xR(t)} are expected to increase and decrease together.
The TF cross spectrum, Sik(t,f), is estimated by using a TF distribution (TFD) belonging to the Cohen’s class[18]:
where Aik(τ ν) is the narrow-band symmetric ambiguity function[19] of signals xi(t) and xk(t), that in (1) is windowed by an elliptical exponential kernel, defined in the ambiguity function domain as[13]:
The kernel function ϕd-D(τ,ν) can be equivalently defined in the TF domain as:
The choice of the kernel function is discussed in Section ‘Time-frequency filtering’.
Time-frequency coherence is estimated as[13]:
Time-frequency coherence quantifies the strength of the local coupling between two non-stationary signals, being γik(t,f) = 1 in the TF regions where the signals are perfectly coupled and γik(t,f) = 0 in the regions where signals are uncorrelated. Note that in[13] it was shown that the local averaging performed in (1) can be used to estimate non-stationary spectra, whose 4 definition includes expectation over different realizations of a given process. This implies that (5) can be used to estimate coherence function from only one pair of signals. Partial coherence is used to assess the coupling of two signals xi(t) and xk(t), after having removed the influence of a third signal xz(t)[20]. Time-frequency partial coherence can be defined as:
In the TF regions where xz(t) is uncorrelated with xi(t) and xk(t), i.e., wherever Szi(t,f) = 0 and Szk(t,f) = 0, TF partial coherence is equal to TF coherence, γik/z(t,f) = γik(t,f). Furthermore, partial coherence vanished, γik/z(t,f) = 0, wherever xi(t) = a xz(t) and xk(t) = b xz(t), with. For any triplet of signals, the interpretation of TF partial coherence is as follows: if around a TF point (t0,f0) γik/z(t,f) < < γik(t,f), then it follows that the TF structure of signal xz(t) matches with that of xi(t) and xk(t) around (t0,f0). Moreover, if in a given TF region γik/z(t,f) < γiz/k(t,f) < γkz/i(t,f), it follows that xz(t) better represents the communality shared by the three signals[21].
Partial spectra in (6) can be obtained as:
Expression (7) shows that Sik(t,f) and Sik/z(t,f) are complex functions characterized by same phase and different magnitude. Time-frequency phase difference (TFPD) spectrum is defined as:
Within this framework, in a given TF region, a change in xi(t) precedes (leads) a correlated change in xk(t) wherever θi k(t,f) ∈[0,Π], while a change in xi(t) lags behind a correlated change in xk(t) wherever θik(t,f) ∈[−Π,0].
Note that according to the open-loop assumption of cross spectral analysis, θik(t,f) = −θki(t,f). Moreover, phase spectra θRS(t,f), θRH(t,f) and θSH(t,f) are related by θSH(t,f) = θRH(t,f)−θRS(t,f). Finally, it is worth mentioning that, given two signals {xi(t),xk(t)}, if xz(t) = −xk(t) ⇒ θiz(t,f) = θik(t,f) ± Π.
Time-frequency filtering
The kernel function determines the degree of TF filtering, and, consequently, the TF resolution of the spectra and the interference terms (ITs) reduction. In this study, time resolution is quantified by Δt, the full width at half maximum of ϕt-f(t,0), while frequency resolution is quantified by Δf, the full width at half maximum of ϕt-f(0,f). These quantities measure the degree of spreading of a line in the TF domain: Δt and Δf are equal to the full width at half maximum of the TFD of a Dirac impulse, evaluated along t for a given frequency, and of a sinusoid, evaluated along f for a given time instant, whose ideal TF representations would be straight lines[13]. The kernel function used in this study gives a TF resolution of {Δt,Δf} = {10.9 s,0.039 Hz}.
The choice of the kernel is especially important in coherence analysis, because to obtain reliable coherence estimates, i.e., γik(t,f) ∈[0,1], the filtering provided by the kernel should be able to completely remove the ITs that characterize the Wigner-Ville distribution[13, 22]. A necessary, but not sufficient, condition to obtain reliable coherence estimates is the positiveness of the auto-spectra[13]. This condition is not sufficient, since residual (oscillating) ITs present in positive TFDs may cause coherence estimates to be higher than 1. Reducing ITs and obtaining reliable coherence estimates are two aspects of the same problem. Thus, our strategy consists in finding a kernel function able to provide γik(t,f) ∈[0,1], which in turn implies ITs canceling.
In this study, we used a kernel function (3) that is a particular case of the multiform, tiltable exponential kernel proposed in[23]. Among all the possibilities given by this kernel, we used a function whose isocontours are, in the ambiguity function domain, ellipses with major and minor axes along τ and ν[24]. The choice of an elliptical shape is motivated by its good concentration around the origin of the ambiguity function domain (where auto-terms are located). Parameters ν0 and τ0 are used to change the length of the ellipse axes aligned along ν (the degree of time filtering) and τ (the degree of frequency filtering), respectively. The parameter λ sets the roll-off of the filter as well as the size of the tails of the kernel. Several simulation studies demonstrated the effectiveness of this kernel in reducing ITs[13, 14, 25, 26]. In particular, in[13], it was shown that this kernel offers the possibility of obtaining TF coherence estimates bounded between 0 and 1, as well as spectra characterized by a better TF resolution than multitaper spectrogram and continuous wavelet transform. The simulation study carried out in Section ‘Simulation study’ confirms these previous results.
The choice of the parameters was made as follows: First, the desired TF resolution {Δt,Δf}, corresponding to the minimum amount of TF filtering, is decided based on a-priori information about the signals and the experimental settings. The set of parameters {τ0, ν0, λ} that provides the desired TF resolution is used as starting point. If using this set of parameters γik(t,f) ∉[0,1], the degree of time (or frequency) filtering is maintained constant, while the frequency (or time) filtering is increased until reaching meaningful estimates over the entire TF domain. If at the end of the process, the frequency (or time) resolution is not satisfactory, the time (or frequency) resolution is decreased, i.e., the corresponding ν0 (or τ0) is increased, and the process iterates. This process allows to adjust the TF filtering to the specific needs of analysis. In this study, we used {τ0, ν0, λ} = {0.05,0.046,0.300}, with 2048 frequency points (corresponding to {Δt,Δf} = {10.9s,0.039Hz}), because among all the explored combination of parameters which gave γik(t,f) ∈[0,1], this was associated to the minimum degree of TF filtering. An example of this scheme was given in[14].
Estimation of synchronization indices
The time course of coherence, partial coherence, and phase difference is extracted by averaging the corresponding TF representation in specific TF regions.
The region where γik(t,f) is significant, i.e., that where the two signals are sharing approximately the same instantaneous frequencies, is defined as:
where γTH(t,f) is a threshold function, estimated by using surrogate data[13], and which depends on the TF resolution of the spectra. Briefly, γTH(t,f) is obtained by estimating the TF coherence between several realizations of uncorrelated white Gaussian noises, and taking at each TF point the 95th percentile of TF coherence estimates[13].
Given that we are interested in assessing the influence of respiration on HPV and SAPV, the TF region where the time course of spectral coherence is estimated is centered around respiratory rate, fR(t), and is defined as:
where Δf is a term related to the frequency resolution. Respiratory rate, fR(t), is estimated from the TF spectrum of the respiratory signal, as the frequency corresponding to the maximum of the instantaneous spectral peak.
The TF region where the time course of phase difference index is estimated is composed of those part of Ω(γ) in which coherence is statistically significant:
In this expression, R(t,f) is a rectangle of sides 2 s Hz and ∘ denotes the opening (processing technique which involves erosion and dilation). The opening excludes from the portions of TF domain that are smaller than R(t,f), thus adding robustness to the final estimates.
The time course of the band coherence, as well as the time course of partial coherence, is then obtained by averaging γik(t,f) and γik/z(t,f) in Ω(γ):
Index θik(t) is estimated (in radians) by averaging the TFPD spectrum in:
The time delay associated to θik(t) is estimated (in seconds) by index, defined as:
A method to reduce the uncertainty of phase difference in non-stationary signals
The cross spectra Sik(t,f) are complex functions, and as such, their phase is θik(t,f) = arg[Sik(t,f)exp(j 2nΠ)], with. The periodicity of 2Π introduces an uncertainty over the actual value of θik(t,f), which may prevent one from drawing conclusions about the temporal sequence of events described by xi(t) and xk(t). For instance, it is not possible to determine whether xi(t) precedes or lags behind xk(t), since values for θik(t) + n 2Π and, where TR(t) = 1/fR(t), are positive for n > 1 and negative for n < 1. Although in cardiovascular applications the range of values for n is usually reduced to n = {−1,0,1} by considering physiological information[10], the uncertainty still remains. In the contest of non-stationary signals, we propose to use TF coherence estimates as control parameters to reduce this uncertainty. The idea is that, the lowest the time delay between non-stationary spectral components, the highest the TF coherence. This is due to the fact that TF coherence is a measure of local correlation. Thus, to determine the actual time delay among, that is associated to phase difference θik(t), one can use the following procedure: (i) Generate a set of pairs of delayed signals {xi(t),xk(t + r t0)}, with and being t0 a small time delay. (ii) Estimate γik(t;r), θik(t;r) and between each pair of signals for each r, as well as their temporal median, and. (iii) Find rm, as the sample for which is maximal. (iv) Among the possible time delays, the closest to rmt0 is the actual one.
Statistical analysis
Statistical analysis of each data set is performed as follows. The time course of each general index from one subject s is denoted as, with,, with {xi(t),xk(t)} ∈ {xH(t),xS(t),−xR(t)}, and s ∈ {1,…,N}, being N the number of subjects.
The median time course of,, estimated across subjects, as well as the interquartile range, is used to describe the pattern of response of the population during a given condition.
The temporal median values of indices, denoted as, are estimated during epochs where stationarity is assumed, and are used to: (i) Assess inter-conditions differences, i.e., whether estimated during a given condition are statistically different from estimated during another condition. (ii) Assess the inter-indices differences, i.e., whether estimated during a given condition are statistically different from estimated during the same condition, with (i,k) ≠ (p,q).
Pairwise comparisons between the same indices evaluated in different epochs are performed by using the Wilcoxon signed rank test, while comparisons between different indices or between the same indices but in different data sets are performed by using the Wilcoxon ranksum test. Statistical significance is assumed for P < 0.05.
Results
Simulation study
A simulation study was carried out to verify that TFDs can be used to correctly estimate the time delay between two non-stationary spectral components of signals xi(t) and xk(t), even when θik(t) > 2Π. The SAPV signal from subjects of TTSR data were used to create pairs of signals characterized by a known time delay:
where, xS(t,s) is the complex analytical signal representation of SAPV of subject s, is the median value of the respiratory rate, and u1(t) and u2(t) are two zero-mean complex white Gaussian noises associated to a SNR equal to 10 dB. The procedure described in the previous section was applied, and coherence, phase difference and time delay between {x1(t),x2(t + r t0)}, with r ∈ {−30,−29,…,30} and t0 = 0.5 s, were estimated.
Results are shown in Figure1. In panel (a), the results from only two signals {x1(t),x2(t)}, obtained from the same xS(t), are shown. The highest coherence estimates (median value) were obtained for a time shift r t0 = 7.5 s, that is very close to the actual time delay s. As expected, for r t0 = 7.5 s, phase differences and time delay were almost zero.
For n ∈ {−1,0,1}, median time delay estimates were equal to {−2.49,2.37,7.22} s, respectively. Given that coherence estimates assessed for the closest r t0 to, {−2.50,2.50,7.00} s, were equal to {0.884,0.936,0.997}, phase difference and time delay were chosen for n = 1, i.e., θ12(t) = θ12(t) + 2Π and s. In Figure1b, global results obtained considering all 14 pairs of signals are shown. Time delays are plotted over their estimates, with n ∈ {−1,0,1}. Owing to the assessment of the coherence estimates, the correct time delay values, those for n = 1, were estimated.
These results confirm those from previous studies, and demonstrate that the methodology described in Sections ‘Cross time-frequency analysis’ offers the possibility of estimating TF coherence and time-varying phase difference accurately[14, 25, 26].
Dynamic interactions during tilt table test
Results from TTSR data set
An illustrative example of subject-based cross TF analysis of signals {xi(t),xk(t)} ∈ {xH(t),xS(t),−xR(t)} by TFD is shown in Figure2. Time-frequency coherence functions γik(t,f), partial coherence γik/z(t,f) and phase difference spectra θik(t,f) are shown in panels a–c, d–f, and g–i, respectively. The time course of the synchronization indices are shown below each TF representation. These indices changed during time, showing that correlations between cardiovascular oscillations are time-varying. Although the respiratory rate fR(t) changed fast during the entire test, the signals were highly correlated inside Ω(γ), which in a–c and d–f is delimited by white and black contour, respectively. The time-course of coherence indicates that the strength of the coupling, although high, was intermittent (see γik(t) around 250 and 400 s). As expected, in (t,f) ∈ Ω(γ), partial coherence γSH/R(t,f) was lower than γSH(t,f). Importantly, the influence of respiration was removed even when the respiratory rate was inside the traditional LF range, fR(t) < 0.15 Hz. The same held true for γRS/H(t,f), coherence between −xR(t) and xS(t) after having removed xH(t), and γRH/S(t,f), coherence between −xR(t) and xH(t) after having removed xS(t). Region, which in panels g–i is delimited by a black contour, only includes portions of Ω(γ), where coherence estimates were statistically significant. In this particular example, differences between Ω(γ) and are visible during the transition from supine to head-up tilt positions and around 400 s. Interestingly, in this subject, the position changes caused fast changes in partial coherence, phase difference and time delay indices.
Indices extracted from population TTSR are given in Figure3, where black lines and shadowed area represent median trends and interquartile ranges across the 14 subjects, respectively. For all indices, the interquartile range is low, indicating that in all subjects cardiovascular signals are related by a common mechanism. In median, spectral coherence, shown in Figure3a–c, was very high during the entire test and decreased after the position changes. The interquartile range of γRS(t,s) was lower than that of γSH(t,s) and γRH(t,s), indicating that the coupling between respiration and SAPV was less affected by inter subject variability. Partial spectral coherence, γik/z(t,s) shown in Figure3d–f, was always lower than spectral coherence γik(t,s). Importantly, during the test, γik/z(t,s) changed regardless of γik(t,s), thus showing that partial coherence reveals new information that is independent from coherence estimates. Orthostatic stress due to position changes caused phase difference and time delay, shown in Figure3g–l, to change. The most pronounced change was that between {−xR(t),xH(t)}. Notably, changes observed in θik(t) were not due to changes in fR(t), since changed as well.
To further portray the comparison between the estimated indices, the median trends are shown in Figure4, while the distributions of the temporal median,, are given in Figure5. Note that in TTSR, median values were estimated in the intervals which are denoted by bold gray lines in Figure4. Numeric values are reported in Table1. Spectral coherence between the three pairs of signals, shown in Figures4a and 5a, followed very similar patterns. After upward and downward movement, from TES to THT and from THT to TLS, coherence abruptly decreased. Previous values were restored in about 1–2 min. During head-up tilt, was statistically higher than during TES and TLS(P < 5·10-4). During head-up tilt, the coupling between −xR(t) and xS(t) was also higher than the coupling between −xR(t) and xH(t), (P < 0.023).
Partial coherence estimates, shown in Figures4b and5b, were always lower than the corresponding coherence estimates. This is not surprising since according to the physiological model, respiratory activity affects both HPV and SAPV. During head-up tilt, was higher than during both TES and TLS(P < 0.043), was lower than during TLS(P < 0.011), and was lower than during TES(P < 0.025). The increase in during head-up tilt indicates that, owing to the position change, the respiratory component in xH(t), was no longer able to explain the coupling between the respiratory components in xR(t) and in xS(t). The fact that was higher during TES indicates that in supine position, the respiratory induced oscillations in SAPV and HPV were not perfectly fitted by −xR(t), and that after removing respiration, the residual coupling represents mutual SAPV-HPV interactions.
The comparison between γik/z(t,s) for different {i,k,z} shows that during (P < 0.011). During THT, (P < 3·10-4), and (P < 0.026). During TLS, no differences were found between partial coherence estimates. The fact that the coupling between {xS(t),xH(t)} after removing −xR(t), γSH/R(t,s), was not lower than both γRS/H(t,s) and γRH/S(t,s) indicates that RSA and the respiratory induced oscillation in SAPV were not merely reflecting −xR(t), and that closed loop interactions between SAPV and HPV had an important role in determining the strength of the coupling between them around the respiratory frequency.
In Figures4c and5c, a comparison between phase difference estimates is shown. Head-up tilt caused to statistically change from both TES and TLS(P < 0.002). Similarly, changes from TES to THT(P < 0.035). During supine position, no differences were found between and, while during head-up tilt the differences between and, for (i,k) ≠ (p,q), were significant. The same changes observed in θik(t) are also observed in.
Comparison between TTSR and TTCR data set
The distribution of the temporal median extracted from TTCR data set is shown in Figure5e–h, while numerical results are reported in Table1. A segment of two representative respiratory signals, xR(t), from TTSR and TTCR data sets, as well as their TF spectra are shown in Figure6a,b. Respiratory signals from TTSR are non-stationary and relatively narrow band, while xR(t) from TTCR are characterized by a wider instantaneous spectra. The distribution of the median respiratory frequency,, reported in Figure6c,d, shows that respiratory frequency had a higher inter-subject variability in TTSR, (respiratory protocol in TTCR was the same for all the participants), while intra-subject variability was higher in TTCR.
Although in the two data set the respiratory pattern differs substantially, results obtained with the same procedure are consistent. Pairwise comparisons between all indices extracted from TTSR and TTCR show that differences were never significant, except for coherence estimates obtained during head-up tilt, that in TTCR were lower than in TTSR (P < 0.001).
Changes during autonomic blockade
The distribution of median indices from the autonomic blockade data set are shown in Figure7, while numerical values are reported in Table2. The comparisons between control and blockade conditions included only those subjects in which a given index was estimated for at least 3 min. As expected, parasympathetic blockade due to atropine injection had a much stronger effect than propranolol sympathetic blockade, especially on the interactions that involve HPV. Atropine had virtually no effect on the interaction between {−xR(t),xS(t)}, while it reduced the strength of coupling between {xS(t),xH(t)} (P < 0.047). Although the coupling between {−xR(t),xH(t)} decreased, this decrease was not statistically significant (P < 0.078). The decrease in partial coherence from 0.623±0.297 to 0.365±0.172 (P < 0.047) indicates that the amount of coupling between {xS(t),xH(t)} that still remained after atropine injection was mainly due to respiration. After atropine injection, increased toward values of. This is likely due to the fact that after parasympathetic blockade, respiratory HPV oscillations decreased; consequently, xH(t) was no longer able to estimate the communality shared by {−xR(t),xS(t)} around respiratory frequency. Due to coherence reduction, and were estimated in only four subjects. In these subjects, parasympathetic blockade caused the latency between respiratory oscillations in {xS(t),xH(t)} to change of about 1.2 s. These changes were not significant, probably due to the small number of subjects involved in the comparison. Interestingly, within this framework, this change may imply a reverse of the causality in HPV-SAPV interactions, that after atropine injection would be characterized by respiratory oscillations in HPV preceding respiratory oscillations in SAPV.
Propranolol caused to slightly decrease (P < 0.078). This may indicate that sympathetic modulation is also involved in the HPV-SAPV interactions. Phase difference also slightly changed (P < 0.078), i.e., after sympathetic blockade phase differences between −xR(t) and xS(t) decreased.
Discussion
The main purpose of this study was to describe a new framework for characterizing the dynamic interactions between non-stationary signals and to demonstrate its usefulness in the analysis of cardiovascular control. The most important contributions of the study are: For the first time, non-parametric time-frequency analysis is used to explore the interactions between the three most relevant cardiovascular signals (heart period, arterial pressure and respiration) simultaneously. Non-parametric time-frequency partial coherence is used, for the first time, in the assessment of the cardiovascular regulation. A new method for reducing the limitation of phase periodicity in the estimation of the latencies between two non-stationary spectral components is proposed. Additionally, in this article we had the unique opportunity to analyze the changes in the interactions between cardiovascular signals after autonomic blockade.
Cross time-frequency analysis
The methodology presented in this article is based on TFD and provides a robust and fast tracking of coherence, partial coherence, phase difference and time delay, by means of a procedure composed of the following steps: (i) TF representations of spectra, coherence, partial coherence and phase differences are estimated through (1)–(9); (ii) A statistical test is performed to localize TF regions, Ωik, where the strength of the local coupling between signals {xi(t),xk(t)} is significant[13]; (iii) A time-varying spectral band, Ω(γ), whose bandwidth is related to the frequency resolution of the TF representation, is defined around the respiratory rate. (iv) Time-varying spectral bands,, which include the portions of Ω(γ) in which coherence estimates are significant are localized. A minimum bandwidth and duration, given by R(t,f), is imposed to by the opening. (v) The time course of coherence and partial coherence is estimated by averaging the corresponding TF representations in Ω(γ), while phase difference changes are estimated by averaging the phase difference spectrum in. (vi) Time delay is estimated to reduce the influence of instantaneous frequency over phase difference estimates.
In this study, the respiratory rate is used to localize Ω(γ) and to convert phase difference into time delay estimates. However, in the case in which the respiratory signal was not available or the interest focused on interactions around another spectral component, the respiratory frequency should be replaced by the central frequency of the instantaneous spectral peak of interest in |Sik(t,f)|.
One of the most important aspects of this procedure is the localization of specific TF regions to estimate any given index. This offers the possibility of solving the problem of estimating the time course of indices that are reliable only around specific time-varying oscillations. Moreover, the statistical assessment of coherence level in the phase difference estimation add further robustness to the estimates.
This TFD has been used because it offers the possibility of independently controlling the time and frequency filtering, and it is characterized by fine resolution. Its capability of accurately estimating the TF structure of physiological and synthetic real-like cardiovascular signals has been demonstrated in different studies[27–31]. In a recent study, TFD has been shown to provide more accurate coherence estimates than multitaper spectrogram and continuous wavelet transform[13]. The estimation of the phase differences in the joint TF domain was used in few studies, where it was performed by wavelet[32, 33], Rihaczek[34], and reduced interference distributions[35]. Among them, no one focuses on the characterization of cardiovascular dynamics. The simulation study performed in this article shows that TFD gives accurate estimates of phase difference and time delay, and confirm previous results[14, 26].
To the extent of our knowledge, this is the first time that TF partial coherence is used. Partial coherence includes information from three signals sharing a similar oscillation[20]. It was shown that in a given triplet, the lowest partial coherence estimate is obtained when the signal which contains the shared oscillation with the better signal to noise ratio is extracted[21]. In other words, the better the extracted signal estimates (explains) the other two, the lower the partial coherence. In this study, partial coherence was the most sensitive index to orthostatic stress, in the sense that in γik/z(t) changes due to head-up tilt were more pronounced than in other indices, at least in spontaneous breathing.
Dependence of phase difference on signal representations
Phase difference estimates depend on how the information carried by a signal is organized on the temporal axis. Thus, the phase of a cross spectrum is sensitive to any transformation that may involve the analyzed signals, such as change of sign, inversion, derivative etc. Continuous cardiovascular signals are often obtained by interpolation of discrete values unevenly distributed in time. The interpolation process directly affects phase difference estimates. Figure8 shows the distribution of and in TTSR data set obtained by using three different representations of HPV. The HPV signal is estimated by interpolating on the first beat (INT1), on the second beat (INT2) or by using the IPFM model as in this study. It is shown that by using INT1 instead of INT2, phase differences are reduced of about 1 rad, which supposing a respiratory rate of 0.2 Hz corresponds to about 800 ms, consistent with an average heart beat. The IPFM model provides estimates in between INT1 and INT2. Importantly, these differences are sufficient to change the sign of θik(t) and consequently the interpretation of the temporal sequence of events in cardiovascular interactions. A reduction of the dependence of phase difference estimates on interpolation may be obtained by modeling the heart period as a point process[36], thus providing a continuous assessment of HPV and cardiovascular interactions[37, 38]. The estimation of SAPV is also critical. For instance, when arterial pressure is estimated in the periphery, its temporal structure also includes the pulse transit time[30]. Taking together, these considerations suggest that caution should be used in interpreting phase difference estimates, especially when they are used to infer causality. However, these issues do not affect the relative changes of the indices over time, which describe how the synchronization of physiological oscillations responds to changes in external conditions.
Limitations
The main limitation of multivariate analysis based on cross-spectral identification is the assumption of an implicit open loop model for the system. This implies that, although causal information may be inferred, to a certain degree, by phase difference estimation[9], our approach is not able to describe feedback and feedforward pathways of closed loops separately, neither to discriminate between direct and indirect interactions[39]. Causal coherence[40–42] or partial directed coherence[43, 44], based on multivariate autoregressive modeling may be used for this purpose. However, multivariate autoregressive analysis is usually applied in stationary conditions and strongly depends on the goodness of fit of the model.
Another limitation is that spectral analysis, contrary to parametric modeling[45] or information based algorithms[46], accounts for linear interactions only. Although cardiovascular variability is known to be affected by non linear phenomena, such as saturation or hysteresis, a linear model characterizes a system by means of indices, such as power, coherence, phase difference, etc. that have a straightforward physical meaning, which in turn can be more useful for the interpretation of the ongoing processes[47].
Finally, cross spectral analysis is not able to assess hemodynamic changes during inspiration and expiration separately, neither to capture morphological information from the analyzed oscillations, thus disregarding potentially valuable information[48, 49].
Physiological results
Major findings were: Intra-subject variability in indices characterizing the dynamic interactions between cardiovascular parameters shows that coupling between physiological rhythms is time-varying and intermittent. In spite of this intra-subject variability, these indices were consistent among subjects, thus suggesting the presence of common underling mechanisms of regulation between HPV, SAPV, and respiration. Systolic arterial pressure changes followed respiratory flow (see Section ‘Interactions between respiration and SAPV’) both in supine and standing positions and even after selective autonomic blockade. During head-up tilt, latencies between inspiratory/expiratory flow and SAPV reduction/increase augmented (see Section ‘Interactions between respiration and SAPV’). Latencies between respiration and HPV were comparable to those between respiration and SAPV during supine position, and significantly increased during standing. As a result, respiratory oscillations in SAPV preceded respiratory oscillations in HPV during standing. In spontaneous breathing, partial coherence was the most sensitive index to orthostatic stress. Phase difference estimates were consistent among spontaneous and controlled breathing patterns, whereas coherence was higher in spontaneous breathing, especially during tilt.
Although in line with other studies[3, 4, 48, 50], caution should be used in drawing conclusions from the results shown in this study, due to the small number of subjects involved in the experiments (14 in TTSR and TTCR, and 7 in ABCR data sets).
Interactions between respiration and SAPV
During inspiration and expiration, intrathoracic pressure changes cause fluctuations in the venous return and cardiac output, which in turn cause respiratory oscillations in arterial pressure[3]. Our results show that xS(t) preceded −xR(t) both in supine and head-up tilt position, both in spontaneous and controlled breathing and even after parasympathetic or sympathetic blockade. This temporal relationship was also documented in other studies[3, 48, 50], and has been explained by suggesting a close mechanical dependence between arterial pressure and the derivative of respiratory volume, i.e., the respiratory flow[3]. Considering the phase shift due to the derivative, −d(xR(t)) / dt, if −xR(t) and xS(t) were synchronous θRS(t) = −Π / 2. Since we found that θRS(t) > −Π / 2, our results suggest that SAPV always follows −d(xR(t)) / dt with a short time delay. This delay increased in standing position[3, 50], and the pattern reported in Figure3c,d show that this change was fast. The local coupling between respirations and SAPV increases after 1–2 min from tilting during spontaneous breathing but decreased during controlled breathing. Interestingly, during head-up tilt, respiratory SAPV better explained the coupling between all the three variables. Autonomic blockade did not affect the interactions between SAPV and respiration significantly, even if it seems that propranolol slightly increased the latencies between them.
Interactions between respiration and HPV
Respiratory sinus arrhythmia, the heart rate variability in synchrony with respiration[2], has both clinical and physiological relevance[1]. Several hypothesis have been made to explain the origin of the respiratory oscillations in HPV. They include a central mechanism[5, 11], the baroreflex[7, 8, 51], the mechanical stretch of the sinus node[3], which is enhanced during exercise[52], and a mixture of them[6, 48, 50]. However, most of the mechanisms underlying RSA are still unclear, and there is no general consensus about its origin[5]. Importantly, the estimation of the time delay between HPV and SAPV around the respiratory frequency is a central point of this debate[5]. Although the discussion of the origin of RSA is beyond the scope of this study, it is important to stress that the presented methodology, which provides robust and reliable phase difference estimates, may be useful to gain some insight into this debate.
Our results show that xH(t) preceded −xR(t), probably due to a similar mechanism as that illustrated in the previous section, i.e., respiratory oscillations in HPV followed changes in the respiratory flow. The fact that changes in HPV anticipated changes in the instantaneous lung volume was also documented in[3, 48, 50]. Change from supine to standing position increased the time delay between respiratory flow and HPV. This increase may be related to an increase of baroreflex contribution to the respiratory oscillations in HPV[3].
Partial coherence analysis shows that during standing, HPV was not able to explain the coupling between respiration and SAPV, likely due to a decrease in the amplitude of respiratory HPV. Parasympathetic blockade reduced the strength of the coupling between HPV and respiration. However, this decrease was not statistically significant and did not reduce coherence estimates below the significance threshold. The amount of HPV not abolished by atropine injection was characterized by a phase difference with respect to −xR(t) close to −Π / 2, which indicates very short time delay with respect to respiratory flow. This supports the hypothesis that non-autonomic, possibly mechanically mediated, mechanisms also contributes to RSA[3].
Interactions between SAPV and HPV
With intact autonomic modulation the removal of respiration influence reduced, but did not abolish, the coupling between respiratory oscillation in SAPV and in HPV. This indicates that some correlation that cannot be perfectly captured by respiration is present in the closed loop between HPV and SAPV around the respiratory frequency. Parasympathetic blockade dampened the local coupling below significance level. Since during parasympathetic blockade vagal baroreflex is abolished, while mechanical influences of respiration on SAPV are intact, the residual interactions between HPV and SAPV are likely due to correlation between residual mechanical oscillations in HPV and respiratory induced oscillations in SAPV. In absence of vagal control, phase differences were consistent to the direction of the feedforward path. During parasympathetic blockade, the amount of residual coupling almost disappeared after removing the influence of respiration. Interestingly, propranolol slightly reduced the strength of the coupling, suggesting a small contribution of sympathetic activity on HPV-SAPV interactions around the respiratory frequency.
Conclusions
In this paper, a new framework for characterizing the dynamic interactions between non-stationary signals was presented and its usefulness in the analysis of cardiovascular control was demonstrated. The most important contributions of this methodology are: For the first time, non-parametric time-frequency analysis is used to explore the interactions between the three most relevant cardiovascular signals (heart period, arterial pressure and respiration) simultaneously. Non-parametric time-frequency partial coherence is used, for the first time, in the assessment of the cardiovascular regulation. A new method for reducing the limitation of phase periodicity in the estimation of the latencies between two non-stationary spectral components is proposed. Additionally, in this article we had the unique opportunity to analyze the changes in the interactions between cardiovascular signals after autonomic blockade.
References
Grossman P, Taylor EW: Toward understanding respiratory sinus arrhythmia: relations to cardiac vagal tone, evolution and biobehavioral functions. Biol. Psychol 2007, 74(2):263-285. 10.1016/j.biopsycho.2005.11.014
Yasuma F, Hayano J: Respiratory sinus arrhythmia∗: why does the heartbeat synchronize with respiratory rhythm? Chest 2004, 125(2):683-690. 10.1378/chest.125.2.683
Saul JP, Berger RD, Albrecht P, Stein SP, Chen MH, Cohen RJ: Transfer function analysis of the circulation: unique insights into cardiovascular regulation. Am. J. Physiol. Heart Circ. Physiol 1991, 261(4):H1231-H1245.
Barbieri R, Bianchi A, Triedman J, Mainardi L, Cerutti S, Saul J: Model dependency of multivariate autoregressive spectral analysis. IEEE Eng. Med. Biol. Mag 1997, 16(5):74-85. 10.1109/51.620498
Eckberg DL: Point:counterpoint: respiratory sinus arrhythmia is due to a central mechanism vs. respiratory sinus arrhythmia is due to the baroreflex mechanism. J. Appl. Physiol 2009, 106(5):1740-1742. discussion 1744 10.1152/japplphysiol.91107.2008
Julien C, Parkes MJ, Tzeng SYC, Sin PYW, Ainslie PN, van de Borne P, Fortrat JO, Custaud MA, Gharib C, Porta A, Vallais F, Baselli G, Pagani M, Lucini D, Hughson RL, Taylor JA, Tan CO, Baekey DM, Dick TE, Paton JFR, Taha B: Comments on point:counterpoint: respiratory sinus arrhythmia is due to a central mechanism vs. respiratory sinus arrhythmia is due to the baroreflex mechanism. J. Appl. Physiol 2009, 106(5):1745-1749. 10.1152/japplphysiol.00196.2009
Karemaker JM: Last word on point:counterpoint: respiratory sinus arrhythmia is due to a central mechanism vs. respiratory sinus arrhythmia is due to the baroreflex mechanism. J. Appl. Physiol 2009, 106(5):1750. 10.1152/japplphysiol.00225.2009
deBoer RW, Karemaker JM, Strackee J: Hemodynamic fluctuations and baroreflex sensitivity in humans: a beat-to-beat model. Am. J. Physiol 1987, 253(3 Pt 2):H680-H689.
Cooke WH, Hoag JB, Crossman AA, Kuusela TA, Tahvanainen KUO, Eckberg DL: Human responses to upright tilt: a window on central autonomic integration. J. Physiol 1999, 517(2):617-628. 10.1111/j.1469-7793.1999.0617t.x
Porta A, Catai AM, Takahashi ACM, Magagnin V, Bassani T, Tobaldini E, van de Borne P, Montano N: Causal relationships between heart period and systolic arterial pressure during graded head-up tilt. Am. J. Physiol. Regul. Integr. Comp. Physiol 2011, 300(2):R378-R386. 10.1152/ajpregu.00553.2010
Badra LJ, Cooke WH, Hoag JB, Crossman AA, Kuusela TA, Tahvanainen KUO, Eckberg DL: Respiratory modulation of human autonomic rhythms. Am. J. Physiol. Heart Circ. Physiol 2001, 280(6):H2674-H2688.
Larsen PD, Lewis CD, Gebber GL, Zhong S: Partial spectral analysis of cardiac-related sympathetic nerve discharge. J. Neurophysiol 2000, 84(3):1168-1179.
Orini M, Bailón R, Mainardi L, Laguna P, Flandrin P: Characterization of the dynamic interactions between cardiovascular signals by time-frequency coherence. IEEE Trans. Biomed. Eng 2012, 59(3):663-673.
Orini M, Laguna P, Mainardi LT, Bailón R: Assessment of the dynamic interactions between heart rate and arterial pressure by the cross time-frequency analysis. Physiol. Meas 2012, 33(3):315-331. 10.1088/0967-3334/33/3/315
Mincholé A, Pueyo E, Rodríguez JF, Zacur E, Doblaré M, Laguna P: Quantification of restitution dispersion from the dynamic changes of the T-wave peak to end, measured at the surface ECG. IEEE Trans. Biomed. Eng 2011, 58(5):1172-1182.
Martinez JP, Almeida R, Olmos S, Rocha AP, Laguna P: A wavelet-based ECG delineator: Evaluation on standard databases. IEEE Trans. Biomed. Eng 2004, 51(4):570-581. 10.1109/TBME.2003.821031
Mateo J, Laguna P: Analysis of heart rate variability in the presence of ectopic beats using the heart timing signal. IEEE Trans. Biomed. Eng 2003, 50: 334-343. 10.1109/TBME.2003.808831
Hlawatsch F, Boudreaux-Bartels GF: Linear and quadratic time-frequency signal representations. IEEE Signal Process. Mag 1992, 9(2):21-67.
Flandrin P: Time-Frequency Signal Analysis and Processing, chap. Ambiguity Function. Philadelphia, Elsevier; 2003.
Bendat J, Piersol A: Random Data: Analysis and Measurement Procedures. New York, John Wiley & Sons; 2010. Wiley Series in Probability and Statistics
Baccalá LA, Sameshima K: Comments on ‘Is partial coherence a viable technique for identifying generators of neural oscillations?’: Why the term ‘Gersch Causality’ is inappropriate: common neural structure inference pitfalls. Biol. Cybern 2006, 95(2):135-141. 10.1007/s00422-006-0075-7
Matz G, Hlawatsch F: Time-frequency coherence analysis of nonstationary random processes. In Proc. Tenth IEEE Workshop Statistical Signal and Array Processing. Pennsylvania, PA, USA; 2000. pp. 554–558
Costa A, Boudreau-Bartels G: Design of time-frequency representations using a multiform, tiltable exponential kernel. IEEE Trans. Signal Process 1995, 43(10):2283-2301. 10.1109/78.469860
Hlawatsch F, Flandrin P: The Wigner Distribution - Theory and Applications in Signal Processing. Philadelphia, Elsevier; 1997. pp. 59–113. Elsevier 1997 chap. The interference structure of the Wigner distribution and related time-frequency signal representations
Orini M, Bailón R, Mainardi L, Mincholé A, Laguna P: Continuous quantification of spectral Coherence using Quadratic Time-Frequency distributions, error analysis and application. In Internat. Conf. Computers in Cardiology. Park City, UT, USA; 2009. pp. 681–684
Orini M, Bailón R, Mainardi LT, Laguna P: Time-frequency phase differences and phase locking to characterize dynamic interactions between cardiovascular signals. In Conf. Proc. IEEE Eng. Med. Biol. Soc.. vol 2011, Boston, MA USA; 2011. pp. 4689–4692
Pola S, Macerata A, Emdin M, Marchesi C: Estimation of the power spectral density in nonstationary cardiovascular time series: assessing the role of the time-frequency representations (TFR). IEEE Trans. Biomed. Eng 1996, 43: 46-59.
Orini M, Bailón R, Enk R, Koelsch S, Mainardi LT, Laguna P: A method for continuously assessing the autonomic response to music-induced emotions through HRV analysis. Med. Biol. Eng. Comput 2010, 48(5):423-433. 10.1007/s11517-010-0592-3
Bailón R, Mainardi L, Orini M, Land Laguna SP: Analysis of heart rate variability during exercise stress testing using respiratory information. Biomed. Signal Process Control 2010, 5: 299-310. 10.1016/j.bspc.2010.05.005
Gil E, Orini M, Bailón R, Vergara JM, Mainardi L, Laguna P: Photoplethysmography pulse rate variability as a surrogate measurement of heart rate variability during non-stationary conditions. Physiol. Meas 2010, 31(9):1271. 10.1088/0967-3334/31/9/015
Orini M, Bailón R, Mainardi L, Laguna P: Synthesis of HRV signals characterized by predetermined time-frequency structure by means of time-varying ARMA models. Biomed. Signal Process Control 2012, 7: 141-150. 10.1016/j.bspc.2011.05.003
Lachaux JP, Rodriguez E, Martinerie J, Varela FJ: Measuring phase synchrony in brain signals. Hum Brain Mapp 1999, 8(4):194-208. 10.1002/(SICI)1097-0193(1999)8:4<194::AID-HBM4>3.0.CO;2-C
Quyen MLV, Foucher J, Lachaux J, Rodriguez E, Lutz A, Martinerie J, Varela FJ: Comparison of Hilbert transform and wavelet methods for the analysis of neuronal synchrony. J. Neurosci. Methods 2001, 111(2):83-98. 10.1016/S0165-0270(01)00372-7
Aviyente S, Bernat EM, Evans WS, Sponheim SR: A phase synchrony measure for quantifying dynamic functional integration in the brain. Hum Brain Mapp 2010, 32(1):80-93.
Shin YJ, Gobert D, Sung SH, Powers EJ, Park JB: Application of cross time-frequency analysis to postural sway behavior: the effects of aging and visual systems. IEEE Trans. Biomed. Eng 2005, 52(5):859-868. 10.1109/TBME.2005.844278
Barbieri R, Matten EC, Alabi AA, Brown EN: A point-process model of human heartbeat intervals: new definitions of heart rate and heart rate variability. Am. J. Physiol. Heart Circ. Physiol 2005, 288: H424-H435.
Chen Z, Brown EN, Barbieri R: Assessment of autonomic control and respiratory sinus arrhythmia using point process models of human heart beat dynamics. IEEE Trans. Biomed. Eng 2009, 56(7):1791-1802.
Chen Z, Purdon P, Harrell G, Pierce E, Walsh J, Brown E, Barbieri R: Dynamic assessment of baroreflex control of heart rate during induction of propofol anesthesia using a point process method. Annals Biomed. Eng 2011, 39: 260-276. 10.1007/s10439-010-0179-z
Barbieri R, Parati G, Saul J: Closed- versus open-loop assessment of heart rate baroreflex. IEEE Eng. Med. Biol. Mag 2001, 20(2):33-42. 10.1109/51.917722
Zhao H, Cupples WA, Ju KH, Chon KH: Time-varying causal coherence function and its application to renal blood pressure and blood flow data. IEEE Trans. Biomed. Eng 2007, 54(12):2142-2150.
Porta A, Furlan R, Rimoldi O, Pagani M, Malliani A, van de Borne P: Quantifying the strength of the linear causal coupling in closed loop interacting cardiovascular variability signals. Biol. Cybern 2002, 86(3):241-251. 10.1007/s00422-001-0292-z
Nollo G, Faes L, Porta A, Antolini R, Ravelli F: Exploring directionality in spontaneous heart period and systolic pressure variability interactions in humans: implications in the evaluation of baroreflex gain. Am. J. Physiol. Heart Circ. Physiol 2005, 288(4):H1777-H1785.
Baccalá LA, Sameshima K: Partial directed coherence: a new concept in neural structure determination. Biol. Cybern 2001, 84: 463-474. 10.1007/PL00007990
Faes L, Porta A, Nollo G: Testing frequency-domain causality in multivariate time series. IEEE Trans. Biomed. Eng 2010, 57(8):1897-1906.
Chen Z, Brown EN, Barbieri R: Characterizing nonlinear heartbeat dynamics within a point process framework. IEEE Trans. Biomed. Eng 2010, 57(6):1335-1347.
Faes L, Nollo G, Porta A: Information domain approach to the investigation of cardio-vascular, cardio-pulmonary and vasculo-pulmonary causal couplings. Frontiers Physiol 2011, 2(80):1-13.
Xiao X, Mullen TJ, Mukkamala R: System identification: a multi-signal approach for probing neural cardiovascular regulation. Physiol. Meas 2005, 26(3):R41-71. 10.1088/0967-3334/26/3/R01
Sin PYW, Galletly DC, Tzeng YC: Influence of breathing frequency on the pattern of respiratory sinus arrhythmia and blood pressure: old questions revisited. Am. J. Physiol. Heart Circ. Physiol 2010, 298(5):H1588-H1599. 10.1152/ajpheart.00036.2010
Tzeng YC, Sin PYW, Galletly DC: Human sinus arrhythmia: inconsistencies of a teleological hypothesis. Am. J. Physiol. Heart Circ. Physiol 2009, 296: H65-H70.
Elstad M, Toska K, Chon KH, Raeder EA, Cohen RJ: Respiratory sinus arrhythmia: opposite effects on systolic and mean arterial pressure in supine humans. J. Physiol 2001, 536(Pt 1):251-259.
Karemaker JM: Counterpoint: respiratory sinus arrhythmia is due to the baroreflex mechanism. J. Appl. Physiol 2009, 106(5):1742-1743. discussion 1744 10.1152/japplphysiol.91107.2008a
Blain G, Meste O, Bermon S: Influences of breathing patterns on respiratory sinus arrhythmia in humans during exercise. Am. J. Physiol. Heart Circ. Physiol 2005, 288(2):H887-95.
Acknowledgements
Thanks to Dr. JP Saul and Dr. A Mincholé for providing the authors with TTSR, TTCB, and ABCR data sets. This study was supported in part by NIH Grants R01-HL084502, Ministerio de Ciencia e Innovación, Spain, under Projects TEC2010-21703-C03-02 and TRA2009-0127, Diputación General de Aragón, Spain, through Grupos Consolidados GTC ref:T30, Instituto de Salud Carlos III, Spain, through CIBER CB06/01/0062, and ARAID and Ibercaja under project “Programa de apoyo a la I+D+i”.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Orini, M., Bailón, R., Laguna, P. et al. A multivariate time-frequency method to characterize the influence of respiration over heart period and arterial pressure. EURASIP J. Adv. Signal Process. 2012, 214 (2012). https://doi.org/10.1186/1687-6180-2012-214
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1687-6180-2012-214