FWE-corrected p values in peak_nii

94 views
Skip to first unread message

Monica Leslie

unread,
Sep 4, 2019, 9:42:04 AM9/4/19
to Aaron's fMRI matlab tools.
Hi Aaron and Donald,

I have just finished using GLM Flex Fast 4 and running peak_nii to identify significant clusters in my data. After perusing JISCMail and this group, I have seen several discussions which refer to both corrected and uncorrected p values in the output from peak_nii; however, I don't seem to be getting any p values in my output at all. 

Peak_nii has been helpful in identifying clusters in my data at the p < .001 uncorrected level and yields T statistics for peaks, but I'm not sure how to generate or find FWE-corrected p values.

I have been using the peak_nii script from the following webpage (https://www.nitrc.org/projects/peak_nii), but I would be very grateful if you could point me in the direction of an updated script or an additional option I should use to generate FWE-corrected p values.

Many thanks,

Monica

Monica Leslie

unread,
Sep 9, 2019, 5:21:13 AM9/9/19
to Aaron's fMRI matlab tools.
Hi All,

After a little more trouble-shooting it seems that whether I use peak_nii directly from the command line calling for [voxels voxelstats clusterstats sigthresh regions map parameters UID] or use FIVE to get Peak Info I get the following error message:

Error using peak_nii
Too many output arguments.

Has anyone else come across this before or know what might be causing the issue? For the moment, it seems that this is preventing me from being able to view FWE-corrected p values for my data.

Many thanks,

Monica

MCLAREN, Donald

unread,
Sep 9, 2019, 9:19:17 PM9/9/19
to fMRI MATLAB Tools
Hi Monica,

I'll try to do some testing this week to identify the issue and the patch.

Best Regards, 
Donald McLaren, PhD



--
You received this message because you are subscribed to the Google Groups "Aaron's fMRI matlab tools." group.
To unsubscribe from this group and stop receiving emails from it, send an email to fmri_matlab_to...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/fmri_matlab_tools/69bb4735-a2ae-40e9-99de-a46cf0a412aa%40googlegroups.com.

Monica Leslie

unread,
Sep 10, 2019, 5:08:26 AM9/10/19
to Aaron's fMRI matlab tools.
Hi Donald,

That would be amazing. Thank you very much.

Monica

On Tuesday, 10 September 2019 02:19:17 UTC+1, Donald McLaren wrote:
Hi Monica,

I'll try to do some testing this week to identify the issue and the patch.

Best Regards, 
Donald McLaren, PhD



On Wed, Sep 4, 2019 at 9:42 AM Monica Leslie <monica...@gmail.com> wrote:
Hi Aaron and Donald,

I have just finished using GLM Flex Fast 4 and running peak_nii to identify significant clusters in my data. After perusing JISCMail and this group, I have seen several discussions which refer to both corrected and uncorrected p values in the output from peak_nii; however, I don't seem to be getting any p values in my output at all. 

Peak_nii has been helpful in identifying clusters in my data at the p < .001 uncorrected level and yields T statistics for peaks, but I'm not sure how to generate or find FWE-corrected p values.

I have been using the peak_nii script from the following webpage (https://www.nitrc.org/projects/peak_nii), but I would be very grateful if you could point me in the direction of an updated script or an additional option I should use to generate FWE-corrected p values.

Many thanks,

Monica

--
You received this message because you are subscribed to the Google Groups "Aaron's fMRI matlab tools." group.
To unsubscribe from this group and stop receiving emails from it, send an email to fmri_mat...@googlegroups.com.

Monica Leslie

unread,
Sep 26, 2019, 10:35:43 AM9/26/19
to Aaron's fMRI matlab tools.
Hi Donald,

I just wanted to follow up to see if you've had a chance to see what might be causing the "Too Many Output Arguments" issue with peak_nii? Unfortunately I still haven't found a way around it on my end.

Kind Regards,

Monica

MCLAREN, Donald

unread,
Sep 28, 2019, 5:39:18 PM9/28/19
to fMRI MATLAB Tools
Hi Monica,

[voxels voxelstats clusterstats sigthresh regions map parameters UID] 

should be

[voxels voxelstats clusterstats sigthresh regions mapparameters UID] 

sigthresh has the threshold values in them.

voxelstats column 1 should have the FWE-corrected values. This requires that you provide the FWHM of the image. This can be found in the SPM.mat file.

I'm also attaching an updated version of peak_nii.m

Let me know if this doesn't work and I can look again.

Best Regards, 
Donald McLaren, PhD


--
You received this message because you are subscribed to the Google Groups "Aaron's fMRI matlab tools." group.
To unsubscribe from this group and stop receiving emails from it, send an email to fmri_matlab_to...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/fmri_matlab_tools/75ef184a-fc48-4941-8ff5-8176c5ca51e3%40googlegroups.com.
peak_nii.m
Reply all
Reply to author
Forward
0 new messages