Biomedical subjects

Riccardo Barbieri

Publications and source records attributed to Riccardo Barbieri.

7 recordsLinked to original sources

PubMed2006

Analysis of heartbeat dynamics by point process adaptive filtering.

Heartbeats are a point process yet, most of the current analysis methods do not model this important characteristic of these data. We describe human heartbeat time series as a history dependent inverse Gaussian model. We present a point process adaptive filter algorithm to estimate the model's time-varying parameters, and use it to compute new measures of heart rate variability. We apply our algorithm to analyze simulated heartbeat data and actual heartbeat data from a tilt table experiment and from healthy subjects and subjects with congestive heart failure during sleep. Our results suggest a new approach for characterizing heartbeat dynamics.

Algorithms↗

PubMed2005

An analysis of hippocampal spatio-temporal representations using a Bayesian algorithm for neural spike train decoding.

Neural spike train decoding algorithms are important tools for characterizing how ensembles of neurons represent biological signals. We present a Bayesian neural spike train decoding algorithm based on a point process model of individual neurons, a linear stochastic state-space model of the biological signal, and a temporal latency parameter. The latency parameter represents the temporal lead or lag between the biological signal and the ensemble spiking activity. We use the algorithm to study whether the representation of position by the ensemble spiking activity of pyramidal neurons in the CA1 region of the rat hippocampus is more consistent with prospective coding, i.e., future position, or retrospective coding, past position. Using 44 simultaneously recorded neurons and an ensemble delay latency of 400 ms, the median decoding error was 5.1 cm during 10 min of foraging in an open circular environment. The true coverage probability for the algorithm's 0.95 confidence regions was 0.71. These results illustrate how the Bayesian neural spike train decoding paradigm may be used to investigate spatio-temporal representations of position by an ensemble of hippocampal neurons.

Algorithms↗

PubMed2004

A point-process model of human heartbeat intervals: new definitions of heart rate and heart rate variability.

Heart rate is a vital sign, whereas heart rate variability is an important quantitative measure of cardiovascular regulation by the autonomic nervous system. Although the design of algorithms to compute heart rate and assess heart rate variability is an active area of research, none of the approaches considers the natural point-process structure of human heartbeats, and none gives instantaneous estimates of heart rate variability. We model the stochastic structure of heartbeat intervals as a history-dependent inverse Gaussian process and derive from it an explicit probability density that gives new definitions of heart rate and heart rate variability: instantaneous R-R interval and heart rate standard deviations. We estimate the time-varying parameters of the inverse Gaussian model by local maximum likelihood and assess model goodness-of-fit by Kolmogorov-Smirnov tests based on the time-rescaling theorem. We illustrate our new definitions in an analysis of human heartbeat intervals from 10 healthy subjects undergoing a tilt-table experiment. Although several studies have identified deterministic, nonlinear dynamical features in human heartbeat intervals, our analysis shows that a highly accurate description of these series at rest and in extreme physiological conditions may be given by an elementary, physiologically based, stochastic model.

Adult↗

PubMed2004

Dynamic analyses of information encoding in neural ensembles.

Neural spike train decoding algorithms and techniques to compute Shannon mutual information are important methods for analyzing how neural systems represent biological signals. Decoding algorithms are also one of several strategies being used to design controls for brain-machine interfaces. Developing optimal strategies to design decoding algorithms and compute mutual information are therefore important problems in computational neuroscience. We present a general recursive filter decoding algorithm based on a point process model of individual neuron spiking activity and a linear stochastic state-space model of the biological signal. We derive from the algorithm new instantaneous estimates of the entropy, entropy rate, and the mutual information between the signal and the ensemble spiking activity. We assess the accuracy of the algorithm by computing, along with the decoding error, the true coverage probability of the approximate 0.95 confidence regions for the individual signal estimates. We illustrate the new algorithm by reanalyzing the position and ensemble neural spiking activity of CA1 hippocampal neurons from two rats foraging in an open circular environment. We compare the performance of this algorithm with a linear filter constructed by the widely used reverse correlation method. The median decoding error for Animal 1 (2) during 10 minutes of open foraging was 5.9 (5.5) cm, the median entropy was 6.9 (7.0) bits, the median information was 9.4 (9.4) bits, and the true coverage probability for 0.95 confidence regions was 0.67 (0.75) using 34 (32) neurons. These findings improve significantly on our previous results and suggest an integrated approach to dynamically reading neural codes, measuring their properties, and quantifying the accuracy with which encoded information is extracted.

