×

通过状态空间建模分解神经多元时间序列。 (英语) Zbl 1221.92020年

摘要:将多元时间序列数据分解为独立的源成分是神经科学中时间分辨数据预处理和分析的重要组成部分。我们简要回顾了用于此目的的可用工具,如因子分析(FA)和独立成分分析(ICA),然后我们展示了如何将线性状态空间建模(一种来自统计时间序列分析的方法)用于相同目的。状态空间建模是经典ARMA建模的推广,非常适合于利用时间序列数据时序中编码的动态信息,而FA和大多数ICA算法仍然无法获取这些信息。因此,可以进行更详细的分解,即使在没有先验信息的情况下,也可以分离出功率谱尖锐的两个分量,例如α分量、正弦伪影或睡眠纺锤波,以及功率谱较宽的两个成分,例如FMRI扫描仪伪影或癫痫尖峰分量。
此外,还讨论了三种推广,第一种是放松独立性假设,第二种是引入驱动动力学的噪声协方差的非平稳性,第三种是通过非线性观测函数允许数据的非高斯性。给出了三个应用实例,一个心电图时间序列和两个脑电图(EEG)时间序列。两个均来自癫痫患者的EEG示例显示了各种伪影的分离和去除,包括嗡嗡声和FMRI扫描仪伪影,以及睡眠纺锤波、癫痫灶和尖峰成分的识别。为了进行比较,显示了两种ICA算法获得的分解。

MSC公司:

