Great.
Sorry I thought it was attached, but I guess txt files don't attach here. I've pasted the beginning and end without all the calculations output below. If you'd like to see everything please let me know and I will send.
So any idea what those warnings are? OK to ignore?
Thanks a lot
Running GLM
*+ WARNING: '-stim_times 2' didn't read any good times from file 'BN100_HappyGo.1D'
*+ WARNING: '-stim_times 3' didn't read any good times from file 'BN100_FearGo.1D'
*+ WARNING: '-stim_times 4' didn't read any good times from file 'BN100_CalmNoGo.1D'
*+ WARNING: '-stim_times 5' didn't read any good times from file 'BN100_HappyNoGo.1D'
++ '-stim_times_AM1 8 BN100_Negative.1D' has 1 auxiliary values per time point
++ '-stim_times_AM1 8': basis function model 'dmUBLOCK' uses 1 parameters,
out of the 1 found in timing file 'BN100_Negative.1D'
++ '-stim_times_AM1 9 BN100_Neutral.1D' has 1 auxiliary values per time point
++ '-stim_times_AM1 9': basis function model 'dmUBLOCK' uses 1 parameters,
out of the 1 found in timing file 'BN100_Neutral.1D'
++ '-stim_times_AM1 10 BN100_Positive.1D' has 1 auxiliary values per time point
++ '-stim_times_AM1 10': basis function model 'dmUBLOCK' uses 1 parameters,
out of the 1 found in timing file 'BN100_Positive.1D'
++ 3dDeconvolve: AFNI version=AFNI_16.0.00 (Jan 1 2016) [64-bit]
++ Authored by: B. Douglas Ward, et al.
++ current memory malloc-ated = 345,015 bytes (about 345 thousand [kilo])
++ loading dataset pb04.BN100.r1.scale+tlrc pb04.BN100.r2.scale+tlrc pb04.BN100.r3.scale+tlrc pb04.BN100.r4.scale+tlrc pb04.BN100.r5.scale+tlrc pb04.BN100.r6.scale+tlrc
++ current memory malloc-ated = 1,239,367,748 bytes (about 1.2 billion [giga])
++ Auto-catenated input datasets treated as multiple imaging runs
++ Auto-catenated datasets start at: 0 190 380 570 760 950
++ STAT automask has 231931 voxels (out of 271633 = 85.4%)
++ Skipping check for initial transients
++ Imaging duration=475.0 s; Automatic polort=4
++ -stim_times using TR=2.5 s for stimulus timing conversion
++ -stim_times using TR=2.5 s for any -iresp output datasets
++ [you can alter the -iresp TR via the -TR_times option]
++ -stim_times 1 using GLOBAL times
++ -stim_times 2 using GLOBAL times
++ -stim_times 3 using GLOBAL times
++ -stim_times 4 using GLOBAL times
++ -stim_times 5 using GLOBAL times
++ -stim_times 6 using GLOBAL times
++ -stim_times 7 using GLOBAL times
++ -stim_times_AM1 8 using GLOBAL times
++ -stim_times_AM1 9 using GLOBAL times
++ -stim_times_AM1 10 using GLOBAL times
++ Number of time points: 1140 (before censor) ; 1087 (after)
+ Number of parameters: 60 [36 baseline ; 24 signal]
++ total shared memory needed = 1,392,934,024 bytes (about 1.4 billion [giga])
++ current memory malloc-ated = 1,239,792,580 bytes (about 1.2 billion [giga])
++ mmap() memory allocated: 1,392,934,024 bytes (about 1.4 billion [giga])
++ Memory required for output bricks = 1,392,934,024 bytes (about 1.4 billion [giga])
++ Wrote matrix image to file BN_CtxEGNG_EmoXRespXCtx_BN100.jpg
++ Wrote matrix values to file CtxEGNG_EmoXResp_Ctx_SPMG3_BN100_bucket.xmat.1D
++ ========= Things you can do with the matrix file =========
++ (a) Linear regression with ARMA(1,1) modeling of serial correlation:
3dREMLfit -matrix CtxEGNG_EmoXResp_Ctx_SPMG3_BN100_bucket.xmat.1D \
-input "pb04.BN100.r1.scale+tlrc pb04.BN100.r2.scale+tlrc pb04.BN100.r3.scale+tlrc pb04.BN100.r4.scale+tlrc pb04.BN100.r5.scale+tlrc pb04.BN100.r6.scale+tlrc" \
-Rbeta CtxEGNG_EmoXResp_Ctx_SPMG3_BN100_cbucket_REML \
-fout -tout -rout -Rbuck CtxEGNG_EmoXResp_Ctx_SPMG3_BN100_bucket_REML -Rvar CtxEGNG_EmoXResp_Ctx_SPMG3_BN100_bucket_REMLvar \
-Rerrts CtxEGNG_EmoXResp_Ctx_SPMG3_BN100_errts_REML -verb
++ N.B.: 3dREMLfit command above written to file CtxEGNG_EmoXResp_Ctx_SPMG3_BN100_bucket.REML_cmd
++ (b) Visualization/analysis of the matrix via ExamineXmat.R
++ (c) Synthesis of sub-model datasets using 3dSynthesize
++ ==========================================================
++ ----- Signal+Baseline matrix condition [X] (1087x60): 7.31609 ++ VERY GOOD ++
++ ----- Signal-only matrix condition [X] (1087x24): 3.341 ++ VERY GOOD ++
++ ----- Baseline-only matrix condition [X] (1087x36): 7.08657 ++ VERY GOOD ++
++ ----- stim_base-only matrix condition [X] (1087x6): 3.74586 ++ VERY GOOD ++
++ ----- polort-only matrix condition [X] (1087x30): 1.11433 ++ VERY GOOD ++
++ +++++ Matrix inverse average error = 2.57612e-15 ++ VERY GOOD ++
++ Matrix setup time = 1.35 s
++ Voxels in dataset: 271633
++ Voxels per job: 135816
++ Job #1: processing voxels 135816 to 271632; elapsed time=11.440
++ Job #0: processing voxels 0 to 135815; elapsed time=11.450
####################
## A LOT CUT OUT###
####################
++ Job #0 waiting for children to finish; elapsed time=473.696
++ Job #0 now finishing up; elapsed time=473.696
++ Wrote cbucket to ./CtxEGNG_EmoXResp_Ctx_SPMG3_BN100_cbucket+tlrc.BRIK
*+ WARNING: Smallest FDR q [27 CalmNoGo#0_Tstat] = 0.207485 ==> very few (if any) true detections
*+ WARNING: Smallest FDR q [29 CalmNoGo#1_Tstat] = 0.84067 ==> very few (if any) true detections
*+ WARNING: Smallest FDR q [31 CalmNoGo#2_Tstat] = 0.50051 ==> very few (if any) true detections
*+ WARNING: Smallest FDR q [37 HappyNoGo#1_Tstat] = 0.336163 ==> very few (if any) true detections
*+ WARNING: Smallest FDR q [39 HappyNoGo#2_Tstat] = 0.817056 ==> very few (if any) true detections
*+ WARNING: Smallest FDR q [47 FearNoGo#2_Tstat] = 0.999814 ==> very few (if any) true detections
++ Wrote bucket dataset into ./CtxEGNG_EmoXResp_Ctx_SPMG3_BN100_bucket+tlrc.BRIK
+ created 46 FDR curves in bucket header
++ Wrote 3D+time dataset into ./CtxEGNG_EmoXResp_Ctx_SPMG3_BN100_errts+tlrc.BRIK
++ Program finished; elapsed time=667.526