plotsig

 Function to plot C where it is higher than a threshold sig
 useful for plotting coherence
 Usage: plotsig(C,sig,t,f)
 Inputs:
 C: input array t x f - also works for a single vector
 sig: significance level
 t: t axis grid for plot
 f: f axis grid for plot.
 c: color to use (default blue)-only meaningful for a line plot

0001 function plotsig(C,sig,t,f,c)
0002 % Function to plot C where it is higher than a threshold sig
0003 % useful for plotting coherence
0004 % Usage: plotsig(C,sig,t,f)
0005 % Inputs:
0006 % C: input array t x f - also works for a single vector
0007 % sig: significance level
0008 % t: t axis grid for plot
0009 % f: f axis grid for plot.
0010 % c: color to use (default blue)-only meaningful for a line plot
0011 if nargin < 4; error('Need at least 4 arguments'); end;
0012 if nargin < 5 | isempty(c); c='b'; end;
0013 [T,F]=size(C);
0014 if F==1; C=C'; [T,F]=size(C);end;
0015 if T~=length(t) | F~=length(f);
0016     error('frequency and/or time axes are incompatible with data'); 
0017 end;
0018 if T==1;
0019     dim=max(T,F);
0020     C=C(:);
0021     indx=find(C>sig);
0022     plot(f,C,c); 
0023 %     plot(f,C,f,mask.*C)
0024     line(get(gca,'xlim'),[sig sig]);
0025     xlabel('f'); ylabel('|C|');
0026 else
0027     mask=zeros(T,F);
0028     for n=1:length(t);
0029         for m=1:length(f);
0030            if C(n,m)>sig
0031               mask(n,m)=1;
0032            end;
0033         end;
0034     end;
0035     imagesc(t,f,(mask.*C)'); axis xy; colorbar
0036     xlabel('t'); ylabel('f');
0037 end;