diff --git a/createTestData/evaluateProFitFittingParameters.m b/createTestData/evaluateProFitFittingParameters.m
index 227baa767be233095425d6fa0bed665d415b2985..cba0a1e0de45c4ae66717290f08b78710643c4d4 100644
--- a/createTestData/evaluateProFitFittingParameters.m
+++ b/createTestData/evaluateProFitFittingParameters.m
@@ -4,8 +4,8 @@ folderName = 'final_simulations';
 exportFolder = {'11 Simulations'};
 pathNameExport = 'ProFit test data';
 pathBaseExportFiles = pathToDataFolder(pathNameExport, exportFolder);
-pathBaseExportFiles = pathBaseExportFiles(1:end-1);
-pathBaseExportFiles = [pathBaseExportFiles ' - Weighted Cost Fun R3 0.25scale m_15 trunc\'];
+% pathBaseExportFiles = pathBaseExportFiles(1:end-1);
+% pathBaseExportFiles = [pathBaseExportFiles ' - Weighted Cost Fun R3 0.25scale m_15 trunc\'];
 dataExportPathBase = strcat(pathBaseExportFiles, exportFolder{1}, '\');
 truncSuffix = '_truncOn';
 folderFigures = [pathBaseExportFiles 'Output\' folderName, truncSuffix, '\'];
@@ -40,7 +40,7 @@ for iParam = 1: numParam
     legendLabels_ = {};
     unitLabel = [titleLabels{iParam}, units{iParam}];
     for iFolder = 1:length(folders)
-        load([dataExportPathBase, folders{iFolder}, '\Diff_Summaries', truncSuffix, '_no_BL_df3.mat'], 'diff_conc_mat', 'diff_em_mat', 'diff_gm_mat', 'diff_em_g_mat', ...
+        load([dataExportPathBase, folders{iFolder}, '\Diff_Summaries', truncSuffix, '_no_BL_df3.mat'], 'diff_conc_mat', 'diff_em_mat', 'diff_gm_mat', ...
             'diff_pc0_mat', 'diff_df2_mat', 'diff_pc12_mat', 'diff_T2_mat');
         size1 = size(diff_conc_mat, 1);
         size2 = size(diff_conc_mat, 2);
@@ -50,8 +50,10 @@ for iParam = 1: numParam
         
         eval(['diff_mat = diff_', parameters{iParam}, '_mat;']);
         
-        diff_mean = mean(mean(abs(diff_mat)));
-        diff_std = std(mean(abs(diff_mat)));
+        diff_mean = mean(mean(abs(diff_mat), 'omitnan'));
+        diff_std = std(mean(abs(diff_mat), 'omitnan'));
+        numberOfNaNs = sum(sum(isnan(diff_mat)))
+        numberOfNaNsV = sum((isnan(diff_mat)),2)'
         
         if useSubPlots
             subplot(3,2,iParam)