wilcoxon_table.cpp

非常著名的曲线拟合程序

// This is -*- C++ -*-
// $Id: wilcoxon_table.cpp,v 1.3 1999/04/08 18:09:47 trow Exp $

#include <math.h>
#include <iostream>
#include <specfns.h>
#include "wilcoxon.h"

double foo[] = {4.9, 6.0, 6.9, 17.6, 4.5, 12.3, 5.7, 5.3, 9.6, 13.5, 15.7, 7.7, -1};

main()
{
  /*
  RealSet rs;
  int i=0;
  while (foo[i] > 0)
    rs.add(foo[i++]);

  cout << wilcoxon_statistic(rs, 12.0) << endl;
  */
  
  for (unsigned n=1; n<51; ++n) {
    double max_err = 0;
    for (unsigned k=0; k<=n*(n+1)/2; ++k) {
      double p = wilcoxon_noties_cdf(n,k);
      double pest = normal_cdf((k+0.5-0.25*n*(n+1))/sqrt((n/24.0)*(n+1)*(2*n+1)));
      double err = fabs(p-pest);
      if (err > max_err) max_err = err;
    }
    cout << n << " " << max_err << endl;
  }
   
}




// $Id: wilcoxon_table.cpp,v 1.3 1999/04/08 18:09:47 trow Exp $