STK::IMixtureComposer Class Reference

Base class for Mixture (composed) model. More...

#include <STK_IMixtureComposer.h>

Inheritance diagram for STK::IMixtureComposer:

Public Member Functions
virtual	~IMixtureComposer ()
	destructor
Clust::modelState	state () const
void	setState (Clust::modelState state)
	set the state of the model : should be used by any strategy
virtual IMixtureComposer *	create () const =0
	create pattern
virtual IMixtureComposer *	clone () const =0
	clone pattern
virtual void	randomInit ()=0
	initialize randomly the parameters of the components of the model
virtual void	paramUpdateStep ()=0
	Compute the proportions and the model parameters given the current tik mixture parameters.
virtual int	cStep ()
	Replace tik by zik.
virtual int	sStep ()
	Simulate zi accordingly to tik and replace tik by zik by calling cStep().
virtual Real	eStep ()
	compute the zi, the lnLikelihood of the current estimates and the next value of the tik.
virtual void	mapStep ()
	Compute zi using the Map estimate.
virtual void	finalizeStep ()
	Finalize the estimation of the model.
virtual void	pStep ()
	Compute proportions using the ML estimates, default implementation.
virtual void	initializeStep ()
	Initialize the model before its first use.
void	randomClassInit ()
	Initialize randomly the labels zi of the model.
void	randomFuzzyInit ()
	Initialize randomly the posterior probabilities tik of the model, then compute the zi values with mapStep, compute the initial parameter values using paramUpdateStep, and compute the tik using eStep().
Public Member Functions inherited from STK::IMixtureStatModel
virtual	~IMixtureStatModel ()
	destructor
int	nbCluster () const
CPointX const &	pk () const
CArrayXX const &	tik () const
CPointX const &	tk () const
CVectorXi const &	zi () const
std::vector< IMixture * > const &	v_mixtures () const
Real	computeLnLikelihood (int i) const
Real	computeLikelihood (int i) const
Real	computeLnLikelihood () const
Real	computeICL () const
IMixture *	getMixture (String const &idData) const
	Utility lookup function allowing to find a Mixture from its idData.
void	registerMixture (IMixture *p_mixture)
	register a mixture to the composer.
void	releaseMixture (String const &idData)
	release a mixture from the composer.
int	computeNbFreeParameters () const
	compute the number of free parameters of the model.
int	computeNbMissingValues () const
	compute the missing values of the model.
virtual Real	lnComponentProbability (int i, int k) const =0
virtual void	imputationStep ()
	Impute the missing values.
virtual void	samplingStep ()
	Simulation of all the latent variables and/or missing data excluding class labels.
virtual void	setParametersStep ()
	Utility method allowing to signal to a mixture to set its parameters.
virtual void	storeIntermediateResults (int iteration)
	This step can be used to signal to the mixtures that they must store results.
virtual void	releaseIntermediateResults ()
	This step can be used to signal to the mixtures that they must release the stored results.
virtual void	writeParameters (ostream &os) const
	write the parameters of the model in the stream os.
template<class Array >
void	setMixtureParameters (Array const &tik)
	set the mixture parameters using an array of posterior probabilities.
template<class Array , class RowVector >
void	setMixtureParameters (Array const &tik, RowVector const &pk)
	set the mixture parameters giving the posterior probabilities and the proportions.
template<class RowVector >
void	setProportions (RowVector const &pk)
	Set proportions of each classes.
template<class Manager , class Parameters >
void	setParameters (IMixtureManager< Manager > const &manager, String const &idData, Parameters const &param)
	Utility method allowing to set the parameters to a specific mixture.
template<class Manager >
void	createMixture (IMixtureManager< Manager > &manager)
	Utility method allowing to create all the mixtures registered in the data handler of a mixture manager and to register them.
template<class Manager >
IMixture *	createMixture (IMixtureManager< Manager > &manager, String const &idData)
	Utility method allowing to create a mixture with a given data set and register it.
template<class Manager >
void	removeMixture (IMixtureManager< Manager > &manager, String const &idData)
	Utility method allowing to release completely a mixture with its data set.
template<class Manager , class Parameters >
void	getParameters (IMixtureManager< Manager > const &manager, String const &idData, Parameters &param) const
	Utility method allowing to get the parameters of a specific mixture.
