help > Nuisance regress via compcor
Showing 1-3 of 3 posts
May 18, 2015 06:05 PM | Frank Hillary - Penn State
Nuisance regress via compcor
Is there a way to access CompCor for regression of CSF and white
matter signal only? That is, for this specific question, all
other preprocessing is handled elsewhere, but we want to use
CompCor to analyze signal from CSF and white matter as regressors
of no interest. How might this limited nuisance regression be
conducted within CONN? Also, can one create an "eroded" mask for
each within CONN?
Much appreciated,
Frank
Much appreciated,
Frank
May 19, 2015 08:05 PM | Alfonso Nieto-Castanon - Boston University
RE: Nuisance regress via compcor
Hi Frank,
Yes, that is perfectly possible. If you have already fully preprocessed your data simply enter your preprocessed functionals in Setup.Functional and your preprocessed structurals (in same space as your functional volumes, e.g. normalized if you have normalized your functionals) in Setup.Structural. If you have already segmented your structural volumes you may enter those gray/white/csf masks (again in the same space as your functional and structural volumes) in Setup.ROIs (or leave those empty otherwise and CONN will run that segmentation). In Setup.Options check the 'create confound-corrected timeseries' checkbox, the 'same space as functionals' option in 'analysis space', the 'implicit masking' option in 'analysis mask', and (if your data is band-pass filtered) select the 'raw units' option in 'analysis units'. Then skip the 'Preprocessing' step and run directly the Setup step. Last, in the Denoising step select the desired denoising options (e.g. enter only the white/csf masks and remove other regressors as well as filtering and detrending if you only want to apply CompCor) and run that step. After this you will find the denoised files in conn_*/results/preprocessing/niftiDATA_Subject*_Condition*.nii (one 4f nifti file per subject/condition). You will also find the eroded masks in the same folders as your original structural or segmentation files (eroded filenames prepend 'e' to the original gray/white/csf mask filenames)
Let me know if that is what you had in mind
Best
Alfonso
EDIT: on second thoughts, if your functional data is already band-pass filtered the 'implicit masking' option is going to have trouble identifying the mask of "in-brain" voxels form the functional data alone (since that mask is computed from the average BOLD signal across your entire timeseries) so I would use instead 'explicit mask' and enter there an appropriately liberal mask (e.g. if your data is normalized to MNI space leaving the default brainmask option will probably suffice since that is quite liberal)
Originally posted by Frank Hillary:
Yes, that is perfectly possible. If you have already fully preprocessed your data simply enter your preprocessed functionals in Setup.Functional and your preprocessed structurals (in same space as your functional volumes, e.g. normalized if you have normalized your functionals) in Setup.Structural. If you have already segmented your structural volumes you may enter those gray/white/csf masks (again in the same space as your functional and structural volumes) in Setup.ROIs (or leave those empty otherwise and CONN will run that segmentation). In Setup.Options check the 'create confound-corrected timeseries' checkbox, the 'same space as functionals' option in 'analysis space', the 'implicit masking' option in 'analysis mask', and (if your data is band-pass filtered) select the 'raw units' option in 'analysis units'. Then skip the 'Preprocessing' step and run directly the Setup step. Last, in the Denoising step select the desired denoising options (e.g. enter only the white/csf masks and remove other regressors as well as filtering and detrending if you only want to apply CompCor) and run that step. After this you will find the denoised files in conn_*/results/preprocessing/niftiDATA_Subject*_Condition*.nii (one 4f nifti file per subject/condition). You will also find the eroded masks in the same folders as your original structural or segmentation files (eroded filenames prepend 'e' to the original gray/white/csf mask filenames)
Let me know if that is what you had in mind
Best
Alfonso
EDIT: on second thoughts, if your functional data is already band-pass filtered the 'implicit masking' option is going to have trouble identifying the mask of "in-brain" voxels form the functional data alone (since that mask is computed from the average BOLD signal across your entire timeseries) so I would use instead 'explicit mask' and enter there an appropriately liberal mask (e.g. if your data is normalized to MNI space leaving the default brainmask option will probably suffice since that is quite liberal)
Originally posted by Frank Hillary:
Is there a way to access CompCor for regression
of CSF and white matter signal only? That is, for this
specific question, all other preprocessing is handled elsewhere,
but we want to use CompCor to analyze signal from CSF and white
matter as regressors of no interest. How might this limited
nuisance regression be conducted within CONN? Also, can one create
an "eroded" mask for each within CONN?
Much appreciated,
Frank
Much appreciated,
Frank
May 15, 2018 02:05 PM | Moa Peter
RE: Nuisance regress via compcor
Hi Alfonso,
I am, as Frank was, interested in the denoising aspects of the CONN toolbox, preferrably using batch scripting. Your reply to Frank has been very helpful for me, but I am still struggling with defining my gm/wm/csf ROIs in a way that makes the toolbox take them as the "Grey Matter", "White Matter" and "CSF" ROIs (which would avoid redoing the segmentation). What I have tried so far is defining the segmented data using batch.Setup.rois.files and naming the files (batch.Setup.rois.name) "Grey Matter" etc, however, this only adds them as additional ROIs and the original "Grey Matter" etc ROIs are still empty.
I am a beginner and would really appreciate any kind of help!
Best regards
Moa
I am, as Frank was, interested in the denoising aspects of the CONN toolbox, preferrably using batch scripting. Your reply to Frank has been very helpful for me, but I am still struggling with defining my gm/wm/csf ROIs in a way that makes the toolbox take them as the "Grey Matter", "White Matter" and "CSF" ROIs (which would avoid redoing the segmentation). What I have tried so far is defining the segmented data using batch.Setup.rois.files and naming the files (batch.Setup.rois.name) "Grey Matter" etc, however, this only adds them as additional ROIs and the original "Grey Matter" etc ROIs are still empty.
I am a beginner and would really appreciate any kind of help!
Best regards
Moa