[Eeglablist] NewCrossf -- bootstrap

[Arnaud Delorme]

Dear Michele and Makoto,

yes, there was a problem with the documentation. The number of values for the bootstrap thresholds is 18 values. Coherence values are complex numbers and the function only return an upper threshold (the absolute value of the coherence should be larger than a specific number). I have updated the help message of the function.

If you want to see the mask of significant values, you may type

mask  = mcoh > repmat(cohboot', [1 size(mcoh,2)]); %or more efficiently mask  = bsxfun(@gt, mcoh, cohboot');
figure; imagesc(mask);

Thanks,

Arno

On Aug 27, 2012, at 10:30 AM, Makoto Miyakoshi wrote:

> Dear Arno,
> 
> She says cohboot does not produce Matrix (nfreqs,2) as explained in
> the help but instead it produces Matrix(freq,1). Would you please help
> us?
> 
> Makoto
> 
>> The function help say that:
>> %       cohboot     = Matrix (nfreqs,2) of [lower;upper] coher signif.
>> limits
>> %                     if 'boottype' is 'trials',  (nfreqs,timesout, 2)
>> 
>> But I obtain always a Matrix(nfreqs,1) and I'm not able to use it.
> 
> 2012/8/23 Michele Barsotti <nomeserio at gmail.com>:
>> Hi dear experts,
>> 
>> I'm using the function newcrossf() to calculate coherence between two
>> electrodes.
>> 
>> As optional INPUT parameters I give the following:
>>         type: 'phasecoher'
>>       amplag: []
>>       subitc: 'off'
>>      shuffle: []
>>      detrend: 'off'
>>        rmerp: 'off'
>>      winsize: 256
>>     timesout: 895
>>     padratio: 1
>>        freqs: [3 40]
>>       nfreqs: 18
>>    freqscale: 'linear'
>>     baseline: [-2.5039e+003 -3.9063]
>>        alpha: 0.0500
>>        naccu: 200
>>     boottype: 'shuffle'
>>     baseboot: 1
>>     condboot: 'abs'
>>      topovec: [5 9]
>>        elocs: [1x13 struct]
>>     chaninfo: [1x1 struct]
>>       maxamp: 1
>>        title: 'coherence between, C3 C4'
>> 
>> And newcrossf returns as OUTPUT:
>>         mcoh: [1x18 single]
>>     timesout: [1x895 double]
>>     freqsout: [1x18 double]
>>      cohboot: [18x1 double]
>>    cohangles: [18x895 single]
>>     allcoher: [18x895x20 single]
>>       alltfX: [18x895x20 single]
>>       alltfY: [18x895x20 single]
>>          coh: [18x895 single]
>> 
>> Now, I would like to work with the coherence bootstraped but I don't know
>> how to use the cohboot:[18x1] with the coherence matrix coh: [18x895].
>> 
>> The function help say that:
>> %       cohboot     = Matrix (nfreqs,2) of [lower;upper] coher signif.
>> limits
>> %                     if 'boottype' is 'trials',  (nfreqs,timesout, 2)
>> 
>> But I obtain always a Matrix(nfreqs,1) and I'm not able to use it.
>> 
>> Any suggestion?
>> 
>> Thank you
>> 
>> Michele
>> 
>> 
>> 
>> 
>> 
>> 
>> 
>> 
>> --
>>        -Michele-
>> 
>> _______________________________________________
>> 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
> 
> 
> 
> -- 
> Makoto Miyakoshi
> JSPS Postdoctral Fellow for Research Abroad
> Swartz Center for Computational Neuroscience
> Institute for Neural Computation, University of California San Diego