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 /Applications/abin_old/matlab (GroupAna.m calls many other functions in the matlab folder) 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