STK_HDMatrixGaussianParameters.h

Go to the documentation of this file.

/*--------------------------------------------------------------------*/
/*     Copyright (C) 2004-2016 Serge Iovleff
 
    This program is free software; you can redistribute it and/or modify
    it under the terms of the GNU Lesser General Public License as
    published by the Free Software Foundation; either version 2 of the
    License, or (at your option) any later version.
 
    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU Lesser General Public License for more details.
 
    You should have received a copy of the GNU Lesser General Public
    License along with this program; if not, write to the
    Free Software Foundation, Inc.,
    59 Temple Place,
    Suite 330,
    Boston, MA 02111-1307
    USA
 
    Contact : S..._DOT_I..._AT_stkpp.org (see copyright for ...)
*/
 
/*
 * Project:  stkpp::Clustering
 * created on: Oct 24, 2013
 * Author:   Serge Iovleff
 **/
 
#ifndef STK_HDMATRIXGAUSSIANPARAMETERS_H
#define STK_HDMATRIXGAUSSIANPARAMETERS_H
 
#include <Arrays/include/STK_Array1D.h>
#include <Arrays/include/STK_CArray.h>
#include <Arrays/include/STK_CArraySquare.h>
#include <Arrays/include/STK_CArrayPoint.h>
#include <Arrays/include/STK_CArrayVector.h>
 
#include <STatistiK/include/STK_Stat_Online.h>
 
namespace STK
{
 
template<class Array_>
struct HDMatrixModelParameters
{
    Array1D<Array_> meank_;
 
    Array1D<CSquareX> rowQk_;
    Array1D<CVectorX> rowAjk_;
    Array1D<Real> rowBk_;
 
    CVectorXi rowDk_;
 
    Array1D<CSquareX> colQk_;
    Array1D<CVectorX> colAjk_;
    Array1D<Real> colBk_;
 
    CVectorXi colDk_;
 
    Array1D< Stat::Online<Array_, Real> > statMeank_;
 
    Array1D< Stat::Online<CSquareX, Real> > statRowQk_;
    Array1D< Stat::Online<CVectorX, Real> > statRowAjk_;
    Array1D< Stat::Online<Real, Real> > statRowBk_;
 
    Array1D< Stat::Online<CSquareX, Real> > statColQk_;
    Array1D< Stat::Online<CVectorX, Real> > statColAjk_;
    Array1D< Stat::Online<Real, Real> > statColBk_;
 
    HDMatrixModelParameters( int nbCluster)
                           : meank_(nbCluster)
                           , rowQk_(nbCluster), rowAjk_(nbCluster)
                           , rowBk_(nbCluster), rowDk_(nbCluster)
                           , colQk_(nbCluster), colAjk_(nbCluster)
                           , colBk_(nbCluster), colDk_(nbCluster)
                           , statMeank_(nbCluster)
                           , statRowQk_(nbCluster), statRowAjk_(nbCluster)
                           , statRowBk_(nbCluster)
                           , statColQk_(nbCluster), statColAjk_(nbCluster)
                           , statColBk_(nbCluster)
    {}
    HDMatrixModelParameters( HDMatrixModelParameters const& param)
                           : meank_(param.meank_)
                           , rowQk_(param.rowQk_), rowAjk_(param.rowAjk_)
                           , rowBk_(param.rowBk_), rowDk_(param.rowDk_)
                           , colQk_(param.colQk_), colAjk_(param.colAjk_)
                           , colBk_(param.colBk_), colDk_(param.colDk_)
                           , statMeank_(param.statMeank_)
                           , statRowQk_(param.statRowQk_), statRowAjk_(param.statRowAjk_)
                           , statRowBk_(param.statRowBk_)
                           , statColQk_(param.statColQk_), statColAjk_(param.statColAjk_)
                           , statColBk_(param.statColBk_)
    {}
    ~HDMatrixModelParameters() {}
    HDMatrixModelParameters& operator=( HDMatrixModelParameters const& param)
    {
      meank_ = param.meank_;
      rowQk_ = param.rowQk_; rowAjk_ = param.rowAjk_;
      rowBk_ = param.rowBk_; rowDk_  = param.rowDk_;
      colQk_ = param.colQk_; colAjk_ = param.colAjk_;
      colBk_ = param.colBk_; colDk_  = param.colDk_;
      statMeank_ = param.statMeank_;
      statRowQk_ = param.statRowQk_; statRowAjk_ = param.statRowAjk_;
      statRowBk_ = param.statRowBk_;
      statColQk_ = param.statColQk_; statColAjk_ = param.statColAjk_;
      statColBk_ = param.statColBk_;
    }
 
    void resize(Range const& rangeRows, Range const& rangeCols);
 
    void updateStatistics()
    {
      for(int k=statMeank_.begin(); k<statMeank_.end(); ++k)
      {
        statMeank_[k].update(meank_[k]);
        statRowQk_[k].update(rowQk_[k]);
        statRowAjk_[k].update(rowAjk_[k]);
        statRowBk_[k].update(rowBk_[k]);
        statColQk_[k].update(colQk_[k]);
        statColAjk_[k].update(colAjk_[k]);
        statColBk_[k].update(colBk_[k]);
      }
    }
    void setStatistics()
    {
      for(int k=statMeank_.begin(); k<statMeank_.end(); ++k)
      {
        meank_[k]  = statMeank_[k].mean();
        rowQk_[k]  = statRowQk_[k].mean();
        rowAjk_[k] = statRowAjk_[k].mean();
        rowBk_[k]  = statRowBk_[k].mean();
        colQk_[k]  = statColQk_[k].mean();
        colAjk_[k] = statColAjk_[k].mean();
        colBk_[k]  = statColBk_[k].mean();
      }
      releaseStatistics();
    }
    void releaseStatistics()
    {
      for(int k=statMeank_.begin(); k<statMeank_.end(); ++k)
      {
        statMeank_[k].release();
        statRowQk_[k].release();
        statRowAjk_[k].release();
        statRowBk_[k].release();
        statColQk_[k].release();
        statColAjk_[k].release();
        statColBk_[k].release();
      }
    }
    template<class Array>
    void setParameters( ExprBase<Array> const& params)
    {}
//    {
//      for(int k=mean_.begin(), kp= params.beginRows(); k<mean_.end(); ++k, kp+=2)
//      {
//        mean_[k] = params.row(kp);
//        sigma_[k] = params.row(kp+1);
//      }
//    }
};
 
} // namespace STK
 
#endif /* STK_HDGAUSSIANPARAMETERS_H */