In the qdec output, there should be a file called Xg.dat and one called
contrast/C.dat (where contrast is the name of the contrast of interest)
X = load('Xg.dat');
C = load('contrast/C.dat');
y = MRIread('input-to-qdec.mgz'); % Not sure what this is called
ymat = fast_vol2mat(y);
yvertex = ymat(:,vertexno+1);
[beta rvar] = fast_glmfit(yvertex,X);
[F p] = fast_fratio(beta,X,rvar,C)
On 6/8/2020 3:55 AM, Jian Ou wrote:
External Email - Use Caution
Hello FreeSurfer Developers, I have a few questions about the p value
in QDEC. I want to use matlab to calculate the p value in a vertex in
the QDEC result. 1. In my design, there is a continuous variable A and
two Nuisance factors B and C. How to calculate the p value in the
correlation analysis of A and thickness. 2. If there are two groups of
subjects F, M. How to calculate the p value.
Thank for your help,
Jian Ou
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer