Hello,
I ran the following analysis to target task-specific voxels at a certain contrast and offset.
func2roi-sess -roidef dACC10_lh -analysis EMerror -anatlabel dACC-lh -maskcontrast ASevfix -maskframe 5 -maskthresh 1 -masktail pos -maskmap sig -sf Subject_files/MTHFRn18-list -d $SUBJECTS_DIR
I would like to create a functional constraint including only the 25% most active voxels (top quartile).
Is there a command that can achieve this?
Thanks,
Dave Brohawn
Sorry, I don't have anything that will do that. You could probably do it in matlab, if you know how to program in matlab.
doug
Dave Brohawn wrote:
Hello,
I ran the following analysis to target task-specific voxels at a certain contrast and offset.
func2roi-sess -roidef dACC10_lh -analysis EMerror -anatlabel dACC-lh -maskcontrast ASevfix -maskframe 5 -maskthresh 1 -masktail pos -maskmap sig -sf Subject_files/MTHFRn18-list -d $SUBJECTS_DIR
I would like to create a functional constraint including only the 25% most active voxels (top quartile).
Is there a command that can achieve this?
Thanks,
Dave Brohawn _______________________________________________ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
On Thursday 08 January 2009 15:29:47 Dave Brohawn wrote:
I would like to create a functional constraint including only the 25% most active voxels (top quartile).
Is there a command that can achieve this?
You can use 'mris_calc' (the development version):
$>mris_calc someVolume.mgz stats
--> will dump min/max/mean/std data
Based on the range max-min determine the 25% threshold, then:
$>mris_calc -o top25pc.mgz someVolume.mgz gte <threshold>
--> top25pc.mgz will contain only the voxels in someVolume.mgz with intensities greater than or equal to <threshold>.
-=R
freesurfer@nmr.mgh.harvard.edu