daten uploaded!

Muthuraman Muthuraman muthuraman10 at HOTMAIL.COM
Wed Jan 10 17:44:25 CET 2007


Hello

discard my question regarding uploading the data in my previous mail
I have uploaded the file in www.yousendit.com
and this is the link

http://download.yousendit.com/543F313F300F6EAF

thanking you

with regards
M.Muthuraman.


>From: Robert Oostenveld <r.oostenveld at FCDONDERS.RU.NL>
>Reply-To: FieldTrip discussion list <FIELDTRIP at NIC.SURFNET.NL>
>To: FIELDTRIP at NIC.SURFNET.NL
>Subject: Re: [FIELDTRIP] Coherence!
>Date: Tue, 9 Jan 2007 17:53:23 +0100
>
>dear Muthuraman
>
>On 9 Jan 2007, at 11:08, Muthuraman Muthuraman wrote:
>>%MTMFFT
>>cfg.channel    = channelselection({'all' 'M1'}, data.label);
>>cfg.channelcmb = channelcombination({'all' 'M1'}, data.label);
>
>for DICS you will also need the cross-spectrum between all EEG  channels.
>Hence you should do cfg.channelcmb = 'all'. The code above  is fine if you
>want to look at coherence on channel level, which is  what you currently
>are looking at.
>
>>cfg.foilim     = [2 30];
>>cfg.taper      = 'hanning';
>>cfg.pad        = 1;
>>cfg.keeptrials = 'yes';
>>cfg.variance   = 'yes';
>>cfg.jackknife  = 'yes';
>>freqmtmfft     = freqanalysis(cfg, data);
>>fdmtmfft       = freqdescriptives([], freqmtmfft);
>
>This uses 1 hanning taper, and you have one-second long data  segments.
>Hence the frequency resolution is 1/1s = 1Hz.
>
>>%MTMWELCH
>>cfg.method     = 'mtmwelch';
>>cfg.output     = 'powandcsd';
>>cfg.foi        = [2:2:30];
>>cfg.t_ftimwin  = ones(1,length(cfg.foi)).*0.5;
>>cfg.taper      = 'hanning';
>>cfg.pad        = 1;
>
>This uses a 500ms window per trial to do the estimation. That gives  you a
>1/0.5 = 2Hz frequency resolution. The lower frequency  resolution implies
>that there is more frequency smoothing, so the  estimate will be different.
>If the true coherence peak is broad-band  (i.e. broader than 1 Hz), then
>the estimate will probably be better  (=higher), if not, then it will be
>worse (=lower).
>
>>cfg.keeptrials = 'yes';
>>cfg.variance   = 'yes';
>>cfg.jackknife  = 'yes';
>>freqmtmwelch   = freqanalysis(cfg, data);
>>fdmtmwelch     = freqdescriptives([], freqmtmwelch);
>
>You probably should do
>   cfg = []
>   cfg.variance   = 'yes';
>   cfg.jackknife  = 'yes';
>and then
>   fdmtmwelch     = freqdescriptives(cfg, freqmtmwelch);
>since the variance options do not apply to freqanalysis, but to
>freqdescriptives instead.
>
>>I have attached a figure with this mail[Coherence between C3 and  M1] in
>>which there is a difference in coherence values
>>for the fieldtrip method and the welch periodogram method with and
>>without laplacian transform
>>
>>Is there any input parameters i can change in the fieldtrip analysis
>>so that i can have the similar values, or is there any major  differences
>>in the methods
>>with fieldtrip and Welch periodogram method.
>
>There is a difference in scalp-current-density (laplacian) data and  EEG
>potential data, so I also expect coherence to be different. It is  better
>to report channel level CSD/laplacian coherence than channel- level
>ordinary potential difference, since SCD/laplacian is reference  free and
>potential is not. Note that you can also use FT to compute  coherence on
>SCD data, assuming that you can export the laplacian  data to a raw file. I
>don't know which software you used to compute  the SCD/laplacian and welsh
>coherence that you have in the figure,  therefore I also cannot comment on
>why the values are different. The  difference between the mtmfft and
>mtmwelsh can be explained (see  above) because the different frequency
>resolutions.
>
>best regards,
>Robert

_________________________________________________________________
MSN cricket features 'Cricketer of the Month'
http://content.msn.co.in/Sports/Cricket/Default.aspx



More information about the fieldtrip mailing list