Dear Martin,
Yes I think so, and I run the mri_glmfit with fsaverage as the subject like:mri_glmfit --glmdir lh.thickness.rmanova --y lh.thickness.sm05.mgh --fsgd rmanova.fsgd doss --C tp1-vs-tp2.mtx --C tp1-vs-tp3.mtx --C tp2-vs-tp3.mtx --C tp-effect.mtx --C mean.mtx --surf fsaverage lh
it reports:
FixVertexAreaFlag = 1UseMaskWithSmoothing 1OneSampleGroupMean 0y /Freesurfer/k/lh.thickness.sm05.mghlogyflag 0usedti 0FSGD rmanova.fsgdlabelmask /Freesurfer//fsaverage/label/lh.cortex.labelmaskinv 0glmdir lh.thickness.rmanovaIllCondOK 0ReScaleX 1DoFFx 0Creating output directory lh.thickness.rmanovaLoading y from /Freesurfer/k/lh.thickness.sm05.mghINFO: gd2mtx_method is dossSaving design matrix to lh.thickness.rmanova/Xg.datNormalized matrix condition is 1.73205Matrix condition is 1.73205Found 149955 points in label.Pruning voxels by thr: 0.000000Found 149955 voxels in maskSaving mask to lh.thickness.rmanova/mask.mghReshaping mriglm->mask...search space = 74612.965197DOF = 0ERROR: DOF = 0
It seems still something wrong was ignored.
Also if I call the preproc or surf2surf commands with the base subject, i.e. k, instead of fsaverage, will it work? or even better?mris_preproc --target k --hemi lh --meas thickness --out lh.thickness.mgh --fsgd rmanova.fsgd
mri_surf2surf --s k --hemi lh --fwhm 5 --sval lh.thickness.mgh --tval lh.thickness.sm05.mgh
Thanks,Kaiming
On Thu, Jan 30, 2014 at 3:18 AM, Martin Reuter <mreuter@nmr.mgh.harvard.edu> wrote:
Your data was probably mapped to fsaverage when you call the preproc or surf2surf commands? If so you should specify fsaverage as the subject nane.Best Martin
Sent via my smartphone, please excuse brevity.
-------- Original message --------
From: Kaiming Yin
Date:01/29/2014 9:07 PM (GMT-05:00)
To: FreeSurfer Mailing List
Subject: [Freesurfer] Longitudinal RepeatedMeasuresAnova mri_glmfit "--surf subject hemi"
Dear FreeSurfer guys,
I am using longitudinal processing stream and with the RepeatedMeasuresAnova.I follow the instruction here from http://surfer.nmr.mgh.harvard.edu/fswiki/RepeatedMeasuresAnova
my fsgd file is:
GroupDescriptorFile 1Title k_longitudinalClass kVariables TP1-vs-TP2 TP1-vs-TP3Input k_1 k 1 1Input k_2 k -1 0Input k_3 k 0 -1
contrasts is:tp1-vs-tp2.mtx tp1-vs-tp3.mtx tp2-vs-tp3.mtx tp-effect.mtx mean.mtx
0 1 0 0 0 1 0 -1 1 0 1 0 1 0 00 0 1
I run the mris_preproc and mri_surf2surf commands as the same for left hemi thickness,
but when I run the mri_glmfit command it reports "ERROR: you must use '--surface subject hemi' with surface data". if i add the flag as:mri_glmfit ...(as the same before)... --surface k lhor --surface k lh thicknessor --surface k_1.long.k lhor --surface k_1.long.k lh thicknessor --surface k_1 lhor --surface k_1 lh thickness
it still goes wrong and even not make an output folder.
Any problems there?
thanks,KaimingThe information in this e-mail is intended only for the person to whom it is
addressed. If you believe this e-mail was sent to you in error and the e-mail
contains patient information, please contact the Partners Compliance HelpLine at
http://www.partners.org/complianceline . If the e-mail was sent to you in error
but does not contain patient information, please contact the sender and properly
dispose of the e-mail.
-- Martin Reuter, Ph.D. Assistant in Neuroscience - Massachusetts General Hospital Instructor in Neurology - Harvard Medical School MGH / HMS / MIT A.A.Martinos Center for Biomedical Imaging 149 Thirteenth Street, Suite 2301 Charlestown, MA 02129 Phone: +1-617-724-5652 Email: mreuter@nmr.mgh.harvard.edu reuter@mit.edu Web : http://reuter.mit.edu