Hi Claudia,<br><br>I think that if you want to do non-parametric tests using fieldtrip the cluster randomization statistics routines would help you a great deal. I think a lot of people using fieldtrip are using this method. You could have a look at this site:<br>
<br><a href="http://fieldtrip.fcdonders.nl/tutorial/cluster_permutation_freq">http://fieldtrip.fcdonders.nl/tutorial/cluster_permutation_freq</a><br><br>and the paper:<br><pre style="font-family: arial,helvetica,sans-serif;">
<i><font size="2">Eric Maris, Robert Oostenveld, Nonparametric statistical testing of EEG- and MEG-data, Journal of Neuroscience Methods, Volume 164, Issue 1, 15 August 2007, Pages 177-190, ISSN 0165-0270, DOI: 10.1016/j.jneumeth.2007.03.024.<br>
(<a href="http://www.sciencedirect.com/science/article/B6T04-4NFXDDW-1/2/ebd219fd45cb52bf39c710fe25635040">http://www.sciencedirect.com/science/article/B6T04-4NFXDDW-1/2/ebd219fd45cb52bf39c710fe25635040</a>)<br>Keywords: Nonparametric statistical testing; Hypothesis testing; EEG; MEG; Multiple comparisons problem</font></i></pre>
<br>Also, this mailing list is not
a help desk but is read by fellow researchers also using fieldtrip. And
I would like to say that although it sometimes takes a little while before you
get your answer, that is because everybody is busy and I am impressed
and happy with the way everybody on this list is nice enough to take
the time to help eachother!<br><br>Best,<br><br>Hanneke<br><br><div class="gmail_quote">On Fri, Sep 18, 2009 at 9:15 PM, Bonin, Claudia (NIH/NINDS) [F] <span dir="ltr"><<a href="mailto:boninc@ninds.nih.gov">boninc@ninds.nih.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">DOES ANYONE KNOW HOW TO USE THE FIELDTRIP FUNCTION FREQSTATISTICS TO RUN A KRUSKALWALLIS TEST?<br>
<div class="im"><br>
-----Original Message-----<br>
From: FieldTrip discussion list [mailto:<a href="mailto:FIELDTRIP@NIC.SURFNET.NL">FIELDTRIP@NIC.SURFNET.NL</a>] On Behalf Of Bonin, Claudia (NIH/NINDS) [F]<br>
Sent: Wednesday, September 16, 2009 1:06 PM<br>
To: <a href="mailto:FIELDTRIP@NIC.SURFNET.NL">FIELDTRIP@NIC.SURFNET.NL</a><br>
</div><div><div></div><div class="h5">Subject: Re: [FIELDTRIP] freqstastics script<br>
<br>
Hi. I am having trouble getting a freqstastics script to run to do a Kruskalwallis test. Thre freqgrandaverage works but I am having trouble with the arguments for the Kruskalwallis. I have pasted in the script and the error message in this email. Has anyone written a script using freqsstatistics to do Kruskalwallis? I would love to see it or get some direction on the arguments that need to be assigned....Also I made up an argument for cfg.design but I don't really understand what this assignment is. Could anyone elaborate?<br>
Thank you!<br>
claudia Bonin<br>
<br>
<br>
<br>
<br>
<br>
here's the error message<br>
<br>
isnan(data.time)<br>
??? Undefined function or method 'isnan' for input arguments of type 'cell'.<br>
<br>
computing average power over 21 subjects<br>
Warning: discarding gradiometer information<br>
> In freqgrandaverage at 121<br>
In tf2_virt_all_plot at 11<br>
<br>
here's the script.....<br>
<br>
<br>
Nsubj = 21;<br>
<br>
for subj = 1:Nsubj<br>
fname = ['/home/boninc/tfresults/subj' 64 + subj];<br>
clear fr<br>
load(fname);<br>
frs(subj) = fr;<br>
end<br>
<br>
cfg = [];<br>
[gr] = freqgrandaverage(cfg, frs(1), frs(2), frs(3), frs(4), frs(5), frs(6), frs(7), frs(8), frs(9), frs(10), frs(11), frs(12), frs(13), frs(14), frs(15), frs(16), frs(17), frs(18), frs(19), frs(20), frs(21));<br>
<br>
%cfg = [];<br>
%cfg.alpha = .05;<br>
%cfg.method = 'stats';<br>
%cfg.design = 1:21;<br>
%cfg.statistic = 'kruskalwallis';<br>
%[stat] = freqstatistics(cfg, frs(1), frs(2), frs(3), frs(4), frs(5), frs(6), frs(7), frs(8), frs(9), frs(10), frs(11), frs(12), frs(13), frs(14), frs(15), frs(16), frs(17), frs(18), frs(19), frs(20), frs(21));<br>
<br>
clf<br>
contourf(gr.time, gr.freq, squeeze(sum(gr.powspctrm))); colorbar<br>
%contourf(gr.time, gr.freq, gr.powspctrm); colorbar<br>
<br>
title '21 Subjects, RandL, RIGHT Central Motor'<br>
xlabel 'time, sec'<br>
ylabel 'frequency, Hz'<br>
<br>
<br>
________________________________________<br>
From: FieldTrip discussion list [<a href="mailto:FIELDTRIP@NIC.SURFNET.NL">FIELDTRIP@NIC.SURFNET.NL</a>] On Behalf Of Hanneke Van Dijk [<a href="mailto:Hanneke.vanDijk@MED.UNI-DUESSELDORF.DE">Hanneke.vanDijk@MED.UNI-DUESSELDORF.DE</a>]<br>
Sent: Thursday, September 10, 2009 7:21 AM<br>
To: <a href="mailto:FIELDTRIP@NIC.SURFNET.NL">FIELDTRIP@NIC.SURFNET.NL</a><br>
Subject: [FIELDTRIP] AW: [FIELDTRIP] freqstastics script<br>
<br>
Dear Claudia,<br>
<br>
Normally when matlab says 'undefined function for input of type ....' it means it cannot find the specific function in the paths you have defined. You can check whether neighbourselection.m is in your fieldtrip folder, if not I have attached it (version 20090428) so you can put it there. If this doesn't fix it please let me know!<br>
<br>
Groetjes Hanneke<br>
<br>
<br>
-----Ursprüngliche Nachricht-----<br>
Von: FieldTrip discussion list im Auftrag von Bonin, Claudia (NIH/NINDS) [F]<br>
Gesendet: Mi 09.09.2009 18:26<br>
An: <a href="mailto:FIELDTRIP@NIC.SURFNET.NL">FIELDTRIP@NIC.SURFNET.NL</a><br>
Betreff: [FIELDTRIP] freqstastics script<br>
<br>
I am trying to do statistics on time-frequency data across 16 subjects, averaging over channels. I received back the following error about 'neighbourselection' which means nothing to me.....Does anyone have any idea what might be wrong here? I have pasted in the exact message below.<br>
thanks!<br>
claudia<br>
<br>
<br>
>> cfg.statistic = 'kruskalwallis';<br>
>> cfg.method = 'stats';<br>
>> [s] = freqstatistics(cfg, fr, fr);<br>
averaging over 3 channels<br>
selected 50 time bins<br>
selected 41 frequency bins<br>
Warning: PACK can only be used from the MATLAB command line.<br>
> In fieldtrip/private/prepare_timefreq_data at 305<br>
In fieldtrip/private/statistics_wrapper at 167<br>
In freqstatistics at 44<br>
??? Undefined function or method 'neighbourselection' for input arguments of type<br>
'struct'.<br>
<br>
Error in ==> statistics_wrapper at 181<br>
cfg = neighbourselection(cfg,varargin{1});<br>
<br>
Error in ==> freqstatistics at 44<br>
[stat] = statistics_wrapper(cfg, varargin{:});<br>
----------------------------------<br>
The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis. See also <a href="http://listserv.surfnet.nl/archives/fieldtrip.html" target="_blank">http://listserv.surfnet.nl/archives/fieldtrip.html</a> and <a href="http://www.ru.nl/neuroimaging/fieldtrip" target="_blank">http://www.ru.nl/neuroimaging/fieldtrip</a>.<br>
<br>
<br>
<br>
----------------------------------<br>
The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis. See also <a href="http://listserv.surfnet.nl/archives/fieldtrip.html" target="_blank">http://listserv.surfnet.nl/archives/fieldtrip.html</a> and <a href="http://www.ru.nl/neuroimaging/fieldtrip" target="_blank">http://www.ru.nl/neuroimaging/fieldtrip</a>.<br>
<br>
----------------------------------<br>
The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis. See also <a href="http://listserv.surfnet.nl/archives/fieldtrip.html" target="_blank">http://listserv.surfnet.nl/archives/fieldtrip.html</a> and <a href="http://www.ru.nl/neuroimaging/fieldtrip" target="_blank">http://www.ru.nl/neuroimaging/fieldtrip</a>.<br>
<br>
----------------------------------<br>
The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis. See also <a href="http://listserv.surfnet.nl/archives/fieldtrip.html" target="_blank">http://listserv.surfnet.nl/archives/fieldtrip.html</a> and <a href="http://www.ru.nl/neuroimaging/fieldtrip" target="_blank">http://www.ru.nl/neuroimaging/fieldtrip</a>.<br>
</div></div></blockquote></div><br>
<p>----------------------------------</p>
<p>The aim of this list is to facilitate the discussion between users of the FieldTrip toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis.</p>
<p> http://listserv.surfnet.nl/archives/fieldtrip.html</p>
<p> http://www.ru.nl/fcdonders/fieldtrip/</p>