Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

We have got a Gaussian distribution covariance matrix(precision matrix) and the potential information, that is, if g is proportional to exp(-X'KX+h'X). However, K here is not positive definite. So we would like to have an approximation of this function.

Vectorizing K and h into one vector, we have want to project the current distribution to a distribution with positive definite K.

There is a conclusion that disregarding the h vector, the best projection will just be to change negative eigenvalues of K to be 0(with least squared sum for each elements on K). But the question is, will the h vector affect this result? Or we may just use the same h vector since they only manipulate the first order terms?

Thank you very much!

share|improve this question

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.