When you run fast_fratio, include a 3rd output (so [F pvalues ces] = ...), then the stderr for the contrast is
stderr = sqrt(F./(ces.^2));
Select the contrast to isolate each group

On 4/21/2020 5:54 PM, Graduate Imaging wrote:

External Email - Use Caution

Hi FreeSurfer experts,

I ran a vertex wise analysis that consisted of six groups with two continuous covariates which yielded a cluster after mri_glm-fit. I've extracted the values from the ocn.dat file. I ran a post hoc analyses in matlab to determine which group is significantly different using the command below:

X = load('Group.Xg.dat'); design matrix from freesurfer
C = load('C.dat'); post hoc contrast(s)
Y = load('Whole.brain.cluster.dat'); values from the ocn.dat file
[beta rvar] = fast_glmfit(Y,X)
[F pvalues] = fast_fratio(beta,X,rvar,C)

This outputs an array of beta values where the first six would be the intercept. How would I calculate the standard error of the intercept for each group? Would it be rvar/(sqroot(n)).
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer