QUESO::InvLogitGaussianVectorRV< V, M > Class Template Reference

A class representing a (transformed) Gaussian vector RV with bounds. More...

#include <InvLogitGaussianVectorRV.h>

Inheritance diagram for QUESO::InvLogitGaussianVectorRV< V, M >:

Collaboration diagram for QUESO::InvLogitGaussianVectorRV< V, M >:

Public Member Functions
Constructor/Destructor methods
	InvLogitGaussianVectorRV (const char *prefix, const BoxSubset< V, M > &imageBoxSubset, const V &lawExpVector, const V &lawVarVector)
	Constructor. More...
	InvLogitGaussianVectorRV (const char *prefix, const BoxSubset< V, M > &imageBoxSubset, const V &lawExpVector, const M &lawCovMatrix)
	Constructor. More...
virtual	~InvLogitGaussianVectorRV ()
	Virtual destructor. More...
Statistical methods
void	updateLawExpVector (const V &newLawExpVector)
	Updates the vector that contains the mean values for the underlying Gaussian. More...
void	updateLawCovMatrix (const M &newLawCovMatrix)
	Updates the covariance matrix. More...
I/O methods
void	print (std::ostream &os) const
	TODO: Prints the vector RV. More...
Public Member Functions inherited from QUESO::BaseVectorRV< V, M >
	BaseVectorRV (const char *prefix, const VectorSet< V, M > &imageSet)
	Constructor. More...
virtual	~BaseVectorRV ()
	Virtual destructor. More...
const BaseEnvironment &	env () const
	QUESO environment; access to private attribute m_env. More...
const VectorSet< V, M > &	imageSet () const
	Image set of the vector RV; access to private attribute m_imageSet. More...
const BaseJointPdf< V, M > &	pdf () const
	Posterior Density Function of the vector RV; access to private attribute m_pdf. More...
const BaseVectorRealizer< V, M > &	realizer () const
	Finds a realization (sample) of the PDF of this vector RV; access to private attribute m_realizer. More...
const BaseVectorCdf< V, M > &	subCdf () const
	Finds the Cumulative Distribution Function of this vector RV, considering only the sub-sequence of data; access to private attribute m_subCdf. More...
const BaseVectorCdf< V, M > &	unifiedCdf () const
	Finds the Cumulative Distribution Function of this vector RV, considering the unified sequence of data; access to private attribute m_unifiedCdf. More...
const BaseVectorMdf< V, M > &	mdf () const
	Finds the Mass Density Function of this vector RV; access to private attribute m_mdf. More...

Additional Inherited Members
Protected Attributes inherited from QUESO::BaseVectorRV< V, M >
const BaseEnvironment &	m_env
std::string	m_prefix
const VectorSet< V, M > &	m_imageSet
BaseJointPdf< V, M > *	m_pdf
BaseVectorRealizer< V, M > *	m_realizer
const BaseVectorCdf< V, M > *	m_subCdf
const BaseVectorCdf< V, M > *	m_unifiedCdf
const BaseVectorMdf< V, M > *	m_mdf

Detailed Description

template<class V = GslVector, class M = GslMatrix>
class QUESO::InvLogitGaussianVectorRV< V, M >

A class representing a (transformed) Gaussian vector RV with bounds.

This class allows the user to compute the value of a (transoformed) Gaussian PDF and to generate realizations (samples) from it.

Definition at line 45 of file InvLogitGaussianVectorRV.h.

Constructor & Destructor Documentation

template<class V , class M >

QUESO::InvLogitGaussianVectorRV< V, M >::InvLogitGaussianVectorRV	(	const char *	prefix,
		const BoxSubset< V, M > &	imageBoxSubset,
		const V &	lawExpVector,
		const V &	lawVarVector
	)

Constructor.

Construct a (transformed) Gaussian vector RV with mean lawExpVector (of the Gaussian, not the transformed Gaussian) and diagonal covariance matrix lawVarVector whose variates live in imageBoxSubset.

Definition at line 35 of file InvLogitGaussianVectorRV.C.

