|
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: [EMAIL PROTECTED]
|
_______________________________________________ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