template<class Manager , class MissingValues >
void	getMissingValues (IMixtureManager< Manager > const &manager, String const &idData, MissingValues &missing) const
	Utility method allowing to get the missing values of a specific mixture.
template<class DataHandler >
void	createMixture (IMixtureManager< DataHandler > &manager)
Public Member Functions inherited from STK::IStatModelBase
int	nbSample () const
Real	lnNbSample () const
int	nbVariable () const
Real	lnLikelihood () const
Real	likelihood () const
int	nbFreeParameter () const
Real	computeBIC () const
Real	computeAIC () const
Real	computeML () const

Protected Member Functions
	IMixtureComposer (int nbSample, int nbCluster)
	Constructor.
	IMixtureComposer (IMixtureComposer const &model)
	copy constructor.
void	sStep (int i)
	Simulate zi accordingly to tik.
void	cStep (int i)
	Replace tik by zik.
Real	eStep (int i)
	compute one zi and the next value of the tik for i fixed
void	mapStep (int i)
	Compute zi using the Map estimate for i fixed.
virtual void	initializeMixtureParameters ()
	Create the mixture model parameters pk_ and tik_.
virtual int	randomTik ()
	generate random tik_
virtual int	randomZi ()
	generate random zi_
Protected Member Functions inherited from STK::IMixtureStatModel
	IMixtureStatModel (int nbSample, int nbCluster)
	Constructor.
	IMixtureStatModel (IMixtureStatModel const &model)
	copy constructor.
void	setNbCluster (int nbCluster)
	set the number of cluster of the model
Protected Member Functions inherited from STK::IStatModelBase
	IStatModelBase ()
	Default constructor.
	IStatModelBase (int nbSample)
	Constructor with specified dimension.
	IStatModelBase (int nbSample, int nbVariable)
	Constructor with specified dimension.
	IStatModelBase (IStatModelBase const &model)
	Copy constructor.
	~IStatModelBase ()
	destructor
void	setNbFreeParameter (int const &nbFreeParameter)
	set the number of free parameters of the model
void	setNbSample (int const &nbSample)
	set the number of samples of the model
void	setNbVariable (int const &nbVariable)
	set the number of variables of the model
void	setLnLikelihood (Real const &lnLikelihood)
	set the log-likelihood of the model
void	initialize (int nbSample, int nbVariable)
	set the dimensions of the parameters of the model

Private Attributes
Clust::modelState	state_
	state of the model
CPointX	lnComp_
	Auxiliary array used in the eStep.

Additional Inherited Members
Public Types inherited from STK::IMixtureStatModel
typedef std::vector< IMixture * >::const_iterator	ConstMixtIterator
typedef std::vector< IMixture * >::iterator	MixtIterator
Protected Attributes inherited from STK::IMixtureStatModel
int	nbCluster_
	number of cluster.
CPointX	pk_
	The proportions of each mixtures.
CArrayXX	tik_
	The tik probabilities.
CPointX	tk_
	The sum of the columns of tik_.
CVectorXi	zi_
	The zi class label.
std::vector< IMixture * >	v_mixtures_
	vector of pointers to the mixtures components

Detailed Description

Base class for Mixture (composed) model.

In this interface we assume there is an underline generative model that will be estimated using either an EM, SEM or CEM algorithm. All mixture parameters are created using the method

void initializeParameters();

in the constructor. They can be accessed from the mixtures using constant accessors.

The pure virtual function to implement in derived class are

virtual IMixtureComposer* create() const = 0;
virtual IMixtureComposer* clone() const = 0;
virtual void randomInit() = 0;
virtual void paramUpdateStep() = 0;
virtual int computeNbFreeParameters() const = 0;

The public virtual function that can be re-implemented in derived class for a specific behavior are:

virtual void initializeStep();
virtual void randomClassInit();
virtual void randomFuzzyInit();
virtual int cStep();
virtual int sStep();
virtual Real eStep();
virtual void mapStep();
virtual void finalizeStep();
virtual void pStep();

See also

IMixture

Note

the virtual method IMixtureComposer::initializeStep is called in all the initialization method. Don't forget to called it in the randomInit implementation.

Definition at line 85 of file STK_IMixtureComposer.h.

Constructor & Destructor Documentation

IMixtureComposer() [1/2]

STK::IMixtureComposer::IMixtureComposer ( int  nbSample, int  nbCluster  )

protected

Constructor.

Parameters

nbCluster,nbSample

number of clusters and samples

Definition at line 49 of file STK_IMixtureComposer.cpp.

                                  : IMixtureStatModel(nbSample, nbCluster)
                                  , state_(Clust::modelCreated_)
#ifndef _OPENMP
                                  , lnComp_(nbCluster)
#endif
{}

IMixtureComposer() [2/2]

STK::IMixtureComposer::IMixtureComposer ( IMixtureComposer const &  model )

protected

copy constructor.

Parameters

model

the model to clone

Definition at line 58 of file STK_IMixtureComposer.cpp.

                                  : IMixtureStatModel(model)
                                  , state_(model.state_)
#ifndef _OPENMP
                                  , lnComp_(model.lnComp_)
#endif
{}

~IMixtureComposer()

STK::IMixtureComposer::~IMixtureComposer ( )

virtual

destructor

Definition at line 66 of file STK_IMixtureComposer.cpp.

{}

Member Function Documentation

clone()

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

pure virtual

clone pattern

Implements STK::IMixtureStatModel.

Implemented in STK::MixtureComposer, STK::MixtureComposerFixedProp, STK::MixtureSemiLearner, and STK::MixtureSemiLearnerFixedProp.

Referenced by STK::FullStrategy::initStep(), and STK::FullStrategy::run().

create()

virtual IMixtureComposer * STK::IMixtureComposer::create ( ) const

pure virtual

create pattern

Implements STK::IMixtureStatModel.

Implemented in STK::MixtureComposer, STK::MixtureComposerFixedProp, STK::MixtureSemiLearner, and STK::MixtureSemiLearnerFixedProp.

Referenced by STK::FullStrategy::initStep(), STK::SimpleStrategy::run(), and STK::XemStrategy::run().

cStep() [1/2]

int STK::IMixtureComposer::cStep ( )

virtual

Replace tik by zik.

Returns

the minimal value of individuals in a class

Reimplemented in STK::MixtureSemiLearner.

Definition at line 149 of file STK_IMixtureComposer.cpp.

{
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("Entering IMixtureComposer::cStep()\n");
#endif
  for (int i=tik_.beginRows(); i < tik_.endRows(); i++)
  { cStep(i);}
  // count the number of individuals in each class
  tk_= Stat::sum(tik_);
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("IMixtureComposer::cStep() tk_ = ") << tk_;
  stk_cout << _T("IMixtureComposer::cStep() done\n");
#endif
  return tk_.minElt();
}

References _T, cStep(), STK::ExprBase< Derived >::minElt(), stk_cout, STK::Stat::sum(), STK::IMixtureStatModel::tik_, and STK::IMixtureStatModel::tk_.

Referenced by cStep(), STK::MixtureSemiLearner::cStep(), randomZi(), STK::CEMAlgo::run(), and sStep().

cStep() [2/2]

void STK::IMixtureComposer::cStep ( int  i )

protected

Replace tik by zik.

Parameters

i	index of the the individual

Definition at line 219 of file STK_IMixtureComposer.cpp.

{ tik_.row(i) = 0.; tik_.elt(i, zi_[i]) = 1.;}

References STK::ICArray< Derived >::row(), STK::IMixtureStatModel::tik_, and STK::IMixtureStatModel::zi_.

eStep() [1/2]

Real STK::IMixtureComposer::eStep ( )

virtual

compute the zi, the lnLikelihood of the current estimates and the next value of the tik.

Returns

the minimal value of tk

Reimplemented in STK::MixtureSemiLearner.

Definition at line 180 of file STK_IMixtureComposer.cpp.

{
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("Entering IMixtureComposer::eStep()\n");
#endif
  Real sum = 0.;
  int i;
#ifdef _OPENMP
#pragma omp parallel for reduction (+:sum)
#endif
  for (i = tik_.beginRows(); i < tik_.endRows(); ++i) { sum += eStep(i);}
  // update ln-likelihood
  setLnLikelihood(sum);
  // compute proportions
  tk_ = Stat::sumByCol(tik_);
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("IMixtureComposer::eStep() lnLikelihood =") << sum << _T("\n");
  stk_cout << _T("IMixtureComposer::eStep() done\n");
  if (!isFinite(sum))
  { stk_cout << _T("tik_ =\n") << tik_.transpose();}
#endif
  return tk_.minElt();
}

References _T, eStep(), STK::isFinite(), STK::ExprBase< Derived >::minElt(), STK::IStatModelBase::setLnLikelihood(), stk_cout, STK::sum(), STK::Stat::sumByCol(), STK::IMixtureStatModel::tik_, STK::IMixtureStatModel::tk_, and STK::ArrayBase< Derived >::transpose().

Referenced by STK::IMixtureAlgoPredict::burnStep(), eStep(), STK::MixtureSemiLearner::eStep(), STK::IMixtureAlgoPredict::predictBayesClassifier(), randomClassInit(), randomFuzzyInit(), STK::MixtureComposer::randomInit(), STK::EMAlgo::run(), STK::CEMAlgo::run(), STK::SEMAlgo::run(), STK::SemiSEMAlgo::run(), STK::EMPredict::run(), and STK::SemiSEMPredict::run().

eStep() [2/2]

Real STK::IMixtureComposer::eStep ( int  i )

protected

compute one zi and the next value of the tik for i fixed

Parameters

i	the individual

Returns

the contribution of the individual i to the log-likelihood

Definition at line 223 of file STK_IMixtureComposer.cpp.

{
#ifdef _OPENMP
  CPointX lnComp_(tik_.cols());
#endif
  // get maximal element of ln(x_i,\theta_k) + ln(p_k)
  for (int k=tik_.beginCols(); k< tik_.endCols(); k++)
  { lnComp_[k] = std::log(pk_[k])+lnComponentProbability(i,k);}
  if (lnComp_.isInfinite().any())
  {
#ifdef STK_MIXTURE_VERY_VERBOSE
    stk_cout << _T("IMixtureComposer::eStep(") << i << _T(")\n");
    stk_cout << _T("pk_ =") << pk_;
    stk_cout << _T("tk_ =") << tk_;
    stk_cout << _T("lnComp_ =") << lnComp_;
    for (int k=tik_.beginCols(); k< tik_.endCols(); k++)
    {
      Real val = lnComponentProbability(i,k);
      stk_cout << _T("lnComponentProbability(") << i << _T(",")<< k << _T(")=") << val << _T("\n");
    }
    writeParameters(stk_cout);
    stk_cout << _T("tik_ =\n") << tik_.transpose();
#endif
    throw(Clust::eStepFail_);
  }
  int kmax;
  Real max = lnComp_.maxElt(kmax);
  // set zi_
  zi_[i] = kmax;
  // return  max + sum_k p_k exp{lnCom_k - lnComp_kmax}
  Real sum =  (tik_.row(i) = (lnComp_ - max).exp()).sum();
  tik_.row(i) /= sum;
  return max + std::log( sum );
}

References _T, STK::Clust::eStepFail_, STK::ExprBase< Derived >::isInfinite(), lnComp_, STK::IMixtureStatModel::lnComponentProbability(), STK::max(), STK::ExprBase< Derived >::maxElt(), STK::IMixtureStatModel::pk_, STK::ICArray< Derived >::row(), stk_cout, STK::sum(), STK::IMixtureStatModel::tik_, STK::IMixtureStatModel::tk_, STK::ArrayBase< Derived >::transpose(), STK::IMixtureStatModel::writeParameters(), and STK::IMixtureStatModel::zi_.

finalizeStep()

void STK::IMixtureComposer::finalizeStep ( )

virtual

Finalize the estimation of the model.

Compute obtained lnLikelihood and set state to finalized.

Reimplemented from STK::IMixtureStatModel.

Reimplemented in STK::MixtureComposer.

Definition at line 281 of file STK_IMixtureComposer.cpp.

{
  setLnLikelihood(computeLnLikelihood());
  setState(Clust::modelFinalized_);
}

References STK::IMixtureStatModel::computeLnLikelihood(), STK::Clust::modelFinalized_, STK::IStatModelBase::setLnLikelihood(), and setState().

Referenced by STK::IMixtureAlgoPredict::predictBayesClassifier(), STK::EMPredict::run(), STK::SemiSEMPredict::run(), and STK::StrategyFacade::run().

initializeMixtureParameters()

void STK::IMixtureComposer::initializeMixtureParameters ( )

