1 / 38

ANALYSIS OF fMRI DATA BASED ON NN-ARx MODELING

ANALYSIS OF fMRI DATA BASED ON NN-ARx MODELING. Biscay-Lirio, R : Inst. of Cybernetics, Mathematics and Physics, Cuba Bosch-Bayard, J. : Cuban Neuroscience Center, Havana, Cuba Riera-Diaz, J. : NICHe, Tohoku University, Sendai, Japan

efuru
Download Presentation

ANALYSIS OF fMRI DATA BASED ON NN-ARx MODELING

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. ANALYSIS OF fMRI DATABASED ON NN-ARx MODELING Biscay-Lirio, R: Inst. of Cybernetics, Mathematics and Physics, Cuba Bosch-Bayard, J.: Cuban Neuroscience Center, Havana, Cuba Riera-Diaz, J.: NICHe, Tohoku University, Sendai, Japan Biscay-Lirio,R.: Inst. of Cybernetics, Mathematics and Physics, Havana, Cuba Galka, A.: Inst. of Experim. and Applied Physics, University of Kiel, Germany Sadato, N.: National Institute for Physiological Science, Okazaki, Japan Valdes-Sosa, P.: Cuban Neuroscience Center, Havana, Cuba Ozaki, Tohru: The Institute of Statistical Mathematics, Tokyo, Japan

  2. THE NN-ARx MODEL

  3. fMRI: Functional Magnetic Resonance Imaging Provides functional information about the state of the brain Measuresthe brain blood oxygenation level at some specific instant of time. fMRI Stimulus Blood vessels dilate bringing more blood highly oxygenated Neuronal activity demands more glucose and O2 Brain Activation fMRI signal increases in this area, detecting the change in the oxygenation level of the blood Metabolism

  4. Standard continuous-time modelfor theBOLD signal U(t) neuronal sinaptic activity X1(t) inducing signal X2(t) blood flow X3(t) blood volume X4(t) de-oxyhemoglobine Friston et al (2000) Riera et al (2004) Buxton et al (1998) BOLD signal

  5. Continuos-time model NN-ARx model fMRI activation maps based on the NN-ARx model. NeuroImage 23 (2004) 680–697 J. Riera, J. Bosch, O. Yamashita, R. Kawashima, N. Sadato, T. Okada,e and T. Ozaki

  6. Equations of the NN-ARx Model eXogenous variable Nearest Neighbors Innovations AR term MAIN FEATURES • Dynamical Model • Spatial dependency x x x x x t

  7. Equation of the NN-ARx Model eXogenous variable Nearest Neighbors Innovations AR term Whiteness Gaussianity Variance MAIN FEATURES • Dynamical Model • Spatial dependency • Innovations analysis • Long-range conectivity analysis Connection between y(v1) & y(v2) ?  v1  v2

  8. ACTIVATION ANALYSIS BASED ON NN-ARx MODELING

  9. Before Starting the Analysis. fMRI preprocessing • Realigning Correcting the fMRI scans for possible head movements, so the time series we see in one voxel over time corresponds approximately to the same site in the brain. • Time slicing • Correcting the time shifting introduced among slices while taking one fMRI scan. We perform these two preprocessing using the SPM software (Statistical Parametric Mapping, by Friston et al).

  10. Some exploratory tools for NN-ARx model fitting

  11. Visual experiment Task :Visual stimulus by black and white shuffled check board Data provided by Prof. N. Sadato, (National Institute of Physiological Sciences) Sampling frequency: 3s Resolution: 64× 64 × 36 # of time points : 60 3 T

  12. HRF activation for a visual experiment

  13. Model fitting for a voxel in the calcarine

  14. HRF in detail for the selected voxel

  15. Map of the innovations variance Double Click

  16. Map of autocorrelations at a selected voxel

  17. Correlation maps for different voxels Lag 0 Cerebelum R = 0.5 Vermix R = 0.7 Lingual R = 0.6

  18. Experiment from Sassa et al, IDAC. Tohoku University 1-Talk to a familiar person 2-Talk to an unfamiliar person 3-Listen from familiar person 4-Listen from unfamiliar person 5-Say an object name

  19. HRF for a voxel at the Hippocampus

  20. HRF for a voxel at the Cunneus Left Click here

  21. Testing for activation Fitted NN-ARx model T2 statistics at each voxel Permutation tests based on all T2 statistics

  22. Difference of Conditions 1 and 2 (1-2)

  23. Difference of Conditions 1 and 3 (1-3)

  24. Difference of Conditions 2 and 3 (2-3)Are these activations significative? T2 test for Condition 1

  25. T2 test for Condition 3

  26. REGIONAL CONNECTIVITY ANALYSIS

  27. Some methodological issues • Functional connectivity (observed dependencies) vs effective connectivity (causal relations). • In general, causal relations can not be inferred from observational data.

  28. Some approaches for connectivity analysis • Standard (zero-lag) correlation analysis • Structural equation modeling • Dynamic causal modeling

  29. v ● w * Correlations between two voxels based on innovations Instantaneous ● * Lagged ● * ● * Notation

  30. Equation of the NN-ARx Model eXogenous variable Nearest Neighbors Innovations AR term Whiteness Gaussianity Variance MAIN FEATURES • Dynamical Model • Spatial dependency • Innovations analysis • Long-range conectivity analysis Connection between y(v1) & y(v2) ?  v1  v2

  31. Regional Correlations • Calculate the vector of all possible correlations between all voxels v in region V vs. all voxels w in region W, pair to pair. • Take the square of the correlations in order to capture both positive and negative correlations. • Summarize the correlation between the two regions by the upper 90th percentile of the values in for all positive and negative lags. • Further, summarize the lagged correlations by: • For statistical significance we use the bootstrap technique.

  32. Results. Significant correlations for a group of subjects under a visual task

  33. Results. Significant correlations for a group of subjects performing a motor task

  34. Results. Significant correlations for a blind subject under a tactile discrimination task

  35. Some concluding remarks • NN-ARx modeling offers a dynamic approach for the analysis of both activation and connectivity from fMRI data. • Connectivity analysis based on innovations permits to clean the data from short-range connections and focus on long-range connections. • Regional connectivity measures that do not involved spatial averaging may be defined to atenue the confounding effects of lack of homogeneity within each region and of errors in brain segmentation.

  36. But…some limitations and cautions • fMRI has low time resolution (in comparison with neural time scale). • Flexibility in defining regional connectivity measures without spatial averaging is achieved at the expense of computer-intensive algorithms for statistical testing. • A high correlation between the past of a region and the future of another region does notimply causal connectivity. • The neurophysiological meaningof innovations in NN-Arx modeling should be further elucidated in the context of fMRI experiments to aid interpretaion of findings.

  37. THANKS

More Related