A joint Bernoulli model is a statistical model of the form: following form. More...

#include <STK_ModelBernoulli_pj.h>

Inheritance diagram for STK::ModelBernoulli_pj< Data_, WColVector_ >:

Public Types
typedef DataBridge< Data_ >	Data
	Type of the container storing the data.

typedef hidden::Traits< Data_ >::Row	RowVector

typedef WColVector_	WColVector
	Type of the array storing the weights of the data.

typedef hidden::Traits< Data_ >::Type	Type
	Type of the data in the container.

typedef IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > >	Base
	Base class.

Public Types inherited from STK::IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > >
typedef hidden::StatModelTraits< ModelBernoulli_pj< Data_, WColVector_ > >::Data	Data
	Type of the container with the data.

typedef hidden::StatModelTraits< ModelBernoulli_pj< Data_, WColVector_ > >::Type	Type
	Type of the data in the container.

typedef hidden::StatModelTraits< ModelBernoulli_pj< Data_, WColVector_ > >::WColVector	WColVector
	Type of the vector with the weights.

typedef hidden::StatModelTraits< ModelBernoulli_pj< Data_, WColVector_ > >::Parameters	Parameters
	Type of the parameters of the Model.

Public Member Functions
	ModelBernoulli_pj ()
	default constructor.

	ModelBernoulli_pj (Data const &data)
	Constructor with data set.

	ModelBernoulli_pj (Data const *p_data)
	Constructor with a ptr on the data set.

	ModelBernoulli_pj (ModelBernoulli_pj const &model)
	Copy constructor.

	~ModelBernoulli_pj ()
	destructor

CPointX const &	prob () const

CPointX const &	lnProb () const
	vector of the log probabilities of the observations

CPointX const &	ln1mProb () const
	vector of the log probabilities of the reversed observations

int	computeNbFreeParameters () const
	compute the number of free parameters

Real	computeLnLikelihood (RowVector const &rowData) const
	compute the log Likelihood of an observation.

void	computeParameters ()
	compute the parameters

void	computeParameters (WColVector const &weights)
	compute the weighted parameters

void	writeParametersImpl (ostream &os) const
	Write the parameters on the output stream os.

Public Member Functions inherited from STK::IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > >
Data const *const	p_data () const

Parameters const &	param () const

String const &	error () const

void	setData (Data const *p_data)
	Set the data set.

bool	run ()
	Estimate the parameters of the model and update the model.

bool	run (WColVector const &weights)
	compute the weighted empirical probability of success based on the observed variables.

void	writeParameters (ostream &os)

void	writeParametersImpl (ostream &os)
	default implementation of the writeParameters method.

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

Public Member Functions inherited from STK::IRecursiveTemplate< Derived >
Derived &	asDerived ()
	static cast : return a reference of this with a cast to the derived class.

Derived const &	asDerived () const
	static cast : return a const reference of this with a cast to the derived class.

Derived *	asPtrDerived ()
	static cast : return a ptr on a `Derived` of this with a cast to the derived class.

Derived const *	asPtrDerived () const
	static cast : return a ptr on a constant `Derived` of this with a cast to the derived class.

Derived *	clone () const
	create a leaf using the copy constructor of the Derived class.

Derived *	clone (bool isRef) const
	create a leaf using the copy constructor of the Derived class and a flag determining if the clone is a reference or not.

Additional Inherited Members
Protected Member Functions inherited from STK::IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > >
	IMultiStatModel ()
	default constructor.

	IMultiStatModel (Data const &data)
	Constructor with data set.

	IMultiStatModel (Data const *p_data)
	Constructor with a ptr on the data set.

	IMultiStatModel (IMultiStatModel const &model)
	Copy constructor.

	~IMultiStatModel ()
	destructor

Parameters &	param ()

Real	computeLnLikelihood () const
	compute the log Likelihood of the statistical model.

void	update ()
	update the model if a new data set is set

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

Protected Member Functions inherited from STK::IRecursiveTemplate< Derived >
	IRecursiveTemplate ()
	constructor.

	~IRecursiveTemplate ()
	destructor.

Protected Attributes inherited from STK::IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > >
Data const *	p_data_
	Pointer on the parameters of the model.

Parameters	param_
	Pointer on the parameters of the model.

String	msg_error_
	String with the last error message.

Detailed Description

template<class Data_, class WColVector_>
class STK::ModelBernoulli_pj< Data_, WColVector_ >

A joint Bernoulli model is a statistical model of the form: following form.