protectedvirtual

Create the mixture model parameters pk_ and tik_.

Default implementation is to set pk_ and tik_ arrays to 1/K value.

Reimplemented in STK::MixtureSemiLearner.

Definition at line 290 of file STK_IMixtureComposer.cpp.

{
#ifdef STK_MIXTURE_VERBOSE
  stk_cout << _T("Entering IMixtureComposer::initializeMixtureParameters\n");
#endif
  pk_  = 1./nbCluster_;
  tik_ = 1./nbCluster_;
  tk_  = sumByCol(tik_);
}

References _T, STK::IMixtureStatModel::nbCluster_, STK::IMixtureStatModel::pk_, stk_cout, STK::sumByCol(), STK::IMixtureStatModel::tik_, and STK::IMixtureStatModel::tk_.

Referenced by initializeStep().

initializeStep()

void STK::IMixtureComposer::initializeStep ( )

virtual

Initialize the model before its first use.

Initialize the values of the mixture parameters pk_ and @© tik_ using virtual method initializeMixtureParameters() and compute tk_ and zk_ using the virtual methods mapStep() and pStep().

Reimplemented from STK::IMixtureStatModel.

Definition at line 74 of file STK_IMixtureComposer.cpp.

{
#ifdef STK_MIXTURE_VERBOSE
  stk_cout << _T("Entering IMixtureComposer::initializeStep\n");
#endif
  // (re)initialize the mixture parameters tik and pk. (virtual method)
  initializeMixtureParameters();
  // compute nk
  tk_  = Stat::sumByCol(tik_);
  // (re) initialize mixtures
  IMixtureStatModel::initializeStep();
  // compute zi (virtual method)
  mapStep();
  // compute proportions pk_ (virtual method)
  pStep();
  // (re)initialize the likelihood
  setLnLikelihood(computeLnLikelihood());
  // compute the number of free parameters
  setNbFreeParameter(computeNbFreeParameters());
  // update state
  setState(Clust::modelInitialized_);
#ifdef STK_MIXTURE_VERBOSE
  stk_cout << _T("IMixtureComposer::initializeStep done\n");
#endif
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("Model parameters\n");
  this->writeParameters(stk_cout);
#endif
}

References _T, STK::IMixtureStatModel::computeLnLikelihood(), STK::IMixtureStatModel::computeNbFreeParameters(), initializeMixtureParameters(), STK::IMixtureStatModel::initializeStep(), mapStep(), STK::Clust::modelInitialized_, pStep(), STK::IStatModelBase::setLnLikelihood(), STK::IStatModelBase::setNbFreeParameter(), setState(), stk_cout, STK::Stat::sumByCol(), STK::IMixtureStatModel::tik_, STK::IMixtureStatModel::tk_, and STK::IMixtureStatModel::writeParameters().

Referenced by STK::IMixtureAlgoPredict::burnStep(), STK::MixtureComposer::createComposer(), STK::IMixtureAlgoPredict::predictBayesClassifier(), randomClassInit(), randomFuzzyInit(), STK::MixtureComposer::randomInit(), STK::EMPredict::run(), STK::SemiSEMPredict::run(), STK::RandomInit::run(), STK::ClassInit::run(), and STK::FuzzyInit::run().

mapStep() [1/2]

void STK::IMixtureComposer::mapStep ( )

virtual

Compute zi using the Map estimate.

Reimplemented in STK::MixtureSemiLearner.

Definition at line 205 of file STK_IMixtureComposer.cpp.

{
  for (int i = zi_.begin(); i< zi_.end(); ++i)
  { mapStep(i);  }
}

References mapStep(), and STK::IMixtureStatModel::zi_.

Referenced by initializeStep(), mapStep(), STK::MixtureSemiLearner::mapStep(), STK::IMixtureAlgoPredict::predictBayesClassifier(), randomFuzzyInit(), STK::MixtureComposer::randomInit(), STK::EMPredict::run(), and STK::SemiSEMPredict::run().

mapStep() [2/2]

void STK::IMixtureComposer::mapStep ( int  i )

protected

Compute zi using the Map estimate for i fixed.

Definition at line 259 of file STK_IMixtureComposer.cpp.

{
  int k;
  tik_.row(i).maxElt(k);
  zi_[i] = k;
}

