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  

|
March 01, 2016 12:07PM
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
Subject Author Posted

stim time warning

3dFoodie February 29, 2016 04:53PM

Re: stim time warning

rick reynolds February 29, 2016 08:31PM

Re: stim time warning

3dFoodie March 01, 2016 11:31AM

Re: stim time warning

rick reynolds March 01, 2016 11:41AM

Re: stim time warning

3dFoodie March 01, 2016 12:07PM

Re: stim time warning

ace327 April 12, 2016 06:12PM

Re: stim time warning

rick reynolds April 12, 2016 08:53PM

Re: stim time warning

cephalopod April 14, 2016 02:08PM

Re: stim time warning

rick reynolds April 14, 2016 04:03PM

Re: stim time warning

cephalopod April 15, 2016 11:56AM

Re: stim time warning

rick reynolds April 15, 2016 12:01PM

Re: stim time warning

rick reynolds July 19, 2017 10:24AM

Re: stim time warning

EXP.2086 October 22, 2017 09:54PM

Re: stim time warning

rick reynolds October 23, 2017 09:10AM

Re: stim time warning

EXP.2086 October 23, 2017 01:43PM

Re: stim time warning

rick reynolds October 23, 2017 09:35PM

Re: stim time warning

EXP.2086 October 23, 2017 09:54PM

Re: stim time warning

rick reynolds October 23, 2017 09:58PM