??? Dynamic Causal Modelling for M/EEG. Electroencephalography (EEG) Dynamic Causal Modelling. M/EEG analysis at sensor level. time.

Similar documents
Dynamic Causal Modelling for EEG and MEG

ST 524 NCSU - Fall 2008 One way Analysis of variance Variances not homogeneous

Analyzing Frequencies

2. Grundlegende Verfahren zur Übertragung digitaler Signale (Zusammenfassung) Informationstechnik Universität Ulm

Soft k-means Clustering. Comp 135 Machine Learning Computer Science Tufts University. Mixture Models. Mixture of Normals in 1D

Outlier-tolerant parameter estimation

The impact of complex network structure on seizure activity induced by depolarization block

The Fourier Transform

Lecture 14. Relic neutrinos Temperature at neutrino decoupling and today Effective degeneracy factor Neutrino mass limits Saha equation

Dynamic Causal Modelling for EEG and MEG. Stefan Kiebel

Grand Canonical Ensemble

Physics of Very High Frequency (VHF) Capacitively Coupled Plasma Discharges

10/7/14. Mixture Models. Comp 135 Introduction to Machine Learning and Data Mining. Maximum likelihood estimation. Mixture of Normals in 1D

Economics 600: August, 2007 Dynamic Part: Problem Set 5. Problems on Differential Equations and Continuous Time Optimization

Review - Probabilistic Classification

A Probabilistic Characterization of Simulation Model Uncertainties

Chapter 6 Student Lecture Notes 6-1

Radial Cataphoresis in Hg-Ar Fluorescent Lamp Discharges at High Power Density

:2;$-$(01*%<*=,-./-*=0;"%/;"-*

Determination of Vibrational and Electronic Parameters From an Electronic Spectrum of I 2 and a Birge-Sponer Plot

Outline. Why speech processing? Speech signal processing. Advanced Multimedia Signal Processing #5:Speech Signal Processing 2 -Processing-

Basic Electrical Engineering for Welding [ ] --- Introduction ---

EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Lecture 12

An Overview of Markov Random Field and Application to Texture Segmentation

Authentication Transmission Overhead Between Entities in Mobile Networks

Electrochemical Equilibrium Electromotive Force. Relation between chemical and electric driving forces

Lecture 3: Phasor notation, Transfer Functions. Context

Phys 774: Nonlinear Spectroscopy: SHG and Raman Scattering

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d)

~ Outline ~ 1. Normalization Quality filtering Rank-Invariant selection for non-differentially expressed genes Fit normalization curve

External Equivalent. EE 521 Analysis of Power Systems. Chen-Ching Liu, Boeing Distinguished Professor Washington State University

The Hyperelastic material is examined in this section.

CHAPTER 33: PARTICLE PHYSICS

MUSIC Based on Uniform Circular Array and Its Direction Finding Efficiency

Introduction to the Fourier transform. Computer Vision & Digital Image Processing. The Fourier transform (continued) The Fourier transform (continued)

UNIT 8 TWO-WAY ANOVA WITH m OBSERVATIONS PER CELL

TOPOLOGY DESIGN OF STRUCTURE LOADED BY EARTHQUAKE. Vienna University of Technology

Random Process Part 1

Today s logistic regression topics. Lecture 15: Effect modification, and confounding in logistic regression. Variables. Example

Ch. 9 Common Emitter Amplifier

EE 570: Location and Navigation: Theory & Practice

Lecture 28 Title: Diatomic Molecule : Vibrational and Rotational spectra

MECH321 Dynamics of Engineering System Week 4 (Chapter 6)

1- Summary of Kinetic Theory of Gases

Studies of Turbulence and Transport in Alcator C-Mod Ohmic Plasmas with Phase Contrast Imaging and Comparisons with GYRO*

Slide 1. Slide 2. Slide 3 DIGITAL SIGNAL PROCESSING CLASSIFICATION OF SIGNALS

Status of LAr TPC R&D (2) 2014/Dec./23 Neutrino frontier workshop 2014 Ryosuke Sasaki (Iwate U.)

GPC From PeakSimple Data Acquisition

Modeling and Simulation Analysis of Power Frequency Electric Field of UHV AC Transmission Line

A Note on Estimability in Linear Models

Logistic Regression I. HRP 261 2/10/ am

Learning Spherical Convolution for Fast Features from 360 Imagery

Econ107 Applied Econometrics Topic 10: Dummy Dependent Variable (Studenmund, Chapter 13)

ACOUSTIC WAVE EQUATION. Contents INTRODUCTION BULK MODULUS AND LAMÉ S PARAMETERS

