Beauchamp:GroupAna

Using GroupAna.m to do ANOVA's with different group sizes http://afni.nimh.nih.gov/sscc/gangc/Interface.html/

1. Make sure you have the most updated version (~2005) http://afni.nimh.nih.gov/afni/community/board/read.php?f=1&i=37940&t=37940#reply_37940

2. Need individual subBRIKs in individual files (you can't just reference a single subBRIK in a bigger dataset)

foreach ec  (CAM CAS CAV CBA CBB CBD CBE CBI CBJ CBK CBL CBM CBN CBO CBP CBQ CBU) 3dbucket -prefix {$ec}_McG {$ec}bl6+tlrc'[2]' 3dbucket -prefix {$ec}_InC {$ec}bl6+tlrc'[5]' 3dbucket -prefix {$ec}_Cong {$ec}bl6+tlrc'[8]' end

3. Run from whatever folder has all of the Matlab files in it (/Applications/abin_old/matlab), since GroupAna.m calls many of the other Matlab functions.

Reference datafiles from the server ( for example, /Volumes/data1/UT/ChildMcGurk/CAM_McG+tlrc.BRIK)

Example of commands

%-- 4/22/11 9:59 AM --%

GroupAna 0 3 3 0 1

PercGroup 2 NonPerc Perc

StimulusType 3 McG Inc Cong

Subject 7 CAM CBM CBP CAV CBK CBB CBQ 10 CBD CAS CBE CBO CBU CBI CBA CBJ CBL CBN 1 1

/Volumes/data1/UT/ChildMcGurk/CAM_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBM_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBP_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CAV_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBK_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBB_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBQ_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CAM_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBM_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBP_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CAV_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBK_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBB_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBQ_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CAM_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBM_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBP_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CAV_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBK_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBB_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBQ_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBD_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CAS_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBE_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBO_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBU_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBI_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBA_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBJ_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBL_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBN_McG+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBD_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CAS_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBE_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBO_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBU_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBI_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBA_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBJ_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBL_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBN_InC+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBD_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CAS_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBE_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBO_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBU_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBI_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBA_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBJ_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBL_Cong+tlrc.BRIK /Volumes/data1/UT/ChildMcGurk/CBN_Cong+tlrc.BRIK GroupAna_Child_6mm 1 3

Perc-NonPerc 2 100 -1 200 1

McG-InC 2 010 1 020 -1

McG-Cong 2 010 1 030 -1 1

PercMcG-NonPercMcG 2 110 -1 210 1 0

4. Add Full F stat (make sure datum type is float!!) 3dbucket -glueto GroupAna_Child_6mm+tlrc -fbuc NewFullF_17_float+tlrc'[0]' 3drefit -sublabel 18 NewFullF GroupAna_Child_6mm+tlrc