AFNI Message Board

Dear AFNI users-

We are very pleased to announce that the new AFNI Message Board framework is up! Please join us at:

https://discuss.afni.nimh.nih.gov

Existing user accounts have been migrated, so returning users can login by requesting a password reset. New users can create accounts, as well, through a standard account creation process. Please note that these setup emails might initially go to spam folders (esp. for NIH users!), so please check those locations in the beginning.

The current Message Board discussion threads have been migrated to the new framework. The current Message Board will remain visible, but read-only, for a little while.

Sincerely, AFNI HQ

History of AFNI updates  

|
April 08, 2022 05:04PM
Hello,

Thank you for replying. I am working on a group level t-test.
The script is below. I have only included the first 5 subjects in each group due to space constraints.


cd /nfs/mint/results/gng/SubjectsFolders/Test_All_Subj

set output_dir = /nfs/mint/results/gng/SubjectsFolders/GroupAnalysis/Take2

#GO_ANGRY

set output = go_angry_gng

3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[1]' \
2 AP1129.results/stats.AP1129+tlrc.'[1]' \
3 AR1342.results/stats.AR1342+tlrc.'[1]' \
4 BK1041.results/stats.BK1041+tlrc.'[1]' \
5 BK1268.results/stats.BK1268+tlrc.'[1]' \
\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[1]' \
50 AD1173.results/stats.AD1173+tlrc.'[1]' \
51 AM1304.results/stats.AM1304+tlrc.'[1]' \
52 AM1319.results/stats.AM1319+tlrc.'[1]' \
53 AP1162.results/stats.AP1162+tlrc.'[1]' \
54 AS1001.results/stats.AS1001+tlrc.'[1]' \


3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#GO_ANGRY_OM


set output = go_angry_om_ttest_gng

3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[4]' \
2 AP1129.results/stats.AP1129+tlrc.'[4]' \
3 AR1342.results/stats.AR1342+tlrc.'[4]' \
4 BK1041.results/stats.BK1041+tlrc.'[4]' \
5 BK1268.results/stats.BK1268+tlrc.'[4]' \
\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[4]' \
50 AD1173.results/stats.AD1173+tlrc.'[4]' \
51 AM1304.results/stats.AM1304+tlrc.'[4]' \
52 AM1319.results/stats.AM1319+tlrc.'[4]' \
53 AP1162.results/stats.AP1162+tlrc.'[4]' \
54 AS1001.results/stats.AS1001+tlrc.'[4]' \


3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#No Go Angry

set output = nogo_angry_ttest_gng


3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[7]' \
2 AP1129.results/stats.AP1129+tlrc.'[7]' \
3 AR1342.results/stats.AR1342+tlrc.'[7]' \
4 BK1041.results/stats.BK1041+tlrc.'[7]' \
5 BK1268.results/stats.BK1268+tlrc.'[7]' \

\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[7]' \
50 AD1173.results/stats.AD1173+tlrc.'[7]' \
51 AM1304.results/stats.AM1304+tlrc.'[7]' \
52 AM1319.results/stats.AM1319+tlrc.'[7]' \
53 AP1162.results/stats.AP1162+tlrc.'[7]' \
54 AS1001.results/stats.AS1001+tlrc.'[7]' \
55 AS1017.results/stats.AS1017+tlrc.'[7]' \


3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#No Go Angry Comm

set output = nogo_angry_comm_ttest_gng

3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[10]' \
2 AP1129.results/stats.AP1129+tlrc.'[10]' \
3 AR1342.results/stats.AR1342+tlrc.'[10]' \
4 BK1041.results/stats.BK1041+tlrc.'[10]' \
5 BK1268.results/stats.BK1268+tlrc.'[10]' \

\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[10]' \
50 AD1173.results/stats.AD1173+tlrc.'[10]' \
51 AM1304.results/stats.AM1304+tlrc.'[10]' \
52 AM1319.results/stats.AM1319+tlrc.'[10]' \
53 AP1162.results/stats.AP1162+tlrc.'[10]' \
54 AS1001.results/stats.AS1001+tlrc.'[10]' \
55 AS1017.results/stats.AS1017+tlrc.'[10]' \


3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#Go Fear

set output = go_fear_ttest_gng


3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 BK1268.results/stats.BK1268+tlrc.'[13]' \
1 AH1346.results/stats.AH1346+tlrc.'[13]'
2 AP1129.results/stats.AP1129+tlrc.'[13]' \
3 AR1342.results/stats.AR1342+tlrc.'[13]' \
4 BK1041.results/stats.BK1041+tlrc.'[13]' \
5 BK1268.results/stats.BK1268+tlrc.'[13]' \
\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[13]' \
50 AD1173.results/stats.AD1173+tlrc.'[13]'
51 AM1304.results/stats.AM1304+tlrc.'[13]' \
52 AM1319.results/stats.AM1319+tlrc.'[13]' \
53 AP1162.results/stats.AP1162+tlrc.'[13]' \
54 AS1001.results/stats.AS1001+tlrc.'[13]' \
55 AS1017.results/stats.AS1017+tlrc.'[13]' \

