Changeset 76

Timestamp:

May 4, 2004, 11:52:46 AM (19 years ago)

Author:

Peter

Message:

updated the test program

File:

• trunk/test/test_roc.cc (modified) (1 diff)

trunk/test/test_roc.cc

@@ -16 +16 @@
 
 {  
-  const double x_tmp = 1.0;
-  const double sigma_tmp = 1.0;
-  double p_tmp = gsl_cdf_gaussian_Q(x_tmp, sigma_tmp);
   bool ok = true;
   theplu::gslapi::vector value(31);
   theplu::gslapi::vector target(31,1);
   for (unsigned int i=0; i<16; i++) 
-    target(2*i) = -1;
+    target(i) = -1;
   for (unsigned int i=0; i<value.size(); i++) 
     value(i)=i;
   theplu::cpptools::ROC roc(target, value);
   double area = roc.area();
-  if (area!=0.5){
-    cerr << "test_roc: area should be 0.5" << endl;
+  if (area!=1.0){
+    cerr << "test_roc: area should be 1.0" << endl;
     ok = false;
   } 
   double p = roc.get_p(area);
-  if (p<0.5 | p>0.51){
-    cerr << "test_roc: p-value not correct" << endl;
+  double p_matlab = 0.00000105;
+  if (p/p_matlab > 1.01 | p/p_matlab < 0.99){
+    cerr << "get_p_approx: p-value not correct" << endl;
     ok = false;
   }
-//  target.set_all(1);
-//  for (unsigned int i=0; i<16; i++) 
-//    target(i) = -1;
-//  roc = theplu::cpptools::ROC(target, value);
-  
+  roc.minimum_size(20);
+  p = roc.get_p(area);
+  if (p > pow(10, -8.0) | p < pow(10, -9.0)){
+    cerr << "get_p_exact: p-value not correct" << endl;
+    ok = false;
+  }
   
   if (ok)