From 91c7bd065174f861bd3e64c289e417d1a990e02f Mon Sep 17 00:00:00 2001
From: Bertrand <bertrand.horel@gmail.com>
Date: Tue, 1 Mar 2016 14:37:15 +0000
Subject: rajout quantile_norm

---
 src/stratified_sampling.cpp | 15 ++++++---------
 1 file changed, 6 insertions(+), 9 deletions(-)

(limited to 'src/stratified_sampling.cpp')

diff --git a/src/stratified_sampling.cpp b/src/stratified_sampling.cpp
index 1c8e670..9586a86 100644
--- a/src/stratified_sampling.cpp
+++ b/src/stratified_sampling.cpp
@@ -1,16 +1,13 @@
 #include "stratified_sampling.hpp"
 
 
-std::pair<double, double> mean_var( std::vector<double> r){
-    std::pair<double, double> p;
-    for(auto &x: r){
-        p.first += x;
-        p.second += x*x;
+//--génération quantiles--
+vector<double> quantile_norm(int n, double sigma){
+    vector<double> q(n);
+    for (int i=0; i<n; i++) {
+        q[i] = gsl_cdf_gaussian_Pinv ((double)(i+1)/n, sigma);
     }
-    p.first /= r.size();
-    p.second /= r.size();
-    p.second -= p.first * p.first;
-    return p;
+    return q;
 }
 
 
-- 
cgit v1.2.3-70-g09d2