More robust pvalue sorting

Author: bogpetre

CanlabCore/Statistics_tools/holm_sidak.m

@@ -34,7 +34,7 @@
 
 pthr = 1 - (1-alpha)^(1/length(sortedP)); % max P-value for sig results. Anything below this is significant.
 
-while sortedP(end) < 1 - (1-alpha)^(1/length(sortedP))
+while ~isempty(sortedP) && sortedP(end) < 1 - (1-alpha)^(1/length(sortedP))
 
     sig(argsort(end)) = true;
     sortedP(end) = [];