<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body bgcolor="#ffffff" text="#000000">
<font face="Georgia">Hi Fieldtrippers,<br>
<br>
I am trying to run a comparison between TFRs from two groups of
subjects and I keep getting some errors that I can't track down
answers for. I seem to recall something like this being discussed
on this list, but I can't for the life of me track down that
thread at the moment. <br>
<br>
I am using data from a neuromag 306 system on which I have
calculated TFRs on the gradiometer data, then ran
ft_combineplanar, followed by ft_freqdescriptives, followed by
ft_freqgrandaverage for each of the two groups with the
cfg.keepindividual = 'yes' option. I am able to plot the grand
averages just fine, so the data seem to be ok (and there are
fairly striking differences between the two). <br>
<br>
When I run a permutation test to quantify the group differences, I
get the following error messages (although the analysis does run):<br>
<br>
<i>computing statistic 500 from 500<br>
Warning: Not all replications are used for the computation of
the statistic.<br>
> > In statfun_indepsamplesT at 57<br>
In statistics_montecarlo at 298<br>
In fieldtrip-20100617/private/statistics_wrapper at 285<br>
In ft_freqstatistics at 105<br>
Warning: Divide by zero.</i><br>
<br>
When I inspect the output file, the .stat contains a single column
of NaNs. Does anyone know what I am doing wrong here? <br>
<br>
Here is the actual code I am using to run the statistics:<br>
<br>
cfg = [];<br>
cfg.channel = {'MEG *3'}; %only combined gradiometer data
<br>
cfg.layout = 'neuromag306cmb_dz.lay';<br>
cfg.latency = [0.2 0.8];<br>
cfg.avgovertime = 'yes';<br>
cfg.avgoverchan = 'no';<br>
cfg.avgoverfreq = 'yes';<br>
cfg.parameter = 'powspctrm';<br>
cfg.frequency = [14 26];<br>
cfg.method = 'montecarlo';<br>
cfg.statistic = 'indepsamplesT';<br>
cfg.correctm = 'no';<br>
cfg.clusteralpha = 0.05;<br>
cfg.clusterstatistic = 'maxsum';<br>
cfg.minnbchan = 2;<br>
cfg.tail = 0;<br>
cfg.clustertail = 0;<br>
cfg.alpha = 0.05;<br>
cfg.numrandomization = 500;<br>
<br>
cfg.design = [1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
18 19 20 21 22 23 24 25 ; % subject number<br>
1 1 1 1 1 1 1 1 1 1 1 1 2 2 2
2 2 2 2 2 2 2 2 2 2]; % 12subjs in grp1, 13subjs in grp2<br>
cfg.uvar = 1; <br>
cfg.ivar = 2; <br>
<br>
GA_low_stat = ft_freqstatistics(cfg, GA_grp1_low, GA_grp2_low)<br>
<br>
Thanks,<br>
David<br>
<br>
<br>
</font>
<div class="moz-signature">-- <br>
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<meta name="ProgId" content="Word.Document">
<meta name="Generator" content="Microsoft Word 11">
<meta name="Originator" content="Microsoft Word 11">
<link rel="File-List" href="email%20signature_files/filelist.xml">
<title>David Ziegler</title>
<o:smarttagtype
namespaceuri="urn:schemas-microsoft-com:office:smarttags"
name="PostalCode">
<o:smarttagtype
namespaceuri="urn:schemas-microsoft-com:office:smarttags"
name="State">
<o:smarttagtype
namespaceuri="urn:schemas-microsoft-com:office:smarttags"
name="City">
<o:smarttagtype
namespaceuri="urn:schemas-microsoft-com:office:smarttags"
name="place">
<o:smarttagtype
namespaceuri="urn:schemas-microsoft-com:office:smarttags"
name="Street">
<o:smarttagtype
namespaceuri="urn:schemas-microsoft-com:office:smarttags"
name="address">
<!--[if gte mso 9]><xml>
<o:DocumentProperties>
<o:Author>daz</o:Author>
<o:Template>Normal</o:Template>
<o:LastAuthor>daz</o:LastAuthor>
<o:Revision>8</o:Revision>
<o:TotalTime>14</o:TotalTime>
<o:Created>2005-05-26T14:55:00Z</o:Created>
<o:LastSaved>2008-09-13T02:21:00Z</o:LastSaved>
<o:Pages>1</o:Pages>
<o:Words>34</o:Words>
<o:Characters>196</o:Characters>
<o:Company>mit</o:Company>
<o:Lines>1</o:Lines>
<o:Paragraphs>1</o:Paragraphs>
<o:CharactersWithSpaces>229</o:CharactersWithSpaces>
<o:Version>11.5606</o:Version>
</o:DocumentProperties>
</xml><![endif]--><!--[if gte mso 9]><xml>
<w:WordDocument>
<w:SpellingState>Clean</w:SpellingState>
<w:GrammarState>Clean</w:GrammarState>
<w:DrawingGridHorizontalSpacing>9.35 pt</w:DrawingGridHorizontalSpacing>
<w:DisplayVerticalDrawingGridEvery>2</w:DisplayVerticalDrawingGridEvery>
<w:ValidateAgainstSchemas/>
<w:SaveIfXMLInvalid>false</w:SaveIfXMLInvalid>
<w:IgnoreMixedContent>false</w:IgnoreMixedContent>
<w:AlwaysShowPlaceholderText>false</w:AlwaysShowPlaceholderText>
<w:Compatibility>
<w:BreakWrappedTables/>
<w:SnapToGridInCell/>
<w:WrapTextWithPunct/>
<w:UseAsianBreakRules/>
<w:DontGrowAutofit/>
<w:UseFELayout/>
<w:SplitPgBreakAndParaMark/>
<w:DontVertAlignCellWithSp/>
<w:DontBreakConstrainedForcedTables/>
<w:DontVertAlignInTxbx/>
<w:Word11KerningPairs/>
<w:CachedColBalance/>
</w:Compatibility>
<w:BrowserLevel>MicrosoftInternetExplorer4</w:BrowserLevel>
<w:TrackMoves>false</w:TrackMoves>
<w:TrackFormatting/>
<w:DoNotPromoteQF/>
<w:LidThemeOther>EN-US</w:LidThemeOther>
<w:LidThemeAsian>X-NONE</w:LidThemeAsian>
<w:LidThemeComplexScript>X-NONE</w:LidThemeComplexScript>
<m:mathPr>
<m:mathFont m:val="Cambria Math"/>
<m:brkBin m:val="before"/>
<m:brkBinSub m:val="--"/>
<m:smallFrac m:val="off"/>
<m:dispDef/>
<m:lMargin m:val="0"/>
<m:rMargin m:val="0"/>
<m:defJc m:val="centerGroup"/>
<m:wrapIndent m:val="1440"/>
<m:intLim m:val="subSup"/>
<m:naryLim m:val="undOvr"/>
</m:mathPr>
</w:WordDocument>
</xml><![endif]--><!--[if gte mso 9]><xml>
<w:LatentStyles DefLockedState="false" LatentStyleCount="156">
</w:LatentStyles>
</xml><![endif]--><!--[if !mso]><object
classid="clsid:38481807-CA0E-42D2-BF39-B33AF135CC4D" id=ieooui></object>
<style>
st1\:*{behavior:url(#ieooui) }
</style>
<![endif]-->
<style>
<!--p.MSONORMAL
{mso-style-unhide:no;
mso-style-qformat:yes;}
li.MSONORMAL
{mso-style-unhide:no;
mso-style-qformat:yes;}
div.MSONORMAL
{mso-style-unhide:no;
mso-style-qformat:yes;}
a:link
{mso-style-priority:99;}
span.MSOHYPERLINK
{mso-style-priority:99;}
a:visited
{mso-style-priority:99;}
span.MSOHYPERLINKFOLLOWED
{mso-style-priority:99;}
.MSOCHPDEFAULT
{mso-default-props:yes;}
table.MSONORMALTABLE
{mso-style-priority:99;
mso-style-qformat:yes;}
/* Font Definitions */
@font-face
{font-family:Garamond;
panose-1:2 2 4 4 3 3 1 1 8 3;
mso-font-charset:0;
mso-generic-font-family:roman;
mso-font-pitch:variable;
mso-font-signature:647 0 0 0 159 0;}
@font-face
{font-family:"Bodoni MT";
panose-1:2 7 6 3 8 6 6 2 2 3;
mso-font-charset:0;
mso-generic-font-family:roman;
mso-font-pitch:variable;
mso-font-signature:3 0 0 0 1 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{mso-style-parent:"";
margin:0in;
margin-bottom:.0001pt;
mso-pagination:widow-orphan;
font-size:12.0pt;
font-family:"Times New Roman";
mso-fareast-font-family:"Times New Roman";}
a:link, span.MsoHyperlink
{mso-style-noshow:yes;
color:blue;
text-decoration:underline;
text-underline:single;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-noshow:yes;
color:purple;
text-decoration:underline;
text-underline:single;}
span.GramE
{mso-style-name:"";
mso-gram-e:yes;}
@page Section1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;
mso-header-margin:.5in;
mso-footer-margin:.5in;
mso-paper-source:0;}
div.Section1
{page:Section1;}
-->
</style><!--[if gte mso 10]>
<style>
/* Style Definitions */
table.MsoNormalTable
{mso-style-name:"Table Normal";
mso-tstyle-rowband-size:0;
mso-tstyle-colband-size:0;
mso-style-noshow:yes;
mso-style-parent:"";
mso-padding-alt:0in 5.4pt 0in 5.4pt;
mso-para-margin:0in;
mso-para-margin-bottom:.0001pt;
mso-pagination:widow-orphan;
font-size:10.0pt;
font-family:"Times New Roman";
mso-fareast-font-family:"Times New Roman";
mso-ansi-language:#0400;
mso-fareast-language:#0400;
mso-bidi-language:#0400;}
</style>
<![endif]-->
<link rel="themeData"
href="email%20signature_files/themedata.thmx">
<link rel="colorSchemeMapping"
href="email%20signature_files/colorschememapping.xml">
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="3074"/>
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1"/>
</o:shapelayout></xml><![endif]-->
<div class="Section1">
<p class="MsoNormal"><span style="font-size: 11pt;
font-family: "Bodoni MT";">David
A. Ziegler<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 11pt;
font-family: "Bodoni MT";">Department
of Brain and Cognitive Sciences <br>
Massachusetts Institute of Technology <br>
<st1:street w:st="on"><st1:address w:st="on">43
Vassar St</st1:address></st1:street><span
class="GramE">,<span style=""> </span>46</span>-5121
<br>
<st1:place w:st="on"><st1:city w:st="on">Cambridge</st1:city>,
<st1:state w:st="on">MA</st1:state> <st1:postalcode
w:st="on">02139</st1:postalcode></st1:place><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 11pt;
font-family: "Bodoni MT";">Tel:
617-258-0765<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 11pt;
font-family: "Bodoni MT";">Fax:
617-253-1504<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 11pt;
font-family: "Bodoni MT";"><a
href="mailto:daz@mit.edu">daz@mit.edu</a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:
Garamond;"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-family:
Garamond;"><o:p> </o:p></span></p>
</div>
</o:smarttagtype></o:smarttagtype></o:smarttagtype></o:smarttagtype></o:smarttagtype></o:smarttagtype></div>
</body>
</html>
<p>---------------------------------------------------------------------------
You are receiving this message because you are subscribed to <br>
the FieldTrip list. The aim of this list is to facilitate the discussion <br>
between users of the FieldTrip toolbox, to share experiences <br>
and to discuss new ideas for MEG and EEG analysis. <br>
See also http://listserv.surfnet.nl/archives/fieldtrip.html <br>
and http://www.ru.nl/neuroimaging/fieldtrip.<br>
---------------------------------------------------------------------------</p>