Aaron,

Hi, can you tell me the version of make_average_surface you are running? Type:

make_average_surface --version

Just in case yours in out-of-date, I've attached the latest make_average_surface script, which should not have the SUBJECTS bug that you describe.

The file is gzipped, so do the following:

- copy the file to $FREESURFER_HOME/bin - change directory to $FREESURFER_HOME/bin - type: gunzip make_average_surface.gz

Then you should be able rerun the make_average_surface utility.

Nick

On Tue, 2005-07-26 at 15:11 -0400, Goldman, Aaron (NIH/NIMH) wrote:

Hi, I’m trying to do a thickness analysis with mris_glm, and I had a couple of questions:

1. I made an average surface using mris_make_average_surface, but I

won’t load it with tksurfer because it “could not find anatomical header information.” Do you know how I can fix this? I also tried “make_average_surface -subjects subject1… subject111,” but it said SUBJECTS wasn’t defined, and wouldn’t run.

2. I’ve been able to successfully run a test mris_glm command using a

design file like:

subject1 0 1

subject2 1 0

etc

and the following commands:

foreach h ( lh rh )

mris_glm    --design $SUBJECTS_DIR/results/test_design.txt \

                    --surfmeas thickness \

                    --hemi $h \

                    --trgsubj Colin \

                    --nsmooth 300 \

                    --beta ${h}_test_beta bfloat \

                    --var ${h}_test_var bfloat

mris_glm    --design $SUBJECTS_DIR/results/test_design.txt \

                    --hemi $h \

                    --trgsubj Colin \

                    --beta_in ${h}_test_beta bfloat \

                    --var_in ${h}_test_var bfloat \

                    --gcv 1 -1 \

                    --sigt $SUBJECTS_DIR/Colin/surf/${h}_test_sigt

paint \

                    --t $SUBJECTS_DIR/$trgsubj/surf/${h}_test_t

paint

end

where a 1 in column 2 indicates one allele of the gene we’re studying, whereas a 1 in column 3 indicates the other. However, I’d like to also covary for age and gender, and I was wondering how I would need to edit my design file and mris_glm commands to do that.

Thanks,

-Aaron-

Aaron L. Goldman

Post-Bacc IRTA

Genes, Cognition and Psychosis Program, NIMH

Building 10, Room 3C108

Phone: (301) 435-0944

Email: aarong@mail.nih.gov

---

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