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