Hello Flora and Paul, I believe we are having a similar issue with step 2 of FABBER, however ours is occurring only at the partial volume correction stage. If this is a different issue I can start a new thread. We have attempted running through about 5 sets of data and we have the same issue with each one. Here is the information I have and I would be willing to share our data if you can tell me how. fsl_ wb1DNR _ox_asl/basil/step1 finalMVN.nii.gz logfile noise_means.nii.gz std_delttiss.nii.gz zstat_delttiss.nii.gz freeEnergy.nii.gz mean_delttiss.nii.gz noise_stdevs.nii.gz std_ftiss.nii.gz zstat_ftiss.nii.gz info.txt mean_ftiss.nii.gz paramnames.txt uname.txt fsl_ wb1DNR _ox_asl/basil/step2 finalMVN.nii.gz mean_delttiss.nii.gz paramnames.txt var_delttiss.nii.gz freeEnergy.nii.gz mean_ftiss.nii.gz std_delttiss.nii.gz var_ftiss.nii.gz info.txt noise_means.nii.gz std_ftiss.nii.gz zstat_delttiss.nii.gz logfile noise_stdevs.nii.gz uname.txt zstat_ftiss.nii.gz fsl_ wb1DNR _ox_asl/pvcorr/basil/step1 finalMVN.nii.gz mean_ftiss.nii.gz paramnames.txt std_pvwm.nii.gz zstat_pvgm.nii.gz freeEnergy.nii.gz mean_fwm.nii.gz std_delttiss.nii.gz uname.txt zstat_pvwm.nii.gz info.txt mean_pvgm.nii.gz std_deltwm.nii.gz zstat_delttiss.nii.gz logfile mean_pvwm.nii.gz std_ftiss.nii.gz zstat_deltwm.nii.gz mean_delttiss.nii.gz noise_means.nii.gz std_fwm.nii.gz zstat_ftiss.nii.gz mean_deltwm.nii.gz noise_stdevs.nii.gz std_pvgm.nii.gz zstat_fwm.nii.gz fsl_ wb1DNR _ox_asl/pvcorr/basil/step2 info.txt logfile paramnames.txt uname.txt STEP 2: Spatial VB Tissue PVE - init with STEP 1 ---------------------- Welcome to FABBER v3.9.2-78-g7e4027c ---------------------- Last commit: Thu Jan 25 13:00:17 2018 Logfile started: /tmp/fsl_wb1DNR_ox_asl/pvcorr/basil/step2/logfile 1%35371370231.389931 47029505458.779755 1.300000 1.600000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 10.014421 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 10.014421 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 76011042523822017122425372672.000000 96576635593542070815122522112.000000 0.000000 0.000000 0.000000 0.000000 96576635593542070815122522112.000000 122706467808867387843872292864.000000 11893025264.716621 103846276214.940628 1.300000 1.600000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 10.014421 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 10.014421 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 8844443584466570525867507712.000000 73798232007530254509063274496.000000 0.000000 0.000000 0.000000 0.000000 73798232007530254509063274496.000000 615774072775176728573502291968.000000 35371370231.389931 47029505458.779755 1.300000 1.600000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 10.014421 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 10.014421 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 76011042523822017122425372672.000000 96576635593542070815122522112.000000 0.000000 0.000000 0.000000 0.000000 96576635593542070815122522112.000000 122706467808867387843872292864.000000 NEWMAT exception caught in fabber: An exception has been thrown Runtime error:- detected by Newmat: matrix is singular MatrixType = Crout # Rows = 6; # Cols = 6 Trace: Crout(lubksb); GeneralSolvI; InvertedMatrix::Evaluate. Final logfile: /tmp/fsl_wb1DNR_ox_asl/pvcorr/basil/step2/logfile End. Cannot open volume /tmp/fsl_wb1DNR_ox_asl/pvcorr/basil/step2/mean_ftiss for reading! Cannot open volume /tmp/fsl_wb1DNR_ox_asl/pvcorr/basil/step2/mean_delttiss for reading! Cannot open volume /tmp/fsl_wb1DNR_ox_asl/pvcorr/basil/step2/mean_fwm for reading! Cannot open volume /tmp/fsl_wb1DNR_ox_asl/pvcorr/basil/step2/mean_deltwm for reading! Thank you for any help you can provide, Anthony On Thu, Jul 12, 2018 at 4:52 AM, Flora Kennedy McConnell < [log in to unmask]> wrote: > Hi Paul, > > > Thanks for this. We can see where the error is arising - for at least one > voxel the number of neighbouring voxels available for consideration in the > spatial VB step is zero even though that voxel should be considered as one > of its own neighbours. > > > We're really not sure what about your data could cause this. Would you be > willing to share one of the datasets with us so we can have a look at it? > > > All the best, > > Flora > > > Flora Kennedy McConnell > ------------------------------ > Postdoctoral Research Associate > > Institute of Biomedical Engineering <http://www.ibme.ox.ac.uk/> > Department of Engineering Science > University of Oxford > Old Road Campus Research Building > Oxford, UK > OX3 7DQ > ------------------------------ > *From:* FSL - FMRIB's Software Library <[log in to unmask]> on behalf of > Paul Mullins <[log in to unmask]> > *Sent:* 12 July 2018 09:27:40 > *To:* [log in to unmask] > *Subject:* Re: [FSL] Oxford_ASL - fabber errors. > > Hi Flora, > > we always use a Mac. Have run the analysis successfully on several > different macs with different data sets (all produced in the same fashion > though, and all from the same study). Ran the command again with debug on > as you requested. > > The files in Step 1 are : > finalMVN.nii.gz info.txt mean_delttiss.nii.gz > paramnames.txt zstat_delttiss.nii.gz > freeEnergy.nii.gz logfile mean_ftiss.nii.gz > uname.txt zstat_ftiss.nii.gz > > and step 2: > info.txt logfile uname.txt > > Hope this helps figure it out. > Paul. > > ######################################################################## > > To unsubscribe from the FSL list, click the following link: > https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=FSL&A=1 > > ------------------------------ > > To unsubscribe from the FSL list, click the following link: > https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=FSL&A=1 > ######################################################################## To unsubscribe from the FSL list, click the following link: https://www.jiscmail.ac.uk/cgi-bin/webadmin?SUBED1=FSL&A=1