function [RMS] = rms_datastack(datastack)
%gives the RMS over all sensors for each time and epoch. 
RMS = sqrt(mean((datastack(:,:,:).^2),2));
RMS = permute(RMS,[1,3,2]);
end