Fixed bug in sort

groupttest.cc

@@ -214,6 +214,7 @@ double pvalue(double tX, double dof) {
 } 
 
 
+
 vector<int> get_sortindex(const Matrix& vals)
 {
   // return the mapping of old indices to new indices in the
@@ -226,12 +227,11 @@ vector<int> get_sortindex(const Matrix& vals)
   sort(sortlist.begin(),sortlist.end());  // O(N.log(N))
   vector<int> idx(length);
   for (int n=0; n<length; n++) {
-    idx[n] = sortlist[n].second;
+    idx[sortlist[n].second-1] = n+1;
   }
   return idx;
 }
 
-
 ////////////////////////////////////////////////////////////////////////////
 
 // Main function - this does all the work