25 using namespace Eigen;
61 "Exact inference method can only use Gaussian likelihood function\n")
63 "Labels must be type of CRegressionLabels\n")
105 float64_t result=(eigen_y-eigen_m).
dot(eigen_alpha)/2.0+
106 eigen_L.diagonal().array().log().sum()+
m_L.
num_rows*
183 a=L.triangularView<Upper>().adjoint().solve(eigen_y-eigen_m);
184 a=L.triangularView<Upper>().solve(a);
217 MatrixXd eigen_V=eigen_L.triangularView<Upper>().adjoint().solve(
223 eigen_V = eigen_V/sigma;
244 eigen_Q=eigen_L.triangularView<Upper>().adjoint().solve(
246 eigen_Q=eigen_L.triangularView<Upper>().solve(eigen_Q);
252 eigen_Q-=eigen_alpha*eigen_alpha.transpose();
258 REQUIRE(!strcmp(param->
m_name,
"scale"),
"Can't compute derivative of "
259 "the nagative log marginal likelihood wrt %s.%s parameter\n",
268 result[0]=(eigen_Q.cwiseProduct(eigen_K)*
m_scale*2.0).
sum()/2.0;
276 REQUIRE(!strcmp(param->
m_name,
"sigma"),
"Can't compute derivative of "
277 "the nagative log marginal likelihood wrt %s.%s parameter\n",
292 result[0]=
CMath::sq(sigma)*eigen_Q.trace();
309 "Length of the parameter %s should not be NULL\n", param->
m_name)
347 "Length of the parameter %s should not be NULL\n", param->
m_name)
368 result[i]=-eigen_dmu.dot(eigen_alpha);
virtual const char * get_name() const =0
virtual void update_deriv()
virtual ELabelType get_label_type() const =0
virtual void update_parameter_hash()
Class that models Gaussian likelihood.
Real Labels are real-valued labels.
SGVector< float64_t > m_alpha
The Inference Method base class.
virtual SGMatrix< float64_t > get_posterior_covariance()
Vector::Scalar dot(Vector a, Vector b)
The class Labels models labels, i.e. class assignments of objects.
real valued labels (e.g. for regression, classifier outputs)
virtual const char * get_name() const
virtual SGVector< float64_t > get_alpha()
static T sum(T *vec, int32_t len)
Return sum(vec)
virtual ELikelihoodModelType get_model_type() const
virtual int32_t get_num_vectors() const =0
virtual void update_mean()
virtual SGVector< float64_t > get_mean_vector(const CFeatures *features) const =0
An abstract class of the mean function.
SGMatrix< float64_t > m_L
virtual void update_chol()
virtual ~CExactInferenceMethod()
virtual SGMatrix< float64_t > get_cholesky()
virtual SGVector< float64_t > get_derivative_wrt_mean(const TParameter *param)
Matrix::Scalar sum(Matrix m, bool no_diag=false)
virtual SGVector< float64_t > get_derivative_wrt_inference_method(const TParameter *param)
virtual SGVector< float64_t > get_diagonal_vector()
static void fill_vector(T *vec, int32_t len, T value)
static CGaussianLikelihood * obtain_from_generic(CLikelihoodModel *lik)
virtual SGVector< float64_t > get_parameter_derivative(const CFeatures *features, const TParameter *param, index_t index=-1)
all of classes and functions are contained in the shogun namespace
virtual float64_t get_negative_log_marginal_likelihood()
The class Features is the base class of all feature objects.
virtual SGVector< float64_t > get_derivative_wrt_kernel(const TParameter *param)
virtual SGMatrix< float64_t > get_parameter_gradient(const TParameter *param, index_t index=-1)
static float64_t log(float64_t v)
virtual void check_members() const
virtual SGVector< float64_t > get_derivative_wrt_likelihood_model(const TParameter *param)
virtual void update_alpha()
virtual void update_cov()
virtual void check_members() const
virtual bool parameter_hash_changed()
The Likelihood model base class.
SGMatrix< float64_t > m_ktrtr
CLikelihoodModel * m_model
virtual SGVector< float64_t > get_posterior_mean()
static const float64_t PI