MIXED EFFECTS MODELS FOR TIME SERIES
|
|
- Laurence Small
- 6 years ago
- Views:
Transcription
1 Outline MIXED EFFECTS MODELS FOR TIME SERIES Cristina Gorrostieta Hakmook Kang Hernando Ombao Brown University Biostatistics Section February 16, 2011
2 Outline OUTLINE OF TALK 1 SCIENTIFIC MOTIVATION 2 BACKGROUND ON MIXED EFFECTS MODELS 3 MIXED EFFECTS VAR 4 SPATIO-SPECTRAL MIXED EFFECTS MODELS
3 Outline OUTLINE OF TALK 1 SCIENTIFIC MOTIVATION 2 BACKGROUND ON MIXED EFFECTS MODELS 3 MIXED EFFECTS VAR 4 SPATIO-SPECTRAL MIXED EFFECTS MODELS
4 Outline OUTLINE OF TALK 1 SCIENTIFIC MOTIVATION 2 BACKGROUND ON MIXED EFFECTS MODELS 3 MIXED EFFECTS VAR 4 SPATIO-SPECTRAL MIXED EFFECTS MODELS
5 Outline OUTLINE OF TALK 1 SCIENTIFIC MOTIVATION 2 BACKGROUND ON MIXED EFFECTS MODELS 3 MIXED EFFECTS VAR 4 SPATIO-SPECTRAL MIXED EFFECTS MODELS
6 ANALYSIS OF BRAIN SIGNALS Electrophysiologic data: multi-channel EEG, local field potentials Hemodynamic data: fmri time series at several ROIs
7 ANALYSIS OF BRAIN SIGNALS Goals of our research Characterize dependence in a brain network Temporal: Y 1 (t) [Y 1 (t 1), Y 2 (t 1),...] Spectral: interactions between oscillatory activities at Y 1, Y 2 Develop estimation and inference methods for connectivity Investigate potential for connectivity as a biomarker Predicting behavior Motor intent (left vs. right movement) State of learning Differentiating patient groups (bipolar vs. healthy children) Connectivity between left DLPFC right STG is greater for bipolar than healthy
8 ANALYSIS OF BRAIN SIGNALS Desirable Components of Model for Connectivity Fixed Effects Condition-specific or group-specific connectivity Test for differences across conditions or groups Random Effects In multiple subjects: subject-specific random deviation from the group effect In multiple trials: trial-specific random deviation from the mean condition effect Model variation in connectivity between subjects/participants in a group; and between trials within a subject
9 MIXED EFFECTS MODELS A Simple Regression Model Data for one subject: {(x t, Y t ), t = 1,...,T} A simple model Y t = f(x t )+ǫ t, ǫ t (0,σ 2 ) f(x t ) = β 0 +β 1 x t parametric f(x t ) = β k ψ k (x t ) non-parametric k β s are fixed effects Estimation and inference via maximum likelihood
10 MIXED EFFECTS MODELS Multi-subject data: {(xt n, Y t n ); t = 1,...,T; n = 1,...,N} Fixed effects only model Y n t = β 0 +β 1 x n t +ǫ n t Y n t Mixed effects model = k β k ψ(x n t )+ǫn t Y n t = β n 0 +βn 1 x n t +ǫ n t Y n t = k β n k ψ(x n t )+ǫ n t Decompose β n k = β k + b n k ; where bn k (0,τ 2 ) Eβk n = β k fixed effect (average across all subjects) bk n subject-specific random deviation from the group effect
11 ILLUSTRATION: ESTIMATING THE HEMODYNAMIC RESPONSE IN FMRI Goal is to estimate the brain hemodynamic response function in olfaction stimuli Standard HRF are estimated from visual stimuli In visual stimuli rise and decay are fast In olfaction, rise is slow, peak is persistent and decay is slow Experiment: 30 fmri trials; each trial 20 seconds; TR = 1 second
12 ILLUSTRATION: ESTIMATING THE HEMODYNAMIC RESPONSE IN FMRI fmri time series for several trials
13 ILLUSTRATION: ESTIMATING THE HEMODYNAMIC RESPONSE IN FMRI {φ k (t)} B-splines basis
14 ILLUSTRATION: ESTIMATING THE HEMODYNAMIC RESPONSE IN FMRI Time series for trial n HRF for trial n Y n (t) = f n (t)+ǫ n (t) f n (t) = f(t)+δ n (t) f(t) = k β k φ k (t) δ n (t) = k b n k φ k(t) f(t) - over-all olfaction HRF δ n trial-specific deviation from the over-all HRF Implemented in R (nlmix) and SAS (proc mixed)
15 ILLUSTRATION: ESTIMATING THE HEMODYNAMIC RESPONSE IN FMRI Trial-specific HRF estimate
16 ILLUSTRATION: ESTIMATING THE HEMODYNAMIC RESPONSE IN FMRI Over-all HRF estimate
17 ESTIMATION Model: Y n = X n β + Z n b n +η n Between-subjects Random Effect: b n N(0, D) Within-subject error: η n N(0, V) Unconditional distribution Y n N(X n β, Z n DZ n + V) Conditional distribution (subject n) Y n b n N(X n β + Z n b n, V) The subject-specific mean function E[Y n b n ] = X n β + Z n b n
18 ESTIMATION Marginal Model Y n = X n β +ǫ n ǫ n = Z n b n +η n N(0,Σ(θ)) ǫ n Y n N(X n β,σ(θ)) Estimates β and θ obtained by REML Random effects estimated by empirical Bayes Prior b n N(0, D) Likelihood Y n N(X n β, Z n DZ n +σ 2 I) Posterior Prior Likelihood 2 Plug in estimates D, β, σ b n = [Z nz n + D 1 σ 2 ] 1 Z n [ Y n X n β ]
19 MIXED-EFFECTS MODELS IN TIME SERIES Experiment N = 15 right-handed college students Experiment: subjects see visual targets and must move joystick Two Conditions Free choice - subject freely chooses any target Instructed - subject must choose the specified target Regions of interest (7 areas that show highest differential activation) PFC, SMA, etc.
20 MIXED-EFFECTS MODELS IN TIME SERIES
21 MIXED-EFFECTS MODELS IN TIME SERIES Total of R ROIs Y n r (t) the fmri time series at the ROI r for subject n Entire network: Y n (t) = [Y n 1 (t),..., Y n R (t)]. General additive model Y n (t) = F n (t) + E n (t) The components F n (t) - mean (deterministic) component E n (t) - stochastic component
22 MIXED-EFFECTS MODELS IN TIME SERIES The mean component F n (t) Decomposition F n (t) = D n (t)+m n (t)+β n 1 X 1(t)+...+β n C X C(t), The mean component includes systematic changes in the BOLD signal that is due to Scanner drift Physiological signals of non-interest (e.g., cardiac and respiratory) Experimental conditions
23 MIXED-EFFECTS MODELS IN TIME SERIES Methods for estimating the mean component Worsley and Friston, 1995; Nichols and Holmes, Implemented in Statistical Parametric Mapping (SPM) FMRIB Statistical Laboratory (FSL) Analysis of Functional NeuroImages (AFNI)
24 MIXED-EFFECTS MODELS IN TIME SERIES The stochastic component E n (t) E n (t) captures between-roi connectivity Cov[Y n (t + h), Y n (t)] = Cov[F n (t + h)+e n (t + h), F n (t)+e n (t)] = Cov[E n (t + h), E n (t)] E n (t) cannot be observed directly; we use the residuals: E n (t) = Y n (t) F n (t) R n (t) = Y n (t) F n (t)
25 MIXED-EFFECTS MODELS IN TIME SERIES Vector Auto-Regressive Model E(t) = [E 1 (t), E 2 (t)] VAR(1) Model E(t) = Φ 1 E(t 1)+η(t) Φ 1 lag-1 connectivity matrix; componentsφ kl VAR(1) Equations E 1 (t) = φ 11 E 1 (t 1)+φ 12 E 2 (t 1)+η 1 (t) E 2 (t) = φ 21 E 1 (t 1)+φ 22 E 2 (t 1)+η 2 (t)
26 MIXED-EFFECTS MODELS IN TIME SERIES Conditional MLE Model E(t) = Φ 1 E(t 1)+η(t), t = 1,...,T Define E(t 1) = {E(t 1), E(t 2),...} Suppose η(t) iid N(0, Σ) Condition on past data E(t 1) E(t) E(t 1) N(Φ 1 X(t 1),Σ) Conditional likelihood function L C [Φ 1 E(1)] = f(e(2) E(1))... f(e(t) E(T 1))) Conditional approach leads to a closed form for the estimator of Φ 1 Full likelihood - no closed form
27 MIXED-EFFECTS MODELS IN TIME SERIES The ME-VAR(1) Model [ ] E (n) (t) = Φ 1,k W 1 (t)+φ 2,k W 2 (t)+b (n) 1 E (n) (t 1)+e (n) (t) W (t) is the indicator function When condition 1 is active then W 1 (t) = 1 and W 2 (t) = 0 When condition 2 is active then W 1 (t) = 0 and W 2 (t) = 1 b (n) 1 models between-subject variation in connectivity b (n) 1 (0,σ 2 1 )
28 MIXED-EFFECTS MODELS IN TIME SERIES Subject-specific connectivity matrix (condition on b (n) 1 ) When W 1 (t) = 1, the connectivity matrix for subject n is Φ 1,1 + b (n) 1 When W 2 (t) = 1, the connectivity matrix for subject n is Φ 1,2 + b (n) 1 The model can be utilized to test for Lagged dependence between each pair of ROIs H 0 : Φ 1,1 = 0,Φ 1,2 = 0 Granger causality in each experimental condition Testing for differences between conditions H 0 : 1 = Φ 1,1 Φ 1,2 = 0
29 CONNECTIVITY ANALYSIS OF FMRI DATA Step 0. Selection of pre-defined ROIs. Step 1. Estimate the mean function at each ROI F n (t). Step 2. Fit the ME-VAR model to the residuals R n (t).
30 CONNECTIVITY ANALYSIS OF FMRI DATA Fitting the ME-VAR Model Optimal lag over was p = 2 using the Bayesian Information Criterion (Pmax = 8) In physical time: TR = 2 seconds; optimal delay/lag is 4 seconds. Computations were carried out in SAS proc mixed (we use existing machinery!)
31 CONNECTIVITY ANALYSIS OF FMRI DATA
32 CONNECTIVITY ANALYSIS OF FMRI DATA
33 CONNECTIVITY ANALYSIS OF FMRI DATA
34 CONCLUSION Statistical vs Clinical significance Connectivity parameter of can be statistically significant. Is it clinically meaningful? How do we measure clinical significance? Impact of increase (or decrease) in connectivity measure on some clinical outcome Should we be testing for H 0 : θ = 0 or H 0 : θ > η?
35 CONCLUSION Granger causality not physiological causality Connectivity in the hemodynamic activity not necessarily neuronal temporal ordering could be switched because of the poor temporal resolution in fmri
36 SPATIO-SPECTRAL MIXED EFFECTS MODELS Spatio-temporal models In fmri, spatial correlation adds computational complications! Current approaches Bowman (2005, 2007, 2008); Worsley (1999); Valdes-Sosa et al. (2004) Spatial - parametric; functional" distance Bayesian hierarchial 4D-Wavelet packets Partial least squares
37 SPATIO-SPECTRAL MIXED EFFECTS MODELS Spatio-Spectral Model With H. Kang (PhD student) [John Van Ryzin Award 2011] Goals: Estimate ROI-specific activation Use all local and global information Approach Spectral domain: Fourier coefficients are approx uncorrelated Spatio-spectral covariance matrix is block diagonal - simplified! Multi-scale correlation - local (within ROI) and global (between ROIs)
38 SPATIO-SPECTRAL MIXED EFFECTS MODELS ROI1 ROI2 ROI3
39 SPATIO-SPECTRAL MIXED EFFECTS MODELS Spatio-temporal Model Y cv (t) = P [βc p + bcv]x p p (t)+d c (t)+ǫ cv (t) p=1 β p c is the ROI-specific fixed effect due to stimulus p; b p cv is a zero-mean voxel-specific random deviation local spatial covariance between voxels within an ROI Cov(b p cv, b p c v ) = δ(c c )ψ b ( v v ) d c (t) is a zero-mean ROI-specific signal global spatial covariation Cov(d c (t), d c (t)) = ψ d (c, c ) ǫ cv (t) is the within-voxel temporal noise
40 SPATIO-SPECTRAL MIXED EFFECTS MODELS Problem: Spatio-temporal covariance structure is complicated! A solution: Fourier transform X(t) stationary time series with covariance matrix Σ d(ω k ) = t X(t) exp( i2πω kt) Fourier coefficient I(ω k ) = d(ω k ) 2 I(ω 1 ),...,I(ω M ) are approx uncorrelated
41 SPATIO-SPECTRAL MIXED EFFECTS MODELS Y cv (t) time series at voxel v in ROI c Y cv (ω) corresponding Fourier coefficient Spatio-Spectral Model P Y cv (ω k ) = [βc p + bcv]x p p (ω k )+d c (ω k )+ǫ cv (ω k ) p=1 β p c is the ROI-specific fixed effect due to stimulus p b p cv is a zero-mean voxel-specific random effect d c (t) is a zero-mean ROI-specific random effect
42 SPATIO-SPECTRAL MIXED EFFECTS MODELS
The General Linear Model (GLM)
he General Linear Model (GLM) Klaas Enno Stephan ranslational Neuromodeling Unit (NU) Institute for Biomedical Engineering University of Zurich & EH Zurich Wellcome rust Centre for Neuroimaging Institute
More informationGranger Mediation Analysis of Functional Magnetic Resonance Imaging Time Series
Granger Mediation Analysis of Functional Magnetic Resonance Imaging Time Series Yi Zhao and Xi Luo Department of Biostatistics Brown University June 8, 2017 Overview 1 Introduction 2 Model and Method 3
More informationCausal modeling of fmri: temporal precedence and spatial exploration
Causal modeling of fmri: temporal precedence and spatial exploration Alard Roebroeck Maastricht Brain Imaging Center (MBIC) Faculty of Psychology & Neuroscience Maastricht University Intro: What is Brain
More informationEffective Connectivity & Dynamic Causal Modelling
Effective Connectivity & Dynamic Causal Modelling Hanneke den Ouden Donders Centre for Cognitive Neuroimaging Radboud University Nijmegen Advanced SPM course Zurich, Februari 13-14, 2014 Functional Specialisation
More informationThe General Linear Model. Guillaume Flandin Wellcome Trust Centre for Neuroimaging University College London
The General Linear Model Guillaume Flandin Wellcome Trust Centre for Neuroimaging University College London SPM Course Lausanne, April 2012 Image time-series Spatial filter Design matrix Statistical Parametric
More informationDynamic Causal Modelling for fmri
Dynamic Causal Modelling for fmri André Marreiros Friday 22 nd Oct. 2 SPM fmri course Wellcome Trust Centre for Neuroimaging London Overview Brain connectivity: types & definitions Anatomical connectivity
More informationMODELING DEPENDENCE IN MULTIVARIATE TIME SERIES WITH APPLICATIONS TO BRAIN SIGNALS
Outline MODELING DEPENDENCE IN MULTIVARIATE TIME SERIES WITH APPLICATIONS TO BRAIN SIGNALS Hernando Ombao University of California at Irvine May 15, 2014 Outline OUTLINE OF TALK 1 SCIENTIFIC MOTIVATION
More informationStatistical Analysis Aspects of Resting State Functional Connectivity
Statistical Analysis Aspects of Resting State Functional Connectivity Biswal s result (1995) Correlations between RS Fluctuations of left and right motor areas Why studying resting state? Human Brain =
More informationPart 2: Multivariate fmri analysis using a sparsifying spatio-temporal prior
Chalmers Machine Learning Summer School Approximate message passing and biomedicine Part 2: Multivariate fmri analysis using a sparsifying spatio-temporal prior Tom Heskes joint work with Marcel van Gerven
More informationOverview of Spatial Statistics with Applications to fmri
with Applications to fmri School of Mathematics & Statistics Newcastle University April 8 th, 2016 Outline Why spatial statistics? Basic results Nonstationary models Inference for large data sets An example
More informationBeyond Univariate Analyses: Multivariate Modeling of Functional Neuroimaging Data
Beyond Univariate Analyses: Multivariate Modeling of Functional Neuroimaging Data F. DuBois Bowman Department of Biostatistics and Bioinformatics Center for Biomedical Imaging Statistics Emory University,
More informationStatistics of Stochastic Processes
Prof. Dr. J. Franke All of Statistics 4.1 Statistics of Stochastic Processes discrete time: sequence of r.v...., X 1, X 0, X 1, X 2,... X t R d in general. Here: d = 1. continuous time: random function
More informationBayesian Analysis. Bayesian Analysis: Bayesian methods concern one s belief about θ. [Current Belief (Posterior)] (Prior Belief) x (Data) Outline
Bayesian Analysis DuBois Bowman, Ph.D. Gordana Derado, M. S. Shuo Chen, M. S. Department of Biostatistics and Bioinformatics Center for Biomedical Imaging Statistics Emory University Outline I. Introduction
More informationModelling temporal structure (in noise and signal)
Modelling temporal structure (in noise and signal) Mark Woolrich, Christian Beckmann*, Salima Makni & Steve Smith FMRIB, Oxford *Imperial/FMRIB temporal noise: modelling temporal autocorrelation temporal
More informationHST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS
HST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS EMERY N. BROWN AND CHRIS LONG NEUROSCIENCE STATISTICS RESEARCH LABORATORY DEPARTMENT
More informationAn Implementation of Dynamic Causal Modelling
An Implementation of Dynamic Causal Modelling Christian Himpe 24.01.2011 Overview Contents: 1 Intro 2 Model 3 Parameter Estimation 4 Examples Motivation Motivation: How are brain regions coupled? Motivation
More informationExperimental design of fmri studies
Experimental design of fmri studies Sandra Iglesias With many thanks for slides & images to: Klaas Enno Stephan, FIL Methods group, Christian Ruff SPM Course 2015 Overview of SPM Image time-series Kernel
More informationExperimental design of fmri studies & Resting-State fmri
Methods & Models for fmri Analysis 2016 Experimental design of fmri studies & Resting-State fmri Sandra Iglesias With many thanks for slides & images to: Klaas Enno Stephan, FIL Methods group, Christian
More informationExperimental design of fmri studies
Methods & Models for fmri Analysis 2017 Experimental design of fmri studies Sara Tomiello With many thanks for slides & images to: Sandra Iglesias, Klaas Enno Stephan, FIL Methods group, Christian Ruff
More informationExperimental design of fmri studies
Experimental design of fmri studies Zurich SPM Course 2016 Sandra Iglesias Translational Neuromodeling Unit (TNU) Institute for Biomedical Engineering (IBT) University and ETH Zürich With many thanks for
More informationDynamic Causal Modelling for EEG/MEG: principles J. Daunizeau
Dynamic Causal Modelling for EEG/MEG: principles J. Daunizeau Motivation, Brain and Behaviour group, ICM, Paris, France Overview 1 DCM: introduction 2 Dynamical systems theory 3 Neural states dynamics
More informationData Analysis I: Single Subject
Data Analysis I: Single Subject ON OFF he General Linear Model (GLM) y= X fmri Signal = Design Matrix our data = what we CAN explain x β x Betas + + how much x of it we CAN + explain ε Residuals what
More informationStochastic Dynamic Causal Modelling for resting-state fmri
Stochastic Dynamic Causal Modelling for resting-state fmri Ged Ridgway, FIL, Wellcome Trust Centre for Neuroimaging, UCL Institute of Neurology, London Overview Connectivity in the brain Introduction to
More informationADAPTIVE SPECTRAL ANALYSIS OF REPLICATED NONSTATIONARY TIME SERIES
ADAPTIVE SPECTRAL ANALYSIS OF REPLICATED NONSTATIONARY TIME SERIES ROBERT KRAFTY Department of Biostatistics University of Pittsburgh Joint work with: Scott Bruce Department of Statistical Science, Temple
More informationThe General Linear Model (GLM)
The General Linear Model (GLM) Dr. Frederike Petzschner Translational Neuromodeling Unit (TNU) Institute for Biomedical Engineering, University of Zurich & ETH Zurich With many thanks for slides & images
More informationGroup analysis. Jean Daunizeau Wellcome Trust Centre for Neuroimaging University College London. SPM Course Edinburgh, April 2010
Group analysis Jean Daunizeau Wellcome Trust Centre for Neuroimaging University College London SPM Course Edinburgh, April 2010 Image time-series Spatial filter Design matrix Statistical Parametric Map
More informationClass: Trend-Cycle Decomposition
Class: Trend-Cycle Decomposition Macroeconometrics - Spring 2011 Jacek Suda, BdF and PSE June 1, 2011 Outline Outline: 1 Unobserved Component Approach 2 Beveridge-Nelson Decomposition 3 Spectral Analysis
More informationWavelet-Based Nonparametric Modeling of Hierarchical Functions in Colon Carcinogenesis
Wavelet-Based Nonparametric Modeling of Hierarchical Functions in Colon Carcinogenesis Jeffrey S. Morris University of Texas, MD Anderson Cancer Center Joint wor with Marina Vannucci, Philip J. Brown,
More informationSpatio-Spectral Mixed-Effects Model for Functional Magnetic Resonance Imaging Data
This article was downloaded by: [Brown University] On: 30 July 2012, At: 10:27 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: Mortimer
More informationMixed effects and Group Modeling for fmri data
Mixed effects and Group Modeling for fmri data Thomas Nichols, Ph.D. Department of Statistics Warwick Manufacturing Group University of Warwick Warwick fmri Reading Group May 19, 2010 1 Outline Mixed effects
More informationOverview of SPM. Overview. Making the group inferences we want. Non-sphericity Beyond Ordinary Least Squares. Model estimation A word on power
Group Inference, Non-sphericity & Covariance Components in SPM Alexa Morcom Edinburgh SPM course, April 011 Centre for Cognitive & Neural Systems/ Department of Psychology University of Edinburgh Overview
More informationSwitching Regime Estimation
Switching Regime Estimation Series de Tiempo BIrkbeck March 2013 Martin Sola (FE) Markov Switching models 01/13 1 / 52 The economy (the time series) often behaves very different in periods such as booms
More informationProblem Set 2 Solution Sketches Time Series Analysis Spring 2010
Problem Set 2 Solution Sketches Time Series Analysis Spring 2010 Forecasting 1. Let X and Y be two random variables such that E(X 2 ) < and E(Y 2 )
More informationGroup Analysis. Lexicon. Hierarchical models Mixed effect models Random effect (RFX) models Components of variance
Group Analysis J. Daunizeau Institute of Empirical Research in Economics, Zurich, Switzerland Brain and Spine Institute, Paris, France SPM Course Edinburgh, April 2011 Image time-series Spatial filter
More informationFunctional Causal Mediation Analysis with an Application to Brain Connectivity. Martin Lindquist Department of Biostatistics Johns Hopkins University
Functional Causal Mediation Analysis with an Application to Brain Connectivity Martin Lindquist Department of Biostatistics Johns Hopkins University Introduction Functional data analysis (FDA) and causal
More informationEvent-related fmri. Christian Ruff. Laboratory for Social and Neural Systems Research Department of Economics University of Zurich
Event-related fmri Christian Ruff Laboratory for Social and Neural Systems Research Department of Economics University of Zurich Institute of Neurology University College London With thanks to the FIL
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationTowards a Regression using Tensors
February 27, 2014 Outline Background 1 Background Linear Regression Tensorial Data Analysis 2 Definition Tensor Operation Tensor Decomposition 3 Model Attention Deficit Hyperactivity Disorder Data Analysis
More informationDetecting fmri activation allowing for unknown latency of the hemodynamic response
Detecting fmri activation allowing for unknown latency of the hemodynamic response K.J. Worsley McGill University J.E. Taylor Stanford University January 7, 006 Abstract Several authors have suggested
More informationDynamic Modeling of Brain Activity
0a Dynamic Modeling of Brain Activity EIN IIN PC Thomas R. Knösche, Leipzig Generative Models for M/EEG 4a Generative Models for M/EEG states x (e.g. dipole strengths) parameters parameters (source positions,
More informationPiotr Majer Risk Patterns and Correlated Brain Activities
Alena My²i ková Piotr Majer Song Song Alena Myšičková Peter N. C. Mohr Peter N. C. Mohr Wolfgang K. Härdle Song Song Hauke R. Heekeren Wolfgang K. Härdle Hauke R. Heekeren C.A.S.E. Centre C.A.S.E. for
More informationDynamic Causal Models
Dynamic Causal Models V1 SPC Will Penny V1 SPC V5 V5 Olivier David, Karl Friston, Lee Harrison, Andrea Mechelli, Klaas Stephan Wellcome Department of Imaging Neuroscience, ION, UCL, UK. Mathematics in
More informationExogeneity and Causality
Università di Pavia Exogeneity and Causality Eduardo Rossi University of Pavia Factorization of the density DGP: D t (x t χ t 1, d t ; Ψ) x t represent all the variables in the economy. The econometric
More informationBayesian modelling of fmri time series
Bayesian modelling of fmri time series Pedro A. d. F. R. Højen-Sørensen, Lars K. Hansen and Carl Edward Rasmussen Department of Mathematical Modelling, Building 321 Technical University of Denmark DK-28
More informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
More informationMachine learning strategies for fmri analysis
Machine learning strategies for fmri analysis DTU Informatics Technical University of Denmark Co-workers: Morten Mørup, Kristoffer Madsen, Peter Mondrup, Daniel Jacobsen, Stephen Strother,. OUTLINE Do
More informationHeteroskedasticity in Time Series
Heteroskedasticity in Time Series Figure: Time Series of Daily NYSE Returns. 206 / 285 Key Fact 1: Stock Returns are Approximately Serially Uncorrelated Figure: Correlogram of Daily Stock Market Returns.
More informationCourse topics (tentative) The role of random effects
Course topics (tentative) random effects linear mixed models analysis of variance frequentist likelihood-based inference (MLE and REML) prediction Bayesian inference The role of random effects Rasmus Waagepetersen
More informationExtracting fmri features
Extracting fmri features PRoNTo course May 2018 Christophe Phillips, GIGA Institute, ULiège, Belgium c.phillips@uliege.be - http://www.giga.ulg.ac.be Overview Introduction Brain decoding problem Subject
More informationA longitudinal Bayesian model for spectral analysis of neuroimaging time series data. Joint work with Ning Dai
A longitudinal Bayesian model for spectral analysis of neuroimaging time series data Joint work with Ning Dai Resting-state fmri is popular for functional connectivity studies. Another analysis approach
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationJean-Baptiste Poline
Edinburgh course Avril 2010 Linear Models Contrasts Variance components Jean-Baptiste Poline Neurospin, I2BM, CEA Saclay, France Credits: Will Penny, G. Flandin, SPM course authors Outline Part I: Linear
More informationLinear Models for Regression
Linear Models for Regression Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr
More information1 Bayesian Linear Regression (BLR)
Statistical Techniques in Robotics (STR, S15) Lecture#10 (Wednesday, February 11) Lecturer: Byron Boots Gaussian Properties, Bayesian Linear Regression 1 Bayesian Linear Regression (BLR) In linear regression,
More informationNonparametric Bayesian Methods (Gaussian Processes)
[70240413 Statistical Machine Learning, Spring, 2015] Nonparametric Bayesian Methods (Gaussian Processes) Jun Zhu dcszj@mail.tsinghua.edu.cn http://bigml.cs.tsinghua.edu.cn/~jun State Key Lab of Intelligent
More informationContents. Introduction The General Linear Model. General Linear Linear Model Model. The General Linear Model, Part I. «Take home» message
DISCOS SPM course, CRC, Liège, 2009 Contents The General Linear Model, Part I Introduction The General Linear Model Data & model Design matrix Parameter estimates & interpretation Simple contrast «Take
More informationStatistical inference for MEG
Statistical inference for MEG Vladimir Litvak Wellcome Trust Centre for Neuroimaging University College London, UK MEG-UK 2014 educational day Talk aims Show main ideas of common methods Explain some of
More informationExperimental Design. Rik Henson. With thanks to: Karl Friston, Andrew Holmes
Experimental Design Rik Henson With thanks to: Karl Friston, Andrew Holmes Overview 1. A Taxonomy of Designs 2. Epoch vs Event-related 3. Mixed Epoch/Event Designs A taxonomy of design Categorical designs
More informationChapter 3 - Temporal processes
STK4150 - Intro 1 Chapter 3 - Temporal processes Odd Kolbjørnsen and Geir Storvik January 23 2017 STK4150 - Intro 2 Temporal processes Data collected over time Past, present, future, change Temporal aspect
More informationSignal Processing for Functional Brain Imaging: General Linear Model (2)
Signal Processing for Functional Brain Imaging: General Linear Model (2) Maria Giulia Preti, Dimitri Van De Ville Medical Image Processing Lab, EPFL/UniGE http://miplab.epfl.ch/teaching/micro-513/ March
More informationA Unified Estimation Framework for State- Related Changes in Effective Brain Connectivity
Universiti Teknologi Malaysia From the SelectedWorks of Chee-Ming Ting PhD 6 A Unified Estimation Framework for State- Related Changes in Effective Brain Connectivity Chee-Ming Ting, PhD S. Balqis Samdin
More informationPart 6: Multivariate Normal and Linear Models
Part 6: Multivariate Normal and Linear Models 1 Multiple measurements Up until now all of our statistical models have been univariate models models for a single measurement on each member of a sample of
More informationDynamic Causal Modelling for EEG and MEG. Stefan Kiebel
Dynamic Causal Modelling for EEG and MEG Stefan Kiebel Overview 1 M/EEG analysis 2 Dynamic Causal Modelling Motivation 3 Dynamic Causal Modelling Generative model 4 Bayesian inference 5 Applications Overview
More informationTheory of Maximum Likelihood Estimation. Konstantin Kashin
Gov 2001 Section 5: Theory of Maximum Likelihood Estimation Konstantin Kashin February 28, 2013 Outline Introduction Likelihood Examples of MLE Variance of MLE Asymptotic Properties What is Statistical
More informationarxiv: v4 [stat.ap] 7 Jul 2017
Estimating Causal Mediation Effects under Correlated Errors arxiv:40.77v4 [stat.ap] 7 Jul 07 Yi Zhao Brown University, Providence RI, USA. Xi Luo Brown University, Providence RI, USA. E-mail: xi.rossi.luo@gmail.com
More information(I AL BL 2 )z t = (I CL)ζ t, where
ECO 513 Fall 2011 MIDTERM EXAM The exam lasts 90 minutes. Answer all three questions. (1 Consider this model: x t = 1.2x t 1.62x t 2 +.2y t 1.23y t 2 + ε t.7ε t 1.9ν t 1 (1 [ εt y t = 1.4y t 1.62y t 2
More informationFunctional Connectivity and Network Methods
18/Sep/2013" Functional Connectivity and Network Methods with functional magnetic resonance imaging" Enrico Glerean (MSc), Brain & Mind Lab, BECS, Aalto University" www.glerean.com @eglerean becs.aalto.fi/bml
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
More informationSpatial Source Filtering. Outline EEG/ERP. ERPs) Event-related Potentials (ERPs( EEG
Integration of /MEG/fMRI Vince D. Calhoun, Ph.D. Director, Image Analysis & MR Research The MIND Institute Outline fmri/ data Three Approaches to integration/fusion Prediction Constraints 2 nd Level Fusion
More informationSTAT Financial Time Series
STAT 6104 - Financial Time Series Chapter 4 - Estimation in the time Domain Chun Yip Yau (CUHK) STAT 6104:Financial Time Series 1 / 46 Agenda 1 Introduction 2 Moment Estimates 3 Autoregressive Models (AR
More informationOptimization of Designs for fmri
Optimization of Designs for fmri UCLA Advanced Neuroimaging Summer School August 2, 2007 Thomas Liu, Ph.D. UCSD Center for Functional MRI Why optimize? Scans are expensive. Subjects can be difficult to
More informationExperimental design of fmri studies
Experimental design of fmri studies Sandra Iglesias Translational Neuromodeling Unit University of Zurich & ETH Zurich With many thanks for slides & images to: Klaas Enno Stephan, FIL Methods group, Christian
More informationM/EEG source analysis
Jérémie Mattout Lyon Neuroscience Research Center Will it ever happen that mathematicians will know enough about the physiology of the brain, and neurophysiologists enough of mathematical discovery, for
More informationGaussian processes for spatial modelling in environmental health: parameterizing for flexibility vs. computational efficiency
Gaussian processes for spatial modelling in environmental health: parameterizing for flexibility vs. computational efficiency Chris Paciorek March 11, 2005 Department of Biostatistics Harvard School of
More informationOverfitting, Bias / Variance Analysis
Overfitting, Bias / Variance Analysis Professor Ameet Talwalkar Professor Ameet Talwalkar CS260 Machine Learning Algorithms February 8, 207 / 40 Outline Administration 2 Review of last lecture 3 Basic
More informationLong-Run Covariability
Long-Run Covariability Ulrich K. Müller and Mark W. Watson Princeton University October 2016 Motivation Study the long-run covariability/relationship between economic variables great ratios, long-run Phillips
More informationDynamic causal modeling for fmri
Dynamic causal modeling for fmri Methods and Models for fmri, HS 2015 Jakob Heinzle Structural, functional & effective connectivity Sporns 2007, Scholarpedia anatomical/structural connectivity - presence
More informationMA 575 Linear Models: Cedric E. Ginestet, Boston University Mixed Effects Estimation, Residuals Diagnostics Week 11, Lecture 1
MA 575 Linear Models: Cedric E Ginestet, Boston University Mixed Effects Estimation, Residuals Diagnostics Week 11, Lecture 1 1 Within-group Correlation Let us recall the simple two-level hierarchical
More informationWill Penny. SPM short course for M/EEG, London 2015
SPM short course for M/EEG, London 2015 Ten Simple Rules Stephan et al. Neuroimage, 2010 Model Structure The model evidence is given by integrating out the dependence on model parameters p(y m) = p(y,
More informationROI analysis of pharmafmri data: an adaptive approach for global testing
ROI analysis of pharmafmri data: an adaptive approach for global testing Giorgos Minas, John A.D. Aston, Thomas E. Nichols and Nigel Stallard Abstract Pharmacological fmri (pharmafmri) is a new highly
More informationFinal Exam November 24, Problem-1: Consider random walk with drift plus a linear time trend: ( t
Problem-1: Consider random walk with drift plus a linear time trend: y t = c + y t 1 + δ t + ϵ t, (1) where {ϵ t } is white noise with E[ϵ 2 t ] = σ 2 >, and y is a non-stochastic initial value. (a) Show
More informationMultivariate Regression Generalized Likelihood Ratio Tests for FMRI Activation
Multivariate Regression Generalized Likelihood Ratio Tests for FMRI Activation Daniel B Rowe Division of Biostatistics Medical College of Wisconsin Technical Report 40 November 00 Division of Biostatistics
More informationUsing Estimating Equations for Spatially Correlated A
Using Estimating Equations for Spatially Correlated Areal Data December 8, 2009 Introduction GEEs Spatial Estimating Equations Implementation Simulation Conclusion Typical Problem Assess the relationship
More informationGeneral linear model: basic
General linear model: basic Introducing General Linear Model (GLM): Start with an example Proper>es of the BOLD signal Linear Time Invariant (LTI) system The hemodynamic response func>on (Briefly) Evalua>ng
More informationWhy experimenters should not randomize, and what they should do instead
Why experimenters should not randomize, and what they should do instead Maximilian Kasy Department of Economics, Harvard University Maximilian Kasy (Harvard) Experimental design 1 / 42 project STAR Introduction
More informationLinear models. Linear models are computationally convenient and remain widely used in. applied econometric research
Linear models Linear models are computationally convenient and remain widely used in applied econometric research Our main focus in these lectures will be on single equation linear models of the form y
More informationBayesian inference J. Daunizeau
Bayesian inference J. Daunizeau Brain and Spine Institute, Paris, France Wellcome Trust Centre for Neuroimaging, London, UK Overview of the talk 1 Probabilistic modelling and representation of uncertainty
More informationFirst Technical Course, European Centre for Soft Computing, Mieres, Spain. 4th July 2011
First Technical Course, European Centre for Soft Computing, Mieres, Spain. 4th July 2011 Linear Given probabilities p(a), p(b), and the joint probability p(a, B), we can write the conditional probabilities
More informationEEG/MEG Inverse Solution Driven by fmri
EEG/MEG Inverse Solution Driven by fmri Yaroslav Halchenko CS @ NJIT 1 Functional Brain Imaging EEG ElectroEncephaloGram MEG MagnetoEncephaloGram fmri Functional Magnetic Resonance Imaging others 2 Functional
More informationStatistical Analysis of fmrl Data
Statistical Analysis of fmrl Data F. Gregory Ashby The MIT Press Cambridge, Massachusetts London, England Preface xi Acronyms xv 1 Introduction 1 What Is fmri? 2 The Scanning Session 4 Experimental Design
More informationEstimation of Dynamic Regression Models
University of Pavia 2007 Estimation of Dynamic Regression Models Eduardo Rossi University of Pavia Factorization of the density DGP: D t (x t χ t 1, d t ; Ψ) x t represent all the variables in the economy.
More informationCourse content (will be adapted to the background knowledge of the class):
Biomedical Signal Processing and Signal Modeling Lucas C Parra, parra@ccny.cuny.edu Departamento the Fisica, UBA Synopsis This course introduces two fundamental concepts of signal processing: linear systems
More informationEstimation and Model Selection in Mixed Effects Models Part I. Adeline Samson 1
Estimation and Model Selection in Mixed Effects Models Part I Adeline Samson 1 1 University Paris Descartes Summer school 2009 - Lipari, Italy These slides are based on Marc Lavielle s slides Outline 1
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationStatistical Inference
Statistical Inference Jean Daunizeau Wellcome rust Centre for Neuroimaging University College London SPM Course Edinburgh, April 2010 Image time-series Spatial filter Design matrix Statistical Parametric
More informationIntroduction to MRI Acquisition
Introduction to MRI Acquisition James Meakin FMRIB Physics Group FSL Course, Bristol, September 2012 1 What are we trying to achieve? 2 What are we trying to achieve? Informed decision making: Protocols
More informationA short introduction to INLA and R-INLA
A short introduction to INLA and R-INLA Integrated Nested Laplace Approximation Thomas Opitz, BioSP, INRA Avignon Workshop: Theory and practice of INLA and SPDE November 7, 2018 2/21 Plan for this talk
More informationHierarchical Modeling for Spatio-temporal Data
Hierarchical Modeling for Spatio-temporal Data Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of
More informationA hierarchical group ICA model for assessing covariate effects on brain functional networks
A hierarchical group ICA model for assessing covariate effects on brain functional networks Ying Guo Joint work with Ran Shi Department of Biostatistics and Bioinformatics Emory University 06/03/2013 Ying
More informationIntroduction to Biomedical Engineering
Introduction to Biomedical Engineering Biosignal processing Kung-Bin Sung 6/11/2007 1 Outline Chapter 10: Biosignal processing Characteristics of biosignals Frequency domain representation and analysis
More information