help > RE: Masking
Jul 29, 2014  02:07 AM | Alfonso Nieto-Castanon - Boston University
RE: Masking
Dear Julia,

First regarding the error message thanks for reporting this, this was actually a bug in rex and I have fixed that for the next release of the toolbox (in the meantime I am attaching a patch for version CONN14i, simply copy the attached file in your conn distribution folder overwriting the file with the same name there). 

Regarding question (1), it is not necessary to orthogonalize your covariate to the 'All' covariate (the results will be identical whether you do or not), and your regression analyses looking at the effect of your covariate are always considered orthogonal to the whole sample average effects, so you do not have any circularity issues there.

Regarding question (2), rather than simply using more liberal thresholds for your new/limited search volume, the more appropriate way to do this would be to use 'SVC' (small volume correction) and selecting your mask image there to obtain statistics that already take into account this smaller search volume (e.g. DMN areas) when applying the corresponding multiple-comparison corrections (e.g. FWE- p-values)

Regarding question (3), if you extract your data from regions that have been defined by the same (or a non-orthogonal) contrast that you are then trying to evaluate (and from your description I take that this was the case here), then yes, you will have a circularity issue there, and you may assume that the resulting measures will be biased (overestimating effect sizes, under-estimating p-values, etc.) Some times you wish to do just that, for example as a descriptive post hoc analysis, and that is perfectly fine (you just need to acknowledge the limitations and potential biases of the resulting measures, and do not use these tests as any form of confirmatory analyses). If your goal, on the other hand, is to obtain accurate effect-size measures within your regions of interest then I would probably suggest using a cross-validation approach to do this. 

Hope this helps
Alfonso



Originally posted by Julia Landsiedel:
Dear Alfonso,

Thank you for your quick reply and the detailed explanation. 

Regarding the usage of REX I encountered the problem that if I start it from the command line that it gives error message when I want to save the output data, presumably because it does not have an output folder defined (and doesn't use the current folder either). This is no problem if I start REX from the seed-to voxel explorer with 'Explore Clusters'

This would be the error message:

Reference to non-existent field 'output_folder'.
Error in rex>rex_do (line 789)
name_dat=fullfile(params.output_folder,[ROInames{rr},'.rex.data.mat']);
Error in rex>rex_gui (line 444)
[data.params.ROIdata,data.params.ROInames,data.params.ROIinfo.basis,data.params.ROIinfo.voxels,data.params.ROIinfo.files,data.params.ROIinfo.select,data.params.ROIinfo.trans]=rex_do(data,~dat
Error while evaluating uicontrol Callback

*****
Besides just three quick follow up questions.

1) I create my DMN mask based on the results of the whole sample. Do I have to orthogonalize my covariate in the Setup.2nd level covariates tab two the 'All' covariate or does this already avoid the circularity problem?

2) After masking is it 'allowed' to use lower height and extend thresholds? Do you have any recommendation for this?

3) With REX, I did not extract the data with the default mode network (DMN) mask but with a mask of only the significant clusters from my DMN-masked covariate analysis in SPM. Would this bias my results?

Best,
Julia
Attachment: rex.m

Threaded View

TitleAuthorDate
Julia Landsiedel Jul 21, 2014
Alfonso Nieto-Castanon Jul 22, 2014
Kaitlin Cassady Sep 5, 2023
Mayron Pereira Picolo Ribeiro Jul 10, 2020
Julia Landsiedel Jul 22, 2014
RE: Masking
Alfonso Nieto-Castanon Jul 29, 2014
Julia Landsiedel Jul 30, 2014
Yifei Zhang Aug 1, 2014
Yifei Zhang Aug 5, 2014
Alfonso Nieto-Castanon Aug 6, 2014