diff --git a/matlab/+SSF/draw_ssf_ncd.m b/matlab/+SSF/draw_ssf_ncd.m
index b58fff8210aef06859c56e7448107592c8e79ed4..788c609b158bcbbfd0647f590688a07073de505d 100644
--- a/matlab/+SSF/draw_ssf_ncd.m
+++ b/matlab/+SSF/draw_ssf_ncd.m
@@ -1,4 +1,4 @@
-function [csf_mean,csf_std] = draw_ssf_ncd(ROOT,classList,A,K)
+function [csf_mean,csf_std,csf] = draw_ssf_ncd(ROOT,classList,A,K)
 
 flist = dir(fullfile(ROOT,'*.LEVER'));
 [classUnique,ia,ic]= unique(classList);
@@ -6,7 +6,7 @@ idxTrue = ic;
 
 [idx,Y] = Cluster.SpectralCluster(A,K);
 
-[csf_mean, csf_std] = CSF.csf_spatial(Y,idxTrue);
+csf = CSF.csf_spatial(Y,idxTrue);
 
 % Y = tsne(Y,'Algorithm','exact','NumDimensions',3);
 clf;hold off;
@@ -43,6 +43,7 @@ if ~isempty(idxWT)
     classUnique{idxWT} = 'WILD TYPE';
 end
 legend(hx(find(hx)),classUnique(find(hx)),'interpreter','none');
-tstr = ['CSF^6_{Erk} = ' jsonencode(round([csf_mean,csf_std],2,'significant'))]; 
+tstr = ['CSF^6_{Erk} = ' jsonencode(round([mean(csf),std(csf)],2,'significant'))]; 
 title(tstr)
 
+4;