= MEG visual tutorial: Group analysis (BIDS) =
''Authors: Francois Tadel, Elizabeth Bock. ''

The aim of this tutorial is to reproduce in the Brainstorm environment the analysis described in the SPM tutorial "[[ftp://ftp.mrc-cbu.cam.ac.uk/personal/rik.henson/wakemandg_hensonrn/Publications/SPM12_manual_chapter.pdf|Multimodal, Multisubject data fusion]]".  It is part of a collective effort to document and standardize MEG/EEG group analysis, see Frontier's research topic: [[https://www.frontiersin.org/research-topics/5158|From raw MEG/EEG to publication: how to perform MEG/EEG group analysis with free academic software]].

The data processed here consists in simultaneous MEG/EEG recordings of  16 subjects performing a simple visual task on a large number of famous,  unfamiliar and scrambled faces. This tutorial follows another page that explains how to process [[Tutorials/VisualSingle|one single subject]] in details.

<<TableOfContents(3,2)>>

== License ==
This dataset was obtained from the OpenNeuro project (https://openneuro.org/),  accession #ds117. It is made available under the Creative Commons  Attribution 4.0 International Public License. Please cite the following  reference if you use these data: <<BR>>Wakeman DG, Henson RN, [[http://www.nature.com/articles/sdata20151|A  multi-subject, multi-modal human neuroimaging dataset]], Scientific Data (2015)<<BR>>Any questions, please contact: rik.henson@mrc-cbu.cam.ac.uk

For citing the analysis, the processing pipeline is published in this  article:<<BR>>Tadel F, Bock E, Niso G, Mosher JC, Cousineau  M, Pantazis D, Leahy RM, Baillet S, [[https://www.frontiersin.org/articles/10.3389/fnins.2019.00076/full|MEG/EEG Group Analysis With Brainstorm]], Frontiers in Neuroscience, Feb 2019

== Download and installation ==
First, make sure you have enough space on your hard drive, at least '''40Gb''':

 * The Brainstorm database with all the data imported, downloaded from this website: '''20Gb'''
 * The same database after this tutorial: '''40Gb'''

You can follow this tutorial after processing the recordings for the 16 good subjects (6 runs per subject) as illustrated in the [[Tutorials/VisualSingle|single subject tutorial]]. Otherwise, we provide a Brainstorm protocol that includes all the imported data, downsampled at 275Hz:

 * Go to the [[http://neuroimage.usc.edu/bst/download.php|Download]] page, download the file '''TutorialGroup.zip''' (20Gb).
 * Unzip this file in your Brainstorm database folder (brainstorm_db).
 * In Brainstorm, menu File > Load protocol > Load from folder > Select '''brainstorm_db/TutorialGroup'''

The database you need in order to follow this tutorial should contain the following:

 * The individual anatomy imported from FreeSurfer for each subject (16 subjects).
 * The sensor level averages (MEG+EEG) for each run (downsampled to 275Hz or not).
 * The forward model, noise covariance and inverse models for each subject and each run. <<BR>><<BR>> {{attachment:db_tutorialgroup.gif||width="422",height="346"}}
 * This protocol TutorialGroup is produced from the single subject protocol TutorialVisual with the script: '''brainstorm3/toolbox/script/tutorial_visual_copy.m''' <<BR>><<BR>><<HTML(<div  style="border:1px solid black;  background-color:#EEEEFF; width:680px;  height:250px; overflow:scroll;  padding:10px; font-family:  Consolas,Menlo,Monaco,Lucida  Console,Liberation Mono,DejaVu Sans  Mono,Bitstream Vera Sans  Mono,Courier New,monospace,sans-serif;  font-size: 13px; white-space: pre;">)>><<EmbedContent("http://neuroimage.usc.edu/bst/viewcode.php?f=tutorial_visual_copy.m")>><<HTML(</div >)>>

== Overview of the analysis ==
=== Coregistration of the acquisition runs ===
For each subject, all the runs have been registered to a common head position with MaxFilter. To verify this, select all the channel files within one subject, right-click > Display sensors > MEG (all). The surfaces representing the MEG helmet are perfectly overlapping for all the runs. When the runs are not aligned, it looks like [[http://neuroimage.usc.edu/brainstorm/Tutorials/ChannelFile#Multiple_runs_and_head_positions|this]].

 . {{attachment:run_coreg.gif||width="486",height="163"}}

This means we can safely average or compare the MEG sensor values across runs within one subject. However, it is not reliable to average MEG recordings across subjects, because of the [[http://neuroimage.usc.edu/brainstorm/Tutorials/Workflows#MEG_recordings|anatomical differences between subjects]].

This doesn't mean we can estimate the sources only once per subject. We have computed different SSP projectors and selected different bad channels for each acquisition run. To be able to use this information efficiently we should estimate the sources for each run separately, then average the sources across runs.

The forward model is the same for all the runs within one subject, therefore it can be computed for the first run and copied to all the other runs.

=== Objectives ===
The objectives for this tutorial are to reproduce the analysis presented in the following documents:

 * Wakeman DG, Henson RN, [[http://www.nature.com/articles/sdata20151|A  multi-subject, multi-modal human neuroimaging dataset]]<<BR>>Scientific Data (2015)
 * Wakeman DG, Henson RN, [[ftp://ftp.mrc-cbu.cam.ac.uk/personal/rik.henson/wakemandg_hensonrn/Publications/WakemanHenson_biomag2010_poster.pdf|Functional and structural connrctivity in face processing: MEG, EEG, fMRI, MRI and DWI data]], Biomag 2016
 * SPM12 manual, chapter 41: [[ftp://ftp.mrc-cbu.cam.ac.uk/personal/rik.henson/wakemandg_hensonrn/Publications/SPM12_manual_chapter.pdf|Multimodal, Multisubject data fusion]]

Summary of the results we will compute in this tutorial:

 * Subject and grand averages for each condition (Famous, Unfamiliar, Scrambled).
 * Normalization of these averages (Z-score for the sources, or ERS/D for the time-frequency maps).
 * Projection of the sources results on a template and spatial smoothing of the source maps.
 * Contrast between faces (famous+unfamiliar) and non-faces (scrambled): difference of averages and significance test.
 * Contrast between famous faces and unfamiliar faces: difference of averages and significance test.
 * Sensors of interest: EEG070 (or EEG060 or EEG065 if EEG070 is marked as bad)

Expected effects:

 * OFA: Bilateral occipital face area
 * FFA: Bilateral fusiform face area
 * STS: Right superior temporal sulcus
 * Around 170ms (N170): Difference between Faces and Scrambled.
 * After 250ms: Difference between Famous and Unfamiliar.

The methodology that we will follow for computing the averages and the other statistics is described in the tutorial "[[Tutorials/Workflows|Workflows]]".

This tutorial page illustrates how to run the analysis using the interface. In practice, with a larger number of subjects, going through all these steps manually with the Brainstorm GUI is not practical, time consuming and prone to manipulation errors. We encourage our users to script their group analyses instead of performing all these steps manually. Typically, one would select a few subjects, prototype the group analysis pipeline with the interactive interface, generate the corresponding Matlab script, and finally run it on all the subjects of the database in reproducible way. Assembling a processing script is discussed in the online tutorial [[https://neuroimage.usc.edu/brainstorm/Tutorials/Scripting|Scripting]] and illustrated with the scripts tutorial_frontiers2018*.m.

== Subject averages: Famous, Unfamiliar, Scrambled ==
We will start by computing the subject-level averages for all the data types we have: sensor-level recordings, source maps and time-frequency results. We will use a weighted average to group the results we have for each run (weighted by the number of good trials). We want to compute the averages for each experimental condition separately (famous, unfamiliar, scrambled).

Since the files are already selected in the interface, we will also compute a grand average across subject. In this context, we will consider that each subject has the same weight in the grand average (option "weighted" not selected). Note that it is not accurate to average MEG recordings across subjects, but we can do it just to get a general idea of the group effects ([[http://neuroimage.usc.edu/brainstorm/Tutorials/Workflows#MEG_recordings|more information]]).

=== MEG/EEG ===
 * Drag and drop all the subjects in Process1. Select button '''[Process recordings]'''.
 * Select process '''Average > Average files''': '''By trial group (subject average)'''<<BR>>Arithmetic average, Weighted
 * Add process  '''Average > Average files''': '''By trial group (grand average)'''<<BR>>Arithmetic average, __Not__ weighted <<BR>><<BR>> {{attachment:average_process_data.gif||width="492",height="500"}}
 * In output, you get three grand averages in the "Group analysis/Inter-subject" folder (top), and three subject averages per subject in the "intra-subject" folders (bottom). {{attachment:average_files_data.gif||width="653",height="273"}}

=== Sources ===
 * Keep all the subjects selected in Process1. Select button '''[Process sources]'''.
 * If you run the averaging process directly with this file selection, it would average together the MEG and the EEG sources (36 files selected per subject instead of 18). To average separately the two modalities, we can use the '''Filter''' box at the bottom-right corner of the Brainstorm window. Enter "EEG" and it would select only the files with a comment including the tag "EEG".
 * Select process '''Average > Average files''': '''By trial group (subject average)'''<<BR>>EEG source files only, Arithmetic average, Weighted
 * Add process '''File > Add tag''': '''EEG'''.<<BR>><<BR>> {{attachment:average_process_sources.gif||width="603",height="504"}}
 * Repeat the same steps with the '''MEG''' source files (select all subjects, filter, average, tag).
 * This time, we cannot compute a grand average directly because the source spaces used to estimate the brain activity do not match between subjects. We have to go through an additional step of projection on a template anatomy. We will also wait to have normalized the source maps before displaying them.<<BR>><<BR>> {{attachment:average_files_sources.gif||width="243",height="208"}}

=== Time-frequency ===
 * Keep all the subjects selected in Process1. Select button '''[Process time-freq]'''.<<BR>>Make sure you reset the Filter search box at the bottom-right of the Brainstorm window.
 * Run process '''Average > Average files''': '''By trial group (subject average)'''<<BR>>Arithmetic average, Weighted, Match signals<<BR>><<BR>> {{attachment:average_process_tf.gif||width="415",height="498"}}
 * Like for the sources, we will first normalize the time-frequency maps with respect with a baseline before computing grand averages. Like the sensor data: it is valid to average EEG recordings across subjects, but it is not accurate to average MEG sensors across subjects.
 * We can average the time-frequency files across runs (it is as valid as averaging the MEG or EEG signals), but we need to select this option "match signals between files using their names" to account for the variable number of bad channels we may have across runs or subjects.  <<BR>><<BR>> {{attachment:average_files_tf.gif||width="633",height="250"}}

== Subject averages: Faces ==
One of the contrast we want to study is faces (Famous and Unfamiliar) vs non-faces (Scrambled). We need to re-average the Famous and Unfamiliar averages together.

=== MEG/EEG ===
 * In Process2: FilesA = all the Famous subject averages, FilesB = all the Unfamiliar subject averages. The subjects must be in the same order in the two lists because the process will process the files by pairs of files (A,B).
 * Select process '''Other > Average A&B''': Weighted.<<BR>>Add process '''File > Set comment''': Comment = "WAvg: Avg: Faces"<<BR>>Add process '''Average > Average files''': '''By trial group (subject average)''': __Not__ weighted<<BR>><<BR>> {{attachment:average_faces_data.gif||width="669",height="531"}}
 * Selecting that many files from the database explorer can be long when done manually, but can be very efficient from a script. Check out the script equivalent at the end of this page.
 * It creates one new average for each pair (A,B) and computes the grand average for condition Faces. <<BR>><<BR>> {{attachment:average_faces_data2.gif||width="299",height="234"}}

=== Sources ===
 * In Process2: FilesA = all the '''Famous '''subject averages (EEG only) [Process sources]
 * In Process2: FilesB = all the '''Unfamiliar''' subject averages (EEG only) [Process sources]
 * Make sure the order of the files matches in the two lists.
 * Select process '''Other > Average A&B''': Weighted.<<BR>>Add process '''File > Set comment''': Comment = "WAvg: Avg: Faces | EEG"<<BR>><<BR>> {{attachment:average_faces_sources.gif||width="614",height="356"}}
 * Repeat the same steps for the '''MEG''' sources.<<BR>><<BR>> {{attachment:average_faces_sources2.gif||width="265",height="270"}}

=== Time-frequency ===
 * In Process2: FilesA = all the '''Famous''' subject averages (EEG only) [Process timefreq]
 * In Process2: FilesB = all the '''Unfamiliar''' subject averages (EEG only) [Process timefreq]
 * Select process '''Other > Average A&B''': Weighted.<<BR>>Add process '''File > Set comment''': Comment = "WAvg: Avg: Faces"<<BR>><<BR>> {{attachment:average_faces_timefreq.gif||width="604",height="331"}}
 * In output: the average for the Face condition, for each subject. <<BR>><<BR>> {{attachment:average_faces_timefreq2.gif||width="284",height="76"}}

== Subject averages: Within-subject differences ==
To detect correctly the differences between two conditions at the source level, we need to estimate the differences of the conditions for each subject, and then normalize the difference (see tutorial [[http://neuroimage.usc.edu/brainstorm/Tutorials/Workflows#Constrained_cortical_sources|Workflows]]).

=== Sources: Faces - Scrambled ===
 * In Process2: FilesA = all the '''Faces''' subject averages (EEG only) [Process sources]
 * In Process2: FilesB = all the '''Scrambled''' subject averages (EEG only) [Process sources]
 * Select process '''Difference > Difference A-B''': Do __'''not'''__ use absolute values.
 * Add process '''File > Set comment''': Comment='''Faces - Scrambled | EEG''' <<BR>><<BR>> {{attachment:diffsubj_faces.gif||width="606",height="331"}}
 * Repeat the same steps for the '''MEG''' sources.

=== Sources: Famous - Unfamiliar ===
 * In Process2: FilesA = all the '''Famous''' subject averages (EEG only) [Process sources]
 * In Process2: FilesB = all the '''Unfamiliar''' subject averages (EEG only) [Process sources]
 * Select process '''Difference > Difference A-B''': Do __'''not'''__ use absolute values.
 * Add process '''File > Set comment''': Comment='''Famous - Unfamiliar | EEG'''
 * Repeat the same steps for the '''MEG''' sources. <<BR>><<BR>> {{attachment:diffsubj_files.gif||width="284",height="240"}}

== Subject averages: Filter and normalize ==
Before comparing the averages across subjects we are going to low-pass filter the signals below 32Hz (to smooth possible latency differences between subjects) and normalize the source and time-frequency values with respect with a baseline (see tutorial [[http://neuroimage.usc.edu/brainstorm/Tutorials/Workflows#Constrained_cortical_sources|Workflows]]).

=== MEG/EEG ===
 * In Process1, select all the '''Intra-subject''' folders from all the subjects, select '''[Process recordings]'''. For a faster selection, you can use the view "Functional data (sorted by conditions)".
 * Select process '''Pre-process > Band-pass filter''': 0Hz-32Hz, MEG,EEG, 60dB, No mirror, Overwrite.
 * Add process '''Extract > Extract time''': Time window='''[-200,900]ms''' <<BR>><<BR>>  {{attachment:filter_data_process.gif||width="622",height="320"}}
 * Two tags are added at the end of the comments of the averaged recordings.  The selection of a smaller time window gets rid of most of the possible edge effects caused by the filter. <<BR>><<BR>> {{attachment:filter_data_files.gif||width="381",height="119"}}

=== Sources ===
 * In Process1, select all the '''Intra-subject''' folders from all the subjects, select '''[Process sources]'''.
 * Select process '''Pre-process > Band-pass filter''': 0Hz-32Hz, MEG,EEG, Mirror, Overwrite
 * Add process '''Extract > Extract time''': Time window='''[-200,900]ms'''
 * Add process '''Standardize > Baseline normalization''':  Baseline=[-200,-5]ms, Z-score, Overwrite <<BR>><<BR>> {{attachment:filter_sources_process.gif||width="671",height="417"}}
 * Three tags are added at the end of the comments of the averaged sources.  <<BR>><<BR>> {{attachment:filter_sources_files.gif||width="403",height="208"}}

=== Time-frequency ===
 * In Process1, select all the '''Intra-subject''' folders from all the subjects, select '''[Process time-freq]'''.
 * Run process '''Standardize > Baseline normalization''':  Baseline=[-200,-5]ms, ERS/ERD, Overwrite <<BR>><<BR>> {{attachment:filter_tf_process.gif||width="434",height="412"}}
 * One tag is added at the end of the comments of the averaged time-frequency maps.  <<BR>><<BR>> {{attachment:filter_tf_files.gif||width="393",height="86"}}

== Subject averages: Screen captures ==
Now we have all the measures ready to be compared across subjects: MEG/EEG, sources, time-frequency. Let's take a few screen captures to make sure the primary visual response look good for all the subjects. The screen captures below represent from left to right:

 * Registration sensors/MRI: Head surface, MEG (inside of the helmet), EEG+headpoints (green points).
 * Faces: EEG recordings
 * Faces: MEG sources at 107ms
 * Faces: Time-frequency for EEG070

sub-01<<BR>> {{attachment:sub002.gif||width="700"}}

sub-02<<BR>> {{attachment:sub003.gif||width="700"}}

sub-03<<BR>> {{attachment:sub004.gif||width="700"}}

sub-04<<BR>> {{attachment:sub006.gif||width="700"}}

sub-05<<BR>> {{attachment:sub007.gif||width="700"}}

sub-06<<BR>> {{attachment:sub008.gif||width="700"}}

sub-07<<BR>> {{attachment:sub009.gif||width="700"}}

sub-08<<BR>> {{attachment:sub010.gif||width="700"}}

sub-09<<BR>> {{attachment:sub011.gif||width="700"}}

sub-10<<BR>> {{attachment:sub012.gif||width="700"}}

sub-11<<BR>> {{attachment:sub013.gif||width="700"}}

sub-12<<BR>> {{attachment:sub014.gif||width="700"}}

sub-13<<BR>> {{attachment:sub015.gif||width="700"}}

sub-14<<BR>> {{attachment:sub017.gif||width="700"}}

sub-15<<BR>> {{attachment:sub018.gif||width="700"}}

sub-16<<BR>> {{attachment:sub019.gif||width="700"}}

== Group analysis: MEG/EEG ==
=== Grand averages ===
We have already computed the group averages for all the conditions. We will not look much further at these results as we are more interested in the contrasts between conditions. Below are screen captures for all group averages (top=MEG MAG, bottom=EEG). Topography at [50,100,150,200,250,300]ms.

Faces: Famous / Unfamiliar<<BR>>

{{attachment:gavg_data_famous.gif||width="350",height="237"}} {{attachment:gavg_data_unfamiliar.gif||width="350",height="236"}}

Scrambled<<BR>> {{attachment:gavg_data_scrambled.gif||width="350",height="237"}}

=== Faces - Scrambled: Differences of averages ===
We could compute the contrasts directly from the grand averages, but we will do it from the subject averages because it will be the same file selection for the statistics.

 * In Process2: FilesA = all the Faces subject averages (from the Intra-subject folders).
 * In Process2: FilesB = all the Scrambled subject averages (from the Intra-subject folders).
 * Run process: '''Test > Difference of means''': Arithmetic average, __Not__ weighted. <<BR>><<BR>> {{attachment:stat_data_diffmean.gif||width="552",height="282"}}
 * Rename the file: '''Faces - Scrambled''' (top=MEG MAG, bottom=EEG). <<BR>><<BR>> {{attachment:diff_data_faces.gif||width="450",height="303"}}

=== Faces - Scrambled: Significance testing ===
We have computed the amplitude of the difference between the two conditions, and we will try now to estimate a p-value indicating the significance of each of these values. In all the screen captures below: top=MEG MAG, bottom=EEG.

 * In Process2: Keep the same file selection.
 * Run process: '''Test > Parametric test: Paired''': All file, All sensors, No average, two-tailed. <<BR>><<BR>> {{attachment:stat_data_param.gif||width="552",height="326"}}
 * Rename the file: '''Faces - Scrambled: Parametric t-test'''. Display with '''<<HTML(&alpha;)>>=0.05, FDR-corrected'''.  <<BR>><<BR>> {{attachment:ttest_data_faces.gif||width="448",height="302"}}
 * You can run other tests in a similar way, with almost identical results.
 * Process: '''Test > Permutation test: Paired''': All file, All sensors, Paired t-test, 1000 randomizations. Display with '''<<HTML(&alpha;)>>=0.05, FDR-corrected'''.  <<BR>><<BR>> {{attachment:ttestperm_data_faces.gif||width="450",height="303"}}
 * Process: '''Test > FieldTrip: ft_timelockstatistics''': All file, EEG, Paired t-test, 1000 randomizations, correction=cluster, cluster alpha=0.05. The cluster-based statistics must be executed on one type of sensors at a time (EEG, MEG MAG or MEG GRAD), because it tries to identify spatio-temporal clusters which groups adjacent sensors. <<BR>><<BR>> {{attachment:cluster_data_faces.gif||width="450",height="151"}}

=== Famous - Unfamiliar: Differences of averages ===
 * In Process2: FilesA = all the Famous subject averages (from the Intra-subject folders).
 * In Process2: FilesB = all the Unfamiliar subject averages (from the Intra-subject folders).
 * Run process: '''Test > Difference of means''': Arithmetic average, __Not__ weighted.
 * Rename the file: '''Famous - Unfamiliar'''. <<BR>><<BR>> {{attachment:diff_data_famous.gif||width="450"}}

=== Famous - Unfamiliar: Significance testing ===
 * In Process2: Keep the same file selection.
 * Run process: '''Test > Parametric test: Paired''': All file, All sensors, No average, two-tailed.
 * Rename the file: '''Faces - Scrambled: Parametric t-test'''. Display with '''<<HTML(&alpha;)>>=0.05, FDR-corrected'''. <<BR>><<BR>> {{attachment:ttest_data_famous.gif||width="450",height="303"}}

 * Run process: '''Test > FieldTrip: ft_timelockstatistics''': All file, EEG, Paired t-test, 1000 randomizations, correction=cluster, cluster alpha=0.05. <<BR>><<BR>> {{attachment:cluster_data_famous.gif||width="450"}}

== Group analysis: Sources ==
=== Project sources on template ===
The sources were estimated on the individual anatomy of each subject, the resulting cortical source maps cannot be averaged directly. We need first to [[http://neuroimage.usc.edu/brainstorm/Tutorials/CoregisterSubjects|re-interpolate]] all the individual results on a common template (the ICBM152 brain, available in the "default anatomy" folder of the protocol). We also need to extract the absolute values for these source maps: the sign of the minimum norm maps are relative to the orientation of the current with respect to the surface normal, which can [[http://neuroimage.usc.edu/brainstorm/Tutorials/Workflows#Constrained_cortical_sources|vary between subjects]].

 * In Process1, select all the '''Intra-subject''' folders from all the subjects, select '''[Process sources]'''. For a faster selection, you can use the view "Functional data (sorted by conditions)".
 * Select process''' Pre-process > Absolute values''': Sensor types=MEG,EEG, Overwrite.
 * Add process '''Sources > Project on default anatomy'''. <<BR>><<BR>> {{attachment:project_sources.gif||width="520",height="324"}}
 * All the 192 subject averages are projected in the same folder "Group analysis / Intra-subject". Re-organize all the projected source files in new folders, one folder per condition.
 * To create a folder, right-click on Group analysis > New folder. To move a file to the new folder: drag and drop it in the database explorer or use the keyboard shortcuts Ctrl+X/Ctrl+V.<<BR>><<BR>> {{attachment:project_sources2.gif||width="296",height="293"}}

=== Spatial smoothing ===
The source maps estimated with constrained orientations can show very focal activity: two adjacent vertices may have very different normals, and therefore very different current values. When averaging multiple subjects, the peaks of activity may not align very well across subjects. Smoothing spatially the source maps may help obtaining better group results.

 * In Process1, select all the source maps in Group analysis.
 * Run process Sources > Spatial smoothing: <<BR>><<BR>> {{attachment:sources_smooth.gif||width="504",height="348"}}

=== MEG: mean(|Faces-Scrambled|) ===
 * In Process1, select all the source files in Group_analysis/Faces-Scrambled_MEG.
 * Run process '''Average > Average files''': By folder (grand average), Not weighted. <<BR>>Display with an arbitrary threshold: Z-score>3. <<BR>><<BR>> {{attachment:sources_absdiff_meg.gif||width="300"}}  {{attachment:sources_absdiff_meg_right.gif||width="300"}}
 * Regions of interest: OFA (Occipital Face Area), FFA (Fusiform Face Area), V1. <<BR>><<BR>> {{attachment:scouts_anat.gif||width="232",height="261"}} {{attachment:scouts_time_absdiff.gif||width="347",height="260"}}
 * This approach identifies correctly the effect size at the subject level. But because the sign is ambiguous across subjects, we need to average the absolute values of the differences. Therefore we are able to identify when and where the two conditions are different, but not which condition leads to a stronger brain response.

=== MEG: Chi2-test |Faces-Scrambled|=0 ===
 * In Process1, select all the source files in Group_analysis/Faces-Scrambled_MEG.
 * Run process '''Test > Parametric test against zero''': All file, One-sample Chi2-test two-tailed.
 * Screen capture: <<HTML(&alpha;)>>'''='''0.05 FDR-corrected<<BR>><<BR>> {{attachment:sources_chi2_meg.gif||width="300"}} {{attachment:sources_chi2_meg_right.gif||width="300"}}
 * This tests the hypotesis H0:[mean_subject(|Faces-Scrambled|) = 0]<<BR>>This  test is appropriate to detect the differences between the two  conditions, but cannot determine which condition is stronger. Its obvious limitation is that it is too sensitive: post-stimulation, everything is shown as significant.

=== MEG: Chi2-test log(|Faces-Scrambled|)=0 ===
 * Right-click on the folder Group_analysis/Faces-Scrambled_MEG > '''File > Duplicate folder'''. <<BR>><<BR>> {{attachment:duplicate_folder.gif||width="351",height="235"}}
 * Rename the new folder: '''Faces-Scrambled_MEG_log'''.
 * In Process1, select all the source files in Group_analysis/Faces-Scrambled_MEG_log.
 * Select process '''Pre-process > Run Matlab command > "Data = log(Data)";''' <<BR>> Add process '''Test > Parametric test against zero''': All file, One-sample Chi2-test two-tailed. <<BR>><<BR>> {{attachment:log_process.gif||width="654",height="378"}}
 * Screen capture: <<HTML(&alpha;)>>=0.05, FDR-corrected<<BR>><<BR>>  {{attachment:sources_chi2log_meg.gif||width="300"}} {{attachment:sources_chi2log_meg_right.gif||width="300"}}

=== MEG: mean(|Faces|)-mean(|Scrambled|) ===
 * In Process2: FilesA = all the source files in Group_analysis/Faces_MEG.
 * In Process2: FilesB = all the source files in Group_analysis/Scrambled_MEG.
 * Run process '''Test > Difference of means''': Arithmetic average, __Not__ weighted. <<BR>>Display with an arbitrary threshold: Z-score>3<<BR>><<BR>> {{attachment:sources_diffmean_meg.gif||width="300"}}   {{attachment:sources_diffmean_meg_right.gif||width="300"}}
 * Regions of interest: OFA (Occipital Face Area), FFA (Fusiform Face Area), V1. <<BR>><<BR>> {{attachment:scouts_anat.gif||width="232",height="261"}}  {{attachment:scouts_time.gif||width="345",height="261"}}
 * This approach is not able to detect correctly all the differences between the two conditions. At the subject level,we discarded the sign of each condition, therefore we are not able to detect differences in the orientation of the sources or phase delays between conditions.

=== MEG: Student t-test |Faces|=|Scrambled| ===
 * In Process2: FilesA = all the source files in Group_analysis/Faces_MEG.
 * In Process2: FilesB = all the source files in Group_analysis/Scrambled_MEG.
 * Run process '''Test > Parametric test: Paired''':  All file, No average, two-tailed.
 * Screen capture: <<HTML(&alpha;)>>=0.05 FDR-corrected <<BR>><<BR>> {{attachment:sources_ttest_meg.gif||width="300"}} {{attachment:sources_ttest_meg_right.gif||width="300"}}
 * This test the hypothesis H0:[mean(abs(Faces))=mean(abs(Scrambled))]. <<BR>>This measure is not appropriate to properly identify the differences between the two conditions because it is insensitive to the sign of the sources within each subject.

=== EEG: Faces-Scrambled ===
 * Repeat the same operations with the EEG sources.
 * Difference of averages: |Faces-Scrambled| (Z-score>3): <<BR>>  {{attachment:sources_diffmean_eeg.gif||width="300"}} {{attachment:sources_diffmean_eeg_right.gif||width="300"}}
 * Parametric Chi2-test: mean_subject(|Faces-Scrambled|) = 0, <<HTML(&alpha;)>>=0.05 FDR-corrected <<BR>> {{attachment:sources_chi2_eeg.gif||width="300"}} {{attachment:sources_chi2_eeg_right.gif||width="300"}}
 * Parametric Chi2-test: '''log'''(mean_subject(|Faces-Scrambled|)) = 0, <<HTML(&alpha;)>>=0.05 FDR-corrected <<BR>> {{attachment:sources_chi2log_eeg.gif||width="300"}} {{attachment:sources_chi2log_eeg_right.gif||width="300"}}
 * Difference of averages: |Faces|-|Scrambled| (Z-score>3) <<BR>> {{attachment:sources_absdiff_eeg.gif||width="300"}} {{attachment:sources_absdiff_eeg_right.gif||width="300"}}
 * Parametric Student t-test: mean(abs(Faces))=mean(abs(Scrambled)), <<HTML(&alpha;)>>=0.05 FDR-corrected <<BR>> {{attachment:sources_ttest_eeg.gif||width="300"}} {{attachment:sources_ttest_eeg_right.gif||width="300"}}

== Scripting ==
The following script from the Brainstorm distribution reproduces the analysis presented in this tutorial page: '''brainstorm3/toolbox/script/tutorial_visual_''''''group.m'''

Execution reports: [[http://neuroimage.usc.edu/bst/examples/report_TutorialGroup_2snapshots.html|Subject averages]], [[http://neuroimage.usc.edu/bst/examples/report_TutorialGroup_3meeg.html|MEG/EEG group results]], [[http://neuroimage.usc.edu/bst/examples/report_TutorialGroup_4sources.html|Source group results]]

<<HTML(<div  style="border:1px solid black;  background-color:#EEEEFF; width:720px;  height:500px; overflow:scroll;  padding:10px; font-family:  Consolas,Menlo,Monaco,Lucida  Console,Liberation Mono,DejaVu Sans  Mono,Bitstream Vera Sans  Mono,Courier New,monospace,sans-serif;  font-size: 13px; white-space: pre;">)>><<EmbedContent("http://neuroimage.usc.edu/bst/viewcode.php?f=tutorial_visual_group.m")>><<HTML(</div >)>>

<<EmbedContent(http://neuroimage.usc.edu/bst/get_feedback.php?Tutorials/VisualGroup)>>