Skip to content
Snippets Groups Projects
Commit 41e115ba authored by ac_fx's avatar ac_fx
Browse files

draw_ssf_ncd with new csf_spatial

parent 632814ba
No related branches found
No related tags found
No related merge requests found
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')); flist = dir(fullfile(ROOT,'*.LEVER'));
[classUnique,ia,ic]= unique(classList); [classUnique,ia,ic]= unique(classList);
...@@ -6,7 +6,7 @@ idxTrue = ic; ...@@ -6,7 +6,7 @@ idxTrue = ic;
[idx,Y] = Cluster.SpectralCluster(A,K); [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); % Y = tsne(Y,'Algorithm','exact','NumDimensions',3);
clf;hold off; clf;hold off;
...@@ -43,6 +43,7 @@ if ~isempty(idxWT) ...@@ -43,6 +43,7 @@ if ~isempty(idxWT)
classUnique{idxWT} = 'WILD TYPE'; classUnique{idxWT} = 'WILD TYPE';
end end
legend(hx(find(hx)),classUnique(find(hx)),'interpreter','none'); 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) title(tstr)
4;
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment