[Eeglablist] How to re-mask ersp output

[Arnaud Delorme]

Dear Mash,

the tftopo is the correct function to use to remask time frequency decompositons. For example,

[ersp,itc,powbase,times,freqs,erspboot,itcboot] = pop_newtimef( EEG, 1, 1, [EEG.xmin EEG.xmax]*1000, [3 0.5], 'alpha', 0.05, 'mcorrect', 'fdr', 'plotersp', 'off', 'plotitc', 'off');
figure; tftopo(ersp, times, freqs, 'signifs', erspboot);

Best,

Arno

On Jul 6, 2012, at 8:18 PM, James Desjardins wrote:

> Hi Mash,
> 
> There may be a more elegant ways of achieving this, but the following works:
> 
> 
> 
> %PERFORM TIMEF CALCULATION (E.G. IC 13 FOR ONE CONDITION)...
> figure;
> [ersp,itc,powbase,times,freqs,erspboot,itcboot,itcphase] = ...
>                  timef(EEG.icaact(13,:,:), ...
> 		       EEG.pnts, ...
> 		       [EEG.xmin*1000 EEG.xmax*1000], ...
> 		       EEG.srate,[3,.5], ...
> 		       'freqs',[3 30], ...
> 		       'alpha',.01);
> 
> 
> %FLAG ERSP VALUES THAT ARE BELOW THE LOWER CONFIDENCE INTERVAL [x(:,:,1)]
> for i=1:length(times);x(:,i,1)=ersp(:,i)'<erspboot(1,:);end
> 
> %FLAG ERSP VALUES THAT ARE ABOVE THE UPPER CONFIDENCE INTERVAL [x(:,:,2)]
> for i=1:length(times);x(:,i,2)=ersp(:,i)'>erspboot(2,:);end
> 
> %CREATE CRIT ARRAY WHERE '1's FLAG THE TF COORDINATES OF SIGNIFICANCE...
> crit=double(any(x,3));
> 
> %PLOT THE ERSP ARRAY...
> figure;
> surf(times,freqs,ersp,'LineStyle','none');
> axis tight;
> view(0,90);
> 
> %PLOT THE CRIT ARRAY...
> figure;
> surf(times,freqs,crit,'LineStyle','none');
> axis tight
> view(0,90)
> 
> %PLOT THE MASKED ERSP ARRAY...
> figure;
> surf(times,freqs,ersp.*crit,'LineStyle','none');
> axis tight
> view(0,90)
> 
> 
> 
> -- 
> James Desjardins, MA
> Technician, Cognitive and Affective Neuroscience Lab
> Department of Psychology, Behavioural Neuroscience
> Brock University
> 500 Glenridge Ave.
> St. Catharines, ON, Canada
> L2S 3A1
> 905-688-5550 x4676
> 
> 
> Quoting "Mash, Clay (NIH/NICHD) [E]" <mashc at mail.nih.gov>:
> 
>> From:   >>timef details,
>> 
>> "If bootstrap statistics are computed (via the 'alpha'   probability  
>> input parameter), (time, freq) points with non-significant     
>> differences from 0 are colored green in the image (but not in the   
>> 'ersp'   output variable - use the output 'erspboot' variable to   
>> re-mask the 'ersp' output if desired)."
>> 
>> How do you re-mask the ersp output?  I would like to report the freq  
>> x time coordinates of the significant differences between conditions.
>> 
>> 
>> Clay Mash
>> 
>> _______________________________________________
>> Eeglablist page: http://sccn.ucsd.edu/eeglab/eeglabmail.html
>> To unsubscribe, send an empty email to eeglablist-unsubscribe at sccn.ucsd.edu
>> For digest mode, send an email with the subject "set digest mime" to  
>> eeglablist-request at sccn.ucsd.edu
>> 
> 
> 
> 
> _______________________________________________
> Eeglablist page: http://sccn.ucsd.edu/eeglab/eeglabmail.html
> To unsubscribe, send an empty email to eeglablist-unsubscribe at sccn.ucsd.edu
> For digest mode, send an email with the subject "set digest mime" to eeglablist-request at sccn.ucsd.edu