Diff for "SpmAnalysis" - Meg Wiki
location: Diff for "SpmAnalysis"
Differences between revisions 1 and 55 (spanning 54 versions)
Revision 1 as of 2007-04-14 13:45:51
Size: 781
Editor: RikHenson
Comment:
Revision 55 as of 2009-01-27 12:13:16
Size: 4445
Editor: RikHenson
Comment:
Deletions are marked like this. Additions are marked like this.
Line 2: Line 2:
SPM5 can perform basic preprocessing of MEG data, and source reconstruction, including ECD and variants of L2-norm imaging solutions. It's weaknesses are: i) ECD is limited (e.g, the Neuromag Xfit program is more flexible), ii) it's graphics are slow. It's strengths are: i) well integrated with MRI and fMRI analysis (eg, segmentation), ii) in Matlab, so relatively easy to understand and extend, iii) includes advanced statistical methods (e.g, for multiple comparisons across space or time), and novel variants of the L2-norm approach to the inverse problem, iv) includes group-based analysis (eg in MNI/Talairach space).
Line 4: Line 3:
For an introduction, see these slides: RiksSlides Many of us here use SPM for EEG and/or MEG analysis. The main reason is that SPM is freeware (requiring only Matlab), and thus modifiable and relatively easily to understand (if you have basic knowledge of a procedural computer language, given that Matlab is a fairly high-level language). Several of us are actively extending SPM for our own purposes, particularly in relation to the Neuromag type of MEG data.
Line 6: Line 5:
First, you will need to convert your *.FIF files into SPM files. SPM also offers some unique features, not available in other packages:

!!TBC!!

 * This page describes the procedure for creating a 3D SensorSpm (topography x time).

== Demos ==

 * For specific demo using data from our Neuromag MEG machine, see SpmDemo
 * For a fuller demo of other EEG/MEG analysis in SPM5 (though from a different MEG machine), including more general features (e.g, time-freq analysis, 3D statistical maps), with proper step-by-step instructions via the GUI, see: http://www.fil.ion.ucl.ac.uk/spm/data/mmfaces.html
 * For a more theoretical introduction to source localisation in SPM5, see these slides: attachment:henson-SPM-Grad08-4-meeg.ppt

== Some relevant papers ==

 * Summary of localisation approach using ReML for evoked and induced responses (mathematical; cites earlier development papers too): attachment:FristonEtAl_hbm_06.pdf
 * Basic considerations for Group Analyses (though using individual meshes): attachment:HensonEtAl_NI_07.pdf
 * Use of inverse-normalised ("canonical") cortical meshes: attachment:MattoutEtAl_JCIN_07.pdf
 * Choice of forward models for MEG (e.g, single-sphere vs BEM), including further validation of canonical meshes: attachment:HensonEtAl_NI_inpress.pdf
 * '''New method of Multiple Sparse Priors (MSP): attachment:FristonEtAl_NI_08_MSP.pdf'''

== The code ==

SPM has been developed at the FIL (in London), with input from people round the world. Here is its home: http://www.fil.ion.ucl.ac.uk/spm/

At the moment, we are using "SPM5" (with any latest updates automatically pulled from the FIL). At the CBU, SPM5 is installed here:

{{{/imaging/local/spm/spm5}}}

We have edited/refined some of these functions. These modifications are stored here:

{{{/imaging/local/spm/spm5/cbu_updates}}}

If you start SPM5 from Linux using Rhodri's wrappers, both these directories will be added to your Matlab path, with the cbu_updates directory higher in the path (so you will be using local CBU versions of any functions duplicated across these directories).

The functions in cbu_updates may be periodically updated, which will be preceded by an email to CBU imagers. More recent changes are under subversion control (SVN): If you are an external collaborator (outside the CBU network), or want access to the most recent changes (at your own risk), consult MeegCodeCbuSvn.
 
In the future, we will move to "SPM8". One exciting reason for this is that SPM8 will share the same data format with Field-Trip and EEGLAB (which are also Matlab packages for analysing EEG+MEG data). This will allow us to take the best features of each (e.g, source localisation and dynamic causal modelling in SPM8; time-frequency analysis and beam-formers in Field-Trip; Independent Component Analysis in EEGLAB). This is another advantage of Matlab-based academic-community-based freeware.

Basic preprocessing and source localisation in SPM8 will not differ much from SPM5. Tne main improvements concern Dynamic Causal Modelling (DCM) for EEG/MEG (together with better graphics and GUI). The reason that we are sticking with SPM5 for the moment is that the data format and forward modelling for Neuromag data is not quite finalised in Field-Trip/SPM8.

