You'll probably have to use matlab, eg, [vertices, label, colortable] = read_annotation(filename) ind = find(label==yourlabel); beta = fast_vol2mat(MRIread('beta.nii.gz')); yourlabelbeta = beta(:,ind);
On 3/6/2023 10:08 AM, Huang, Yujing wrote:
Hi Andy,
Can you try the mask option?
https://surfer.nmr.mgh.harvard.edu/fswiki/mri_glmfit
Best,
Yujing
*From:* freesurfer-bounces@nmr.mgh.harvard.edu freesurfer-bounces@nmr.mgh.harvard.edu *On Behalf Of *Andy Kim *Sent:* Tuesday, February 28, 2023 10:47 PM *To:* freesurfer freesurfer@nmr.mgh.harvard.edu *Subject:* [Freesurfer] How to get individual betas (parameter estimates) in certain ROI?
* External Email - Use Caution *
Dear FreeSurfer Developers,
Hi. Thanks to FreeSurfer, I finished a group activation analysis with FS-FAST.
Now I want to draw a ROC curve to check how well my task conditions can distinguish M1 and S1. (M1 is more activated under certain condition of my task and S1 is more activated under another condition.)
To do so, I need to get every betas of individual nodes in certain ROIs (which I created my own, but let's just say M1 and S1).
Betas are in beta.mgh file which is a product of mri_glmfit. Then how do I get betas of only M1 nodes or only S1 nodes?
Kind Regards,
Andy (Junghyun) Kim
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer