<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 12 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0in;
margin-bottom:.0001pt;
font-size:10.5pt;
font-family:Consolas;}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
{mso-style-priority:99;
mso-style-link:"Balloon Text Char";
margin:0in;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Tahoma","sans-serif";}
span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:Consolas;}
span.BalloonTextChar
{mso-style-name:"Balloon Text Char";
mso-style-priority:99;
mso-style-link:"Balloon Text";
font-family:"Tahoma","sans-serif";}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoPlainText>Greetings,<o:p></o:p></p><p class=MsoPlainText>Based on FT tutorials and user advice I have managed to create a stat structure using ft_freqstatistics. <o:p></o:p></p><p class=MsoPlainText>stat = <o:p></o:p></p><p class=MsoPlainText>prob [2x158x1800 double]<o:p></o:p></p><p class=MsoPlainText>posclusters: [1x43 struct]<o:p></o:p></p><p class=MsoPlainText>posclusterslabelmat: [2x158x1800 double]<o:p></o:p></p><p class=MsoPlainText>posdistribution: [1x500 double]<o:p></o:p></p><p class=MsoPlainText>negclusters: [1x63 struct]<o:p></o:p></p><p class=MsoPlainText>negclusterslabelmat: [2x158x1800 double]<o:p></o:p></p><p class=MsoPlainText>negdistribution: [1x500 double]<o:p></o:p></p><p class=MsoPlainText>mask: [2x158x1800 logical]<o:p></o:p></p><p class=MsoPlainText>stat: [2x158x1800 double]<o:p></o:p></p><p class=MsoPlainText>ref: [2x158x1800 double]<o:p></o:p></p><p class=MsoPlainText>dimord: 'chan_freq_time'<o:p></o:p></p><p class=MsoPlainText>label: {2x1 cell} (i.e., {['MLT'];['MRT']}<o:p></o:p></p><p class=MsoPlainText>freq: [1x158 double]<o:p></o:p></p><p class=MsoPlainText>time: [1x1800 double]<o:p></o:p></p><p class=MsoPlainText>cfg: [1x1 struct]<o:p></o:p></p><p class=MsoPlainText><o:p> </o:p></p><p class=MsoPlainText>Please note that the cluster permuation tests were performed on source space time-frequency representations, thus data is only available for two channels. As such the statistics was done without cfg.minnbchan and cfg.neighbours=[ ]. and clustering was done on the time-frequency bins throuought the entire epoch for wide-band data (i.e., no averaging over frequencies). The initial conditional power spectra used in the statistics were 10 x 2 x 158 x 2880 (subjects x channels x frequencies x time). I did manage to use ft_singleplotTFR(cfg, stat) to plot the results for single channel -after masking time-frequency points by clusterlabels information belonging to a specific cluster (see attached). But I cannot really make much sense out of the results. Am I doing something wrong here? <o:p></o:p></p><p class=MsoPlainText><o:p> </o:p></p><p class=MsoPlainText>When I use <o:p></o:p></p><p class=MsoPlainText>cfg = [];<o:p></o:p></p><p class=MsoPlainText>cfg.alpha = 0.025;<o:p></o:p></p><p class=MsoPlainText>cfg.zparam = 'stat';<o:p></o:p></p><p class=MsoPlainText>cfg.zlim = [-4 4];<o:p></o:p></p><p class=MsoPlainText>ft_clusterplot(cfg, stat); <o:p></o:p></p><p class=MsoPlainText><o:p> </o:p></p><p class=MsoPlainText>I get the following:<o:p></o:p></p><p class=MsoPlainText>??? Error using ==> ft_clusterplot at 75<o:p></o:p></p><p class=MsoPlainText>stat contains multiple frequencies which is not allowed because it should be<o:p></o:p></p><p class=MsoPlainText>averaged over frequencies<o:p></o:p></p><p class=MsoPlainText><o:p> </o:p></p><p class=MsoPlainText>Technically I don't have any apriori frequency band(s) of interest, thus I used the wide-band data. How does ft_clusterplot accept frequency averaging? And is it possible to work around this criterion? Is it possible to customize cfg.layout?<o:p></o:p></p><p class=MsoPlainText>Thanks in advance -Kambiz<o:p></o:p></p><p class=MsoPlainText><o:p> </o:p></p><p class=MsoPlainText><o:p> </o:p></p><p class=MsoPlainText>-K.<o:p></o:p></p><p class=MsoPlainText>-------------------------------------------------------------<o:p></o:p></p><p class=MsoPlainText>Kambiz Tavabi PhD<o:p></o:p></p><p class=MsoPlainText>Biomedical Imaging Laboratory<o:p></o:p></p><p class=MsoPlainText>The Children's Hospital of Philadelphia<o:p></o:p></p><p class=MsoPlainText>34th Street and Civic Center Boulevard<o:p></o:p></p><p class=MsoPlainText>Philadelphia, Pa. 19104<o:p></o:p></p><p class=MsoPlainText>Tel: 267.426.0302<o:p></o:p></p><p class=MsoPlainText>email: <a href="mailto:tavabik@email.chop.edu">tavabik@email.chop.edu</a><o:p></o:p></p><p class=MsoPlainText>-------------------------------------------------------------<o:p></o:p></p></div></body></html>