SHOGUN
4.1.0
|
Class that models dual variational logit likelihood.
This likelihood model is described in the reference paper Mohammad Emtiyaz Khan, Aleksandr Y. Aravkin, Michael P. Friedlander, Matthias Seeger Fast Dual Variational Inference for Non-Conjugate Latent Gaussian Models. ICML2013
The mathematically definition (equation 19 in the paper) is as below
\[ Fenchel_i(\alpha_i,\lambda_i) = max_{h_i,\rho_i}{\alpha_i h_i+\lambda_i \rho_i /2 - E_{q(f_i|h_i,\rho_i)}(-log(p(y_i|f_i)))} \]
where \(\alpha_i\), \(\lambda_i\) are Lagrange multipliers with respective to constraints \(h_i=\mu_i\) and \(\rho_i=\sigma_i^2\) respectively, \(\mu\) and \(\sigma_i\) are variational Gaussian parameters, \(y_i\) is data label, \(q(f_i)\) is the variational Gaussian distribution, and \(p(y_i)\) is the data distribution to be specified. In this setting, \(\alpha\) and \(\lambda\) are called dual parameters for \(\mu\) and \(\sigma^2\) respectively.
Note that \(p(y_i)\) is Logistic distribution and a local variational bound defined as below is used to approximate -E_{q(f_i|h_i,)}(-log(p(y_i|f_i)))
The local variational bound used here is
\[ log(x) \leq t^{-1}x+log(t)-1 \]
, where t is a local variable and the inequality holds for every t>0. See Bernoulli-logit in Table 2 of the paper for detailed information
Definition at line 75 of file LogitDVGLikelihood.h.
Public Attributes | |
SGIO * | io |
Parallel * | parallel |
Version * | version |
Parameter * | m_parameters |
Parameter * | m_model_selection_parameters |
Parameter * | m_gradient_parameters |
uint32_t | m_hash |
Protected Member Functions | |
virtual void | init_likelihood () |
virtual void | precompute () |
virtual CVariationalGaussianLikelihood * | get_variational_likelihood () const |
virtual void | set_likelihood (CLikelihoodModel *lik) |
virtual void | load_serializable_pre () throw (ShogunException) |
virtual void | load_serializable_post () throw (ShogunException) |
virtual void | save_serializable_pre () throw (ShogunException) |
virtual void | save_serializable_post () throw (ShogunException) |
Protected Attributes | |
SGVector< float64_t > | m_lambda |
float64_t | m_strict_scale |
bool | m_is_valid |
SGVector< float64_t > | m_mu |
SGVector< float64_t > | m_s2 |
SGVector< float64_t > | m_lab |
CLikelihoodModel * | m_likelihood |
default constructor
Definition at line 47 of file LogitDVGLikelihood.cpp.
|
virtual |
Definition at line 53 of file LogitDVGLikelihood.cpp.
|
virtualinherited |
this method is used for adjusting step size to ensure the updated value satisfied lower/upper bound constrain
The updated value is defined as below. lambda_new = m_lambda + direction * step
direction | direction for m_lambda update |
step | original step size (non-negative) |
Definition at line 111 of file DualVariationalGaussianLikelihood.cpp.
|
inherited |
Builds a dictionary of all parameters in SGObject as well of those of SGObjects that are parameters of this object. Dictionary maps parameters to the objects that own them.
dict | dictionary of parameters to be built. |
Definition at line 597 of file SGObject.cpp.
|
virtualinherited |
Creates a clone of the current object. This is done via recursively traversing all parameters, which corresponds to a deep copy. Calling equals on the cloned object always returns true although none of the memory of both objects overlaps.
Definition at line 714 of file SGObject.cpp.
|
virtualinherited |
A deep copy. All the instance variables will also be copied.
Definition at line 198 of file SGObject.cpp.
|
virtual |
whether the lower bound is strict
Implements CDualVariationalGaussianLikelihood.
Definition at line 125 of file LogitDVGLikelihood.h.
|
virtualinherited |
check whether the dual parameters are valid or not.
Definition at line 182 of file DualVariationalGaussianLikelihood.cpp.
|
virtual |
whether the upper bound is strict
Implements CDualVariationalGaussianLikelihood.
Definition at line 119 of file LogitDVGLikelihood.h.
Recursively compares the current SGObject to another one. Compares all registered numerical parameters, recursion upon complex (SGObject) parameters. Does not compare pointers!
May be overwritten but please do with care! Should not be necessary in most cases.
other | object to compare with |
accuracy | accuracy to use for comparison (optional) |
tolerant | allows linient check on float equality (within accuracy) |
Definition at line 618 of file SGObject.cpp.
|
virtual |
get the derivative of the dual objective function with respect to param
param | parameter |
Implements CDualVariationalGaussianLikelihood.
Definition at line 93 of file LogitDVGLikelihood.cpp.
|
virtual |
get the lower bound for dual parameter (lambda)
Implements CDualVariationalGaussianLikelihood.
Definition at line 113 of file LogitDVGLikelihood.h.
evaluate the dual objective function
Implements CDualVariationalGaussianLikelihood.
Definition at line 74 of file LogitDVGLikelihood.cpp.
|
virtual |
get the upper bound for dual parameter (lambda)
Implements CDualVariationalGaussianLikelihood.
Definition at line 107 of file LogitDVGLikelihood.h.
|
virtualinherited |
get derivative of log likelihood \(log(p(y|f))\) with respect to given parameter
lab | labels used |
func | function location |
param | parameter |
Reimplemented from CLikelihoodModel.
Definition at line 88 of file VariationalLikelihood.cpp.
|
virtualinherited |
get derivative of log likelihood \(log(p(y|f))\) with respect to given hyperparameter Note that variational parameters are NOT considered as hyperparameters
param | parameter |
Implements CVariationalLikelihood.
Definition at line 90 of file DualVariationalGaussianLikelihood.cpp.
|
virtualinherited |
returns the first moment of a given (unnormalized) probability distribution \(q(f_i) = Z_i^-1 p(y_i|f_i)\mathcal{N}(f_i|\mu,\sigma^2)\), where \( Z_i=\int p(y_i|f_i)\mathcal{N}(f_i|\mu,\sigma^2) df_i\).
This method is useful for EP local likelihood approximation.
mu | mean of the \(\mathcal{N}(f_i|\mu,\sigma^2)\) |
s2 | variance of the \(\mathcal{N}(f_i|\mu,\sigma^2)\) |
lab | labels \(y_i\) |
i | index i |
Implements CLikelihoodModel.
Definition at line 140 of file VariationalLikelihood.cpp.
|
virtualinherited |
returns the first moment of a given (unnormalized) probability distribution \(q(f_i) = Z_i^-1 p(y_i|f_i)\mathcal{N}(f_i|\mu,\sigma^2)\) for each \(f_i\), where \( Z_i=\int p(y_i|f_i)\mathcal{N}(f_i|\mu,\sigma^2) df_i\).
Wrapper method which calls get_first_moment multiple times.
mu | mean of the \(\mathcal{N}(f_i|\mu,\sigma^2)\) |
s2 | variance of the \(\mathcal{N}(f_i|\mu,\sigma^2)\) |
lab | labels \(y_i\) |
Definition at line 72 of file LikelihoodModel.cpp.
|
inherited |
|
inherited |
|
inherited |
|
virtualinherited |
get derivative of log likelihood \(log(p(y|f))\) with respect to location function \(f\)
lab | labels used |
func | function location |
i | index, choices are 1, 2, and 3 for first, second, and third derivatives respectively |
Implements CLikelihoodModel.
Definition at line 125 of file VariationalLikelihood.cpp.
|
virtualinherited |
Returns the logarithm of the point-wise likelihood \(log(p(y_i|f_i))\) for each label \(y_i\).
One can evaluate log-likelihood like: \( log(p(y|f)) = \sum_{i=1}^{n} log(p(y_i|f_i))\)
lab | labels \(y_i\) |
func | values of the function \(f_i\) |
Implements CLikelihoodModel.
Definition at line 118 of file VariationalLikelihood.cpp.
|
virtualinherited |
Returns the log-likelihood \(log(p(y|f)) = \sum_{i=1}^{n} log(p(y_i|f_i))\) for each of the provided functions \( f \) in the given matrix.
Wrapper method which calls get_log_probability_f multiple times.
lab | labels \(y_i\) |
F | values of the function \(f_i\) where each column of the matrix is one function \( f \). |
Definition at line 51 of file LikelihoodModel.cpp.
|
virtualinherited |
returns the zeroth moment of a given (unnormalized) probability distribution:
\[ log(Z_i) = log\left(\int p(y_i|f_i) \mathcal{N}(f_i|\mu,\sigma^2) df_i\right) \]
for each \(f_i\).
mu | mean of the \(\mathcal{N}(f_i|\mu,\sigma^2)\) |
s2 | variance of the \(\mathcal{N}(f_i|\mu,\sigma^2)\) |
lab | labels \(y_i\) |
Implements CLikelihoodModel.
Definition at line 132 of file VariationalLikelihood.cpp.
|
virtualinherited |
get model type
Reimplemented from CLikelihoodModel.
Definition at line 112 of file VariationalLikelihood.cpp.
|
inherited |
Definition at line 498 of file SGObject.cpp.
|
inherited |
Returns description of a given parameter string, if it exists. SG_ERROR otherwise
param_name | name of the parameter |
Definition at line 522 of file SGObject.cpp.
|
inherited |
Returns index of model selection parameter with provided index
param_name | name of model selection parameter |
Definition at line 535 of file SGObject.cpp.
get the dual parameter (alpha) for variational mu
Note that alpha = m_lambda - label For detailed information, please refer to the paper.
Implements CDualVariationalGaussianLikelihood.
Definition at line 64 of file LogitDVGLikelihood.cpp.
|
virtual |
returns the name of the likelihood model
Reimplemented from CDualVariationalGaussianLikelihood.
Definition at line 87 of file LogitDVGLikelihood.h.
|
virtualinherited |
returns the logarithm of the predictive density of \(y_*\):
\[ log(p(y_*|X,y,x_*)) = log\left(\int p(y_*|f_*) p(f_*|X,y,x_*) df_*\right) \]
which approximately equals to
\[ log\left(\int p(y_*|f_*) \mathcal{N}(f_*|\mu,\sigma^2) df_*\right) \]
where normal distribution \(\mathcal{N}(\mu,\sigma^2)\) is an approximation to the posterior marginal \(p(f_*|X,y,x_*)\).
NOTE: if lab equals to NULL, then each \(y_*\) equals to one.
mu | posterior mean of a Gaussian distribution \(\mathcal{N}(\mu,\sigma^2)\), which is an approximation to the posterior marginal \(p(f_*|X,y,x_*)\) |
s2 | posterior variance of a Gaussian distribution \(\mathcal{N}(\mu,\sigma^2)\), which is an approximation to the posterior marginal \(p(f_*|X,y,x_*)\) |
lab | labels \(y_*\) |
Reimplemented in CSoftMaxLikelihood.
Definition at line 45 of file LikelihoodModel.cpp.
|
virtualinherited |
returns mean of the predictive marginal \(p(y_*|X,y,x_*)\)
NOTE: if lab equals to NULL, then each \(y_*\) equals to one.
mu | posterior mean of a Gaussian distribution \(\mathcal{N}(\mu,\sigma^2)\), which is an approximation to the posterior marginal \(p(f_*|X,y,x_*)\) |
s2 | posterior variance of a Gaussian distribution \(\mathcal{N}(\mu,\sigma^2)\), which is an approximation to the posterior marginal \(p(f_*|X,y,x_*)\) |
lab | labels \(y_*\) |
Implements CLikelihoodModel.
Definition at line 72 of file VariationalLikelihood.cpp.
|
virtualinherited |
returns variance of the predictive marginal \(p(y_*|X,y,x_*)\)
NOTE: if lab equals to NULL, then each \(y_*\) equals to one.
mu | posterior mean of a Gaussian distribution \(\mathcal{N}(\mu,\sigma^2)\), which is an approximation to the posterior marginal \(p(f_*|X,y,x_*)\) |
s2 | posterior variance of a Gaussian distribution \(\mathcal{N}(\mu,\sigma^2)\), which is an approximation to the posterior marginal \(p(f_*|X,y,x_*)\) |
lab | labels \(y_*\) |
Implements CLikelihoodModel.
Definition at line 80 of file VariationalLikelihood.cpp.
|
virtualinherited |
get derivative of the first derivative of log likelihood with respect to function location, i.e. \(\frac{\partial log(p(y|f))}{\partial f}\) with respect to given parameter
lab | labels used |
func | function location |
param | parameter |
Reimplemented from CLikelihoodModel.
Definition at line 96 of file VariationalLikelihood.cpp.
|
virtualinherited |
returns the second moment of a given (unnormalized) probability distribution \(q(f_i) = Z_i^-1 p(y_i|f_i)\mathcal{N}(f_i|\mu,\sigma^2)\), where \( Z_i=\int p(y_i|f_i)\mathcal{N}(f_i|\mu,\sigma^2) df_i\).
This method is useful for EP local likelihood approximation.
mu | mean of the \(\mathcal{N}(f_i|\mu,\sigma^2)\) |
s2 | variance of the \(\mathcal{N}(f_i|\mu,\sigma^2)\) |
lab | labels \(y_i\) |
i | index i |
Implements CLikelihoodModel.
Definition at line 148 of file VariationalLikelihood.cpp.
|
virtualinherited |
returns the second moment of a given (unnormalized) probability distribution \(q(f_i) = Z_i^-1 p(y_i|f_i)\mathcal{N}(f_i|\mu,\sigma^2)\) for each \(f_i\), where \( Z_i=\int p(y_i|f_i)\mathcal{N}(f_i|\mu,\sigma^2) df_i\).
Wrapper method which calls get_second_moment multiple times.
mu | mean of the \(\mathcal{N}(f_i|\mu,\sigma^2)\) |
s2 | variance of the \(\mathcal{N}(f_i|\mu,\sigma^2)\) |
lab | labels \(y_i\) |
Definition at line 89 of file LikelihoodModel.cpp.
|
virtualinherited |
get derivative of the second derivative of log likelihood with respect to function location, i.e. \(\frac{\partial^{2} log(p(y|f))}{\partial f^{2}}\) with respect to given parameter
lab | labels used |
func | function location |
param | parameter |
Reimplemented from CLikelihoodModel.
Definition at line 104 of file VariationalLikelihood.cpp.
get the dual parameter (lambda) for variational s2
Implements CDualVariationalGaussianLikelihood.
Definition at line 57 of file LogitDVGLikelihood.cpp.
returns the expection of the logarithm of a given probability distribution wrt the variational distribution given m_mu and m_s2
Implements CVariationalLikelihood.
Definition at line 66 of file DualVariationalGaussianLikelihood.cpp.
|
virtualinherited |
get derivative of the variational expection of log likelihood with respect to given parameter
param | parameter |
Implements CVariationalLikelihood.
Definition at line 78 of file DualVariationalGaussianLikelihood.cpp.
|
protectedvirtualinherited |
this method is used to dynamic-cast the likelihood model, m_likelihood, to variational likelihood model.
Definition at line 55 of file DualVariationalGaussianLikelihood.cpp.
|
protectedvirtual |
this method is called to initialize m_likelihood in init()
Implements CVariationalGaussianLikelihood.
Definition at line 121 of file LogitDVGLikelihood.cpp.
|
virtualinherited |
If the SGSerializable is a class template then TRUE will be returned and GENERIC is set to the type of the generic.
generic | set to the type of the generic if returning TRUE |
Definition at line 296 of file SGObject.cpp.
|
virtualinherited |
Load this object from file. If it will fail (returning FALSE) then this object will contain inconsistent data and should not be used!
file | where to load from |
prefix | prefix for members |
Definition at line 369 of file SGObject.cpp.
|
protectedvirtualinherited |
Can (optionally) be overridden to post-initialize some member variables which are not PARAMETER::ADD'ed. Make sure that at first the overridden method BASE_CLASS::LOAD_SERIALIZABLE_POST is called.
ShogunException | will be thrown if an error occurs. |
Reimplemented in CKernel, CWeightedDegreePositionStringKernel, CList, CAlphabet, CLinearHMM, CGaussianKernel, CInverseMultiQuadricKernel, CCircularKernel, and CExponentialKernel.
Definition at line 426 of file SGObject.cpp.
|
protectedvirtualinherited |
Can (optionally) be overridden to pre-initialize some member variables which are not PARAMETER::ADD'ed. Make sure that at first the overridden method BASE_CLASS::LOAD_SERIALIZABLE_PRE is called.
ShogunException | will be thrown if an error occurs. |
Reimplemented in CDynamicArray< T >, CDynamicArray< float64_t >, CDynamicArray< float32_t >, CDynamicArray< int32_t >, CDynamicArray< char >, CDynamicArray< bool >, and CDynamicObjectArray.
Definition at line 421 of file SGObject.cpp.
|
virtualinherited |
Definition at line 262 of file SGObject.cpp.
|
protectedvirtualinherited |
compute common variables later used in get_variational_expection and get_variational_first_derivative. Note that this method will automatically be called when set_variational_distribution is called
Definition at line 213 of file DualVariationalGaussianLikelihood.cpp.
|
inherited |
prints all parameter registered for model selection and their type
Definition at line 474 of file SGObject.cpp.
|
virtualinherited |
prints registered parameters out
prefix | prefix for members |
Definition at line 308 of file SGObject.cpp.
|
virtualinherited |
Save this object to file.
file | where to save the object; will be closed during returning if PREFIX is an empty string. |
prefix | prefix for members |
Definition at line 314 of file SGObject.cpp.
|
protectedvirtualinherited |
Can (optionally) be overridden to post-initialize some member variables which are not PARAMETER::ADD'ed. Make sure that at first the overridden method BASE_CLASS::SAVE_SERIALIZABLE_POST is called.
ShogunException | will be thrown if an error occurs. |
Reimplemented in CKernel.
Definition at line 436 of file SGObject.cpp.
|
protectedvirtualinherited |
Can (optionally) be overridden to pre-initialize some member variables which are not PARAMETER::ADD'ed. Make sure that at first the overridden method BASE_CLASS::SAVE_SERIALIZABLE_PRE is called.
ShogunException | will be thrown if an error occurs. |
Reimplemented in CKernel, CDynamicArray< T >, CDynamicArray< float64_t >, CDynamicArray< float32_t >, CDynamicArray< int32_t >, CDynamicArray< char >, CDynamicArray< bool >, and CDynamicObjectArray.
Definition at line 431 of file SGObject.cpp.
set dual parameters for variational parameters
the_lambda | dual parameter for variational mean |
lab | labels/data used |
Note that dual parameter (alpha) for the variational variance is implicitly set based on lambda
Definition at line 157 of file DualVariationalGaussianLikelihood.cpp.
|
inherited |
Definition at line 41 of file SGObject.cpp.
|
inherited |
Definition at line 46 of file SGObject.cpp.
|
inherited |
Definition at line 51 of file SGObject.cpp.
|
inherited |
Definition at line 56 of file SGObject.cpp.
|
inherited |
Definition at line 61 of file SGObject.cpp.
|
inherited |
Definition at line 66 of file SGObject.cpp.
|
inherited |
Definition at line 71 of file SGObject.cpp.
|
inherited |
Definition at line 76 of file SGObject.cpp.
|
inherited |
Definition at line 81 of file SGObject.cpp.
|
inherited |
Definition at line 86 of file SGObject.cpp.
|
inherited |
Definition at line 91 of file SGObject.cpp.
|
inherited |
Definition at line 96 of file SGObject.cpp.
|
inherited |
Definition at line 101 of file SGObject.cpp.
|
inherited |
Definition at line 106 of file SGObject.cpp.
|
inherited |
Definition at line 111 of file SGObject.cpp.
|
inherited |
set generic type to T
|
inherited |
|
inherited |
set the parallel object
parallel | parallel object to use |
Definition at line 241 of file SGObject.cpp.
|
inherited |
set the version object
version | version object to use |
Definition at line 283 of file SGObject.cpp.
|
protectedvirtualinherited |
this method used to set m_likelihood
Definition at line 49 of file VariationalLikelihood.cpp.
|
virtualinherited |
set a non-negative noise factor in order to correct the variance if variance is close to zero or negative setting 0 means correction is not applied
noise_factor | noise factor |
The default value is 1e-6.
Reimplemented from CVariationalGaussianLikelihood.
Definition at line 72 of file DualVariationalGaussianLikelihood.cpp.
|
virtualinherited |
set the m_strict_scale
strict_scale | must be between 0 and 1 exclusively |
Definition at line 103 of file DualVariationalGaussianLikelihood.cpp.
|
virtualinherited |
set the variational distribution given data and parameters
mu | mean of the variational distribution |
s2 | variance of the variational distribution |
lab | labels/data used |
Note that the variational distribution is Gaussian
Reimplemented from CVariationalGaussianLikelihood.
Definition at line 96 of file DualVariationalGaussianLikelihood.cpp.
|
virtualinherited |
A shallow copy. All the SGObject instance variables will be simply assigned and SG_REF-ed.
Reimplemented in CGaussianKernel.
Definition at line 192 of file SGObject.cpp.
|
virtualinherited |
return whether likelihood function supports binary classification
Reimplemented from CLikelihoodModel.
Definition at line 162 of file VariationalLikelihood.cpp.
|
virtualinherited |
return whether likelihood function supports computing the derivative wrt hyperparameter Note that variational parameters are NOT considered as hyperparameters
Implements CVariationalLikelihood.
Definition at line 84 of file DualVariationalGaussianLikelihood.cpp.
|
virtualinherited |
return whether likelihood function supports multiclass classification
Reimplemented from CLikelihoodModel.
Definition at line 168 of file VariationalLikelihood.cpp.
|
virtualinherited |
return whether likelihood function supports regression
Reimplemented from CLikelihoodModel.
Definition at line 156 of file VariationalLikelihood.cpp.
|
inherited |
unset generic type
this has to be called in classes specializing a template class
Definition at line 303 of file SGObject.cpp.
|
virtualinherited |
Updates the hash of current parameter combination
Definition at line 248 of file SGObject.cpp.
|
inherited |
io
Definition at line 369 of file SGObject.h.
|
inherited |
parameters wrt which we can compute gradients
Definition at line 384 of file SGObject.h.
|
inherited |
Hash of parameter values
Definition at line 387 of file SGObject.h.
|
protectedinherited |
whether m_lambda is satisfied lower bound and/or upper bound condition.
Definition at line 238 of file DualVariationalGaussianLikelihood.h.
the label of data
Definition at line 277 of file VariationalLikelihood.h.
The dual variables (lambda) for the variational parameter s2.
Note that in variational Gaussian inference, there is a relationship between lambda and alpha, where alpha is the dual parameter for variational parameter mu
Therefore, the dual variables (alpha) for variational parameter mu is not explicitly saved.
Definition at line 227 of file DualVariationalGaussianLikelihood.h.
|
protectedinherited |
the distribution used to model data
Definition at line 280 of file VariationalLikelihood.h.
|
inherited |
model selection parameters
Definition at line 381 of file SGObject.h.
The mean of variational Gaussian distribution
Definition at line 79 of file VariationalGaussianLikelihood.h.
|
inherited |
parameters
Definition at line 378 of file SGObject.h.
The variance of variational Gaussian distribution
Definition at line 82 of file VariationalGaussianLikelihood.h.
|
protectedinherited |
The value used to ensure strict bound(s) for m_lambda in adjust_step_wrt_dual_parameter()
Note that the value should be between 0 and 1 exclusively.
The default value is 1e-5.
Definition at line 235 of file DualVariationalGaussianLikelihood.h.
|
inherited |
parallel
Definition at line 372 of file SGObject.h.
|
inherited |
version
Definition at line 375 of file SGObject.h.