References STK::ICArray< Derived >::row(), STK::IMixtureStatModel::tik_, and STK::IMixtureStatModel::zi_.

paramUpdateStep()

void STK::IMixtureComposer::paramUpdateStep ( )

pure virtual

Compute the proportions and the model parameters given the current tik mixture parameters.

Implemented in STK::MixtureComposer.

Definition at line 269 of file STK_IMixtureComposer.cpp.

{ pStep();
  /* implement specific parameters estimation in concrete class. */
}

References pStep().

Referenced by randomClassInit(), randomFuzzyInit(), STK::EMAlgo::run(), STK::CEMAlgo::run(), STK::SEMAlgo::run(), and STK::SemiSEMAlgo::run().

pStep()

void STK::IMixtureComposer::pStep ( )

virtual

Compute proportions using the ML estimates, default implementation.

Set as virtual in case we impose fixed proportions in derived Composer.

Reimplemented in STK::MixtureComposerFixedProp, and STK::MixtureSemiLearnerFixedProp.

Definition at line 275 of file STK_IMixtureComposer.cpp.

{ pk_ = Stat::meanByCol(tik_);}

References STK::Stat::meanByCol(), STK::IMixtureStatModel::pk_, and STK::IMixtureStatModel::tik_.

Referenced by initializeStep(), paramUpdateStep(), STK::EMAlgo::run(), STK::CEMAlgo::run(), STK::SEMAlgo::run(), and STK::SemiSEMAlgo::run().

randomClassInit()

void STK::IMixtureComposer::randomClassInit

(

)

Initialize randomly the labels zi of the model.

Initialize the model parameters using initializeStep() if it has not been already called. Simulate the zi, compute tik using eStep(), update the parameters using paramUpdateStep() and terminate using eStep().

Definition at line 105 of file STK_IMixtureComposer.cpp.

{
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("Entering IMixtureComposer::randomClassInit()\n");
#endif
  // initialize mixture model if necessary
  if (state() < Clust::modelInitialized_) { initializeStep();}
  // generate random zi, compute tik and nk
  if (randomZi()<2) throw(Clust::randomClassInitFail_);
  // update parameters
  paramUpdateStep();
  // compute eStep()
  eStep();
  // model initialized
  setState(Clust::modelParamInitialized_);
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("IMixtureComposer::randomClassInit() done\n");
#endif
}

References _T, eStep(), initializeStep(), STK::Clust::modelInitialized_, STK::Clust::modelParamInitialized_, paramUpdateStep(), STK::Clust::randomClassInitFail_, randomZi(), setState(), state(), and stk_cout.

Referenced by STK::ClassInit::run().

randomFuzzyInit()

void STK::IMixtureComposer::randomFuzzyInit

(

)

Initialize randomly the posterior probabilities tik of the model, then compute the zi values with mapStep, compute the initial parameter values using paramUpdateStep, and compute the tik using eStep().

Definition at line 126 of file STK_IMixtureComposer.cpp.

{
#ifdef STK_MIXTURE_VERBOSE
  stk_cout << _T("Entering IMixtureComposer::randomFuzzyInit(). state= ") << state() << _T("\n");
#endif
  // initialize mixture model if necessary
  if (state() < Clust::modelInitialized_) { initializeStep();}
  // create random tik and compute nk
  if (randomTik()<2) throw(Clust::randomFuzzyInitFail_);
  // compute zi
  mapStep();
  // update paramters values
  paramUpdateStep();
  // eStep
  eStep();
  // model intialized
  setState(Clust::modelParamInitialized_);
#ifdef STK_MIXTURE_VERBOSE
  stk_cout << _T("IMixtureComposer::randomFuzzyInit() done\n");
#endif
}

References _T, eStep(), initializeStep(), mapStep(), STK::Clust::modelInitialized_, STK::Clust::modelParamInitialized_, paramUpdateStep(), STK::Clust::randomFuzzyInitFail_, randomTik(), setState(), state(), and stk_cout.

Referenced by STK::FuzzyInit::run(), STK::FullStrategy::run(), STK::SimpleStrategy::run(), and STK::XemStrategy::run().

randomInit()

virtual void STK::IMixtureComposer::randomInit ( )

pure virtual

initialize randomly the parameters of the components of the model

Implements STK::IMixtureStatModel.

