RSCHMM_DetermineCommunities_Experimental.m
1f50b9a012abmaster

RSCHMM_DetermineCommunities_Experimental.m
View Options

			function [T,Z,Cutoff] = RSCHMM_DetermineCommunities_Experimental(G)

			R = size(G,1);
			G2 = NaN(R-1,R);

			for i = 1:size(G,2)
			tmp = G(:,i);
			tmp(i) = [];
			G2(:,i) = tmp;
			end

			Z = linkage(G2','ward');

			for n = 1:10000

			G_null = G2;

			for r = 1:R
			G_null(:,r) = G_null(randperm(R-1),r);
			end

			Z_null = linkage(G_null','ward');
			D_null(n) = max(Z_null(:,3));
			end

			Cutoff = prctile(D_null,95);

			T = cluster(Z,'Cutoff',Cutoff,'Criterion','distance');
			end