Closed RobinGuillard closed 4 years ago
something like :
amplitude = rms(epoch)
power_num = np.mean(amplitude[labels]**2)
power_denum = np.mean(amplitude[np.invert(labels)]**2)
power_ratio = power_num / power_denum
?
option 2
amplitude = rms(epoch)
power_num = np.mean(amplitude[labels]**2)
power_denum = np.mean(amplitude**2)
power_ratio = power_num / power_denum
Better option 1
It could be cool to report the mean ratio of power between events and (adaptive) baseline values: in fact, one hypothesis to test is "are MEMA activity for tinnitus patients more powerful than for controls". When I inspect data of "pure" middle ear events, I think it's possible this hypothesis is true.