CDbw.h

#ifndef CDbw_H
#define CDbw_H

#include "Som.h"
#include "defs.h"

class CDbw {

public:
        CDbw( SOM * som, TIntVector cluster, int num_clusters ) : som(som), v(cluster) , c(num_clusters)
        {
                d = som->getData()->getDimension();
                m = som->getMapcode()->getMapSize();
                n = som->getData()->getDataSize();
                entries                 = som->getData()->getEntries();
                codebook                = som->getMapcode()->getCodebook();
                entries_label   = som->getData()->getId(); // getLabel();
                codebook_label  = som->getMapcode()->getLabel();
                Init();
        };
        Value_Type Index();

private:
        // Init variables
        void Init();
        void getStdev_v();
        TMatrix getClose_rep(int c_i, int c_j);
        Value_Type getDensity( TMatrix close_rep, int c_i, int c_j );
        Value_Type Inter_dens();
        Value_Type Sep();
        Value_Type getDensity_v( int c_i );
        Value_Type Intra_dens();
        
        // Overal average standard deviation
        Value_Type stdev();

        TIntVector v;
        TMatrix stdev_v;
        TVector stdev_vv;
        TMatrix mean_v;
        TVector num_v;  
        int c;
        int d;
        int m;
        int n;
        TMatrix entries;
        TMatrix codebook;
        TvecLabel entries_label;
        TLabel  codebook_label;
        SOM * som;
};

void 
CDbw::Init() 
{
        getStdev_v();
};

void 
CDbw::getStdev_v() 
{
        int i, j, l, k;
        stdev_v = create_matrix(1, c, 1, d ); 
        mean_v = create_matrix(1, c, 1, d ); 
        stdev_vv = create_vector( 1, c );
        
        for (i=1; i<=c; i++)  {
                stdev_vv[i] = 0;
                for (j=1; j<=d; j++) { 
                        stdev_v[i][j] = 0; 
                        mean_v[i][j] = 0; 
                };
        };

        num_v = create_vector(1, c ); for (i=1; i<=c; i++) num_v[i] = 0;

        for (i=1; i<=m; i++ ) {
                if ( v[i] > 0 ) {
                        for (l=0; l< codebook_label[i-1].size(); l++ )
                                for (j=1; j<=n; j++) {
                                        if ( codebook_label[i-1][l] == entries_label[j-1] ) {
                                                for (k=1; k<=d; k++ ) {
                                                        mean_v[v[i]][k] += entries[j][k];
                                                }
                                                ++num_v[v[i]];
                                                break;
                                        }
                        }
                }
        }
        for (i=1; i<=c; i++) for (j=1; j<=d; j++)
                mean_v[i][j] /= num_v[i];
        for (i=1; i<=m; i++ ) {
                if ( v[i] > 0 ) {
                        for (l=0; l< codebook_label[i-1].size(); l++ )
                                for (j=1; j<=n; j++) {
                                        if ( codebook_label[i-1][l] == entries_label[j-1] ) {
                                                for (k=1; k<=d; k++ ) {
                                                        stdev_v[v[i]][k] += pow((entries[j][k] - mean_v[v[i]][k]),2 );
                                                        stdev_vv[v[i]] += pow((entries[j][k] - mean_v[v[i]][k]),2 );
                                                }                                               
                                                break;
                                        }
                        }
                }
        }
        for (i=1; i<=c; i++) {
                if ( (num_v[i] - 1) > 0 )
                        stdev_vv[i] = sqrt( stdev_vv[i]/(num_v[i] - 1) );
                else
                        stdev_vv[i] = 0;
                for (j=1; j<=d; j++) {
                        if ((num_v[i] - 1) > 0 )
                                stdev_v[i][j] = sqrt( stdev_v[i][j]/(num_v[i] - 1) );
                        else
                                stdev_v[i][j] = 0;
                }
        }
};

Value_Type 
CDbw::stdev() {
        Value_Type sum = 0;
        int i;
        for (i=1; i<= c; i++ ) 
                sum += stdev_vv[i]*stdev_vv[i]; 
        return sqrt( sum / c );
};

