STK::SemiSEMAlgo Class Reference

Implementation of the SemiSEM algorithm. More...

#include <STK_MixtureAlgo.h>

Inheritance diagram for STK::SemiSEMAlgo:

Public Member Functions
	SemiSEMAlgo ()
	default constructor
	SemiSEMAlgo (SemiSEMAlgo const &algo)
	Copy constructor.
virtual	~SemiSEMAlgo ()
	destructor
virtual SemiSEMAlgo *	clone () const
	clone pattern
virtual bool	run ()
	run the algorithm on the model calling sStep, mStep and eStep of the model until the maximal number of iteration is reached.
Public Member Functions inherited from STK::IMixtureAlgo
virtual	~IMixtureAlgo ()
	destructor
int	nbIterMax () const
int	epsilon () const
Real	threshold () const
void	setModel (IMixtureComposer *p_model)
	set model
void	setNbIterMax (int nbIterMax)
	set maximal number of iterations
void	setEpsilon (Real epsilon)
	set tolerance value
void	setThreshold (Real threshold)
	set threshold value
Public Member Functions inherited from STK::IRunnerBase
String const &	error () const
	get the last error message.

Additional Inherited Members
Protected Member Functions inherited from STK::IMixtureAlgo
	IMixtureAlgo ()
	default constructor
	IMixtureAlgo (IMixtureAlgo const &algo)
	Copy constructor.
Protected Member Functions inherited from STK::IRunnerBase
	IRunnerBase ()
	default constructor
	IRunnerBase (IRunnerBase const &runner)
	copy constructor
virtual	~IRunnerBase ()
	destructor
virtual void	update ()
	update the runner.
Protected Attributes inherited from STK::IMixtureAlgo
IMixtureComposer *	p_model_
	pointer on the mixture model
int	nbIterMax_
	number of iterations of the algorithm
Real	epsilon_
	tolerance of the algorithm.
Real	threshold_
	Minimal number of individuals.
Protected Attributes inherited from STK::IRunnerBase
String	msg_error_
	String with the last error message.
bool	hasRun_
	true if run has been used, false otherwise

Detailed Description

Implementation of the SemiSEM algorithm.

The SemiSEM algorithm calls alternatively the steps:

• samplingStep()
• pStep()
• paramUpdateStep()
• eStep()
• storeIntermediateResults(iter) until the maximum number of iterations is reached.

Definition at line 140 of file STK_MixtureAlgo.h.

Constructor & Destructor Documentation

SemiSEMAlgo() [1/2]

STK::SemiSEMAlgo::SemiSEMAlgo ( )

inline

default constructor

Definition at line 144 of file STK_MixtureAlgo.h.

: IMixtureAlgo() {}

Referenced by clone().

SemiSEMAlgo() [2/2]

STK::SemiSEMAlgo::SemiSEMAlgo ( SemiSEMAlgo const &  algo )

inline

Copy constructor.

Parameters

algo	the algorithm to copy

Definition at line 147 of file STK_MixtureAlgo.h.

: IMixtureAlgo(algo) {}

~SemiSEMAlgo()

virtual STK::SemiSEMAlgo::~SemiSEMAlgo ( )

inlinevirtual

destructor

Definition at line 149 of file STK_MixtureAlgo.h.

{}

Member Function Documentation

clone()

virtual SemiSEMAlgo * STK::SemiSEMAlgo::clone ( ) const

inlinevirtual

clone pattern

Definition at line 151 of file STK_MixtureAlgo.h.

{ return new SemiSEMAlgo(*this);}

References SemiSEMAlgo().

run()

bool STK::SemiSEMAlgo::run ( )

virtual

run the algorithm on the model calling sStep, mStep and eStep of the model until the maximal number of iteration is reached.

Returns

true if no error occur, false otherwise.

Implements STK::IRunnerBase.

Definition at line 231 of file STK_MixtureAlgo.cpp.

{
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("---------------------------------\n");
  stk_cout << _T("Entering SemiSEMAlgo::run() with:\n")
           << _T("nbIterMax_ = ") << nbIterMax_ << _T("\n")
           << _T("p_model_->lnLikelihood = ") << p_model_->lnLikelihood() << _T("\n");
#endif
  bool result = true;
  try
  {
    int iter;
    for (iter = 0; iter < nbIterMax_; ++iter)
    {
      p_model_->samplingStep();
      p_model_->pStep();
      p_model_->paramUpdateStep();
      Real nb = p_model_->eStep();
      if (nb<threshold_)
      {
        msg_error_ = STKERROR_1ARG(SemiSEMAlgo::run,nb,Not enough individuals after eStep\n);
#ifdef STK_MIXTURE_VERBOSE
  stk_cout << _T("An exception occur in SemiSEMAlgo::run(): ") << msg_error_ << _T("\n");
#endif
        result = false;
        break;
      }
      p_model_->storeIntermediateResults(iter+1); // store current parameters
    }
#ifdef STK_MIXTURE_VERBOSE
    stk_cout << _T("In SemiSEMAlgo::run() iteration ") << iter << _T("terminated.\n")
             << _T("p_model_->lnLikelihood = ") << p_model_->lnLikelihood() << _T("\n");
#endif
  }
  catch (Clust::exceptions const& error)
  {
    msg_error_ = Clust::exceptionToString(error);
#ifdef STK_MIXTURE_VERBOSE
  stk_cout << _T("An exception occur in SemiSEMAlgo::run(): ") << msg_error_ << _T("\n");
#endif
    result = false;
  }
  // set averaged parameters
  if (result)
  {
    p_model_->setParametersStep();
#ifdef STK_MIXTURE_VERY_VERBOSE
    stk_cout << _T("\nIn SemiSEMAlgo::run(), setParameters done.\n")
             << _T("p_model_->lnLikelihood = ") << p_model_->lnLikelihood() << _T("\n");
#endif
  }
  else // no result
  { p_model_->releaseIntermediateResults();}
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("Terminating SemiSEMAlgo::run()\n");
  stk_cout << _T("------------------------------\n");
#endif
  return result;
}

References _T, STK::IRunnerBase::error(), STK::IMixtureComposer::eStep(), STK::Clust::exceptionToString(), STK::IStatModelBase::lnLikelihood(), STK::IRunnerBase::msg_error_, STK::IMixtureAlgo::nbIterMax_, STK::IMixtureAlgo::p_model_, STK::IMixtureComposer::paramUpdateStep(), STK::IMixtureComposer::pStep(), STK::IMixtureStatModel::releaseIntermediateResults(), run(), STK::IMixtureStatModel::samplingStep(), STK::IMixtureStatModel::setParametersStep(), stk_cout, STKERROR_1ARG, STK::IMixtureStatModel::storeIntermediateResults(), and STK::IMixtureAlgo::threshold_.

Referenced by run().

---

The documentation for this class was generated from the following files:

• STK_MixtureAlgo.h
• STK_MixtureAlgo.cpp