$f(\mathbf{x}_i = x|\theta) = \prod_{j=1}^p p_{j}^{x_i^j} (1-p_{j})^{1-x_i^j}, \quad x_i^j\in\{0,1\}, \quad j=1,\ldots,p, \quad i=1,\ldots,n.$

Definition at line 142 of file STK_ModelBernoulli_pj.h.

Member Typedef Documentation

◆ Base

template<class Data_ , class WColVector_ >

typedef IMultiStatModel< ModelBernoulli_pj<Data_, WColVector_> > STK::ModelBernoulli_pj< Data_, WColVector_ >::Base

Base class.

Definition at line 153 of file STK_ModelBernoulli_pj.h.

◆ Data

template<class Data_ , class WColVector_ >

typedef DataBridge<Data_> STK::ModelBernoulli_pj< Data_, WColVector_ >::Data

Type of the container storing the data.

Definition at line 146 of file STK_ModelBernoulli_pj.h.

◆ RowVector

template<class Data_ , class WColVector_ >

typedef hidden::Traits<Data_>::Row STK::ModelBernoulli_pj< Data_, WColVector_ >::RowVector

Definition at line 147 of file STK_ModelBernoulli_pj.h.

◆ Type

template<class Data_ , class WColVector_ >

typedef hidden::Traits<Data_>::Type STK::ModelBernoulli_pj< Data_, WColVector_ >::Type

Type of the data in the container.

Definition at line 151 of file STK_ModelBernoulli_pj.h.

◆ WColVector

template<class Data_ , class WColVector_ >

typedef WColVector_ STK::ModelBernoulli_pj< Data_, WColVector_ >::WColVector

Type of the array storing the weights of the data.

Definition at line 149 of file STK_ModelBernoulli_pj.h.

Constructor & Destructor Documentation

◆ ModelBernoulli_pj() [1/4]

template<class Data_ , class WColVector_ >

STK::ModelBernoulli_pj< Data_, WColVector_ >::ModelBernoulli_pj ( )

inline

default constructor.

Definition at line 158 of file STK_ModelBernoulli_pj.h.

158: Base() {}

STK::ModelBernoulli_pj::Base

IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > > Base

Base class.

Definition STK_ModelBernoulli_pj.h:153

◆ ModelBernoulli_pj() [2/4]

template<class Data_ , class WColVector_ >

STK::ModelBernoulli_pj< Data_, WColVector_ >::ModelBernoulli_pj ( Data const & data )

inline

Constructor with data set.

Definition at line 160 of file STK_ModelBernoulli_pj.h.

160: Base(data) {}

◆ ModelBernoulli_pj() [3/4]

template<class Data_ , class WColVector_ >

STK::ModelBernoulli_pj< Data_, WColVector_ >::ModelBernoulli_pj ( Data const * p_data )

inline

Constructor with a ptr on the data set.

Definition at line 162 of file STK_ModelBernoulli_pj.h.

162: Base(p_data) {}

STK::IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > >::p_data

Data const *const p_data() const

Definition STK_IMultiStatModel.h:152

◆ ModelBernoulli_pj() [4/4]

template<class Data_ , class WColVector_ >

STK::ModelBernoulli_pj< Data_, WColVector_ >::ModelBernoulli_pj ( ModelBernoulli_pj< Data_, WColVector_ > const & model )

inline

Copy constructor.

Definition at line 164 of file STK_ModelBernoulli_pj.h.

164: Base(model) {}

◆ ~ModelBernoulli_pj()

template<class Data_ , class WColVector_ >

STK::ModelBernoulli_pj< Data_, WColVector_ >::~ModelBernoulli_pj ( )

inline

destructor

Definition at line 166 of file STK_ModelBernoulli_pj.h.

166{}

Member Function Documentation

◆ computeLnLikelihood()

template<class Data_ , class WColVector_ >

Real STK::ModelBernoulli_pj< Data_, WColVector_ >::computeLnLikelihood ( RowVector const & rowData ) const

compute the log Likelihood of an observation.

Definition at line 191 of file STK_ModelBernoulli_pj.h.

{
  Real sum =0.;
  for (Integer j= rowData.begin(); j < rowData.end(); ++j)
  {
    sum += rowData[j] * lnProb()[j] + (1-rowData[j] * ln1mProb()[j] );
  }
  return sum;
}

References STK::sum().

◆ computeNbFreeParameters()

template<class Data_ , class WColVector_ >

int STK::ModelBernoulli_pj< Data_, WColVector_ >::computeNbFreeParameters ( ) const

inline

compute the number of free parameters

Definition at line 176 of file STK_ModelBernoulli_pj.h.

176{ return p_data()->dataij().sizeCols();}

References STK::IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > >::p_data().

◆ computeParameters() [1/2]

template<class Data_ , class WColVector_ >

void STK::ModelBernoulli_pj< Data_, WColVector_ >::computeParameters ( )

compute the parameters

Definition at line 203 of file STK_ModelBernoulli_pj.h.

{
  for (int j=p_data()->dataij().beginCols(); j < p_data()->dataij().endCols(); ++j)
  {
    Real sum=0.;
    int nbObs=p_data()->dataij().sizeRows();
 
    for (int i=p_data()->dataij().beginRows(); i<p_data()->dataij().endRows(); ++i)
    { (p_data()->dataij().elt(i,j) == binaryNA_) ? --nbObs : sum += p_data()->dataij().elt(i,j);}
    if (nbObs != 0) { param().prob_[j] = sum/nbObs;}
               else { param().prob_[j] = Arithmetic<Real>::NA();}
  }
}

References STK::binaryNA_, STK::Arithmetic< Type >::NA(), and STK::sum().

◆ computeParameters() [2/2]

template<class Data_ , class WColVector_ >

void STK::ModelBernoulli_pj< Data_, WColVector_ >::computeParameters ( WColVector const & weights )

compute the weighted parameters

Parameters

weights the weights of the samples

Definition at line 218 of file STK_ModelBernoulli_pj.h.

{
  // compute
  for (int j=p_data()->dataij().beginCols(); j < p_data()->dataij().endCols(); ++j)
  {
    Real sum=0., wsum = 0.;
    for (int i=p_data()->dataij().beginRows(); i<p_data()->dataij().endRows(); ++i)
    { if (p_data()->dataij().elt(i,j) != binaryNA_)
      { sum  += weights[i]*p_data()->dataij().elt(i,j);
        wsum += weights[i];
      }
    }
    if (wsum != 0) { param().prob_[j] = sum/wsum;}
              else { param().prob_[j] = Arithmetic<Real>::NA();}
  }
}

References STK::binaryNA_, STK::Arithmetic< Type >::NA(), and STK::sum().

◆ ln1mProb()

template<class Data_ , class WColVector_ >

CPointX const & STK::ModelBernoulli_pj< Data_, WColVector_ >::ln1mProb ( ) const

inline

vector of the log probabilities of the reversed observations

Definition at line 173 of file STK_ModelBernoulli_pj.h.

173{ return param().ln1mProb();}

References STK::IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > >::param().

◆ lnProb()

template<class Data_ , class WColVector_ >

CPointX const & STK::ModelBernoulli_pj< Data_, WColVector_ >::lnProb ( ) const

inline

vector of the log probabilities of the observations

Definition at line 171 of file STK_ModelBernoulli_pj.h.

171{ return param().lnProb();}

References STK::IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > >::param().

◆ prob()

template<class Data_ , class WColVector_ >

CPointX const & STK::ModelBernoulli_pj< Data_, WColVector_ >::prob ( ) const

inline

Returns: the vector of the probabilities of the observations

Definition at line 169 of file STK_ModelBernoulli_pj.h.

169{ return param().prob();}

References STK::IMultiStatModel< ModelBernoulli_pj< Data_, WColVector_ > >::param().

◆ writeParametersImpl()

template<class Data_ , class WColVector_ >

void STK::ModelBernoulli_pj< Data_, WColVector_ >::writeParametersImpl ( ostream & os ) const

Write the parameters on the output stream os.

Definition at line 237 of file STK_ModelBernoulli_pj.h.

{
  os << _T("prob = ") << prob();
  os << _T("lnProb = ") << lnProb();
  os << _T("ln1mProb = ") << ln1mProb();
}

References _T.

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

STK_ModelBernoulli_pj.h

Public Types

Public Member Functions

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ Base

◆ Data

◆ RowVector

◆ Type

◆ WColVector

Constructor & Destructor Documentation

◆ ModelBernoulli_pj() [1/4]

◆ ModelBernoulli_pj() [2/4]

◆ ModelBernoulli_pj() [3/4]

◆ ModelBernoulli_pj() [4/4]

◆ ~ModelBernoulli_pj()

Member Function Documentation

◆ computeLnLikelihood()

◆ computeNbFreeParameters()

◆ computeParameters() [1/2]

◆ computeParameters() [2/2]

◆ ln1mProb()

◆ lnProb()

◆ prob()

◆ writeParametersImpl()