SPM5 can read raw and averaged FIF files, though you will probably first want to run your raw data through the [:Maxfilter:Maxfilter utility], particularly if you 1) used Active Shielding during acquisition, 2) if you want to apply (temporal) SSS to remove noise, 3) if you used continuous HPI and/or 4) if you want to transform all subjects to a common (device) space. Max Filter can also downsample (eg from 1000Hz to 200Hz) and convert the data into different datatypes (e.g, short), which will help reduce filesize and processing time.

SPM5 uses Brainstorm (another freeware Matlab package) for creating forward models. You do not need to know how to use Brainstorm, but if you are interested, here is a link:
http://neuroimage.usc.edu/brainstorm/
 

Analysis of MEG Data in SPM5

Many of us here use SPM for EEG and/or MEG analysis. The main reason is that SPM is freeware (requiring only Matlab), and thus modifiable and relatively easily to understand (if you have basic knowledge of a procedural computer language, given that Matlab is a fairly high-level language). Several of us are actively extending SPM for our own purposes, particularly in relation to the Neuromag type of MEG data.

SPM also offers some unique features, not available in other packages:

!!TBC!!

  • This page describes the procedure for creating a 3D SensorSpm (topography x time).

Demos

  • For specific demo using data from our Neuromag MEG machine, see SpmDemo

  • For a fuller demo of other EEG/MEG analysis in SPM5 (though from a different MEG machine), including more general features (e.g, time-freq analysis, 3D statistical maps), with proper step-by-step instructions via the GUI, see: http://www.fil.ion.ucl.ac.uk/spm/data/mmfaces.html

  • For a more theoretical introduction to source localisation in SPM5, see these slides: attachment:henson-SPM-Grad08-4-meeg.ppt

Some relevant papers

  • Summary of localisation approach using ReML for evoked and induced responses (mathematical; cites earlier development papers too): attachment:FristonEtAl_hbm_06.pdf

  • Basic considerations for Group Analyses (though using individual meshes): attachment:HensonEtAl_NI_07.pdf

  • Use of inverse-normalised ("canonical") cortical meshes: attachment:MattoutEtAl_JCIN_07.pdf

  • Choice of forward models for MEG (e.g, single-sphere vs BEM), including further validation of canonical meshes: attachment:HensonEtAl_NI_inpress.pdf

  • New method of Multiple Sparse Priors (MSP): attachment:FristonEtAl_NI_08_MSP.pdf

The code

SPM has been developed at the FIL (in London), with input from people round the world. Here is its home: http://www.fil.ion.ucl.ac.uk/spm/

At the moment, we are using "SPM5" (with any latest updates automatically pulled from the FIL). At the CBU, SPM5 is installed here:

/imaging/local/spm/spm5

We have edited/refined some of these functions. These modifications are stored here:

/imaging/local/spm/spm5/cbu_updates

If you start SPM5 from Linux using Rhodri's wrappers, both these directories will be added to your Matlab path, with the cbu_updates directory higher in the path (so you will be using local CBU versions of any functions duplicated across these directories).

The functions in cbu_updates may be periodically updated, which will be preceded by an email to CBU imagers. More recent changes are under subversion control (SVN): If you are an external collaborator (outside the CBU network), or want access to the most recent changes (at your own risk), consult MeegCodeCbuSvn.

In the future, we will move to "SPM8". One exciting reason for this is that SPM8 will share the same data format with Field-Trip and EEGLAB (which are also Matlab packages for analysing EEG+MEG data). This will allow us to take the best features of each (e.g, source localisation and dynamic causal modelling in SPM8; time-frequency analysis and beam-formers in Field-Trip; Independent Component Analysis in EEGLAB). This is another advantage of Matlab-based academic-community-based freeware.

Basic preprocessing and source localisation in SPM8 will not differ much from SPM5. Tne main improvements concern Dynamic Causal Modelling (DCM) for EEG/MEG (together with better graphics and GUI). The reason that we are sticking with SPM5 for the moment is that the data format and forward modelling for Neuromag data is not quite finalised in Field-Trip/SPM8.

SPM5 can read raw and averaged FIF files, though you will probably first want to run your raw data through the [:Maxfilter:Maxfilter utility], particularly if you 1) used Active Shielding during acquisition, 2) if you want to apply (temporal) SSS to remove noise, 3) if you used continuous HPI and/or 4) if you want to transform all subjects to a common (device) space. Max Filter can also downsample (eg from 1000Hz to 200Hz) and convert the data into different datatypes (e.g, short), which will help reduce filesize and processing time.

SPM5 uses Brainstorm (another freeware Matlab package) for creating forward models. You do not need to know how to use Brainstorm, but if you are interested, here is a link: http://neuroimage.usc.edu/brainstorm/

CbuMeg: SpmAnalysis (last edited 2013-03-08 10:02:27 by localhost)