Yes, that looks correct. If you are interested in the offset, you should probably demean the thickness. You can do this by running

mri_glmfit --y rh.thickness.10B.mgh --osgm --surface fsaverage rh --o glm.thickness --save-eres

This will produce a file called eres.mgh in the output. This is the demeaned thicknesses. You can then pass that to your command below instead of the raw thicknesses

On 4/16/18 2:08 PM, Piero Chiacchiaretta wrote:

Hi there,

I wonder if the fsgd file and the command lines listed below are corrected.

My sample included 4 groups with different diagnosis (g1, g2, g3, g4).

I'm interested to map the vertices that show a relationship between the FC (obtained by FS-FAST) and cortical thickeness, regressing out the diagnosis.

To this aim I used the --pvr option:

mri_glmfit --y ces.nii.gz --wls cesvar.nii.gz --surface fsaverage rh --glmdir

Thickness.wls --nii.gz --fsgd g4v1_thickness.fsgd --C intercept.mtx --C

slope.mtx --pvr $SUBJECTS_DIR/rh.thickness.10B.mgh

The fsgd file lists:

---

GroupDescriptorFile 1

Title Relationship FC-thick reg out dignosis

Class g1

Class g2

Class g3

Class g4

Input XX1 g1

Input XX2 g1

...

Input YY1 g2

Input XX2 g2

...

Input XY1 g3

Input XY2 g3

...

Input YX1 g4

Input YX2 g4

...

---

intercept.mtx 0.25 0.25 0.25 0.25 0

slope.mtx 0 0 0 0 1

thanks piero

Piero Chiacchiaretta, PhD Institute for Advanced Biomedical Technologies - ITAB - University “G. d’Annunzio” of Chieti - Pescara Department of Neuroscience, Imaging and Clinical Sciences Via dei Vestini, 31 66013 Chieti, Italy

tel: 39-0871-3556919 tel://tel:%2039-0871-3556919 fax: 39-0871-3556930 tel://fax:%2039-0871-3556930 e-mail: p.chiacchiaretta@unich.it mailto:p.chiacchiaretta@unich.it

---

Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer