Biomedical subjects

William D Penny

Publications and source records attributed to William D Penny.

5 recordsLinked to original sources

PubMed2005

MEG source localization under multiple constraints: an extended Bayesian framework.

To use Electroencephalography (EEG) and Magnetoencephalography (MEG) as functional brain 3D imaging techniques, identifiable distributed source models are required. The reconstruction of EEG/MEG sources rests on inverting these models and is ill-posed because the solution does not depend continuously on the data and there is no unique solution in the absence of prior information or constraints. We have described a general framework that can account for several priors in a common inverse solution. An empirical Bayesian framework based on hierarchical linear models was proposed for the analysis of functional neuroimaging data [Friston, K., Penny, W., Phillips, C., Kiebel, S., Hinton, G., Ashburner, J., 2002. Classical and Bayesian inference in neuroimaging: theory. NeuroImage 16, 465-483] and was evaluated recently in the context of EEG [Phillips, C., Mattout, J., Rugg, M.D., Maquet, P., Friston, K., 2005. An empirical Bayesian solution to the source reconstruction problem in EEG. NeuroImage 24, 997-1011]. The approach consists of estimating the expected source distribution and its conditional variance that is constrained by an empirically determined mixture of prior variance components. Estimation uses Expectation-Maximization (EM) to give the Restricted Maximum Likelihood (ReML) estimate of the variance components (in terms of hyperparameters) and the Maximum A Posteriori (MAP) estimate of the source parameters. In this paper, we extend the framework to compare different combinations of priors, using a second level of inference based on Bayesian model selection. Using Monte-Carlo simulations, ReML is first compared to a classic Weighted Minimum Norm (WMN) solution under a single constraint. Then, the ReML estimates are evaluated using various combinations of priors. Both standard criterion and ROC-based measures were used to assess localization and detection performance. The empirical Bayes approach proved useful as: (1) ReML was significantly better than WMN for single priors; (2) valid location priors improved ReML source localization; (3) invalid location priors did not significantly impair performance. Finally, we show how model selection, using the log-evidence, can be used to select the best combination of priors. This enables a global strategy for multiple prior-based regularization of the MEG/EEG source reconstruction.

Bayes Theorem↗

PubMed2005

Conjunction revisited.

The aim of this note is to revisit the analysis of conjunctions in imaging data. We review some conceptual issues that have emerged from recent discussion (Nichols, T., Brett, M., Andersson, J., Wager, T., Poline, J.-B., 2004. Valid Conjunction Inference with the Minimum Statistic.) and reformulate the conjunction of null hypotheses as a conjunction of k or more effects. Analyses based on minimum statistics have typically used the null hypothesis that k = 0. This enables inferences about one or more effects (k > 0). However, this does not provide control over false-positive rates (FPR) for inferences about a conjunction of k = n effects, over n tests. This is the key point made by Nichols et al., who suggest a procedure based on supremum P values that provides an upper bound on FPR for k = n. Although valid, this is a very conservative procedure, particularly in the context of multiple comparisons. We suggest that an inference on a conjunction of k = n effects is generally unnecessary and distinguish between congruent contrasts that test for the same treatment and incongruent contrasts of the sort used in cognitive conjunctions. For congruent contrasts, the usual inference, k > 0, is sufficient. With incongruent contrasts it is sufficient to infer a conjunction of k >u effects, where u is the number of contrasts that share some uninteresting effect. The issues highlighted by Nichols et al., have important implications for the design and analysis of cognitive conjunction studies and have motivated a change to the SPM software, that affords a test for the more general hypothesis k >u. This more general conjunction test is described.

Arousal↗

PubMed2005

Bayesian fMRI time series analysis with spatial priors.

We describe a Bayesian estimation and inference procedure for fMRI time series based on the use of General Linear Models (GLMs). Importantly, we use a spatial prior on regression coefficients which embodies our prior knowledge that evoked responses are spatially contiguous and locally homogeneous. Further, using a computationally efficient Variational Bayes framework, we are able to let the data determine the optimal amount of smoothing. We assume an arbitrary order Auto-Regressive (AR) model for the errors. Our model generalizes earlier work on voxel-wise estimation of GLM-AR models and inference in GLMs using Posterior Probability Maps (PPMs). Results are shown on simulated data and on data from an event-related fMRI experiment.

Bayes Theorem↗

PubMed2003

Modeling regional and psychophysiologic interactions in fMRI: the importance of hemodynamic deconvolution.

The analysis of functional magnetic resonance imaging (fMRI) time-series data can provide information not only about task-related activity, but also about the connectivity (functional or effective) among regions and the influences of behavioral or physiologic states on that connectivity. Similar analyses have been performed in other imaging modalities, such as positron emission tomography. However, fMRI is unique because the information about the underlying neuronal activity is filtered or convolved with a hemodynamic response function. Previous studies of regional connectivity in fMRI have overlooked this convolution and have assumed that the observed hemodynamic response approximates the neuronal response. In this article, this assumption is revisited using estimates of underlying neuronal activity. These estimates use a parametric empirical Bayes formulation for hemodynamic deconvolution.

Bayes Theorem↗

PubMed2002

Bayesian nonstationary autoregressive models for biomedical signal analysis.

We describe a variational Bayesian algorithm for the estimation of a multivariate autoregressive model with time-varying coefficients that adapt according to a linear dynamical system. The algorithm allows for time and frequency domain characterization of nonstationary multivariate signals and is especially suited to the analysis of event-related data. Results are presented on synthetic data and real electroencephalogram data recorded in event-related desynchronization and photic synchronization scenarios.

Algorithms↗