Probabilistic Modelling of Probabilistic Modelling of Brain Imaging Data Brain Imaging Data Will Penny Will Penny The Wellcome Department of Imaging The Wellcome Department of Imaging Neuroscience, UCL Neuroscience, UCL http//:www.fil.ion.ucl.ac.uk/~wpenny http//:www.fil.ion.ucl.ac.uk/~wpenny
Probabilistic Modelling of Brain Imaging Data. Will Penny The Wellcome Department of Imaging Neuroscience, UCL http//:www.fil.ion.ucl.ac.uk/~wpenny. Overview. Multiple levels of Bayesian Inference 2. A model of fMRI time series: The Noise 3. A model of fMRI time series: The Signal. - PowerPoint PPT Presentation
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Probabilistic Modelling of Probabilistic Modelling of Brain Imaging DataBrain Imaging Data
Probabilistic Modelling of Probabilistic Modelling of Brain Imaging DataBrain Imaging Data
Will PennyWill Penny
The Wellcome Department of Imaging Neuroscience, The Wellcome Department of Imaging Neuroscience, UCLUCL
fMRI time series modelfMRI time series modelfMRI time series modelfMRI time series model
• Use a General Linear Model at each voxel:Use a General Linear Model at each voxel:
y = X y = X + e + e
where X contains task-related regressors.where X contains task-related regressors.
• Use a General Linear Model at each voxel:Use a General Linear Model at each voxel:
y = X y = X + e + e
where X contains task-related regressors.where X contains task-related regressors.
fMRI time series modelfMRI time series modelfMRI time series modelfMRI time series model
= +
y X e
Time-seriesat one spatial
location
Putative effects ofexperimentalmanipulation
Sizeof effects
Residuals
fMRI time series modelfMRI time series modelfMRI time series modelfMRI time series model
• Use a General Linear Model at each voxel:Use a General Linear Model at each voxel:
y = X y = X + e + e
where X contains task-related regressors.where X contains task-related regressors.
• The errors are modelled as an AR(p) process.The errors are modelled as an AR(p) process. ((Parametric spectral estimationParametric spectral estimation))
• The order can be selected using Bayesian evidenceThe order can be selected using Bayesian evidence
• Use a General Linear Model at each voxel:Use a General Linear Model at each voxel:
y = X y = X + e + e
where X contains task-related regressors.where X contains task-related regressors.
• The errors are modelled as an AR(p) process.The errors are modelled as an AR(p) process. ((Parametric spectral estimationParametric spectral estimation))
• The order can be selected using Bayesian evidenceThe order can be selected using Bayesian evidence
Synthetic GLM-AR(3) DataSynthetic GLM-AR(3) DataSynthetic GLM-AR(3) DataSynthetic GLM-AR(3) Data
Map of AR model order, pMap of AR model order, pMap of AR model order, pMap of AR model order, p
p=0,1,2,3FaceData
AngiogramsAngiogramsAngiogramsAngiograms
Other subjects, aOther subjects, a11Other subjects, aOther subjects, a11
Ring ofvoxels with
highly correlatederror
Other subjects, aOther subjects, a11Other subjects, aOther subjects, a11
Unmodelledsignal
orincreasedcardiac
artifact due to increasedblood flow?
OverviewOverviewOverviewOverview
1.1. Multiple levels of Bayesian InferenceMultiple levels of Bayesian Inference
2. A model of fMRI time series: The Noise2. A model of fMRI time series: The Noise
3. 3. A model of fMRI time series: The SignalA model of fMRI time series: The Signal
1.1. Multiple levels of Bayesian InferenceMultiple levels of Bayesian Inference
2. A model of fMRI time series: The Noise2. A model of fMRI time series: The Noise
3. 3. A model of fMRI time series: The SignalA model of fMRI time series: The Signal
fMRI time series modelfMRI time series modelfMRI time series modelfMRI time series model
• Use a General Linear Model Use a General Linear Model
y = X y = X + e + e
• Priors factorise into groups:Priors factorise into groups:
p(p() = p() = p(11) p() p(22) p() p(33))
• Priors in each group may be smoothness Priors in each group may be smoothness priors or Gaussianspriors or Gaussians
• Use a General Linear Model Use a General Linear Model
y = X y = X + e + e
• Priors factorise into groups:Priors factorise into groups:
p(p() = p() = p(11) p() p(22) p() p(33))
• Priors in each group may be smoothness Priors in each group may be smoothness priors or Gaussianspriors or Gaussians
Rik’s dataRik’s dataRik’s dataRik’s data
24 Transverse Slices acquired with TR=2s
Press left key if famous, right key if not
Time series of 351 images
Part of larger study lookingat factors influencing repetition suppresion
Every face presented twice
Modelling the SignalModelling the SignalModelling the SignalModelling the Signal
Assumption: Neuronal Event Stream is Identical to the Experimental Event Stream
Convolve event-stream with basis functions to account for the HRF
FIR modelsFIR modelsFIR modelsFIR models
Timeafterevent
Sizeof
signal
FIR modelFIR modelFIR modelFIR model
Separate smoothness priors for each event type
Design matrixfor FIR model with
8 time bins in a 20-second window
Q. Is this a good prior ?
FIR basis setFIR basis setFIR basis setFIR basis set
Left occipital cortex (x=-33, y=-81, z=-24)
FIR model average responses
FIR basis setFIR basis setFIR basis setFIR basis set
Right fusiform cortex (x=45, y=-60, z=-18)
FIR model average responses
RFX-Event modelRFX-Event modelRFX-Event modelRFX-Event model
Design Matrix
97 parameters ! But only 24 effective parameters
Responses to each event of type A are randomly distributed about some typical “type A” response