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  

|
February 11, 2020 09:14AM
I posted a few days ago about a new problem Im encountering using 3dMVM. Rick Reynolds has been so kind to help troubleshoot (thank you, Rick), but unfortunately I am still having the same problem. I wanted to post a new topic with more information in hopes of increasing visibility to others that may have experience with my situation.

I have updated afni and r packages; when I run afni_system_check.py, there are no problems found. I also deleted everything in /Library/Frameworks/R.framework/ and re-installed R. I am on a 2018 iMac, quad core i7, 32 gb DDR4. I am re-running a 2x2x2 RM ANOVA using 3dMVM with 35 participant data sets (I successfully ran this 1 year ago). The only difference in re-running this ANOVA is the -resid option. This is the top of my 3dmvm txt file:


#!/bin/tcsh

3dMVM -mask clippedfullmaskCV20.nii -prefix Accuracy_CV20 -jobs 4 \
-model 1 \
-resid /Users/jessicadefenderfer/Desktop/CV15_WAV72/ACC_CV20_RESID.nii \
-wsVars "oxy*acc*noiseType" \
-num_glt 3 \
-gltLabel 1 OxyEff -gltCode 1 'oxy : 1*HbO -1*HbR' \
-gltLabel 2 AccEff -gltCode 2 'acc : 1*right -1*wrong' \
-gltLabel 3 NoiseTypeEff -gltCode 3 'noiseType : 1*vocod -1*noise' \
-dataTable \
Subj oxy acc noiseType InputFile \
S1001 HbO wrong vocod /Users/jessicadefenderfer/Desktop/CV15_WAV72/S1001/viewer/Subject/1001_C7_unmaskO_CV20.nii \
S1001 HbO wrong noise /Users/jessicadefenderfer/Desktop/CV15_WAV72/S1001/viewer/Subject/1001_C6_unmaskO_CV20.nii \
S1001 HbO right vocod /Users/jessicadefenderfer/Desktop/CV15_WAV72/S1001/viewer/Subject/1001_C5_unmaskO_CV20.nii \
S1001 HbO right noise /Users/jessicadefenderfer/Desktop/CV15_WAV72/S1001/viewer/Subject/1001_C4_unmaskO_CV20.nii \
S1001 HbR wrong vocod /Users/jessicadefenderfer/Desktop/CV15_WAV72/S1001/viewer/Subject/1001_C7_unmaskR_CV20.nii \
S1001 HbR wrong noise /Users/jessicadefenderfer/Desktop/CV15_WAV72/S1001/viewer/Subject/1001_C6_unmaskR_CV20.nii \
S1001 HbR right vocod /Users/jessicadefenderfer/Desktop/CV15_WAV72/S1001/viewer/Subject/1001_C5_unmaskR_CV20.nii \
S1001 HbR right noise /Users/jessicadefenderfer/Desktop/CV15_WAV72/S1001/viewer/Subject/1001_C4_unmaskR_CV20.nii \

etc...


3dMVM is able to read all my input files, and processes all 256 slices... once it gets to the 256th slice, it throws the error below. I can only guess its a memory problem, so I started with 6 participants in the ANOVA and through painstaking trial and error, have found that I can include up to 14 subject datasets before I get some version of this error:

*** caught segfault ***
address 0xfffffffd2524d030, cause 'memory not mapped'
Traceback:
1: write.c.AFNI(filename, dset = brk, label = label, space = space, note = note, origin = origin, delta = delta, orient = orient, idcode = idcode, defhead = defhead, verb = verb, maskinf = maskinf, scale = scale, addFDR = addFDR, statsym = statsym, view = view, com_hist = com_hist, type = type, TR = TR, overwrite = overwrite)
2: write.AFNI(lop$outFN, out, brickNames, defhead = head, idcode = newid.AFNI(), com_hist = lop$com_history, statsym = statsym, addFDR = 1, type = "MRI_short", overwrite = lop$overwrite)
An irrecoverable exception occurred. R is aborting now ...
Segmentation fault

Without that -resid option, the ANOVA finishes without error. However, my memory pressure monitor shows no sign of running out memory and CPU load hardly reaches 15% usage running it with or without -resid.

During one of the successful runs with fewer participants, I noticed something odd in the output (see attachment 1, 'congrats'). It warns me that the dataset name "Accuracy_cv20" and the resid filename conflicts with an existing file, so they werent written to disk. But there were no files of that name in this directory. Further, when Im monitoring the CPU/memory usage in the Mac's activity monitor, R appears up to 7 or 8 times in the process column. Is it possible I have some weird program settings in R that has it trying to write (or re-write over) the same ANOVA output before its even finished? I can't make any sense of this.

Thank you for any help. -Jess
Attachments:
open | download - congrats.png (378.1 KB)
open | download - R_multiprocess.png (727.8 KB)
Subject Author Posted

3dMVM R error *** caught segfault *** cause 'memory not mapped' Attachments

defenderfer February 11, 2020 09:14AM

Re: 3dMVM R error *** caught segfault *** cause 'memory not mapped'

gang February 12, 2020 04:28PM

Re: 3dMVM R error *** caught segfault *** cause 'memory not mapped'

defenderfer February 13, 2020 08:50AM