<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=utf-8">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@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:"Segoe UI";
panose-1:2 11 5 2 4 2 4 2 2 3;}
@font-face
{font-family:"\@SimSun";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","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;}
span.EmailStyle17
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
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="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Dear Ronny,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">What error do you get when running this code?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">You can perform a (relative) baseline correction using ft_freqbaseline with cfg.baselinetype = ‘relative’ or ‘relchange’ and cfg.baseline = [your baseline window];.
If your code works on absolute values it should also work on baseline-corrected values using ft_freqbaseline.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Best,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Jim<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> fieldtrip-bounces@science.ru.nl [mailto:fieldtrip-bounces@science.ru.nl]
<b>On Behalf Of </b>Ronny Ibrahim<br>
<b>Sent:</b> maandag 7 december 2015 5:36<br>
<b>To:</b> fieldtrip@science.ru.nl<br>
<b>Subject:</b> [FieldTrip] Clusterstatistics on relative change (%) time frequency data<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Dear fellow fieldtrip users,<o:p></o:p></p>
<div>
<p class="MsoNormal">I would like to implement a clusterstatistics independentT test on a time-frequency datasets that I've acquired using the ft_freqanalysis function. When I run ft_freqstatistics on the data it would give me stats result which is based on
comparing the absolute values of power spectrum instead the relative change with respect to the baseline. I have tried to perform to do my relative change power spectrum data manually but it seems that when I run the ft_freqstatistics it produces an error
of which I could not figure out on why.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Here is my code:<span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black"> </span><o:p></o:p></p>
</div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black"> <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">C_1 = {outP10{1} outP10{2} outP10{3} outP10{4} outP10{5} outP10{6} outP10{7} outP10{8} outP10{9} outP10{10} outP10{11} outP10{12} outP10{13} outP10{14} outP10{15}
outP10{16} outP10{17} outP10{18}}; <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black"> <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">C_2 = {outM10{1} outM10{2} outM10{3} outM10{4} outM10{5} outM10{6} outM10{7} outM10{8} outM10{9} outM10{10} outM10{11} outM10{12} outM10{13} outM10{14} outM10{15}
outM10{16} outM10{17} outM10{18}}; <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black"> <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">%% Performing cluster statistics for the paradigm<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg = [];<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.layout = 'EEG_KR.lay';<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">lay = ft_prepare_layout(cfg);<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">% Needs to change this<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.latency = [0 0.6]; %Need to look this up according to what <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.frequency = [8 12];% Need to do it for different frequency bands (theta, alpha, beta, gamma).<o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2601"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.channel = 'all';<o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2599"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.avgovertime = 'no';<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.avgoverfreq = 'yes';<o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2597"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.avgoverchan = 'no';<o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2596"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.statistic = 'ft_statfun_indepsamplesT';<o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2605"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.numrandomization = 1000;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.correctm = 'cluster';<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.method = 'montecarlo';<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.computeprob = 'yes';<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.computecritval = 'yes';<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.clusteralpha = 0.05;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.clusterstatistic = 'maxsum';<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.minnbchan = 3;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.tail = 0;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.clustertail = 0;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.alpha = 0.05;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg_neighb.layout = lay;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg_neighb.method = 'triangulation';<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.neighbours = ft_prepare_neighbours(cfg_neighb); <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.ivar = 1;<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.design = [ones(1,18) 2*ones(1,18)];<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">stat = ft_freqstatistics(cfg, C_2{:}, C_1{:});<o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2616"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black"> <o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2614"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">%% Plotting significant clusters<o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2613"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black"> <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg = [];<o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2611"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.alpha = 0.05;<o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2610"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.parameter = 'stat';<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">cfg.layout = lay;<o:p></o:p></span></p>
<p class="MsoNormal" id="yui_3_16_0_1_1449462105095_2621"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black">ft_clusterplot(cfg, stat);<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif";color:black"> </span> <br clear="all">
<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">where the variable out is essentially my time-frequency data output structure. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I would greatly appreciate the input and assistance. Thank you.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Kind Regards,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Ronny<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</div>
</body>
</html>