GaussianMixture.h

Go to the documentation of this file.

// ----------------------------------------------------------------------------
//
// $Id$
//
// Copyright 2008, 2009, 2010, 2011, 2012  Antonio Franchi and Paolo Stegagno    
//
// This file is part of MIP.
//
// MIP is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// MIP 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 General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with MIP. If not, see <http://www.gnu.org/licenses/>.
//
// Contact info: antonio.franchi@tuebingen.mpg.de stegagno@diag.uniroma1.it
//
// ----------------------------------------------------------------------------


#ifndef __PDF_H
#define __PDF_H

#include "Gaussian.h"
#include "Discrete.h"

using namespace MipBaselib;
using namespace arma;

namespace MipBaselib{
 
 /* @{ */
 
 
 class GaussianComponent:public Gaussian{
  public:
   Decimal weight;
   
   GaussianComponent();
   
   GaussianComponent(const GaussianComponent& a);
   
   GaussianComponent(Gaussian g, Decimal w);
   
   GaussianComponent(DCol m, DMat &c, Decimal w=0.0);
   
   GaussianComponent& operator=(const GaussianComponent& a);
   
   virtual string print();
 };
 
 
 class GaussianMixture:public vector<GaussianComponent>{
  private:
   
   bool _expectedValueComputed;
   bool _covarianceComputed;
   bool _normalized;
   
   DCol _expectedValue;
   DMat _covariance;
   
   void _computeExpectedValue();
   
   void _computeCovariance();
   
  public:
   
   bool normalized();
   
   void normalize();
   
   GaussianMixture();
   
   DCol expectedValue();
   
   Decimal covariance(int row, int col);
   
   Decimal evaluate(DCol p);
   
   string print();
 };
 
 
 
 
 class GaussianComponentID:public GaussianComponent{
  public:
   int id;
   
   GaussianComponentID();
   
   GaussianComponentID(const GaussianComponentID& a);
   
   GaussianComponentID(Gaussian &g, Decimal w, int ide);
   
   GaussianComponentID(DCol &m, DMat &c, Decimal w=0.0, int ide=0);
   
   GaussianComponentID(stringstream &ss);
   
   GaussianComponentID& operator=(const GaussianComponentID& a);
   
   virtual string print() const;
   
   string exportAsString() const;
 };
 
 
 class GaussianMixtureID:public vector<GaussianComponentID>{
  private:
   
   bool _expectedValueComputed;
   bool _covarianceComputed;
   bool _normalized;
   
   DCol _expectedValue;
   DMat _covariance;
   
   void _computeExpectedValue();
   
   void _computeCovariance();
   
  public:
   
   bool normalized();
   
   void normalize();
   
   GaussianMixtureID();
   
   GaussianMixtureID(stringstream& ss);
   
   GaussianMixtureID(string& s);
   
   DCol expectedValue();
   
   Decimal covariance(int row, int col);
   
   Decimal evaluate(DCol p);
   
   string print() const;
   
   string exportAsString() const;
 };
 
 
 
 
 
 class GaussianComponentJointDiscrete: public GaussianComponent, public Discrete{
  
  public:
  inline GaussianComponent* _thisGC(){ return (GaussianComponent*)this; }
  inline Discrete* _thisD(){ return (Discrete*)this; }
  
  public:
   
   GaussianComponentJointDiscrete();
   
   GaussianComponentJointDiscrete(GaussianComponent gc, Discrete d);
   
   GaussianComponentJointDiscrete(const GaussianComponentJointDiscrete &gcjd);
   
   GaussianComponentJointDiscrete(stringstream &ss);
   
   GaussianComponentJointDiscrete& operator=(const GaussianComponentJointDiscrete &gcjd);
   
   string exportAsString();
 };
 
 
 
 class GaussianMixtureJointId:public vector<GaussianComponentJointDiscrete>{
  private:
   
   bool _expectedValueComputed;
   bool _covarianceComputed;
   bool _normalized;
   
   DCol _expectedValue;
   DMat _covariance;
   
   void _computeExpectedValue();
   
   void _computeCovariance();
   
  public:
   
   
   bool normalized();
   
   void normalize();
   
   GaussianMixtureJointId();
   
   GaussianMixtureJointId(stringstream &ss);
   
   GaussianMixtureJointId(string &s);
   
   DCol expectedValue();
   
   Decimal covariance(int row, int col);
   
   Decimal evaluate(DCol p);
   
   string print();
   
   string exportAsString();
 };
 
 /* @} */
}


#endif