Added NaNs for low concentration (=0) metabolites when plotted as difference spectrum.

createTestData/compareFitInput_Output.m

@@ -3,14 +3,14 @@ function [diff_conc, diff_em, diff_gm, diff_pc0, diff_df2, diff_pc12, diff_T2] =
 saveFigures = false;
 
 fileNameOff = strrep(filenameData,'On', 'Off');
-ls([dataExportPath, fileNameOff '*']);
+xx = ls([dataExportPath, fileNameOff '*']);
 
 load([dataExportPath, fileNameOff, '.mat'], 'summedFidWithNoise', 'current_em', 'currentConcentrationRm', ...
     'current_pc0', 'current_pc1', 'current_df2', 'current_gm', 'current_em_std_factor', 'conc_std_factor');
 
 fileNameProFit = [filenameData '_profit'];
 % close all;
-load(strcat(dataExportPath,fileNameProFit, '.mat'),'fitresult', 'data');
+load(strcat(dataExportPath,fileNameProFit, '.mat'),'fitresult');
 
 currentFitresult = fitresult{1,iteration};
 
@@ -27,6 +27,7 @@ end
 
 currentConcentrationRm(end) = currentConcentrationRm(end)*1e-8;
 diff_conc = currentFitresult.fitted_values.conc * currentFitresult.sub_mask.conc.full - currentConcentrationRm(sortedIndeces);
+diff_conc(currentFitresult.fitted_values.conc==0)=NaN;
 diff_em   = currentFitresult.fitted_values.em * currentFitresult.sub_mask.em.full - current_em(sortedIndeces);
 diff_gm   = currentFitresult.fitted_values.gm * currentFitresult.sub_mask.gm.full- current_gm;
 diff_pc0  = currentFitresult.fitted_values.pc0 * currentFitresult.sub_mask.pc0.full - current_pc0;