hist.cpp

Go to the documentation of this file.

#include "StdAfx.h"

//#define WANT_STREAM
#define WANT_MATH

#include "include.h"
#include "extreal.h"
#include "newran.h"

#ifdef use_namespace
using namespace NEWRAN;
#endif

#define al_trace        printf

void Histogram(Random* rx, int n)          // draw histogram with n obsv
{
   int i,j; int count[20];
   Real* a = new Real[n];
   for (i = 0; i < n; i++) a[i] = rx->Next();
   Real amax = a[0]; Real amin = a[0]; Real mean = a[0]; Real sd = 0;
   for (i = 1; i < n; i++)
   {
      if (amin > a[i]) amin = a[i]; else if (amax < a[i]) amax = a[i];
      mean += a[i];
   }
   mean /= n;
   for (i = 0; i < 20; i++) count[i]=0;
   for (i = 0; i < n; i++)
   {
     Real rat= (amax != amin) ? (a[i] - amin)/(amax - amin) : 1.0;
     j = (int)( 19.999 * rat ); count[j]++;
     Real diff = a[i] - mean; sd += diff*diff;
   }
   sd = sqrt(sd/(n-1));
   j = 0;
   for (i = 0; i < 20; i++) { if (j < count[i]) j = count[i]; }
   if (j > 70) { for (i = 0; i < 20; i++) count[i] = (int)((70L*count[i])/j); }

   al_trace( "\n" );
   for (i = 0; i < 20; i++)
        {
                al_trace( "\n|" );
                for (j = 1; j < count[i]; j = j+1) 
                        al_trace( "*" );
        }
/*   al_trace( "\n%s\np. mean = %9f, p. var = %9f\ns. mean = %9f, s. var = %9f\n" , 
                rx->Name() , rx->Mean() , rx->Variance() , mean , sd * sd );*/

   delete a;
}

void Hist(Random* rx, int n , const int cnt )          // draw histogram with n obsv
{
   int i,j; int *count = new int[cnt];
   Real* a = new Real[n];
   for (i = 0; i < n; i++) a[i] = rx->Next();
   Real amax = a[0]; Real amin = a[0]; Real mean = a[0]; Real sd = 0;
   for (i = 1; i < n; i++)
   {
      if (amin > a[i]) amin = a[i]; else if (amax < a[i]) amax = a[i];
      mean += a[i];
   }
   mean /= n;
   for (i = 0; i < cnt ;  i++) count[i]=0;
   for (i = 0; i < n; i++)
   {
     Real rat= (amax != amin) ? (a[i] - amin)/(amax - amin) : 1.0;
     j = (int)( (cnt - 0.001) * rat ); count[j]++;
     Real diff = a[i] - mean; sd += diff*diff;
   }
   sd = sqrt(sd/(n-1));
   j = 0;
   for (i = 0; i < cnt; i++) { if (j < count[i]) j = count[i]; }
   if (j > 70) { for (i = 0; i < cnt; i++) if( count[i] )       count[i] = (int)((70L*count[i])/j); }

   al_trace( "\n" );
   for (i = 0; i < cnt; i++)
        {
                al_trace( "\n%02d|" , i );
                for (j = 1; j < count[i] ; j = j+1) 
                        al_trace( "*" );
        }
//   al_trace( "\n%s\np. mean = %9f, p. var = %9f\ns. mean = %9f, s. var = %9f\n" , 
//              rx->Name() , rx->Mean() , rx->Variance() , mean , sd * sd );

        delete count;
   delete a;
}
#ifdef PLATFORM_WIN32
void alHistogram( Random* rx , int n , Real *a , const int ncount )          // draw histogram with n obsv
{
   int i,j; int *count = new int[ncount];

   Real amax = a[0]; Real amin = a[0]; Real mean = a[0]; Real sd = 0;
   for (i = 1; i < n; i++)
   {
      if (amin > a[i]) amin = a[i]; else if (amax < a[i]) amax = a[i];
      mean += a[i];
   }
   mean /= n;
   for (i = 0; i < ncount; i++) count[i]=0;
   for (i = 0; i < n; i++)
   {
     Real rat= (amax != amin) ? (a[i] - amin)/(amax - amin) : 1.0;
     j = (int)( (double( ncount ) - 0.001) * rat ); count[j]++;
     Real diff = a[i] - mean; sd += diff*diff;
   }
   sd = sqrt(sd/(n-1));
/*
   j = 0;
   for (i = 0; i < ncount; i++) { if (j < count[i]) j = count[i]; }
   if (j > 70) { for (i = 0; i < ncount; i++) count[i] = (int)((70L*count[i])/j); }
*/
   al_trace( "\n" );
   for (i = 0; i < ncount; i++)
        {
                al_trace( "\n%04d|" , i );
                for (j = 1; j < count[i]; j = j+1) 
                        al_trace( "*" );
        }
   al_trace( "\n%s\np. mean = %9f, p. var = %9f\ns. mean = %9f, s. var = %9f\n" , 
                rx->Name() , rx->Mean() , rx->Variance() , mean , sd * sd );

        delete [] count;
}
#endif

---