<!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">
    Hi all.<br>
    <br>
    I've been trying to define a region of interest (ROI) in order to
    restrict my analyses to a particular region on the temporal lobe
    during sourcestatistics.<br>
    <br>
    I have two grandaverage sources like:<br>
    <small>pos: [33480x3 double]<br>
      dim: [31 36 30]<br>
      avg: [1x1 struct]<br>
      var: [1x1 struct]<br>
      dimord: 'voxel'<br>
      trial: [1x13 struct]<br>
      inside: [12773x1 double]<br>
      outside: [20707x1 double]<br>
      df: [33480x1 double]<br>
      cfg: [1x1 struct]</small><br>
    <br>
    And I apply ft_sourcestatistics with the following configuration:<br>
    <small>cfg             =[];<br>
      cfg.dim         = GA_dev_source.dim;<br>
      cfg.method      = 'montecarlo';<br>
      cfg.statistic   = 'depsamplesT';<br>
      cfg.parameter   = 'pow';<br>
      cfg.correctm    = 'bonferoni';<br>
      cfg.numrandomization = 1000;<br>
      cfg.alpha       = 0.05;<br>
      cfg.tail        = 0;<br>
      nsubj           = length(GA_dev_source.trial);<br>
      cfg.design(1,:) = [1:nsubj 1:nsubj];<br>
      cfg.design(2,:) = [ones(1,nsubj) ones(1,nsubj)*2];<br>
      cfg.uvar        = 1;<br>
      cfg.ivar        = 2;</small><br>
    <small>stat = ft_sourcestatistics(cfg, GA_dev_source,
      GA_stad_source)</small><br>
    <br>
    However it's not very straightforward to me to apply a ROI
    restriction to this analysis.<br>
    I downloaded the files from <a moz-do-not-send="true"
      class="moz-txt-link-freetext" href="http://fmri.wfubmc.edu">http://fmri.wfubmc.edu</a>
    (TD_lobe.img, TD_lobe.mat, TD_lobe_Border.mat, TD_lobe_List.mat)
    which I am suposed to use for the ROI.<br>
    <br>
    I tried to add this to the previous configuration:<br>
    <small>cfg.atlas = '.../TD_lobe.img';        (the template brain
      containing the different lobes)<br>
      cfg.inputcoord = 'mni';<br>
      cfg.roi = ROI(12).Nom_C;            The names of the roi i want to
      extract ('Temporal Lobe')</small><br>
    <br>
    But then...<br>
    <small>stat_roi = <br>
      prob: [31x36x30 double]    <---- All zeros.<br>
      mask: [31x36x30 logical]<br>
      stat: [31x36x30 double]<br>
      ref: [31x36x30 double]<br>
      dim: [31 36 30]<br>
      inside: [0x1 double]       <---- which means there's nothing
      inside the volume<br>
      outside: [33480x1 double]<br>
      pos: [33480x3 double]<br>
      cfg: [1x1 struct]<br>
      <big><br>
        Amb I defining the ROI incorrectly??? Is it because the
        dimensions of the atlas and the source are not equal? Could
        anyone give a good example of how to do it?<br>
        <br>
        Any suggestion will be appreciated. Thanks in advance.</big></small>
    <pre class="moz-signature" cols="72">-- 
Marc Recasens
Tel.: +34 639 24 15 98</pre>
  </body>
</html>