gaussComparer.h

/*  gaussComparer.h

    Mark Woolrich, FMRIB Image Analysis Group

    Copyright (C) 1999-2000 University of Oxford  */

/*  CCOPYRIGHT  */

#include <iostream>
#include <fstream>
#define WANT_STREAM
#define WANT_MATH

#include "histogram.h"
#include "newmatap.h"
#include "newmatio.h"

extern "C"
{
#include "libprob.h"
};

#ifndef NO_NAMESPACE
using namespace NEWMAT;
using namespace SIGPROC;
namespace TACO{
#endif

#if !defined(__gaussComparer_h)
#define __gaussComparer_h
 
  class GaussComparer
    {
    public:
      GaussComparer(ColumnVector& pzscores) 
	: zscores(pzscores), histogram(zscores, -10.0F, 10.0F, 2000) {}
    
      void setup(){histogram.generate();} // sets up histogram
     
      float computeRatio(const float prob, ofstream& out); // returns ratio of tail size in zscores compared to a Gaussian distribution

    private:
      GaussComparer();
      const GaussComparer& operator=(GaussComparer&);
      GaussComparer(GaussComparer&);
     
      const ColumnVector& zscores;

      Histogram histogram;
    };


#ifndef NO_NAMESPACE
       }
#endif

#endif