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  

Previous Message Next Message
Message List New Topic Print View
|
pshuk
April 28, 2014 09:04AM
Hi,
I am using 3dMVM because I want to run ANOVA with unequal subjects in groups. As a sanity check, I am comparing my t-test results from 3dMVM and 3dttest++ but they are not exactly same. They look similar but not same. Is it expected (like they will be different? if So, why). Or maybe I am using some wrong options. Can you guys please check and guide. Here are my two scripts:

1. for 3dMVM and 2 for 3dttest++

*****************3dMVM script**********************************************************

#!/bin/tcsh
set grp1 = PD
set grp2 = ET
set data1 = /mnt/Processing1/Priyank/RESULTS_motion1/PD
set data2 = /mnt/Processing1/Priyank/RESULTS_motion1/ET
3dMVM -prefix test_mvm -jobs 6 \
-model Group \
-wsVars 'task' \
-num_glt 5 \
-gltLabel 1 PD_Mean -gltCode 1 'Group : 1*PD' \
-gltLabel 2 ET_Mean -gltCode 2 'Group : 1*ET' \
-gltLabel 3 PD-ET_inter -gltCode 3 'Group : 1*PD -1*ET' \
-gltLabel 4 PD-ET_same_vis -gltCode 4 'task : 1*same_vis Group : 1*PD -1*ET' \
-gltLabel 5 PD-ET_same_nv -gltCode 5 'task : 1*same_nv Group : 1*PD -1*ET' \
-dataTable \
Subj Group task InputFile \
18803 PD same_vis /mnt/Processing1/Priyank/RESULTS_motion1/PD/18803/stats_blur4.18803+tlrc'[same_vis#0_Coef]' \
18803 PD same_nv /mnt/Processing1/Priyank/RESULTS_motion1/PD/18803/stats_blur4.18803+tlrc'[same_nv#0_Coef]' \
20286 PD same_vis /mnt/Processing1/Priyank/RESULTS_motion1/PD/20286/stats_blur4.20286+tlrc'[same_vis#0_Coef]' \
........
.......
.......
ET_022_1 ET same_nv /mnt/Processing1/Priyank/RESULTS_motion1/ET/ET_022_1/stats_blur4.ET_022_1+tlrc'[same_nv#0_Coef]' \
ET_023_1 ET same_vis /mnt/Processing1/Priyank/RESULTS_motion1/ET/ET_023_1/stats_blur4.ET_023_1+tlrc'[same_vis#0_Coef]' \
ET_023_1 ET same_nv /mnt/Processing1/Priyank/RESULTS_motion1/ET/ET_023_1/stats_blur4.ET_023_1+tlrc'[same_nv#0_Coef]' \

**********************************************END**********************************************************************************


2. ***************************** 3dttest++ script **********************************************************
#!/bin/tcsh
3dttest++ \
-prefix /mnt/Processing1/Priyank/ANALYSIS/MVM/WB.t++.PD-ET.same_vis \
-setA PD \
18803 '/mnt/Processing1/Priyank/RESULTS_motion1/PD/18803/stats_blur4.18803+tlrc.HEAD[same_vis#0_Coef]' \
20286 '/mnt/Processing1/Priyank/RESULTS_motion1/PD/20286/stats_blur4.20286+tlrc.HEAD[same_vis#0_Coef]' \
.......
.......
-setB ET \
16119 '/mnt/Processing1/Priyank/RESULTS_motion1/ET/16119/stats_blur4.16119+tlrc.HEAD[same_vis#0_Coef]' \
ET_022_1 '/mnt/Processing1/Priyank/RESULTS_motion1/ET/ET_022_1/stats_blur4.ET_022_1+tlrc.HEAD[same_vis#0_Coef]' \
ET_023_1 '/mnt/Processing1/Priyank/RESULTS_motion1/ET/ET_023_1/stats_blur4.ET_023_1+tlrc.HEAD[same_vis#0_Coef]' \
.......
.......
-unpooled
**********************************************END**********************************************************************************
Reply Quote
RSS
Subject Author Posted

3dMVM and t-test++ results different!!!

pshuk April 28, 2014 09:04AM

Re: 3dMVM and t-test++ results different!!!

gang April 28, 2014 11:20AM