ComputePValuesNodalFeatures.m
257e25cb59ccmaster

ComputePValuesNodalFeatures.m
View Options

			function [P_final] = ComputePValuesNodalFeatures(Feature,Group_Label,N_null,n_TP,n_ROI)

			% Bonferroni correction (across time windows and regions)
			for t = 1:n_TP
			for r = 1:n_ROI
			p_eig(t,r) = AssessGroupDifference(squeeze(Feature(r,:,t)),Group_Label,N_null);
			end
			end

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

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

			for delta = 80:-1:40

			% Construction of the index vector for time
			IV = ConstructIndexVector(n_TP,delta) + 1i*ones(1,n_TP);

			% Construction of the index vector for space
			IV2 = ones(1,n_ROI) + 1i*ConstructIndexVector(n_ROI,delta);

			% Generation of the final index vector of interest (size n_TP x
			% n_ROI)
			IV3 = IV' * IV2;
			IV4 = nan(n_TP*n_ROI,1);

			tmp = unique(IV3(:));

			for i = 1:length(tmp)
			IV4(IV3(:) == tmp(i)) = i;
			end

			IV4 = reshape(IV4,n_TP,n_ROI);

			% On a concatenated version, I find my p-values
			tmp_p = relaxPvalues(p_eig_BON(:),IV4(:),0.05,0.05,'Bonferroni',false);

			P_final = P_final + MakePVector2(tmp_p,IV4(:),n_TP,n_ROI);

			end

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

			end