title: STATE-SPACE SOLUTIONS TO THE DYNAMIC MAGNETOENCEPHALOGRAPHY INVERSE PROBLEM USING HIGH PERFORMANCE COMPUTING.
abstract: Determining the magnitude and location of neural sources within the brain that are responsible for generating magnetoencephalography (MEG) signals measured on the surface of the head is a challenging problem in functional neuroimaging. The number of potential sources within the brain exceeds by an order of magnitude the number of recording sites. As a consequence, the estimates for the magnitude and location of the neural sources will be ill-conditioned because of the underdetermined nature of the problem. One well-known technique designed to address this imbalance is the minimum norm estimator (MNE). This approach imposes an L(2) regularization constraint that serves to stabilize and condition the source parameter estimates. However, these classes of regularizer are static in time and do not consider the temporal constraints inherent to the biophysics of the MEG experiment. In this paper we propose a dynamic state-space model that accounts for both spatial and temporal correlations within and across candidate intra-cortical sources. In our model, the observation model is derived from the steady-state solution to Maxwell's equations while the latent model representing neural dynamics is given by a random walk process. We show that the Kalman filter (KF) and the Kalman smoother [also known as the fixed-interval smoother (FIS)] may be used to solve the ensuing high-dimensional state-estimation problem. Using a well-known relationship between Bayesian estimation and Kalman filtering, we show that the MNE estimates carry a significant zero bias. Calculating these high-dimensional state estimates is a computationally challenging task that requires High Performance Computing (HPC) resources. To this end, we employ the NSF Teragrid Supercomputing Network to compute the source estimates. We demonstrate improvement in performance of the state-space algorithm relative to MNE in analyses of simulated and actual somatosensory MEG experiments. Our findings establish the benefits of high-dimensional state-space modeling as an effective means to solve the MEG source localization problem.
authors.affiliation: Imperial College London.
authors.forename: Christopher J
authors.initials: CJ
authors.lastname: Long
authors.forename: Patrick L
authors.initials: PL
authors.lastname: Purdon
authors.forename: Simona
authors.initials: S
authors.lastname: Temereanca
authors.forename: Neil U
authors.initials: NU
authors.lastname: Desai
authors.forename: Matti S
authors.initials: MS
authors.lastname: Hämäläinen
authors.forename: Emery N
authors.initials: EN
authors.lastname: Brown
references.citation: Neuroimage. 2008 Jul 1;41(3):924-40
references.pmid: 18455439
references.citation: Neuroimage. 2008 Feb 1;39(3):1104-20
references.pmid: 17997111
references.citation: IEEE Trans Biomed Eng. 1992 Jun;39(6):541-57
references.pmid: 1601435
references.citation: IEEE Trans Biomed Eng. 1997 Sep;44(9):867-80
references.pmid: 9282479
references.citation: Neuroimage. 2004 Oct;23(2):435-53
references.pmid: 15488394
references.citation: Int J Psychophysiol. 1994 Oct;18(1):49-65
references.pmid: 7876038
references.citation: Neuroimage. 2007 Aug 1;37(1):102-15
references.pmid: 17574444
references.citation: Neuroimage. 2005 Apr 15;25(3):734-45
references.pmid: 15808975
references.citation: Int J Psychophysiol. 2009 Jul;73(1):43-52
references.pmid: 19233235
references.citation: Electroencephalogr Clin Neurophysiol. 1994 Oct;91(4):237-48
references.pmid: 7523073
references.citation: Neuroimage. 2006 Apr 15;30(3):753-67
references.pmid: 16368248
references.citation: Neuroimage. 1999 Aug;10(2):173-80
references.pmid: 10417249
references.citation: IEEE Trans Biomed Eng. 1998 Nov;45(11):1342-54
references.pmid: 9805833
references.citation: Neuroimage. 1999 Feb;9(2):179-94
references.pmid: 9931268
references.citation: Neuroimage. 2005 Feb 15;24(4):997-1011
references.pmid: 15670677
references.citation: IEEE Trans Biomed Eng. 1998 Jun;45(6):716-23
references.pmid: 9609936
references.citation: Med Biol Eng Comput. 1994 Jan;32(1):35-42
references.pmid: 8182960
references.citation: Neuroimage. 2009 Feb 1;44(3):932-46
references.pmid: 18603008
references.citation: Proc Natl Acad Sci U S A. 1996 Aug 20;93(17):8809-15
references.pmid: 8799107
references.citation: Neuroimage. 2007 Oct 15;38(1):67-81
references.pmid: 17761440
references.citation: Neuroimage. 2010 Jan 1;49(1):641-55
references.pmid: 19596072
references.citation: Proc IEEE Int Symp Biomed Imaging. 2007 May 15;2007:1092-1095
references.pmid: 20407591
references.citation: Electroencephalogr Clin Neurophysiol. 1995 Dec;95(6):444-52
references.pmid: 8536573
references.citation: Hum Brain Mapp. 2004 Apr;21(4):221-35
references.pmid: 15038004
references.citation: IEEE Trans Biomed Eng. 1989 Feb;36(2):165-71
references.pmid: 2917762