checkElectrodeSignificance.m

function [p,h] = checkElectrodeSignificance(base_ch, signal_ch, zeroi, varargin)
% function [p,h] = checkElectrodeSignificance(base_ch, signal_ch, zeroi, varargin)
%
% Tests the time-averaged signal (in 1 sec window centered at time 0)
% versus the time-averaged baseline period. Tests via a simple Wilcoxon rank sum test 
% Vargargin: 1 = signal time window to test, 2 = histogram axis handle
% Output: self explanatory

if nargin > 3
    signali = varargin{1};
else
    signali = (zeroi-500):(zeroi+500);
end
if nargin > 4
    hist_ah = varargin{2};
else
    hist_ah = [];
end

mean_signal = mean(signal_ch(signali,:), 1);
mean_baseline = mean(base_ch,1);

[p, h] = ranksum(mean_signal, mean_baseline);

if ~isempty(hist_ah)
    %let's plot the Z'd data, just so it's more consistent
    s_base = std(mean(base_ch,1));
    mu_base = mean(mean_baseline);
    h1 = histogram(hist_ah, (mean_baseline-mu_base)./s_base, 'Normalization', 'probability');
    hold(hist_ah,'on'); 
    h2 = histogram(hist_ah, (mean_signal-mu_base)./s_base, 'Normalization', 'probability'); 
    h1.BinWidth = .1; h2.BinWidth = .1;
    h1.EdgeColor = 'none';  h2.EdgeColor = 'none'; 
    text(hist_ah, 0, .05, ['p=', num2str(p, 4)]);
end