Hi Gang-
What I'm trying to do is (I think) super-simple. I am only interested in how a single covariate modifies the effect of Subj.
This is the output of the 3dttest++ script I ran.
---
Dataset File: Pessig_Ttest_PE_New_Patients_CovBPRS+tlrc
Identifier Code: XYZ_LtPNovy7stOZkRiDL-Winw Creation Date: Mon May 8 14:53:56 2017
Template Space: TLRC
Dataset Type: Func-Bucket (-fbuc)
Byte Order: LSB_FIRST [this CPU native = LSB_FIRST]
Storage Mode: BRIK
Storage Space: 21,959,824 (22 million [mega]) bytes
Geometry String: "MATRIX(1.5,0,0,-79.5,0,1.5,0,-79.5,0,0,1.5,-65):107,127,101"
Data Axes Tilt: Plumb
Data Axes Orientation:
first (x) = Right-to-Left
second (y) = Anterior-to-Posterior
third (z) = Inferior-to-Superior [-orient RAI]
R-to-L extent: -79.500 [R] -to- 79.500 [L] -step- 1.500 mm [107 voxels]
A-to-P extent: -79.500 [A] -to- 109.500 [P] -step- 1.500 mm [127 voxels]
I-to-S extent: -65.000
-to- 85.000 [S] -step- 1.500 mm [101 voxels]
Number of values stored at each pixel = 4
-- At sub-brick #0 'Patients_mean' datum type is float: -0.160391 to 0.160074
-- At sub-brick #1 'Patients_Tstat' datum type is float: -5.70935 to 7.06039
statcode = fitt; statpar = 23
-- At sub-brick #2 'Patients_bprs_tot' datum type is float: -0.2191 to 0.257849
-- At sub-brick #3 'Patients_bprs_tot_Tstat' datum type is float: -5.669 to 5.80341
statcode = fitt; statpar = 23
----- HISTORY -----
[jwaltz@admins: Mon May 8 14:53:56 2017] 3dttest++ -prefix Pessig_Ttest_PE_New_Patients_CovBPRS -mask ./mask_group+tlrc -setA Patients 12380 'statsGN25_new_PEonly.12380+tlrc.[All_PE#1_Coef]' 12390 'statsGN25_new_PEonly.12390+tlrc.[All_PE#1_Coef]' 12425 'statsGN25_new_PEonly.12425+tlrc.[All_PE#1_Coef]' 12463 'statsGN25_new_PEonly.12463+tlrc.[All_PE#1_Coef]' 12468 'statsGN25_new_PEonly.12468+tlrc.[All_PE#1_Coef]' 12481 'statsGN25_new_PEonly.12481+tlrc.[All_PE#1_Coef]' 12505 'statsGN25_new_PEonly.12505+tlrc.[All_PE#1_Coef]' 12609 'statsGN25_new_PEonly.12609+tlrc.[All_PE#1_Coef]' 12614 'statsGN25_new_PEonly.12614+tlrc.[All_PE#1_Coef]' 12616 'statsGN25_new_PEonly.12616+tlrc.[All_PE#1_Coef]' 12645 'statsGN25_new_PEonly.12645+tlrc.[All_PE#1_Coef]' 12679 'statsGN25_new_PEonly.12679+tlrc.[All_PE#1_Coef]' 12787 'statsGN25_new_PEonly.12787+tlrc.[All_PE#1_Coef]' 12790 'statsGN25_new_PEonly.12790+tlrc.[All_PE#1_Coef]' 12833 'statsGN25_new_PEonly.12833+tlrc.[All_PE#1_Coef]' 12834 'statsGN25_new_PEonly.12834+tlrc.[All_PE#1_Coef]' 12863 'statsGN25_new_PEonly.12863+tlrc.[All_PE#1_Coef]' 12864 'statsGN25_new_PEonly.12864+tlrc.[All_PE#1_Coef]' 12874 'statsGN25_new_PEonly.12874+tlrc.[All_PE#1_Coef]' 12882 'statsGN25_new_PEonly.12882+tlrc.[All_PE#1_Coef]' 12906 'statsGN25_new_PEonly.12906+tlrc.[All_PE#1_Coef]' 12953 'statsGN25_new_PEonly.12953+tlrc.[All_PE#1_Coef]' 12962 'statsGN25_new_PEonly.12962+tlrc.[All_PE#1_Coef]' 13063 'statsGN25_new_PEonly.13063+tlrc.[All_PE#1_Coef]' 13125 'statsGN25_new_PEonly.13125+tlrc.[All_PE#1_Coef]' -covariates /data/brutus_data12/jwaltz/pessig/scripts/group/BPRS_Totals.txt
---
I tried to do the same thing with 3dMVM, but I keep getting the following error:
---
Error in seq.default(wd + ii, len, wd) : wrong sign in 'by' argument
Calls: process.MVM.opts -> data.frame -> cbind -> seq -> seq.default
Execution halted
Test_Pessig_3dMVM_AllPE_BPRS.sh: line 13: 12380: command not found
---
This is that script:
---
#! bin/sh
cd /Users/jameswaltz/Analyses/AFNIana/PessigMRI/group.results.parametric
name='[All_PE#1_Coef]'
3dMVM -prefix Pessig_AllPE_BPRS_Param.3dMVM -mask ./mask_group+tlrc. \
-bsVars "bprstot" \
-qVars "bprstot" \
-jobs 8 \
-dataTable \
Subj bprstot InputFile \
12380 150 ./statsGN25_new_PEonly.12380+tlrc.${name} \
12390 325 ./statsGN25_new_PEonly.12390+tlrc.${name} \
12425 200 ./statsGN25_new_PEonly.12425+tlrc.${name} \
12463 100 ./statsGN25_new_PEonly.12463+tlrc.${name} \
12468 125 ./statsGN25_new_PEonly.12468+tlrc.${name} \
12481 100 ./statsGN25_new_PEonly.12481+tlrc.${name} \
12505 100 ./statsGN25_new_PEonly.12505+tlrc.${name} \
12609 150 ./statsGN25_new_PEonly.12609+tlrc.${name} \
12614 100 ./statsGN25_new_PEonly.12614+tlrc.${name} \
12616 400 ./statsGN25_new_PEonly.12616+tlrc.${name} \
12645 100 ./statsGN25_new_PEonly.12645+tlrc.${name} \
12679 150 ./statsGN25_new_PEonly.12679+tlrc.${name} \
12787 100 ./statsGN25_new_PEonly.12787+tlrc.${name} \
12790 100 ./statsGN25_new_PEonly.12790+tlrc.${name} \
12833 175 ./statsGN25_new_PEonly.12833+tlrc.${name} \
12834 150 ./statsGN25_new_PEonly.12834+tlrc.${name} \
12863 125 ./statsGN25_new_PEonly.12863+tlrc.${name} \
12864 150 ./statsGN25_new_PEonly.12864+tlrc.${name} \
12874 225 ./statsGN25_new_PEonly.12874+tlrc.${name} \
12882 125 ./statsGN25_new_PEonly.12882+tlrc.${name} \
12906 275 ./statsGN25_new_PEonly.12906+tlrc.${name} \
12953 150 ./statsGN25_new_PEonly.12953+tlrc.${name} \
12962 300 ./statsGN25_new_PEonly.12962+tlrc.${name} \
13063 100 ./statsGN25_new_PEonly.13063+tlrc.${name} \
13125 275 ./statsGN25_new_PEonly.13125+tlrc.${name}
