ComputePValuesDissimilarity.m
257e25cb59ccmaster

ComputePValuesDissimilarity.m
View Options

			function [P_final] = ComputePValuesDissimilarity(IDX_2D,Group_Label,N_null,n_TP,n_CTR,n_ASD)

			for t = 1:n_TP
			DG(t) = AssessGroupSim(IDX_2D(t,Group_Label==1),IDX_2D(t,Group_Label==2),n_CTR,n_ASD);
			end

			for null = 1:N_null

			null

			idd = randperm(n_CTR+n_ASD);
			tmp_idx2 = IDX_2D(:,idd);

			for t = 1:n_TP
			DG_null(t,null) = AssessGroupSim(tmp_idx2(t,Group_Label==1),tmp_idx2(t,Group_Label==2),n_CTR,n_ASD);
			end
			end

			for t = 1:n_TP
			p_DIS(t) = sum(abs(DG_null(t,:)) > abs(DG(t)))/N_null;
			end

			p_eig_BON = p_DIS * 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