Muscle detection with inverse filtering
Artifacts are indicated with NaN's.
[INI,S,E] = detectmuscle(S [, iter [,1]])
[INI,S,E] = detectmuscle(S [, Fs, Mode]) with Mode>1
[INI,S,E] = detectmuscle(S,arg2,Mode)
iter number of iterations [default:1]
INI.MU mean of S
INI.InvFilter coefficients of inverse filter
S outlier replaced by NaN
E isnan(E) indicates muscle artifact
Mode 1: [default] inverse filtering
2: based on gradient, range and amplitude [BrainVision method]
3: slope > 11 uV/sample [1]
4: beta2 > 0.9 uV^2/Hz [1]
References:
[1] Van de Velde, M., Van Erp, G., Cluitmans, P., 1998.
Detection of muscle artefact in the normal human awake EEG.
Electroencephalography and Clinical Neurophysiology 107 (2), 149-158.