Transcript Coupling measures derived from Multivariate Autoregressive

COST B27 ENOC Joint WGs Meeting
Swansea UK, 16-18 September 2006
Analyzing EEG Coupling
Alois Schlögl
University of Technology Graz, Austria
<[email protected]>
1
• Offer:
– Methods for analysing EEG coupling using
Multivariate autoregressive modelling
– Coherency, PDC, DTF, phase, etc.
– BioSig http://biosig.sf.net
• Asking for:
– EEG data on interesting research topics
– collaboration
2
Multivariate AutoRegressive
(MVAR) models
• Multivariate: spatio-temporal correlation
• Estimators:
– Levinson-Wiggins-Robinson (LWR):
Multivariate Yule-Walker
– Nutall-Strand method (multivariate burg method)
– Vieira-Morf (multivariate geometric lattice)
• Software:
– TSA-toolbox: MVAR.M
– BioSig http://biosig.sf.net
Schlögl (2006), Comparison of MVAR estimators, Signal Processing.
3
4
Time-Frequency Analysis
PDC (Hypothesis: PDC>0)
Subject K3
Left hand
1
0
Coupling almost all the time in all frequencies !?
5
Event-related PDC
(Hypothesis: PDC != ref)
Subject K3
Left hand
ref=pdc(0-3s)
pdc>ref
pdc<ref
Increases and decreases of coupling can be observed ! 6
7
MVAR estimators
N=400
5
N=100
Comparison of Multivariate
Autoregressive Estimators.
Signal processing
N=50
N=60
ARFIT
multichannel Yule-Walker
4
ARFIT-0
3
Nuttall-Strand with biased covariance
2
Schlögl, A. (2006)
N=70
N=40
ranking
Comparison of
• ARFIT
• Multichannel YuleWalker
• Multichannel Burg
(Nutall-Strand)
Nuttall-Strand with unbiased covariance
1
1
2
10
5
mean square predicition error
20
8
Signal processing (II)
•
•
•
•
•
•
•
•
Spectrum
Coherence (absCOH, imagCOH)
Partial coherence (pCOH)
Partial directed coherence (PDC)
Directed transfer function (DTF)
Full-Frequency DTF (ffDTF)
Directed coupling (DC)
Phase
9
Auto/Crossspectra & Coherence
COH
logS
20
30
40
50
0
0
10
20
30
40
50
0
10
20
30
40
50
0
10
20
30
40
50
0
10
10
20
30
40
50
0
0
10
20
30
40
10
20
30
40
50
0
10
20
30
40
50
50
0
10
20
30
40
10
20
30
40
50
10
20
30
40
20
30
40
50
30
40
50
10
20
30
40
50
20
30
40
10
20
30
40
20
30
40
50
40
10
20
30
40
0
10
20
30
40
0
0
10
20
30
40
50
0
10
20
30
40
0
10
0
10
20
30
40
50
0
10
20
30
40
50
#05
1
1
0.5
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
0
50 0
0
50 0
0
50 0
0.5
0.5
0.5
0
0
1
50
#04
1
0
0
1
50
#03
1
0
0
1
50
#02
1
0
0
1
50
10
0
10
0
50
10
0
30
0
0
50
20
10
0
10
10
0
10
0
0
10
0
50
0
10
20
0
10
0
10
10
0
50
0
0
0
10
0
40
10
0
10
30
10
10
0
50
10
0
20
10
0
10
10
0
10
0
0
10
0
0
10
0
10
10
0
#01
0
10
#01
#05
#02
#02
#03
10
#04
#03
0
10
0
#04
#03
#04
0
10
0
#05
#02
#05
#01
#01
0.5
0
0
50
50
50
50
50
10
Imaginare und partiel Coherence
iCOH
1
1
0
0
0
0
0
0
0
0
-1
0
1
0
-1
50 0
1
0
-1
0
1
0
-1
-1
50 0
1
-1
50 0
1
0
0
-1
50 0
-1
50 0
1
0
-1
50 0
1
0
-1
50 0
1
0
-1
50 0
1
0
-1
50 0
-1
50 0
1
0
-1
50 0
1
0
-1
50 0
1
0
-1
50 0
1
0
-1
50 0
-1
50 0
1
0
-1
50 0
1
0
-1
50 0
1
0
-1
50 0
1
0
-1
50 0
50
#05
1
1
0.5
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
0
50 0
0
50 0
0
50 0
0.5
0.5
0.5
0
0
1
50
#04
1
0
0
1
50
#03
1
0
0
1
50
#02
1
0
0
1
50
#03
0
-1
50 0
1
#02
1
#01
#01
#05
#04
#02
#03
#04
1
-1
0
1
#04
pCOH
#03
1
-1
0
1
#05
#02
#05
#01
#01
0.5
0
0
50
50
50
50
50
11
Directed transfer function DTF &
partial directed coherence (PDC)
DTF
1
0.5
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
0
50 0
0
50 0
0
50 0
0.5
0.5
0
0
1
0.5
0
0
1
0.5
0
0
50
#05
1
1
0.5
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
1
0
50 0
1
0
50 0
1
0
50 0
1
0.5
0.5
0.5
0.5
0
50 0
0
50 0
0
50 0
0
50 0
0.5
0.5
0.5
0
0
1
50
#04
1
0
0
1
50
#03
1
0
0
1
50
#02
1
0
0
1
50
#02
1
#03
1
#01
#01
#05
#04
#02
#03
#04
1
0
0
1
#04
PDC
#03
1
0
0
1
#05
#02
#05
#01
#01
0.5
0
0
50
50
50
50
50
12
Thank You for the Attention
• References:
– Schlögl, 2006, Comparison of Multivariate
Autoregressive Estimators.
Signal processing
– Schlögl and Supp (in press), Progress In Brain
Research
Contact:
Alois Schlögl
[email protected]
http://biosig.sf.net
13
DTF (Hypothesis: DTF>0)
Subject K3
Left hand
1
0
Coupling almost all the time in all frequencies !?
14
Event-related PDC
(Hypothesis: PDC != ref)
Subject K3
Right hand
ref=pdc(0-3s)
pdc>ref
pdc<ref
Increases and decreases of coupling can be observed !15