Heisenberg Model. Sayed Mohammad Mahdi Sadrnezhaad. Supervisor: Prof. Abdollah Langari

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let

Performance assessment of the window-wall interface. PhD-meeting Nathan Van Den Bossche - 04/06/2010 Department of Architecture Ghent University

Stress-Based Finite Element Methods for Dynamics Analysis of Euler-Bernoulli Beams with Various Boundary Conditions

Filter Design Techniques

Chapter 1 Late 1800 s Several failures of classical (Newtonian) physics discovered

Optimal Topology Design for Replaceable of Reticulated Shell Based on Sensitivity Analysis

Chapter 13 Aggregate Supply

Finite Element Based Implementation of Fiala s Thermal Manikin in THESEUS-FE

Study of Dynamic Aperture for PETRA III Ring K. Balewski, W. Brefeld, W. Decking, Y. Li DESY

2. Laser physics - basics

Folding of Regular CW-Complexes

Chapter 3 Lecture 14 Longitudinal stick free static stability and control 3 Topics

COMPLEX NUMBER PAIRWISE COMPARISON AND COMPLEX NUMBER AHP

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

IMGS-261 Solutions to Homework #9

Lecture 2: Discrete-Time Signals & Systems. Reza Mohammadkhani, Digital Signal Processing, 2015 University of Kurdistan eng.uok.ac.

FEFF and Related Codes

Modelling of new generation plasma optical devices

Full Waveform Inversion Using an Energy-Based Objective Function with Efficient Calculation of the Gradient

Decision-making with Distance-based Operators in Fuzzy Logic Control

Answer Homework 5 PHA5127 Fall 1999 Jeff Stark

Characteristics of Gliding Arc Discharge Plasma

Motion Perception Under Uncertainty. Hongjing Lu Department of Psychology University of Hong Kong

SENSOR networks are wireless ad hoc networks used for. Minimum Energy Fault Tolerant Sensor Networks

A Study on Nonlinear Forced Vibration of an Axial Moving Viscoelasticity Beam using the Multi-Scale Approach


The general linear model for fmri

Lecture 21. Boltzmann Statistics (Ch. 6)

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)

Chp6. pn Junction Diode: I-V Characteristics II

CHAPTER 4. The First Law of Thermodynamics for Control Volumes

Problem Set #2 Due: Friday April 20, 2018 at 5 PM.

Introduction to Condensed Matter Physics

te Finance (4th Edition), July 2017.

Using Markov Chain Monte Carlo for Modeling Correct Enumeration and Match Rate Variability

Abstract Interpretation: concrete and abstract semantics

Physics 178/278 - David Kleinfeld - Fall checked Winter 2014

VII. Quantum Entanglement

From Structural Analysis to FEM. Dhiman Basu

Lecture Outline. Skin Depth Power Flow 8/7/2018. EE 4347 Applied Electromagnetics. Topic 3e

Gamma-ray burst spectral evolution in the internal shock model

IYPT 2000 Problem No. 3 PLASMA

NON-SYMMETRY POWER IN THREE-PHASE SYSTEMS

167 T componnt oftforc on atom B can b drvd as: F B =, E =,K (, ) (.2) wr w av usd 2 = ( ) =2 (.3) T scond drvatv: 2 E = K (, ) = K (1, ) + 3 (.4).2.2

Naresuan University Journal: Science and Technology 2018; (26)1

Transcription:

Elctroncphalography EEG Dynamc Causal Modllng for M/EEG ampltud μv tm ms tral typ 1 tm channls channls tral typ 2 C. Phllps, Cntr d Rchrchs du Cyclotron, ULg, Blgum Basd on slds from: S. Kbl M/EEG analyss at snsor lvl Dynamc Causal Modllng tm channls tral typ 1 Approach: Rduc vokd rspons to a fw varabls,.g.: Th avrag ovr a fw channls n pr-stmulus tm.??? Buld a modl for spatotmporal data: Assum that both ERPs ar gnratd by tmporal dynamcs of a ntwork of a fw sourcs A2 channls tral typ 2 Dffrnt approach that tlls us mor about th nuronal dynamcs of localzd bran sourcs? Dynamc Causal Modllng Dscrb tmporal dynamcs by dffrntal quatons Each sourc projcts to th snsors, followng physcal laws Solv for th modl paramtrs usng Baysan modl nvrson x & = f x, u, θ p θ y, m p y m

