Xu Haojie, Lu Yunfeng, Zhu Shanan, He Bin
IEEE Trans Biomed Eng. 2014 Jul;61(7):1979-88. doi: 10.1109/TBME.2014.2311034.
It is of significance to assess the dynamic spectral causality among physiological signals. Several practical estimators adapted from spectral Granger causality have been exploited to track dynamic causality based on the framework of time-varying multivariate autoregressive (tvMVAR) models. The nonzero covariance of the model's residuals has been used to describe the instantaneous effect phenomenon in some causality estimators. However, for the situations with Gaussian residuals in some autoregressive models, it is challenging to distinguish the directed instantaneous causality if the sufficient prior information about the "causal ordering" is missing. Here, we propose a new algorithm to assess the time-varying causal ordering of tvMVAR model under the assumption that the signals follow the same acyclic causal ordering for all time lags and to estimate the instantaneous effect factor (IEF) value in order to track the dynamic directed instantaneous connectivity. The time-lagged adaptive directed transfer function (ADTF) is also estimated to assess the lagged causality after removing the instantaneous effect. In this study, we first investigated the performance of the causal-ordering estimation algorithm and the accuracy of IEF value. Then, we presented the results of IEF and time-lagged ADTF method by comparing with the conventional ADTF method through simulations of various propagation models. Statistical analysis results suggest that the new algorithm could accurately estimate the causal ordering and give a good estimation of the IEF values in the Gaussian residual conditions. Meanwhile, the time-lagged ADTF approach is also more accurate in estimating the time-lagged dynamic interactions in a complex nervous system after extracting the instantaneous effect. In addition to the simulation studies, we applied the proposed method to estimate the dynamic spectral causality on real visual evoked potential (VEP) data in a human subject. Its usefulness in time-variant spectral causality assessment was demonstrated through the mutual causality investigation of brain activity during the VEP experiments.
评估生理信号之间的动态频谱因果关系具有重要意义。基于时变多元自回归(tvMVAR)模型框架,已经开发了几种从频谱格兰杰因果关系改编而来的实用估计器来跟踪动态因果关系。在一些因果关系估计器中,模型残差的非零协方差已被用于描述瞬时效应现象。然而,对于某些自回归模型中具有高斯残差的情况,如果缺少关于“因果顺序”的足够先验信息,区分有向瞬时因果关系具有挑战性。在此,我们提出一种新算法,在信号在所有时间滞后都遵循相同无环因果顺序的假设下,评估tvMVAR模型的时变因果顺序,并估计瞬时效应因子(IEF)值,以跟踪动态有向瞬时连通性。还估计了时滞自适应定向传递函数(ADTF),以在去除瞬时效应后评估滞后因果关系。在本研究中,我们首先研究了因果顺序估计算法的性能和IEF值的准确性。然后,通过各种传播模型的模拟,将IEF和时滞ADTF方法的结果与传统ADTF方法进行比较。统计分析结果表明,新算法能够在高斯残差条件下准确估计因果顺序并对IEF值进行良好估计。同时,在提取瞬时效应后,时滞ADTF方法在估计复杂神经系统中的时滞动态相互作用方面也更准确。除了模拟研究,我们还将所提出的方法应用于估计人类受试者真实视觉诱发电位(VEP)数据上的动态频谱因果关系。通过VEP实验期间大脑活动的相互因果关系研究,证明了其在时变频谱因果关系评估中的有用性。