Implemented in STK::MixtureComposer.

Referenced by STK::RandomInit::run().

randomTik()

int STK::IMixtureComposer::randomTik ( )

protectedvirtual

generate random tik_

Reimplemented in STK::MixtureSemiLearner.

Definition at line 301 of file STK_IMixtureComposer.cpp.

{
  tk_ = 0.;
  tik_.randUnif();
  for (int i = tik_.beginRows(); i < tik_.endRows(); ++i)
  {
    // create a reference on the i-th row
    CPointX tikRowi(tik_.row(i), true);
    tikRowi = tikRowi * pk_;
    tikRowi /= tikRowi.sum();
    tk_ += tikRowi;
  }
  return tk_.minElt();
}

References STK::ExprBase< Derived >::minElt(), STK::IMixtureStatModel::pk_, STK::ArrayBase< Derived >::randUnif(), STK::ICArray< Derived >::row(), STK::IMixtureStatModel::tik_, and STK::IMixtureStatModel::tk_.

Referenced by randomFuzzyInit(), and STK::MixtureComposer::randomInit().

randomZi()

int STK::IMixtureComposer::randomZi ( )

protectedvirtual

generate random zi_

Reimplemented in STK::MixtureSemiLearner.

Definition at line 317 of file STK_IMixtureComposer.cpp.

{
  Law::Categorical law(pk_);
  zi_.rand(law);
  return cStep();
}

References cStep(), STK::IMixtureStatModel::pk_, STK::ArrayBase< Derived >::rand(), and STK::IMixtureStatModel::zi_.

Referenced by randomClassInit().

setState()

void STK::IMixtureComposer::setState ( Clust::modelState  state )

inline

set the state of the model : should be used by any strategy

Definition at line 105 of file STK_IMixtureComposer.h.

{ state_ = state;}

References state(), and state_.

Referenced by finalizeStep(), initializeStep(), randomClassInit(), randomFuzzyInit(), STK::MixtureComposer::randomInit(), and STK::FullStrategy::run().

sStep() [1/2]

int STK::IMixtureComposer::sStep ( )

virtual

Simulate zi accordingly to tik and replace tik by zik by calling cStep().

Returns

the minimal value of individuals in a class

Reimplemented in STK::MixtureSemiLearner.

Definition at line 167 of file STK_IMixtureComposer.cpp.

{
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("Entering IMixtureComposer::sStep()\n");
#endif
  // simulate zi
  for (int i = zi_.begin(); i< zi_.end(); ++i) { sStep(i);}
#ifdef STK_MIXTURE_VERY_VERBOSE
  stk_cout << _T("IMixtureComposer::sStep() done\n");
#endif
  return cStep();
}

References _T, cStep(), sStep(), stk_cout, and STK::IMixtureStatModel::zi_.

Referenced by STK::IMixtureAlgoPredict::burnStep(), STK::SEMAlgo::run(), sStep(), and STK::MixtureSemiLearner::sStep().

sStep() [2/2]

void STK::IMixtureComposer::sStep ( int  i )

protected

Simulate zi accordingly to tik.

Parameters

i	index of the the individual

Definition at line 214 of file STK_IMixtureComposer.cpp.

{ zi_.elt(i) = Law::Categorical::rand(tik_.row(i));}

References STK::Law::Categorical::rand(), STK::ICArray< Derived >::row(), STK::IMixtureStatModel::tik_, and STK::IMixtureStatModel::zi_.

state()

Clust::modelState STK::IMixtureComposer::state ( ) const

inline

Returns

the state of the model

Definition at line 102 of file STK_IMixtureComposer.h.

{ return state_;}

References state_.

Referenced by randomClassInit(), randomFuzzyInit(), STK::MixtureComposer::randomInit(), STK::SimpleStrategy::run(), STK::XemStrategy::run(), and setState().

Member Data Documentation

lnComp_

CPointX STK::IMixtureComposer::lnComp_

private

Auxiliary array used in the eStep.

Definition at line 194 of file STK_IMixtureComposer.h.

Referenced by eStep().

state_

Clust::modelState STK::IMixtureComposer::state_

private

state of the model

Definition at line 191 of file STK_IMixtureComposer.h.

Referenced by setState(), and state().

---

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

• STK_IMixtureComposer.h
• STK_IMixtureComposer.cpp