connectivity of amri and fmri data
DESCRIPTION
Connectivity of aMRI and fMRI data. Keith Worsley Arnaud Charil Jason Lerch Francesco Tomaiuolo Department of Mathematics and Statistics, McConnell Brain Imaging Centre, Montreal Neurological Institute, McGill University. Effective connectivity. - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/1.jpg)
Connectivity of aMRI and fMRI data
Keith WorsleyArnaud CharilJason Lerch
Francesco Tomaiuolo
Department of Mathematics and Statistics, McConnell Brain Imaging Centre, Montreal Neurological Institute,
McGill University
![Page 2: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/2.jpg)
Effective connectivity• Measured by the correlation between residuals at
pairs of voxels:
Voxel 2
Voxel 1
++ +++ +
Activation onlyVoxel 2
Voxel 1++
+
+
+
+
Correlation only
![Page 3: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/3.jpg)
Types of connectivity
• Focal
• Extensive
![Page 4: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/4.jpg)
-3
-2
-1
0
1
2
3
8
4
0
9
5
1
10
6
2
11
7
3
cor=0.58
Focal correlation
n = 120frames
-2 0 2-3
-2
-1
0
1
2
3
![Page 5: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/5.jpg)
-3
-2
-1
0
1
2
3
8
4
0
9
5
1
10
6
2
11
7
3
-2 0 2-3
-2
-1
0
1
2
3
cor=0.13
Extensive correlation
![Page 6: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/6.jpg)
Methods
1. Seed
2. Iterated seed
3. Thresholding correlations
4. PCA
![Page 7: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/7.jpg)
Method 1: ‘Seed’
Friston et al. (19??): Pick one voxel, then find all others that are correlated with it:
Problem: how to pick the ‘seed’ voxel?
![Page 8: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/8.jpg)
Method 2: Iterated ‘seed’
• Problem: how to find the rest of the connectivity network?
• Hampson et al., (2002): Find significant correlations, use them as new seeds, iterate.
![Page 9: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/9.jpg)
Method 3: All correlations
• Problem: how to find isolated parts of the connectivity network?
• Cao & Worsley (1998): find all correlations (!)
• 6D data, need higher threshold to compensate
![Page 10: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/10.jpg)
Thresholds are not as high as you might think:
E.g. 1000cc search region, 10mm smoothing, 100 df, P=0.05:
dimensions D1 D2 Cor T
Voxel1 - Voxel2 0 0 0.165 1.66
One seed voxel - volume 0 3 0.448 4.99
Volume – volume (auto-correlation) 3 3 0.609 7.64
Volume1 – volume2 (cross-correlation) 3 3 0.617 7.81
![Page 11: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/11.jpg)
Practical details
• Find threshold first, then keep only correlations > threshold
• Then keep only local maxima i.e.cor(voxel1, voxel2)
> cor(voxel1, 6 neighbours of voxel2),
> cor(6 neighbours of voxel1, voxel2),
![Page 12: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/12.jpg)
Method 4: Principal Components Analysis (PCA)
• Friston et al: (1991): find spatial and temporal components that capture as much as possible of the variability of the data.
• Singular Value Decomposition of time x space matrix:
Y = U D V’ (U’U = I, V’V = I, D = diag)
• Regions with high score on a spatial component (column of V) are correlated or ‘connected’
![Page 13: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/13.jpg)
Which is better:
thresholding correlations,
or
PCA?
![Page 14: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/14.jpg)
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
8
4
0
9
5
1
10
6
2
11
7
3
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
8
4
0
9
5
1
10
6
2
11
7
3
-6
-4
-2
0
2
4
6
8
4
0
9
5
1
10
6
2
11
7
3
-6
-4
-2
0
2
4
6
8
4
0
9
5
1
10
6
2
11
7
3
Summary
Extensive correlationFocal correlation
Thresholding T statistic
(=correlations)
PCA
![Page 15: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/15.jpg)
0
500
1000First scan of fMRI data
-5
0
5
T statistic for hot - warm effect
0 100 200 300
870880890 hot
restwarm
Highly significant effect, T=6.59
0 100 200 300
800
820hotrestwarm
No significant effect, T=-0.74
0 100 200 300
790800810
Drift
Time, seconds
fMRI data: 120 scans, 3 scans each of hot, rest, warm, rest, hot, rest, …
T = (hot – warm effect) / S.d. ~ t110 if no effect
![Page 16: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/16.jpg)
0 20 40 60 80 100 1205
4
3
2
1
0
Co
mp
on
en
t
Temporal components (sd, % variance explained)
0.68, 46.9%
0.29, 8.6%
0.17, 2.9%
0.15, 2.4%
-1
-0.5
0
0.5
1
Slice (0 based)
Co
mp
on
en
t
Spatial components
0 2 4 6 8 10 12
1
2
3
4
PCA of time space:
1: excludefirst frames
2: drift
3: long-range correlationor anatomicaleffect: removeby converting to % of brain
4: signal?
Frame
![Page 17: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/17.jpg)
MS lesions and cortical thickness(Arnaud et al., 2004)
• n = 425 mild MS patients
• Lesion density, smoothed 10mm
• Cortical thickness, smoothed 20mm
• Find connectivity i.e. find voxels in 3D, nodes in 2D with high
cor(lesion density, cortical thickness)
![Page 18: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/18.jpg)
0 10 20 30 40 50 60 70 801.5
2
2.5
3
3.5
4
4.5
5
5.5
Average lesion volume
Ave
rag
e co
rtic
al t
hic
kne
ss
n=425 subjects, correlation = -0.56826
![Page 19: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/19.jpg)
Normalization
• Simple correlation:
Cor( LD, CT )
• Subtracting global mean thickness:
Cor( LD, CT – avsurf(CT) )
• And removing overall lesion effect:
Cor( LD – avWM(LD), CT – avsurf(CT) )
![Page 20: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/20.jpg)
0
0.5
1
1.5
2
2.5x 10
5
distance (mm)
corr
elat
ion
Same hemisphere
0 50 100 150-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
0
0.2
0.4
0.6
0.8
1
distance (mm)
corr
elat
ion
Correlation = 0.091943
0 50 100 150-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
0
0.5
1
1.5
2
2.5
x 105
distance (mm)
corr
elat
ion
Different hemisphere
0 50 100 150-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
0
0.2
0.4
0.6
0.8
1
distance (mm)
corr
elat
ion
Correlation = -0.1257
0 50 100 150-0.5
-0.4
-0.3
-0.2
-0.1
0
0.1
threshold
thresholdthreshold
threshold
![Page 21: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/21.jpg)
Deformation Based Morphometry (DBM) (Tomaiuolo et al., 2004)
• n1 = 19 non-missile brain trauma patients, 3-14 days in coma,
• n2 = 17 age and gender matched controls
• Data: non-linear vector deformations needed to warp each MRI to an atlas standard
• Locate damage: find regions where deformations are different, hence shape change
• Is damage connected? Find pairs of regions with high canonical correlation.
![Page 22: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/22.jpg)
-6
-4
-2
0
2
4
6
8
4
0
9
5
1
10
6
2
11
7
3Seed
T = sqrt(df) cor / sqrt (1 - cor2)
T max = 7.81P=0.00000004
![Page 23: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/23.jpg)
PCA, component 1
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
8
4
0
9
5
1
10
6
2
11
7
3
![Page 24: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/24.jpg)
-6
-4
-2
0
2
4
6
8
4
0
9
5
1
10
6
2
11
7
3Seed
T max = 4.17P = 0.59
T, extensive correlation
![Page 25: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/25.jpg)
PCA, focal correlation
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
8
4
0
9
5
1
10
6
2
11
7
3
![Page 26: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/26.jpg)
Modulated connectivity
• Looking for correlations not very interesting – ‘resting state networks’
• More intersting: how does connectivity change with- task or condition (external)- response at another voxel (internal)
• Friston et al., (1995): add interaction to the linear model:
Data ~ task + seed + task*seed Data ~ seed1 + seed2 + seed1*seed2
![Page 27: Connectivity of aMRI and fMRI data](https://reader036.vdocument.in/reader036/viewer/2022062421/56812b12550346895d8f05b0/html5/thumbnails/27.jpg)
Fit a linear model for fMRI time series with AR(p) errors
• Linear model: ? ? Yt = (stimulust * HRF) b + driftt c + errort
• AR(p) errors: ? ? ? errort = a1 errort-1 + … + ap errort-p + s WNt
• Subtract linear model to get residuals.• Look for connectivity.
unknown parameters