...
In Gonzalez-Moreno et al., 2014, a SNR was computed to compare the quality of artifact correction by MaxFilter SSS and Maxfiter tSSS. Maybe in our case we could compute a SNR before and after Maxfilter (and also before and after the bad channel detection)? It would be better if this step is an App. This App would take in input the .fif file before MaxFilter and the .fif file after MaxFilter and would return a report.
To compute this SNR:
Select only MEG channels and exclude bad channels
Create events
Create epochs based on the events
Compute for each epoch its mean amplitude across all electrodes and times
Compute the mean of all the mean amplitudes of each epoch
Compute the standard error of that mean
SNR = result step 5. / result step 6.
Compute SNR on a subset of channels?
Select only magnetometers or gradiometers?
Maybe when the data has events, it’s best to create the epochs based on these existing events
In several Apps there is the function to compute the SNR the way it is described here. For now, this function is not used
It would be interesting to add this function to the
helper.pyonce we agreed on the way to compute the SNR.