3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#Go Fear Omission

set output = go_fear_om_ttest_gng


3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[16]' \
2 AP1129.results/stats.AP1129+tlrc.'[16]' \
3 AR1342.results/stats.AR1342+tlrc.'[16]' \
4 BK1041.results/stats.BK1041+tlrc.'[16]' \
5 BK1268.results/stats.BK1268+tlrc.'[16]' \
\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[16]' \
50 AD1173.results/stats.AD1173+tlrc.'[16]' \
51 AM1304.results/stats.AM1304+tlrc.'[16]' \
52 AM1319.results/stats.AM1319+tlrc.'[16]' \
53 AP1162.results/stats.AP1162+tlrc.'[16]' \
54 AS1001.results/stats.AS1001+tlrc.'[16]' \
55 AS1017.results/stats.AS1017+tlrc.'[16]' \


3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#No Go Fear

set output = nogo_fear_ttest_gng


3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[19]' \
2 AP1129.results/stats.AP1129+tlrc.'[19]' \
3 AR1342.results/stats.AR1342+tlrc.'[19]' \
4 BK1041.results/stats.BK1041+tlrc.'[19]' \
5 BK1268.results/stats.BK1268+tlrc.'[19]' \

\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[19]' \
50 AD1173.results/stats.AD1173+tlrc.'[19]' \
51 AM1304.results/stats.AM1304+tlrc.'[19]' \
52 AM1319.results/stats.AM1319+tlrc.'[19]' \
53 AP1162.results/stats.AP1162+tlrc.'[19]' \
54 AS1001.results/stats.AS1001+tlrc.'[19]' \
55 AS1017.results/stats.AS1017+tlrc.'[19]' \

3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#No Go Fear Comm

set output = nogo_fear_comm_ttest_gng


3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[22]' \
2 AP1129.results/stats.AP1129+tlrc.'[22]' \
3 AR1342.results/stats.AR1342+tlrc.'[22]' \
4 BK1041.results/stats.BK1041+tlrc.'[22]' \
5 BK1268.results/stats.BK1268+tlrc.'[22]' \

\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[22]' \
50 AD1173.results/stats.AD1173+tlrc.'[22]' \
51 AM1304.results/stats.AM1304+tlrc.'[22]' \
52 AM1319.results/stats.AM1319+tlrc.'[22]' \
53 AP1162.results/stats.AP1162+tlrc.'[22]' \
54 AS1001.results/stats.AS1001+tlrc.'[22]' \
55 AS1017.results/stats.AS1017+tlrc.'[22]' \



3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#Go Happy

set output = go_happy_ttest_gng


3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[25]' \
2 AP1129.results/stats.AP1129+tlrc.'[25]' \
3 AR1342.results/stats.AR1342+tlrc.'[25]' \
4 BK1041.results/stats.BK1041+tlrc.'[25]' \
5 BK1268.results/stats.BK1268+tlrc.'[25]' \

\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[25]' \
50 AD1173.results/stats.AD1173+tlrc.'[25]' \
51 AM1304.results/stats.AM1304+tlrc.'[25]' \
52 AM1319.results/stats.AM1319+tlrc.'[25]' \
53 AP1162.results/stats.AP1162+tlrc.'[25]' \
54 AS1001.results/stats.AS1001+tlrc.'[25]' \
55 AS1017.results/stats.AS1017+tlrc.'[25]' \

3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#Go Happy Omission

set output = go_happy_om_ttest_gng


3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 BK1268.results/stats.BK1268+tlrc.'[28]' \
1 AH1346.results/stats.AH1346+tlrc.'[28]' \
2 AP1129.results/stats.AP1129+tlrc.'[28]' \
3 AR1342.results/stats.AR1342+tlrc.'[28]' \
4 BK1041.results/stats.BK1041+tlrc.'[28]' \
5 BK1268.results/stats.BK1268+tlrc.'[28]' \

\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[28]' \
50 AD1173.results/stats.AD1173+tlrc.'[28]' \
51 AM1304.results/stats.AM1304+tlrc.'[28]' \
52 AM1319.results/stats.AM1319+tlrc.'[28]' \
53 AP1162.results/stats.AP1162+tlrc.'[28]' \
54 AS1001.results/stats.AS1001+tlrc.'[28]' \
55 AS1017.results/stats.AS1017+tlrc.'[28]' \



3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#No Go Happy

set output = nogo_happy_ttest_gng


3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[31]' \
2 AP1129.results/stats.AP1129+tlrc.'[31]' \
3 AR1342.results/stats.AR1342+tlrc.'[31]' \
4 BK1041.results/stats.BK1041+tlrc.'[31]' \
5 BK1268.results/stats.BK1268+tlrc.'[31]' \
\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[31]' \
50 AD1173.results/stats.AD1173+tlrc.'[31]' \
51 AM1304.results/stats.AM1304+tlrc.'[31]' \
52 AM1319.results/stats.AM1319+tlrc.'[31]' \
53 AP1162.results/stats.AP1162+tlrc.'[31]' \
54 AS1001.results/stats.AS1001+tlrc.'[31]' \
55 AS1017.results/stats.AS1017+tlrc.'[31]' \

