Assessing the performance of an artifact detection/correction implemented in our Apps is a very important step. To have an idea how to do it, we sent the following message to the mailing list of different tools:
Dear XXX,
I performed an artifact detection/correction procedure and I am looking for a good way to visualize its results but also to assess its performance. I would imagine something akin to a script displaying a summary plot per participant showing all trials x channels and the scores of one or several measurements done before and after artifact detection/correction. More precisely, I would want to find a way to assess the quality of the artifact detection/correction or the quality of the corrected data, thanks to an index for instance.
Do you know a good application/function/script/favorite procedure that would allow me to do this? A paper reference would be of course most welcome!
Many thanks in advance,
Best
Message sent to | Date | Commentary | Answer |
---|---|---|---|
| |||
| I had to subscribe to the mailing list |
| |
Maybe ask directly Guio | |||
| Max sent an email (not exactly the same as above), but maybe ask on the forum instead |
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
Add Comment