25 research outputs found
Algorithms for the analysis of ensemble neural spiking activity using simultaneous-event multivariate point-process models
Understanding how ensembles of neurons represent and transmit information in the patterns of their joint spiking activity is a fundamental question in computational neuroscience. At present, analyses of spiking activity from neuronal ensembles are limited because multivariate point process (MPP) models cannot represent simultaneous occurrences of spike events at an arbitrarily small time resolution. Solo recently reported a simultaneous-event multivariate point process (SEMPP) model to correct this key limitation. In this paper, we show how Solo's discrete-time formulation of the SEMPP model can be efficiently fit to ensemble neural spiking activity using a multinomial generalized linear model (mGLM). Unlike existing approximate procedures for fitting the discrete-time SEMPP model, the mGLM is an exact algorithm. The MPP time-rescaling theorem can be used to assess model goodness-of-fit. We also derive a new marked point-process (MkPP) representation of the SEMPP model that leads to new thinning and time-rescaling algorithms for simulating an SEMPP stochastic process. These algorithms are much simpler than multivariate extensions of algorithms for simulating a univariate point process, and could not be arrived at without the MkPP representation. We illustrate the versatility of the SEMPP model by analyzing neural spiking activity from pairs of simultaneously-recorded rat thalamic neurons stimulated by periodic whisker deflections, and by simulating SEMPP data. In the data analysis example, the SEMPP model demonstrates that whisker motion significantly modulates simultaneous spiking activity at the 1 ms time scale and that the stimulus effect is more than one order of magnitude greater for simultaneous activity compared with non-simultaneous activity. Together, the mGLM, the MPP time-rescaling theorem and the MkPP representation of the SEMPP model offer a theoretically sound, practical tool for measuring joint spiking propensity in a neuronal ensemble
State-space solutions to the dynamic magnetoencephalography inverse problem using high performance computing
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 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 intracortical 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.Comment: Published in at http://dx.doi.org/10.1214/11-AOAS483 the Annals of
Applied Statistics (http://www.imstat.org/aoas/) by the Institute of
Mathematical Statistics (http://www.imstat.org
Local Field Potentials and the Encoding of Whisker Deflections by Population Firing Synchrony in Thalamic Barreloids
Left-Lateralized Contributions of Saccades to Cortical Activity During a One-Back Word Recognition Task
Saccadic eye movements are an inherent component of natural reading, yet their contribution to information processing at subsequent fixation remains elusive. Here we use anatomically-constrained magnetoencephalography (MEG) to examine cortical activity following saccades as healthy human subjects engaged in a one-back word recognition task. This activity was compared with activity following external visual stimulation that mimicked saccades. A combination of procedures was employed to eliminate saccadic ocular artifacts from the MEG signal. Both saccades and saccade-like external visual stimulation produced early-latency responses beginning ~70 ms after onset in occipital cortex and spreading through the ventral and dorsal visual streams to temporal, parietal and frontal cortices. Robust differential activity following the onset of saccades vs. similar external visual stimulation emerged during 150â350 ms in a left-lateralized cortical network. This network included: (i) left lateral occipitotemporal (LOT) and nearby inferotemporal (IT) cortex; (ii) left posterior Sylvian fissure (PSF) and nearby multimodal cortex; and (iii) medial parietooccipital (PO), posterior cingulate and retrosplenial cortices. Moreover, this left-lateralized network colocalized with word repetition priming effects. Together, results suggest that central saccadic mechanisms influence a left-lateralized language network in occipitotemporal and temporal cortex above and beyond saccadic influences at preceding stages of information processing during visual word recognition