References QUESO::BaseEnvironment::displayVerbosity(), QUESO::BaseVectorRV< V, M >::m_env, QUESO::BaseVectorRV< V, M >::m_imageSet, QUESO::BaseVectorRV< V, M >::m_mdf, QUESO::BaseVectorRV< V, M >::m_pdf, QUESO::BaseVectorRV< V, M >::m_prefix, QUESO::BaseVectorRV< V, M >::m_realizer, QUESO::BaseVectorRV< V, M >::m_subCdf, QUESO::BaseVectorRV< V, M >::m_unifiedCdf, queso_require_greater_msg, and QUESO::BaseEnvironment::subDisplayFile().

  : BaseVectorRV<V, M>(((std::string)(prefix)+"invlogit_gau").c_str(),
      imageBoxSubset)
{
  if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) {
    *m_env.subDisplayFile() << "Entering InvLogitGaussianVectorRV<V,M>::constructor() [1]"
                            << ": prefix = " << m_prefix
                            << std::endl;
  }

  queso_require_greater_msg(lawVarVector.getMinValue(), 0.0, "Covariance matrix is not symmetric positive definite.");

  m_pdf = new InvLogitGaussianJointPdf<V,M>(m_prefix.c_str(),
      dynamic_cast<const BoxSubset<V, M> & >(m_imageSet), lawExpVector,
      lawVarVector);

  V cholDiag(lawVarVector);
  cholDiag.cwSqrt();
  M lowerCholLawCovMatrix(cholDiag);
  lowerCholLawCovMatrix.zeroUpper(false);

  m_realizer = new InvLogitGaussianVectorRealizer<V,M>(m_prefix.c_str(),
      dynamic_cast<const BoxSubset<V, M> & >(m_imageSet), lawExpVector,
      lowerCholLawCovMatrix);

  m_subCdf     = NULL; // FIX ME: complete code
  m_unifiedCdf = NULL; // FIX ME: complete code
  m_mdf        = NULL; // FIX ME: complete code

  if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) {
    *m_env.subDisplayFile() << "Leaving InvLogitGaussianVectorRV<V,M>::constructor() [1]"
                            << ": prefix = " << m_prefix
                            << std::endl;
  }
}

template<class V , class M >

QUESO::InvLogitGaussianVectorRV< V, M >::InvLogitGaussianVectorRV	(	const char *	prefix,
		const BoxSubset< V, M > &	imageBoxSubset,
		const V &	lawExpVector,
		const M &	lawCovMatrix
	)

Constructor.

Construct a (transformed) Gaussian vector RV with mean lawExpVector (of the Gaussian, not the transformed Gaussian) and covariance matrix lawCovMatrix whose variates live in imageBoxSubset.

Definition at line 76 of file InvLogitGaussianVectorRV.C.

References QUESO::BaseEnvironment::displayVerbosity(), QUESO::BaseVectorRV< V, M >::m_env, QUESO::BaseVectorRV< V, M >::m_imageSet, QUESO::BaseVectorRV< V, M >::m_mdf, QUESO::BaseVectorRV< V, M >::m_pdf, QUESO::BaseVectorRV< V, M >::m_prefix, QUESO::BaseVectorRV< V, M >::m_realizer, QUESO::BaseVectorRV< V, M >::m_subCdf, QUESO::BaseVectorRV< V, M >::m_unifiedCdf, queso_require_msg, and QUESO::BaseEnvironment::subDisplayFile().

  : BaseVectorRV<V, M>(((std::string)(prefix)+"invlogit_gau").c_str(),
      imageBoxSubset)
{
  if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) {
    *m_env.subDisplayFile() << "Entering InvLogitGaussianVectorRV<V,M>::constructor() [2]"
                            << ": prefix = " << m_prefix
                            << std::endl;
  }

  m_pdf = new InvLogitGaussianJointPdf<V, M>(m_prefix.c_str(),
      dynamic_cast<const BoxSubset<V, M> & >(m_imageSet),
      lawExpVector, lawCovMatrix);

  M lowerCholLawCovMatrix(lawCovMatrix);
  int iRC = lowerCholLawCovMatrix.chol();
  lowerCholLawCovMatrix.zeroUpper(false);
  if (iRC) {
    std::cerr << "In InvLogitGaussianVectorRV<V,M>::constructor() [2]: chol failed, will use svd\n";
    if (m_env.subDisplayFile()) {
      *m_env.subDisplayFile() << "In InvLogitGaussianVectorRV<V,M>::constructor() [2]: chol failed; will use svd; lawCovMatrix contents are\n";
      *m_env.subDisplayFile() << lawCovMatrix; // FIX ME: might demand parallelism
      *m_env.subDisplayFile() << std::endl;
    }
    M matU (lawCovMatrix);
    M matVt(m_imageSet.vectorSpace().zeroVector());
    V vecS (m_imageSet.vectorSpace().zeroVector());
    iRC = lawCovMatrix.svd(matU,vecS,matVt);
    queso_require_msg(!(iRC), "Cholesky decomposition of covariance matrix failed.");

    vecS.cwSqrt();
    m_realizer = new InvLogitGaussianVectorRealizer<V,M>(m_prefix.c_str(),
        dynamic_cast<const BoxSubset<V, M> & >(m_imageSet), lawExpVector, matU,
        vecS, matVt);
  }
  else {
    m_realizer = new InvLogitGaussianVectorRealizer<V, M>(m_prefix.c_str(),
        dynamic_cast<const BoxSubset<V, M> & >(m_imageSet), lawExpVector, lowerCholLawCovMatrix);
  }

  m_subCdf     = NULL; // FIX ME: complete code
  m_unifiedCdf = NULL; // FIX ME: complete code
  m_mdf        = NULL; // FIX ME: complete code

  if ((m_env.subDisplayFile()) && (m_env.displayVerbosity() >= 54)) {
    *m_env.subDisplayFile() << "Leaving InvLogitGaussianVectorRV<V,M>::constructor() [2]"
                            << ": prefix = " << m_prefix
                            << std::endl;
  }
}

