<div dir="ltr">hi all,<div><br></div><div>I'm worried that something is not quite right with how the channel neighborhood configuration is used in the formation of spatial clusters.</div><div><br></div><div>I first noticed it when I got a significant cluster like below. Plotted are t statistics for a between-group comparison, with the significant "cluster" indicated with white electrodes. Notice F5 and P8 here which are disconnected from the rest. I should also mention that I have only one "time point" (and no frequencies), so it's not possible that these channels are somehow connected via a "time bridge" (if that would even be possible).<img src="cid:ii_150d8b6f5baecf47" alt="Inline image 3" width="545" height="436"> </div><div><br></div><div>I've inspected my channel layout and that looks good:</div><div><br></div><div><img src="cid:ii_150d91fa6ed2bfcf" alt="Inline image 4" width="545" height="461"><br></div><div>I've looked at my neighborhood plot and that looks good too, although I notice that frontal channels are plotted to the right here (not sure if that's simply by default):</div><div><br></div><div><img src="cid:ii_150d9222b043ef2c" alt="Inline image 5" width="545" height="464"><br></div><div><br></div><div>Looking into the disconnected channel F5, the neighborhood structure seems in order:</div><div><br></div><div><div><i>cfg.neighbours(15)</i></div><div><br></div><div>ans = </div><div><br></div><div> label: 'F5'</div><div> neighblabel: {'Af7' 'F3' 'F7' 'Fc5'}</div></div><div><br></div><div><br></div><div>I call ft_timelockstatistics, which then calls ft_statistics_montecarlo, which, on line 165 calls </div><div><br></div><div><i>cfg.connectivity = channelconnectivity(cfg);</i></div><div><br></div><div>Stepping into this function, on line 20 it says:</div><div><br></div><div><i>chans=cfg.channel;</i></div><div><br></div><div>Checking what this variable contains, it has all my channels in <b>alphabetical order</b>, but my electrode location order is most definitely not. I've tried to force the orders to be the same by setting cfg.channel = cfg.label before calling ft_timelockstatistics, but somewhere down the line my channel order is made alphabetical again.</div><div><br></div><div>So I suspect that the neighborhood structure, reflecting my actual channel order, is applied to the alphabetical channel order. That would explain why a cluster could be scattered across the brain (and also why there's still some spatial continuity given that alphabetically close channels are usually close together in space).</div><div><br></div><div>Has anyone ever seen this? Any suggestions how to solve this (other than making my entire channel location file alphabetical)? I'm also posting this because it may be of relevance to others.</div><div><br></div><div>Roy</div>
<div><br></div><div><br></div></div>