Action Potentials↗

PubMed2004

Dynamic analysis of neural encoding by point process adaptive filtering.

Neural receptive fields are dynamic in that with experience, neurons change their spiking responses to relevant stimuli. To understand how neural systems adapt their representations of biological information, analyses of receptive field plasticity from experimental measurements are crucial. Adaptive signal processing, the well-established engineering discipline for characterizing the temporal evolution of system parameters, suggests a framework for studying the plasticity of receptive fields. We use the Bayes' rule Chapman-Kolmogorov paradigm with a linear state equation and point process observation models to derive adaptive filters appropriate for estimation from neural spike trains. We derive point process filter analogues of the Kalman filter, recursive least squares, and steepest-descent algorithms and describe the properties of these new filters. We illustrate our algorithms in two simulated data examples. The first is a study of slow and rapid evolution of spatial receptive fields in hippocampal neurons. The second is an adaptive decoding study in which a signal is decoded from ensemble neural spiking activity as the receptive fields of the neurons in the ensemble evolve. Our results provide a paradigm for adaptive estimation for point process observations and suggest a practical approach for constructing filtering algorithms to track neural receptive field dynamics on a millisecond timescale.

Action Potentials↗

PubMed2002

Heart rate control and mechanical cardiopulmonary coupling to assess central volume: a systems analysis.

Small negative changes of central volume reduce cardiac output without significant alterations of arterial blood pressure (ABP), suggesting an adequate regulatory response. Furthermore, evidence has arisen supporting a Bainbridge reflex (tachycardia with hypervolemia) in humans. To investigate these phenomena, multivariate autoregressive techniques were used to evaluate the beat-to-beat interactions between respiration, R-R interval, and ABP at six levels of decreased and increased central volume. With reductions of central volume below control, baroreflex and respiratory sinus arrhythmia gains were reduced, while with increases of volume above control, gains increased for the first two levels but decreased again at the highest volume level, suggesting the presence of a Bainbridge reflex in healthy human subjects. The mechanical influence of respiration on central venous pressure (CVP) had an unexpected shift in phase at the point of mild central hypervolemia, with the expected negative relation at lower volumes (inspiration lowers CVP) but a positive relation at higher volumes (inspiration raises CVP). We conclude that multivariate techniques can quantify the relations between a variety of respiratory and hemodynamic parameters, allowing for the in vivo assessment of complex cardiorespiratory interactions during manipulations of central volume. The results identify the presence of a Bainbridge reflex in humans and suggest that short-term cardiovascular control is optimized at mild hypervolemia.

Adult↗

PubMed2002

The time-rescaling theorem and its application to neural spike train data analysis.

Measuring agreement between a statistical model and a spike train data series, that is, evaluating goodness of fit, is crucial for establishing the model's validity prior to using it to make inferences about a particular neural system. Assessing goodness-of-fit is a challenging problem for point process neural spike train models, especially for histogram-based models such as perstimulus time histograms (PSTH) and rate functions estimated by spike train smoothing. The time-rescaling theorem is a well-known result in probability theory, which states that any point process with an integrable conditional intensity function may be transformed into a Poisson process with unit rate. We describe how the theorem may be used to develop goodness-of-fit tests for both parametric and histogram-based point process models of neural spike trains. We apply these tests in two examples: a comparison of PSTH, inhomogeneous Poisson, and inhomogeneous Markov interval models of neural spike trains from the supplementary eye field of a macque monkey and a comparison of temporal and spatial smoothers, inhomogeneous Poisson, inhomogeneous gamma, and inhomogeneous inverse gaussian models of rat hippocampal place cell spiking activity. To help make the logic behind the time-rescaling theorem more accessible to researchers in neuroscience, we present a proof using only elementary probability theory arguments. We also show how the theorem may be used to simulate a general point process model of a spike train. Our paradigm makes it possible to compare parametric and histogram-based neural spike train models directly. These results suggest that the time-rescaling theorem can be a valuable tool for neural spike train data analysis.

Action Potentials↗