Reduced-order modeling and analysis of dynamic cerebral autoregulation via diffusion maps

Objective.A data-driven technique for parsimonious modeling and analysis of dynamic cerebral autoregulation (DCA) is developed based on the concept of diffusion maps. Specifically, first, a state-space description of DCA dynamics is considered based on arterial blood pressure, cerebral blood flow velocity, and their time derivatives. Next, an eigenvalue analysis of the Markov matrix of a random walk on a graph over the dataset domain yields a low-dimensional representation of the intrinsic dynamics. Further dimension reduction is made possible by accounting only for the two most significant eigenvalues. The value of their ratio indicates whether the underlying system is governed by active or hypoactive dynamics, indicating healthy or impaired DCA function, respectively. We assessed the reliability of the technique by considering healthy individuals and patients with unilateral internal carotid artery (ICA) stenosis or occlusion. We computed the sensitivity of the technique to detect the presumed side-to-side difference in the DCA function of the second group (assuming hypoactive dynamics on the occluded or stenotic side), using McNemar's chi square test. The results were compared with transfer function analysis (TFA). The performance of the two methods was also compared under the assumption of missing data.Main results.Both diffusion maps and TFA suggested a physiological side-to-side difference in the DCA of ICA stenosis or occlusion patients with a sensitivity of 81% and 71%, respectively. Further, both two methods suggested the difference between the occluded or stenotic side and any two sides of the healthy group. However, the diffusion maps captured additional difference between the unoccluded side and the healthy group, that TFA did not. Furthermore, compared to TFA, diffusion maps exhibited superior performance when subject to missing data.Significance.The eigenvalues ratio derived using the diffusion maps technique can be potentially used as a reliable and robust biomarker for assessing how active the intrinsic dynamics of the autoregulation is and for indicating healthy versus impaired DCA function.