3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#No Go Happy Comm

set output = nogo_happy_comm_ttest_gng


3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[34]' \
2 AP1129.results/stats.AP1129+tlrc.'[34]' \
3 AR1342.results/stats.AR1342+tlrc.'[34]' \
4 BK1041.results/stats.BK1041+tlrc.'[34]' \
5 BK1268.results/stats.BK1268+tlrc.'[34]' \
\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[34]' \
50 AD1173.results/stats.AD1173+tlrc.'[34]' \
51 AM1304.results/stats.AM1304+tlrc.'[34]' \
52 AM1319.results/stats.AM1319+tlrc.'[34]' \
53 AP1162.results/stats.AP1162+tlrc.'[34]' \
54 AS1001.results/stats.AS1001+tlrc.'[34]' \
55 AS1017.results/stats.AS1017+tlrc.'[34]' \

3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#Go Control

set output = go_control_ttest_gng

3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[37]' \
2 AP1129.results/stats.AP1129+tlrc.'[37]' \
3 AR1342.results/stats.AR1342+tlrc.'[37]' \
4 BK1041.results/stats.BK1041+tlrc.'[37]' \
5 BK1268.results/stats.BK1268+tlrc.'[37]' \
\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[37]' \
50 AD1173.results/stats.AD1173+tlrc.'[37]' \
51 AM1304.results/stats.AM1304+tlrc.'[37]' \
52 AM1319.results/stats.AM1319+tlrc.'[37]' \
53 AP1162.results/stats.AP1162+tlrc.'[37]' \
54 AS1001.results/stats.AS1001+tlrc.'[37]' \
55 AS1017.results/stats.AS1017+tlrc.'[37]' \

3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#Go Control Omission

set output = go_control_om_ttest_gng

3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[40]' \
2 AP1129.results/stats.AP1129+tlrc.'[40]' \
3 AR1342.results/stats.AR1342+tlrc.'[40]' \
4 BK1041.results/stats.BK1041+tlrc.'[40]' \
5 BK1268.results/stats.BK1268+tlrc.'[40]' \
\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[40]' \
50 AD1173.results/stats.AD1173+tlrc.'[40]' \
51 AM1304.results/stats.AM1304+tlrc.'[40]' \
52 AM1319.results/stats.AM1319+tlrc.'[40]'
53 AP1162.results/stats.AP1162+tlrc.'[40]' \
54 AS1001.results/stats.AS1001+tlrc.'[40]' \
55 AS1017.results/stats.AS1017+tlrc.'[40]' \

3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#No Go Control

set output = nogo_control_ttest_gng

3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[43]' \
2 AP1129.results/stats.AP1129+tlrc.'[43]' \
3 AR1342.results/stats.AR1342+tlrc.'[43]' \
4 BK1041.results/stats.BK1041+tlrc.'[43]' \
5 BK1268.results/stats.BK1268+tlrc.'[43]' \
\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[43]' \
50 AD1173.results/stats.AD1173+tlrc.'[43]' \
51 AM1304.results/stats.AM1304+tlrc.'[43]' \
52 AM1319.results/stats.AM1319+tlrc.'[43]' \
53 AP1162.results/stats.AP1162+tlrc.'[43]' \
54 AS1001.results/stats.AS1001+tlrc.'[43]' \
55 AS1017.results/stats.AS1017+tlrc.'[43]' \

3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

#No Go Control Comm

set output = nogo_control_comm_ttest_gng

3dttest++ -prefix ${output_dir}/${output} \
\
-setA TD \
\
1 AH1346.results/stats.AH1346+tlrc.'[46]' \
2 AP1129.results/stats.AP1129+tlrc.'[46]' \
3 AR1342.results/stats.AR1342+tlrc.'[46]' \
4 BK1041.results/stats.BK1041+tlrc.'[46]' \
5 BK1268.results/stats.BK1268+tlrc.'[46]' \
\
-setB ADHD \
\
49 AA1299.results/stats.AA1299+tlrc.'[46]' \
50 AD1173.results/stats.AD1173+tlrc.'[46]' \
51 AM1304.results/stats.AM1304+tlrc.'[46]' \
52 AM1319.results/stats.AM1319+tlrc.'[46]' \
53 AP1162.results/stats.AP1162+tlrc.'[46]' \
54 AS1001.results/stats.AS1001+tlrc.'[46]' \
55 AS1017.results/stats.AS1017+tlrc.'[46]' \

3dcalc -a ${output_dir}/rsttest+tlrc. -b ${output_dir}/${output}+tlrc. -expr 'ispositive(a)*b' -prefix ${output_dir}/${output}_masked

echo "end"
Subject Author Posted

issue with running t-test in afni

ananru April 05, 2022 04:32PM

Re: issue with running t-test in afni

gang April 05, 2022 10:39PM

Re: issue with running t-test in afni

ananru April 08, 2022 05:04PM

Re: issue with running t-test in afni

gang April 09, 2022 12:10PM