ComputePValuesNetworkFeature.m
View Options

			function [P_final] = ComputePValuesNetworkFeature(Medeig_plus,Group_Label,N_null,n_TP)

			% Bonferroni correction
			for t = 1:n_TP
			p_eig(t) = AssessGroupDifference(Medeig_plus(:,t),Group_Label,N_null);
			end

			p_eig_BON = p_eig * n_TP;
			p_eig_BON(p_eig_BON > 1) = 1;

			% Will contain the final p values
			P_final = zeros(1,n_TP);

			for delta = 80:-1:40

			% Construction of the index vector
			IV = ConstructIndexVector(n_TP,delta);

			tmp_p = relaxPvalues(p_eig_BON,IV,0.05,0.05,'Bonferroni',false);

			P_final = P_final + MakePVector(tmp_p,IV);

			end

			P_final = P_final/length(80:-1:40);

			end

ISFCClustering/Matlab_Scripts/Functions/ComputePValuesNetworkFeature.m257e25cb59ccmaster