TMatrix 
CDbw::getClose_rep(int c_i, int c_j) 
{
        int i, j;
        Value_Type tmp, min_dist = 999;
        TMatrix close_rep = create_matrix( 1, 2, 1, d );
        for (i=1; i<=2; i++) for (j=1; j<=d; j++) close_rep[i][j] = 0;

        for (i=1; i<=m; i++)
                if ( v[i] == c_i ) {
                        for (j=1; j<=m; j++)
                                if ( v[j] == c_j ) {
                                        tmp = dist( codebook[i], codebook[j], d );
                                        if ( tmp < min_dist ) {
                                                min_dist = tmp;
                                                close_rep[1] = codebook[i];
                                                close_rep[2] = codebook[j];
                                        }
                                }
                }
        return close_rep;
};

Value_Type 
CDbw::getDensity( TMatrix close_rep, int c_i, int c_j ) 
{
        Value_Type tmp, min_dist, sum = 0;
        TVector u = create_vector( 1, d );
        Value_Type stdev1 = 0, stdev2 = 0;
        int i, l, j;

        for (i=1; i<=d; i++ ) {
                u[i] = (close_rep[1][i] + close_rep[2][i])/2;
                //stdev1 += stdev_v[ c_i ][i]*stdev_v[ c_i ][i];
                //stdev2 += stdev_v[ c_j ][i]*stdev_v[ c_j ][i];
        }
        
        min_dist = ( stdev_vv[c_i] + stdev_vv[ c_j ] ) / 2;

        for (i=1; i<= m; i++ ) {
                if ((v[i] == c_i) || (v[i] == c_j)) {
                        for (l=0; l< codebook_label[i-1].size(); l++ ) {
                                for ( j=1; j<= n; j++) {
                                        if ( codebook_label[i-1][l] == entries_label[j-1] ) {
                                                tmp = dist( entries[j], u , d );
                                                if ( tmp <= min_dist )
                                                        ++sum;
                                                break;
                                        }
                                }
                        }
                }
        }
        return sum/( num_v[c_i] + num_v[c_j] );
};

Value_Type 
CDbw::Inter_dens() 
{
        int i, j;
        TMatrix close_rep;
        Value_Type sum = 0;
        Value_Type stdev1 = 0, stdev2 = 0;

        for (i=1; i<= c; i++)
                for (j=1; j<= c; j++ )
                        if ( i != j ) {
                                close_rep = getClose_rep(i, j);                                         
                                if ( ( stdev_vv[i] + stdev_vv[j] ) > 0 )
                                        sum += ( dist( close_rep[1], close_rep[2], d ) * getDensity( close_rep, i, j ) )/( stdev_vv[i] + stdev_vv[j] );                         
                        }
        return  sum;
};

Value_Type 
CDbw::Sep() 
{
        int i, j;
        TMatrix close_rep;
        Value_Type  sum = 0;

        for (i=1; i<= c; i++)           
                for (j=1; j<= c; j++ )
                        if ( i != j ) {
                                close_rep = getClose_rep(i, j);
                                sum += dist( close_rep[1], close_rep[2], d );
                        }
        
        return sum / (1 + Inter_dens() ) ;
};

Value_Type 
CDbw::getDensity_v( int c_i )
{
        int i, j, l, k, r_j;
        Value_Type tmp, sum = 0;
        Value_Type stdev_value = stdev();
        int count = 0;

        for (k=1; k<= m; k++ ) {
                if ( v[k] == c_i ) {
                        r_j = k; ++count;
                        for (i=1; i<=m; i++ )
                                if ( v[i] == c_i ) 
                                        for (l=0; l < codebook_label[i-1].size(); l++ )
                                                for (j=1; j<=n; j++ )
                                                        if ( codebook_label[i-1][l] == entries_label[j-1] ) {
                                                                tmp = dist( entries[j], codebook[r_j], d );
                                                                if ( tmp <= stdev_value ) 
                                                                        ++sum;
                                                                break;
                                                        }
                }
        }
        return ( sum / count );
};


Value_Type 
CDbw::Intra_dens() {
        int i;
        Value_Type sum = 0;
        for (i=1; i<= c; i++ ) {
                sum += getDensity_v( i );
                //cout << "Density_v(" << i << ") = " << getDensity_v( i ) << endl;
        }
        return ( sum / c );
};


Value_Type 
CDbw::Index() 
{
        Value_Type intra = Intra_dens();
        Value_Type sep   = Sep();
        //cout << "Inter_dens = " << Inter_dens() << endl;
        //cout << "Intra_dens = " << intra << endl;
        //cout << "Sep = " << sep << endl;
        //cout << "Stdev = " << stdev();
        //printmatrix( stdev_v, 1, c, 1, d );
        
        return intra*sep;
};

#endif

---