92C20美元 神经生物学
62M10个 统计学中的时间序列、自相关、回归等(GARCH)
62小时25分 因子分析和主成分;对应分析
62页第10页 统计学在生物学和医学中的应用;元分析
92 C55 生物医学成像和信号处理
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] 阿伊特·萨哈利亚,Y.,&;Kimmel,R.(2007)。随机波动率模型的最大似然估计。J.财务。经济。,83, 413–452. ·doi:10.1016/j.fineco.2005.10.006
[2] Akaike,H.(1974年a)。随机过程的马尔可夫表示及其在自回归滑动平均过程分析中的应用。安.Inst.Stat.Math。,26, 363–387. ·Zbl 0335.62058号 ·doi:10.1007/BF02479833
[3] Akaike,H.(1974年b)。统计模型识别的新视角。IEEE传输。自动。控制,19,716–723·Zbl 0314.62039号 ·doi:10.1109/TAC.1974.1100705
[4] Akaike,H.和;Nakagawa,T.(1988)。动态系统的统计分析和控制。多德雷赫特:克鲁沃学院·Zbl 0691.93022号
[5] Allen,P.J.、Josephs,O.和;Turner,R.(2000)。一种从功能MRI期间记录的连续EEG中去除图像伪影的方法。神经影像,12,230–239·doi:10.1006/nimg.2000.0599
[6] 奥斯特罗姆,K.J.(1980)。最大似然法和预测误差法。Automatica,16551–574·Zbl 0441.93036号 ·doi:10.1016/0005-1098(80)90078-3
[7] 阿提亚斯·H·;Schreiner,C.E.(1998)。盲源分离与反褶积:动态分量分析算法。神经计算。,10, 1373–1424. ·doi:10.1162/neco.1998.10.6.1373
[8] Baldick,R.(编辑)(2006年)。应用优化:工程系统的公式和算法。剑桥:剑桥大学出版社·兹比尔1127.65032
[9] Bar-Shalom,Y.,&;Fortmann,T.(1988)。跟踪和数据关联。圣地亚哥:学术出版社·Zbl 0634.93001号
[10] Barros,A.K.,&;Cichocki,A.(2001年)。提取具有时间结构的特定信号。神经计算。,1995年至2000年·Zbl 1003.94505号 ·doi:10.1162/089976601750399272
[11] 巴兹列夫斯基(1994)。统计因素分析及相关方法:理论与应用。纽约:Wiley-Interscience·兹比尔1130.62341
[12] Beckmann,C.和;Smith,S.(2004)。功能磁共振成像的概率独立成分分析。IEEE传输。医学成像,23137-152·doi:10.1109/TMI.2003.822821
[13] Beckmann,C.和;Smith,S.(2005)。多主体FMRI分析中独立成分分析的张量扩展。《神经影像》,25,294–311·doi:10.1016/j.neuroimage.2004.10.043
[14] Belouchrani,A.、Abed-Meraim,K.、Cardoso,J.-F.和;Moulines,E.(1997年)。一种使用二阶统计量的盲源分离技术。IEEE传输。信号处理。,45, 434–444. ·数字对象标识代码:10.1109/78.554307
[15] Bollerslev,T.(1986)。广义自回归条件异方差。《经济学杂志》。,31, 307–327. ·Zbl 0616.62119号 ·doi:10.1016/0304-4076(86)90063-1
[16] 盒子、G.E.P.和;Jenkins,G.M.(1970年)。时间序列分析、预测和控制。旧金山:霍尔顿日·Zbl 0249.62009号
[17] Brockwell,P.J.和;Davis,R.A.(1987年)。时间序列:理论和方法。柏林:斯普林格·Zbl 0604.62083号
[18] Cheung,Y.M.和;Xu,L.(2003)。用于时间信号分离的状态空间双多元自回归建模。IEEE传输。系统。人类网络。,33, 386–398.
[19] Choi,S.、Cichocki,A.、Park,H.和;Lee,S.(2005)。盲源分离和独立成分分析:综述。神经信息处理。莱特。第6版,第1-57页·Zbl 1119.81021号 ·doi:10.1007/s11128-006-0037-y
[20] Chui,C.K.和;Chen,G.(1999年)。信息科学史普林格系列:第17卷。卡尔曼滤波:实时应用(第三版)。柏林:斯普林格。
[21] Cichocki,A.和;Amari,S.(2002年)。自适应盲信号和图像处理。奇切斯特:威利·Zbl 0999.93013号
[22] Comon,P.(1994年)。独立成分分析,一个新概念?信号处理。,36, 287–314. ·Zbl 0791.62004号 ·doi:10.1016/0165-1684(94)90029-9
[23] Delorme,A.、Sejnowski,T.和;Makeig,S.(2007)。使用高阶统计和独立分量分析增强EEG数据中伪影的检测。神经影像,341443-1449·doi:10.1016/j.neuroimage.2006.11.004
[24] Durbin,J.和;Koopman,S.J.(2001)。用状态空间方法进行时间序列分析。牛津:牛津大学出版社·Zbl 0995.62504号
[25] Dyrholm,M.、Makeig,S.和;Hansen,L.K.(2007)。卷积ICA的模型选择及其在脑电时空分析中的应用。神经计算。,19, 934–955. ·Zbl 1210.94033号 ·doi:10.1162/neco.2007.19.4.934
[26] 恩格尔、R.F.和;Watson,M.(1981)。大都市工资率的单因素多元时间序列模型。美国统计协会杂志,76,774–781·doi:10.1080/01621459.1981.1047720
[27] 加尔卡,A.,山下,O.,&Ozaki,T.(2004)。逆解动态估计中协方差的GARCH建模。物理学。莱特。A、 333261-268·Zbl 1123.62317号 ·doi:10.1016/j.physleta.2004.10.045
[28] Galka,A.、Ozaki,T.、Bosch Bayard,J.和;Yamashita,O.(2006年)。白化作为估计时空数据集中相互信息的工具。《统计物理学杂志》。,124, 1275–1315. ·Zbl 1108.62090号 ·doi:10.1007/s10955-006-9131-x
[29] Galka,A.、Wong,K.和;Ozaki,T.(2010)。非平稳EEG时间序列建模的广义状态空间模型。在A.Steyn-Ross中;M.Steyn-Ross(编辑),计算神经科学中的Springer系列。建模大脑中的相变(第27-52页)。柏林:斯普林格。
[30] Gevers,M.(2006)。系统识别发展的个人观点。IEEE控制系统。杂志,26,93–105·doi:10.1109/MCS.2006.252834
[31] Gnedenko,B.V.(1969年)。概率论。莫斯科:和平号出版社·Zbl 0191.46702号
[32] Grewal,M.S.和;Andrews,A.P.(2001)。卡尔曼滤波:使用MATLAB的理论和实践。纽约:Wiley-Interscience·Zbl 1322.93001号
[33] Harman,H.H.(1976)。现代因素分析(第三版)。芝加哥:芝加哥大学出版社·Zbl 0161.39805号
[34] Harvey,A.、Koopman,S.J.和;Shephard,N.(编辑)(2004年)。状态空间和未观察到的组件模型。剑桥:剑桥大学出版社·Zbl 1053.62006年
[35] Hyvärinen,A.(1999年)。独立分量分析的快速稳健定点算法。IEEE传输。神经网络。,10, 626–634. ·数字对象标识代码:10.1109/72.761722
[36] Hyvärinen,A.、Karhunen,J.和;Oja,E.(2001)。独立成分分析。纽约:Wiley。
[37] C·詹姆斯;Hesse,C.(2005)。生物医学信号的独立成分分析。生理学。测量。,26,R15–R39·doi:10.1088/0967-3334/26/1/R02
[38] Jung,A.和;Kaiser,A.(2003)。在独立成分分析中考虑时间结构。In:程序。第四个内部症状。ICA BSS,ICA 2003(第95-100页)。日本奈良,2003年4月。
[39] Jung,T.-P.,Makeig,S.,McKeown,M.,Bell,A.,Lee,T.-W.,&Sejnowski,T.(2001)。使用独立成分分析成像大脑动力学。IEEE程序。,88, 1107–1122. ·数字对象标识代码:10.1109/5.939827
[40] Kailath,T.(1968年)。最小二乘估计的创新方法——第一部分:加性白噪声中的线性滤波。IEEE传输。自动。对照组,13646–655·doi:10.1109/TAC.1968.1099025
[41] Kailath,T.(1980)。信息与系统科学系列。线性系统。恩格尔伍德悬崖:普伦蒂斯·霍尔·Zbl 0454.93001号
[42] Kallenberg,O.(2002)。现代概率的基础。柏林:斯普林格·Zbl 0996.60001号
[43] Kalman,R.E.(1960年)。一种解决线性滤波和预测问题的新方法。《基础工程杂志》,82,35–45·数字对象标识代码:10.1115/1.3662552
[44] Kalman,R.E.、Falb,P.L.和;Arbib,M.A.(1969年)。纯数学和应用数学国际系列。数学系统理论的主题。纽约:McGraw-Hill·Zbl 0231.49001号
[45] Ljung,L.(1999)。系统识别:用户理论(第二版)。恩格尔伍德悬崖:普伦蒂斯·霍尔。
[46] Mehra,R.K.(1971)。使用卡尔曼滤波器表示法识别随机线性系统。AIAA J.,第9、28–31页·Zbl 0249.93047号 ·doi:10.2514/3.6120
[47] Mehra,R.K.(1974年)。控制和计量经济学中的识别:异同。《经济学年鉴》。Soc.Meas.公司。,3, 21–47.
[48] Meinecke,F.、Ziehe,A.、Kawanabe,M.和;Müller,K.-R.(2002)。一种重采样方法,用于估计单维或多维独立分量的稳定性。IEEE传输。生物识别。工程师,49,1514-1525·doi:10.1109/TBME.2002.805480
[49] Miwakeichi,F.、Martínez-Montes,E.、Valdés-Sosa,P.、Nishiyama,N.、Mizuhara,H.和;Yamaguchi,Y.(2004)。使用并行因子分析将EEG数据分解为空间-时间-频率分量。神经影像,221035-1045·doi:10.1016/j.neuroimage.2004.03.039
[50] Molenaar,P.C.(1985)。多元时间序列分析的动态因子模型。《心理测量学》,50181-202·Zbl 0603.62099号 ·doi:10.1007/BF02294246
[51] Molgedey,L.和;舒斯特,H.G.(1994)。使用延时相关分离独立信号的混合物。物理学。修订稿。,72, 3634–3637. ·doi:10.1103/PhysRevLett.72.3634
[52] 根岸,M.、阿比尔德加德,M.,尼克松,T;Constable,R.(2004)。从连续功能磁共振成像期间获得的脑电图数据中去除时变梯度伪影。临床。神经生理学。,115, 2181–2192. ·doi:10.1016/j.clinph.2004.04.005
[53] Neumaier,A.和;Schneider,T.(2001)。多元自回归模型的参数和特征模式的估计。ACM事务处理。数学。软质。,27, 27–57. ·Zbl 1070.65503号 ·doi:10.1145/382043.382304
[54] Niazy,R.、Beckmann,C.、Iannetti,D.、Brady,J.和;Smith,S.(2005)。使用最优基集从EEG数据中去除FMRI环境伪影。《神经影像》,28720-737·doi:10.1016/j.neuroimage.2005.06.067
[55] Otter,P.(1986)。间接观测下的动态结构系统:从系统论角度看的可识别性和估计方面。《心理测量学》,51,415–428·Zbl 0621.62107号 ·doi:10.1007/BF02294064
[56] Ozaki,T.和;Iino,M.(2001)。非高斯时间序列分析的创新方法。J.应用。可能性。,第38页,78–92页·Zbl 1008.62086号 ·doi:10.1239/jap/1085496593
[57] Pagan,A.R.(1975年)。关于从时间序列中提取成分的说明。《计量经济学》,43,163–168·Zbl 0291.62124号 ·doi:10.2307/1913421
[58] Pearlmutter,B.A.,&;Parra,L.C.(1997年)。最大似然盲源分离:ICA的上下文敏感推广。M.C.Mozer、M.I.Jordan和;T.Petsche(编辑),《神经信息处理系统的进展》(第9卷,第613-619页)。剑桥:麻省理工学院出版社。
[59] 普洛特(1990)。随机积分和微分方程。柏林:斯普林格·Zbl 0694.60047号
[60] 劳赫·H·E、董·G和;Striebel,C.T.(1965)。线性动力系统的最大似然估计。AIAA期刊,31445–1450·数字对象标识代码:10.2514/3.3166
[61] Schwarz,G.(1978年)。估算模型的维度。Ann.Stat.,第6461–464页·Zbl 0379.62005年 ·doi:10.1214/aos/1176344136
[62] Schweppe,F.(1965年)。高斯信号的似然函数评估。IEEE传输。信息理论,11,61-70·Zbl 0127.10805号 ·doi:10.1109/TIT.1965.1053737号文件
[63] Sorenson,H.W.(1970年)。最小二乘估计:从高斯到卡尔曼。IEEE规范。,7, 63–68. ·doi:10.1109/MSPEC.1970.5213471
[64] Stögbauer,H.、Kraskov,A.、Astakhov,S.A.和;Grassberger,P.(2004)。基于相互信息的租赁相关组件分析。物理学。版本E,70066123·doi:10.1103/PhysRevE.70.066123
[65] Tong,L.,Liu,R.,Soon,V.C.和;Huang,Y.(1991)。盲分离的不确定性和可识别性。IEEE传输。电路系统。,38, 499–509. ·Zbl 0737.93016号 ·数字对象标识代码:10.1109/31.76486
[66] Vigário,R.、Sarela,J.、Jousmiki,V.、Hamalainen,M.和;Oja,E.(2000年)。EEG和MEG记录分析的独立分量方法。IEEE传输。生物识别。工程师,47、589–593·doi:10.10109/10.841330
[67] Waheed,K.,&;Salem,F.M.(2005)。动态环境中盲源恢复的线性状态空间前馈和反馈结构。神经过程。莱特。,22, 325–344. ·doi:10.1007/s11063-005-1000-0
[68] Wong,K.F.K.,Galka,A.,Yamashita,O.,&Ozaki,T.(2006)。利用状态空间GARCH模型对脑电时间序列中的非平稳方差进行建模。计算。生物医学,361327-1335·doi:10.1016/j.compbiomed.2005.10.001
[69] Zhang,L.和;Cichocki,A.(2000年)。动力学系统的盲反褶积:状态空间方法。J.信号处理。,4, 111–130.
[70] Ziehe,A.和;Müller,K.-R.(1998)。TDSEP–一种利用时间结构进行盲分离的有效算法。L.Niklasson、M.Bodén和;T.Ziemke(编辑),Proc。第8期国际比较人工神经网络,ICANN’98(第675-680页)。柏林:斯普林格。
此参考列表基于出版商或数字数学图书馆提供的信息。它的项目与zbMATH标识符启发式匹配,并且可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。