Changeset 511 for trunk/lib/statistics/utility.cc

Timestamp:

Feb 18, 2006, 3:46:46 PM (16 years ago)

Author:

Peter

Message:

added percentile function for gslapi::vector

File:

• trunk/lib/statistics/utility.cc (modified) (1 diff)

trunk/lib/statistics/utility.cc

@@ -32 +32 @@
   }
 
+  double percentile(const gslapi::vector& vec, const double p, 
+                    const bool sorted)
+  {
+    if (!sorted){
+      gslapi::vector vec_c(vec);
+      vec_c.sort();
+      return gsl_stats_quantile_from_sorted_data(vec_c.gsl_vector_p()->data,
+                                                 vec_c.gsl_vector_p()->stride,
+                                                 vec_c.gsl_vector_p()->size,
+                                                 p);
+    }
+    return gsl_stats_quantile_from_sorted_data (vec.gsl_vector_p()->data, 
+                                              vec.gsl_vector_p()->stride,
+                                              vec.gsl_vector_p()->size,
+                                              p);
+  }
+
 }} // of namespace statistics and namespace theplu