template<class V , class M >

QUESO::InvLogitGaussianVectorRV< V, M >::~InvLogitGaussianVectorRV ( )

virtual

Virtual destructor.

Definition at line 132 of file InvLogitGaussianVectorRV.C.

{
  delete m_mdf;
  delete m_unifiedCdf;
  delete m_subCdf;
  delete m_realizer;
  delete m_pdf;
}

Member Function Documentation

template<class V , class M >

void QUESO::InvLogitGaussianVectorRV< V, M >::print ( std::ostream &  os ) const

virtual

TODO: Prints the vector RV.

Todo:

: implement me!

Implements QUESO::BaseVectorRV< V, M >.

Definition at line 192 of file InvLogitGaussianVectorRV.C.

{
  os << "InvLogitGaussianVectorRV<V,M>::print() says, 'Please implement me.'" << std::endl;
  return;
}

template<class V , class M >

void QUESO::InvLogitGaussianVectorRV< V, M >::updateLawCovMatrix

(

const M &

newLawCovMatrix

)

Updates the covariance matrix.

This method tries to use Cholesky decomposition; and if it fails, the method then calls a SVD decomposition.

Definition at line 156 of file InvLogitGaussianVectorRV.C.

References queso_require_msg.

Referenced by QUESO::TransformedScaledCovMatrixTKGroup< V, M >::updateLawCovMatrix().

{
  // We are sure that m_pdf (and m_realizer, etc) point to associated Gaussian
  // classes, so all is well
  (dynamic_cast<InvLogitGaussianJointPdf<V,M> * >(m_pdf))->updateLawCovMatrix(
      newLawCovMatrix);

  M newLowerCholLawCovMatrix(newLawCovMatrix);
  int iRC = newLowerCholLawCovMatrix.chol();
  newLowerCholLawCovMatrix.zeroUpper(false);
  if (iRC) {
    std::cerr << "In InvLogitGaussianVectorRV<V,M>::updateLawCovMatrix(): chol failed, will use svd\n";
    if (m_env.subDisplayFile()) {
      *m_env.subDisplayFile() << "In InvLogitGaussianVectorRV<V,M>::updateLawCovMatrix(): chol failed; will use svd; newLawCovMatrix contents are\n";
      *m_env.subDisplayFile() << newLawCovMatrix; // FIX ME: might demand parallelism
      *m_env.subDisplayFile() << std::endl;
    }
    M matU (newLawCovMatrix);
    M matVt(m_imageSet.vectorSpace().zeroVector());
    V vecS (m_imageSet.vectorSpace().zeroVector());
    iRC = newLawCovMatrix.svd(matU,vecS,matVt);
    queso_require_msg(!(iRC), "Cholesky decomposition of covariance matrix failed.");

    vecS.cwSqrt();
    (dynamic_cast<InvLogitGaussianVectorRealizer<V, M> * >(m_realizer))->
      updateLowerCholLawCovMatrix(matU, vecS, matVt);
  }
  else {
    (dynamic_cast<InvLogitGaussianVectorRealizer<V, M> * >(m_realizer))->
      updateLowerCholLawCovMatrix(newLowerCholLawCovMatrix);
  }
  return;
}

template<class V , class M >

void QUESO::InvLogitGaussianVectorRV< V, M >::updateLawExpVector

(

const V &

newLawExpVector

)

Updates the vector that contains the mean values for the underlying Gaussian.

Definition at line 143 of file InvLogitGaussianVectorRV.C.

Referenced by QUESO::TransformedScaledCovMatrixTKGroup< V, M >::rv().

{
  // We are sure that m_pdf (and m_realizer, etc) point to associated Gaussian
  // classes, so all is well
  (dynamic_cast<InvLogitGaussianJointPdf<V, M> * >(m_pdf))->updateLawExpVector(
      newLawExpVector);
  (dynamic_cast<InvLogitGaussianVectorRealizer<V, M> * >(m_realizer))->
    updateLawExpVector(newLawExpVector);
  return;
}

---

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

• src/stats/inc/InvLogitGaussianVectorRV.h
• src/stats/src/InvLogitGaussianVectorRV.C