Msmatch ngatvty MMN mod 1 Modl for msmatch ngatvty Oddball paradgm standards dvants mod 2 psudo-random audtory squnc 8% standard tons Hz 2% dvant tons 5 Hz tm prprocssng mod 3 raw data 128 EEG scalp lctrods convrt to matlab fl fltr poch down sampl artfact corrcton avrag ERPs / ERFs data rducton to prncpal spatal mods xplanng most of th varanc tm ms Garrdo t al., PNAS, 8 Th gnratv modl Nural mass quatons and connctvty Sourc dynamcs f Spatal forward modl g y = g x, θ x & = f x, u, θ stats x paramtrs θ Evokd rspons data y Input u Davd t al., NuroImag, 5.

Nural mass quatons and connctvty Nural mass quatons and connctvty nhbtory ntrnurons spny stllat clls Stat quatons x & = f x,u,θ Extrnsc forward connctons A F S x 7 = x8 H B L 8 = A + A + 2x8 x7 3I S x 2 4 3 1 = x4 H F L 2x4 x1 4 = A + A + 1I S x + Cu 2 Intrnsc 1 connctons 2 Extrnsc latral A L S x connctons pyramdal clls x = x5 x6 2 = x5 H B L 5 = A + A S x + 3 = x6 H 2x6 x3 6 = 4S x7 2 2x5 x2 2S x1 2 Extrnsc backward A B S x connctons x man trans-mmbran potntals and currnts of th dffrnt nuron populatons. nuronal sourc modl Spatal modl Baysan modl nvrson x L L θ Dpolarsaton of pyramdal clls Spatal modl Snsor data y Masurd data Expctaton-Maxmzaton algorthm Spcfy gnratv forward modl wth pror dstrbutons of paramtrs Itratv procdur: 1. Comput modl rspons usng currnt st of paramtrs 2. Compar modl rspons wth data 3. Improv paramtrs, f possbl 1. Postror dstrbutons of paramtrs p θ y, m 2. Modl vdnc p y m

Modl comparson: Whch modl s th bst? MMN xampl IFG IFG IFG data y Modl 1 Modl 2 p y m1 p θ y, m 1 p y m 2 p θ y, m 2 Modl comparson: Slct modl wth hghst modl vdnc Forward - F Backward - B Forward and Backward - FB bst?... Modl n p y mn p θ y, m n p y m Forward Backward Latral Forward Backward Latral Forward Backward Latral modulaton of ffctv connctvty Group modl comparson Evokd and nducd rsponss Baysan Modl Comparson Group lvl log-vdnc Forward F subjcts Backward B Forward and Backward FB Garrdo t al., 7, NuroImag Trnds Cogn Sc. 1999 Apr;34:151-162

Modllng of nducd rsponss Fac data EEG: Ntwork of four sourcs Invrson of lctromagntc modl L x = 36mm y = 33mm z = 8mm LF RF x = 36mm y = 31mm z = 18mm Tm-srs data n channl spac Dynamc powr data n sourc spac Am: Explan dynamc powr spctrum of ach sourc as functon of powr from othr sourcs. x = 4mm y = 8mm z = 8mm LV RV x = 4mm y = 8mm z = 1mm Chn t al., Nuromag, 8 Obsrvd powr spctra Sngl subjct rsults: Couplng functons LV RV LF RF tral 1: sourc 1 tral 1: sourc 2 tral 1: sourc 3 tral 1: sourc 4 LF RF α β α RV - 1-1 - 1-1 Tm ms frquncy frquncy frquncy frquncy Frquncy Hz LV RV RF Chn t al., Nuromag, 8

Obsrvd and fttd powr spctra Summary LV RV LF RF tral 1: sourc 1 tral 1: sourc 2 tral 1: sourc 3 tral 1: sourc 4-1 - 1-1 - 1 Tm ms DCM combns stat-quatons for nural mass dynamcs wth spatal forward modl. Dffrncs btwn rsponss acqurd undr dffrnt condtons ar modlld as modulaton of connctvty wthn and btwn sourcs. frquncy frquncy frquncy frquncy Frquncy Hz Baysan modl comparson allows on to compar many dffrnt modls and dntfy th bst on. - tral 1: sourc 1 prdctd - tral 1: sourc 2 prdctd - tral 1: sourc 3 prdctd - tral 1: sourc 4 prdctd Mak nfrnc about postror dstrbuton of paramtrs.g., ffctv connctvty, locaton of dpols, tc.. fttd 1 1 1 1 Many xtnsons to DCM for M/EEG ar avalabl n SPM8. channl channl channl channl Thanks to Karl Frston Marta Garrdo CC Chn Jan Daunzau