1 / 40

Detecting Sparse Connectivity: MS Lesions, Cortical Thickness, and the ‘Bubbles’ Task in an fMRI Experiment

This study examines the use of the 'bubbles' task in fMRI experiments to detect spatial point data, such as MS lesions and cortical thickness. The results show that the task can accurately identify and analyze these features, but further statistical analysis is needed to determine if they are real or noise.

bare
Télécharger la présentation

Detecting Sparse Connectivity: MS Lesions, Cortical Thickness, and the ‘Bubbles’ Task in an fMRI Experiment

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Detecting Sparse Connectivity: MS Lesions, Cortical Thickness, and the ‘Bubbles’ Task in an fMRI Experiment Keith Worsley, Nicholas Chamandy,McGill Jonathan Taylor, Stanford and Université de Montréal Robert Adler, Technion Philippe Schyns, Fraser Smith, Glasgow Frédéric Gosselin, Université de Montréal Arnaud Charil, Alan Evans, Montreal Neurological Institute

  2. Three examples of spatial point data Multiple Sclerosis (MS) lesions Galaxies ‘Bubbles’ I hope to show the ‘connections’ …

  3. Astrophysics 우리의 은하

  4. Sloan Digital Sky Survey, data release 6, Aug. ‘07

  5. What is ‘bubbles’?

  6. Nature (2005)

  7. Subject is shown one of 40 faces chosen at random … Happy Sad Fearful Neutral

  8. … but face is only revealed through random ‘bubbles’ • First trial: “Sad” expression • Subject is asked the expression: “Neutral” • Response: Incorrect 75 random bubble centres Smoothed by a Gaussian ‘bubble’ What the subject sees Sad

  9. Your turn … • Trial 2 Subject response: “Fearful” CORRECT

  10. Your turn … • Trial 3 Subject response: “Happy” INCORRECT (Fearful)

  11. Your turn … • Trial 4 Subject response: “Happy” CORRECT

  12. Your turn … • Trial 5 Subject response: “Fearful” CORRECT

  13. Your turn … • Trial 6 Subject response: “Sad” CORRECT

  14. Your turn … • Trial 7 Subject response: “Happy” CORRECT

  15. Your turn … • Trial 8 Subject response: “Neutral” CORRECT

  16. Your turn … • Trial 9 Subject response: “Happy” CORRECT

  17. Your turn … • Trial 3000 Subject response: “Happy” INCORRECT (Fearful)

  18. Bubbles analysis • E.g. Fearful (3000/4=750 trials): Trial 1 + 2 + 3 + 4 + 5 + 6 + 7 + … + 750 = Sum Correct trials Thresholded at proportion of correct trials=0.68, scaled to [0,1] Use this as a bubble mask Proportion of correct bubbles =(sum correct bubbles) /(sum all bubbles)

  19. Results • Mask average face • But are these features real or just noise? • Need statistics … Happy Sad Fearful Neutral

  20. Statistical analysis • Correlate bubbles with response (correct = 1, incorrect = 0), separately for each expression • Equivalent to 2-sample Z-statistic for correct vs. incorrect bubbles, e.g. Fearful: • Very similar to the proportion of correct bubbles: Z~N(0,1) statistic Trial 1 2 3 4 5 6 7 … 750 Response 0 1 1 0 1 1 1 … 1

  21. Results • Thresholded at Z=1.64 (P=0.05) • Multiple comparisons correction? • Need random field theory … Z~N(0,1) statistic Average face Happy Sad Fearful Neutral

  22. Euler Characteristic Heuristic Euler characteristic (EC) = #blobs - #holes (in 2D) Excursion set Xt = {s: Z(s)≥t}, e.g. for neutral face: EC = 0 0 -7 -11 13 14 9 1 0 30 Heuristic: At high thresholds t, the holes disappear, EC ~ 1 or 0, E(EC) ~ P(max Z≥t). Observed Expected 20 10 EC(Xt) 0 • Exact expression for E(EC) for all thresholds, • E(EC) ~ P(max Z ≥ t) is extremely accurate. -10 -20 -4 -3 -2 -1 0 1 2 3 4 Threshold, t

  23. 2 ( ( ( ) ) ) ( ) f d ¯ l d < L I Z Z S N G i i i i 0 1 ¸ t t 2 ½ s s a n s o r o p c a u s s a n r a n o m e s » 0 0 ; , , @ Z ¡ ¢ 2 h ¸ V I i t w = 2 2 £ @ , s ( ( ) ) L Z S ¸ t ½ 1 1 µ ¶ ( ) ( ( f ( ) g ) ) P E Z E C S Z ¸ ¸ t t \ ¼ m a x s s : s ( ( ) ) L Z S ¸ t S 2 ½ s 2 2 1 Z 1 2 = 2 ¡ ( ) z d E C S £ e z = = ( ) 1 2 2 ¼ t 1 2 = 2 ¡ 1 t ( ) ¸ P S i t + £ e r m e e r e 2 2 ¼ 1 2 = 2 2 ¡ t ( ) ¸ A S t + £ r e a e = ( ) 3 2 2 ¼ ( ) f h l d I Z i i i t s s w e n o s e c o n v o v e h G i i i i t t w a n s o r o p c a u s s a n ¯ l f l l d h l f F W H i t t t e r o u a a h M F W H M i t a x m u m e n p l 4 2 o g ¸ : = F W H M The result Lipschitz-Killing curvatures of S (=Resels(S)×c) EC densities of Z above t Z(s) white noise filter = * FWHM

  24. Results, corrected for search • Random field theory threshold: Z=3.92 (P=0.05) 3.82 3.80 3.81 3.80 • Saddle-point approx (2007): Z=↑ (P=0.05) • Bonferroni: Z=4.87 (P=0.05) – nothing Z~N(0,1) statistic Average face Happy Sad Fearful Neutral

  25. Scale space: smooth Z(s) with range of filter widths w = continuous wavelet transform adds an extra dimension to the random field: Z(s,w) Scale space, no signal 34 8 22.7 6 4 15.2 2 10.2 0 -2 6.8 -60 -40 -20 0 20 40 60 w = FWHM (mm, on log scale) One 15mm signal 34 8 22.7 6 4 15.2 2 10.2 0 -2 6.8 Z(s,w) -60 -40 -20 0 20 40 60 s (mm) 15mm signal is best detected with a 15mm smoothing filter

  26. Matched Filter Theorem (= Gauss-Markov Theorem): “to best detect signal + white noise, filter should match signal” 10mm and 23mm signals 34 8 22.7 6 4 15.2 2 10.2 0 -2 6.8 -60 -40 -20 0 20 40 60 w = FWHM (mm, on log scale) Two 10mm signals 20mm apart 34 8 22.7 6 4 15.2 2 10.2 0 -2 6.8 Z(s,w) -60 -40 -20 0 20 40 60 s (mm) But if the signals are too close together they are detected as a single signal half way between them

  27. Scale space can even separate two signals at the same location! 8mm and 150mm signals at the same location 10 5 0 -60 -40 -20 0 20 40 60 170 20 76 15 34 w = FWHM (mm, on log scale) 10 15.2 5 6.8 Z(s,w) -60 -40 -20 0 20 40 60 s (mm)

  28. 2 ( ( ( ) ) ) ( ) f d ¯ l d < L I Z Z S N G i i i i 0 1 ¸ t t 2 ½ s s a n s o r o p c a u s s a n r a n o m e s » 0 0 ; , , @ Z ¡ ¢ 2 h ¸ V I i t w = 2 2 £ @ , s ( ( ) ) L Z S ¸ t ½ 1 1 µ ¶ ( ) ( ( f ( ) g ) ) P E Z E C S Z ¸ ¸ t t \ ¼ m a x s s : s ( ( ) ) L Z S ¸ t S 2 ½ s 2 2 1 Z 1 2 = 2 ¡ ( ) z d E C S £ e z = = ( ) 1 2 2 ¼ t 1 2 = 2 ¡ 1 t ( ) ¸ P S i t + £ e r m e e r e 2 2 ¼ 1 2 = 2 2 ¡ t ( ) ¸ A S t + £ r e a e = ( ) 3 2 2 ¼ ( ) f h l d I Z i i i t s s w e n o s e c o n v o v e h G i i i i t t w a n s o r o p c a u s s a n ¯ l f l l d h l f F W H i t t t e r o u a a h M F W H M i t a x m u m e n p l 4 2 o g ¸ : = F W H M The result Lipschitz-Killing curvatures of S (=Resels(S)×c) EC densities of Z above t Z(s) white noise filter = * FWHM

  29. ( ) l h b h l l l F L K i i i i i t t t o r s c a e - s p a c e s e a r c n g e w e e n p s c z - n g c u r v a u r e s r e s e s ? ? ? ? ? ? ? ? ? L L L L L L ¸ 0 1 2 0 1 2 ; ; ; ; ; ( ) d d d d l h l l l b L K i i i i i i t t t a a n e x r a m e n s o n a n r e p a c e p s c z - n g c u r v a u r e s r e s e s y ? L L = 0 0 ? ? ? ? ? ? ? L L L L L + ¡ 1 1 1 2 2 ? l L L + + o g = 1 0 L ? ? 2 4 ¼ 1 ? ? ? L L + 2 2 ? ? ? L L L + ¡ = 2 1 1 2 ? ? ? L L ¡ 2 2 L = 3 2 Random field theory for scale-space

  30. Unsmoothed fMRI data: T stat for visual stimulus Rotation space:Try all rotated elliptical filters Threshold Z=5.25 (P=0.05) Maximum filter

  31. 2 ( ( ) ) ( ) ( ) f d ¯ l d < L I Z S N G i i i i 0 1 Z t ¸ t 2 s s a n s o r o p c a u s s a n r a n o m e s ½ » 0 0 ; , , @ Z ¡ ¢ 2 h ¸ V I i t w = 2 2 £ @ , s ( ) L ( S ) Z ¸ t ½ 1 1 µ ¶ ( ) ( ( f ( ) g ) ) P E Z E C S Z ¸ ¸ t t \ ¼ m a x s s : s ( ) L S ( ) Z ¸ S t 2 s ½ 2 2 1 Z 1 2 = 2 ¡ ( ) z d E C S £ e z = = ( ) 1 2 2 ¼ t 1 2 = 2 ¡ 1 t ( ) ¸ P S i t + £ e r m e e r e 2 2 ¼ 1 2 = 2 2 ¡ t ( ) ¸ A S t + £ r e a e = ( ) 3 2 2 ¼ ( ) f h l d I Z i i i t s s w e n o s e c o n v o v e h G i i i i t t w a n s o r o p c a u s s a n ¯ l f l l d h l f F W H i t t t e r o u a a h M F W H M i t a x m u m e n p l 4 2 o g ¸ : = F W H M The result Lipschitz-Killing curvatures of S (=Resels(S)×c) EC densities of Z above t Z(s) white noise filter = * FWHM

  32. D ( ) < b h d ¯ l d L T S i i t t t 2 ½ e s s e a s m o o s o r o p c r a n o m e , , . f ( ) g b h h d L X T S i i i ¸ t t t t t e s : s e e e e x c u r s o n s e n s e = t . h T e n D X ( ( ) ) ( ) ( ) E E C S X L S T ¸ t \ ½ = d d t : d 0 = ( ) ( ( ) ) h f f f d d d d N T Z i i i i i t t t t t o w s u p p o s e a s s s a u n c o n o n e p e n e n a n e n - = ( ) ( ( ) ( ) ) l l d b d d ¯ l d h h G Z Z Z i i i i t t t c a y s r u e a u s s a n r a n o m e s s s s e a c w = 1 n ; : : : ; , ( ) @ Z 2 s ( ) ( ) ( ) i d ¸ V Z N I 0 1 s a n » = D D i £ @ ; . s f ( ) g f b h f L R T j i i ¸ t t t t e z : z e e r e e c o n r e g o n o = t . Theorem (1981, 1995) Example: the chi-bar random field, a special case of a random field of test statistics for the magnitude of an fMRI response in the presence of unknown delay of the hemodynamic response function.

  33. f f g g µ µ R X Z Z Z i ¸ ¸ t t + ¹ ¹ ¹  m a x s c : : o   s s n = = = 1 2 t t = µ · · 0 2 ¼ Example: chi-bar random field Z1~N(0,1) Z2~N(0,1) s2 s1 Excursion sets, Rejection regions, Threshold t Z2 Search Region, S Z1

  34. µ ¶ @ Z ¸ d S = @ s ( ) ( ) d E C R D h l l L i K L S i i i i t t t e n s y ½ p s c z - n g c u r v a u r e D d d t d = 1 X d 2 ( ) 2 ¼ ( ( ) ) ( ) ( ) X E E C S X L S R ¼ \ d X ½ j ( ) j ( ) = d d d b ¸ T S L S t t ( ( ) ) ( ) b P T R R u e r r = d D ¡ u e r ½ r = d ; t t ( = ) d ¡ 2 1 + ; d ! d 0 = d 0 = d 0 = Adler & Taylor (2007), Ann. Math, (submitted) Beautiful symmetry: Steiner-Weyl Tube Formula (1930) Taylor Gaussian Tube Formula (2003) • Put a tube of radius r about the search region λS and rejection region Rt: Z2~N(0,1) Rt r Tube(λS,r) Tube(Rt,r) r λS Z1~N(0,1) t t-r • Find volume or probability, expand as a power series in r, pull off coefficients:

  35. Bubbles task in fMRI scanner • Correlate bubbles with BOLD at every voxel: • Calculate Z for each pair (bubble pixel, fMRI voxel) • a 5D “image” of Z statistics … Trial 1 2 3 4 5 6 7 … 3000 fMRI

  36. µ ¶ ( ) ( f ( ) g ) P E C E C S T C ¸ ¸ t t t 2 2 ¼ m a x s c s : s c ; ; ; S T t 2 2 s ; ( ) ( ) d d S T i i m m X X ( ) ( ) ( ) L S L T C ¸ ½ c = i j i j i j 0 0 = = b ( ) = c h 1 2 ¡ k k h 2 ¡ ¡ ( ) ! ! n i j 2 1 ¡ X X X ( ) = k h k h k 1 2 2 1 2 ¡ ¡ ¡ ¡ + ( ) ( ) ( ) n C 1 1 ¸ ¡ ¡ ½ c c c = i j = h 2 1 + ¼ k l 0 0 0 m = = = j ¡ i ¡ ( ) ( n ) n l ¡ ¡ + + m 2 2 ( ) ( ) ( ) ( ) l ! ! k l ! h l k ! k l ! k l ! i j 1 1 ¡ ¡ ¡ ¡ + + + ¡ ¡ ¡ + ¡ ¡ + m m n m m m Thresholding? Cross correlation random field • Correlation between 2 fields at 2 different locations, searched over all pairs of locations, one in S, one in T: • Bubbles data: P=0.05, n=3000, c=0.113, T=6.22 Cao & Worsley, Annals of Applied Probability (1999)

  37. MS lesions and cortical thickness • Idea: MS lesions interrupt neuronal signals, causing thinning in down-stream cortex • Data: n = 425 mild MS patients 5.5 5 4.5 4 Average cortical thickness (mm) 3.5 3 2.5 Correlation = -0.568, T = -14.20 (423 df) 2 Charil et al, NeuroImage (2007) 1.5 0 10 20 30 40 50 60 70 80 Total lesion volume (cc)

  38. MS lesions and cortical thickness at all pairs of points • Dominated by total lesions and average cortical thickness, so remove these effects as follows: • CT = cortical thickness, smoothed 20mm • ACT = average cortical thickness • LD = lesion density, smoothed 10mm • TLV = total lesion volume • Find partial correlation(LD, CT-ACT) removing TLV via linear model: • CT-ACT ~ 1 + TLV + LD • test for LD • Repeat for all voxels in 3D, nodes in 2D • ~1 billion correlations, so thresholding essential! • Look for high negative correlations … • Threshold: P=0.05, c=0.300, T=6.48

  39. Cluster extent rather than peak height (Friston, 1994) • Choose a lower level, e.g. t=3.11 (P=0.001) • Find clusters i.e. connected components of excursion set • Measure cluster extent by resels • Distribution: • fit a quadratic to the peak: • Distribution of maximum cluster extent: • Find distribution for a single cluster • Bonferroni on N = #clusters ~ E(EC). Z D=1 extent t Peak height s Cao and Worsley, Advances in Applied Probability (1999)

  40. Three examples of spatial point data Multiple Sclerosis (MS) lesions Galaxies ‘Bubbles’ I hope I have shown the connections ... 너를 감사하십시요 !

More Related