def calc_performance(evoked_raw, evoked_clean): ''' Gives a measure of the performance of the artifact reduction. Percentage value returned as output. ''' from jumeg import jumeg_math as jmath diff = evoked_raw.data - evoked_clean.data rms_diff = jmath.calc_rms(diff, average=1) rms_meg = jmath.calc_rms(evoked_raw.data, average=1) arp = (rms_diff / rms_meg) * 100.0 return np.round(arp)