Dear all, I recently have come across an error in the ‘fitting water reference’ stage of the Ospreyfit analysis pipeline.
The error seems to occur only when I have earlier selected a custom metabolite basis set made up of 26 metabolites (including macromolecules) in the job file. The error does not occur if I select the custom metabolite basis set and then include all metabolites available (36) in the set or if I run it with the metabolites included in the default model plus one extra (glucose)(28 metabolites included). It is drawing from the Philips PRESS, TE 40 set and I am running it through the batch processing commands (Ospreyjob, Ospreyload, Ospreyprocess, Ospreyfit etc), not the gui.
Please see below for the error. Any help/suggestions of what I may be doing, or what is going wrong are appreciated.
Fitting water reference from dataset 1 out of 30 total datasets…
Index in position 2 exceeds array bounds.
Error in fit_resampleBasis (line 48)
specs_interp(:,ll,rr) = interp1(ppmRangeBasis(ppmIsInDataRange), basisSet.specs(ppmIsInDataRange,ll,rr), ppmRangeData, ‘pchip’,
Error in fit_runFitWater (line 38)
resBasisSet = fit_resampleBasis(dataToFit, basisSet);
Error in osp_fitWater (line 86)
[fitParamsWater, resBasisSetWater] = fit_runFitWater(dataToFit, basisSet, fitModel, fitOpts);
Error in OspreyFit (line 102)
[MRSCont] = osp_fitWater(MRSCont, kk, ‘ref’);
All the best,