3Dttest++ -ETAC Fatal Signal 11 ERROR: nifti_image_write_hdr_img: NBL does not match nim

Trying to implement the new ETAC correction for analysis, got this error at the end of the 3dXClustSim step.

Input:
#!/bin/csh
3dttest++ -prefix ALLFACEvshape_3dttestETAC
-mask group_mask_AUDanxiety+tlrc
-ETAC
-setA ‘/raid10/studies/CII/results/FaceMatching_glt/NV*/stats.NV._REML+tlrc.BRIK[Faces-Base#0_Coef]’
-setB '/raid10/studies/CII/results/FaceMatching_glt/alc
/stats.alc.*_REML+tlrc.BRIK[Faces-Base#0_Coef]’
-labelA HC
-labelB AUD

++ 3dttest++: AFNI version=AFNI_17.1.12 (Jun 23 2017) [64-bit]
++ Authored by: Zhark++
++ 34990 voxels in -mask dataset
++ Number of 3dXClustSim threads set to 60

  • Default -ETAC prefix set to ‘TT.ewJRndftwVF’
    ++ option -setA :: processing as SHORT form (all values are datasets)
    ++ option -setB :: processing as SHORT form (all values are datasets)
    ++ 2-sample test: ‘-AminusB’ option is assumed
    ++ 2-sample test: results will be HC - AUD
    ++ loading -setA datasets
    ++ loading -setB datasets
    ++ Memory usage now = 64,872,878 (65 million [mega])
    ++ t-testing:0123456789.0123456789.0123456789.0123456789.0123456789.!
    ++ ---------- End of analyses – freeing workspaces ----------
    ++ Memory usage now = 88,271,122 (88 million [mega])
    ++ Creating FDR curves in output dataset
    ++ Smallest FDR q [1 HC-AUD_Zscr] = 0.0358173
    ++ Smallest FDR q [3 HC_Zscr] = 4.21587e-14
    ++ Smallest FDR q [5 AUD_Zscr] = 8.9581e-13
  • Added 3 FDR curves to dataset
    [7m** ERROR:[0m output dataset name ‘ALLFACEvshape_3dttestETAC’ conflicts with existing file
    [7m** ERROR:[0m dataset NOT written to disk!
    ++ Output dataset ./TT.ewJRndftwVF.resid.nii
  • 2-sample test: results are HC - AUD
    ++ ================ Starting -ETAC calculations ================
  • === temporary files will have prefix TT.ewJRndftwVF ===
  • === running 60 -randomsign jobs (667 iterations per job) ===
  • === creating 2,799,200,000 (2.8 billion [giga]) bytes of pseudo-data in .sdat files ===
  • — 3dXClustSim reads .sdat files to compute cluster-threshold statistics —
    ++ 3dttest++: AFNI version=AFNI_17.1.12 (Jun 23 2017) [64-bit]
    ++ Authored by: Zhark++
    ++ 34990 voxels in -mask dataset
    ++ option -setA :: processing as SHORT form (all values are datasets)
    ++ option -setB :: processing as SHORT form (all values are datasets)
    ++ random seeds are 1275651492 1589810757
    ++ 2-sample test: ‘-AminusB’ option is assumed
    ++ 2-sample test: results will be SetA - SetB
    ++ opened file ./TT.ewJRndftwVF.0000.sdat for output
    ++ loading -setA datasets
    ++ loading -setB datasets
    ++ Memory usage now = 21,356,811 (21 million [mega])
    ++ t-test randomsign:0123456789.0123456789.0123456789.0123456789.0123456789.0!
    ++ saving main effect t-stat MIN/MAX values in ./TT.ewJRndftwVF.0000.minmax.1D
    ++ output short-ized file ./TT.ewJRndftwVF.0000.sdat
  • ===== simulation jobs have finished (438.7 s elapsed) =====
    ++ Global 5% FPR points for simulated z-stats:
    4.290 = 1-sided 5% FPR
    4.460 = 2-sided 5% FPR
  • [above results also in file TT.ewJRndftwVF.5percent.txt]
    
  • ===== starting 3dXClustSim : elapsed = 467.0 s =====
    ++ 3dXClustSim: AFNI version=AFNI_17.1.12 (Jun 23 2017) [64-bit]
    ++ Authored by: Lamont Cranston
    ++ FPR goal of 5.00% will be used
    ++ Loading -insdat datasets
    ++ 3dXClustSim: Using 60 OpenMP threads
    ++ STEP 1a: start 2-sided clustering with NN=2
  • STEP 1b: merge cluster lists
  • STEP 1c: compute minimum thresholds
    ++ STEP 2: start cluster dilations
    ++ STEP 3: re-loading & sorting FOM vectors after dilations
    ++ STEP 4: adjusting per-voxel FOM thresholds to reach FPR=5.00%
  •  #1: Testing threshold images at 0.0006 ==> 24
    
  •      FPR = 1.14%
    
  •  #2: Testing threshold images at 0.0013332 ==> 54
    
  •      FPR = 3.65%
    
  •  #3: Testing threshold images at 0.00182568 ==> 73
    
  •      FPR = 5.57%
    
  •  #4: Testing threshold images at 0.00167963 ==> 68
    
  •      FPR = 5.02%
    

** ERROR: nifti_image_write_hdr_img: NBL does not match nim
++ Output dataset ./TT.ewJRndftwVF.default.ETAC.mthresh.A.nii
++ === 3dXClustSim ends: Elapsed time = 236.0 s
++ — merging 1 blur cases to make pos 1-sided activation mask —

Fatal Signal 11 (SIGSEGV) received
3dttest++ main
Bottom of Debug Stack
** Command line was:
3dttest++ -prefix ALLFACEvshape_3dttestETAC -mask group_mask_AUDanxiety+tlrc -ETAC -setA ‘/raid10/studies/CII/results/FaceMatching_glt/NV*/stats.NV._REML+tlrc.BRIK[Faces-Base#0_Coef]’ -setB '/raid10/studies/CII/results/FaceMatching_glt/alc/stats.alc.*_REML+tlrc.BRIK[Faces-Base#0_Coef]’ -labelA HC -labelB AUD
** AFNI version = AFNI_17.1.12 Compile date = Jun 23 2017
** [[Precompiled binary linux_fedora_21_64: Jun 23 2017]]
** Program Death **
** If you report this crash to the AFNI message board,
** please copy the error messages EXACTLY, and give
** the command line you used to run the program, and
** any other information needed to repeat the problem.
** You may later be asked to upload data to help debug.
** Memory usage: chunks=84771 bytes=24014835
** Crash log is appended to file /home/fedesj/.afni.crashlog

Hello,

That is an impressive form of -setA and -setB! :slight_smile:

To be sure, what does this show?

3dinfo group_mask_AUDanxiety+tlrc | head -n 25

And how about the same for one of the inputs?

There is also an error about an existing output file.

It looks like you are at the NIH. If you are in
building 10 and would like me to take a peek at this,
just let me know.

  • rick