Expectation of Univariate Quadratic Form under Multivariate Gaussian

459 Views Asked by At

Is there an obvious trick I am missing for solving the following integral:

$$ \int_x P(y|x) W(x) (-x^TMx+2x^Tm -c)dx$$

Distributions are Gaussians and $M$ is symmetric.

I know how to do the expectation for a univariate Gaussian but I'm not sure of this multivariate case since $$ P(y|x) W(x) = V(y,x) $$

$W$ is the PDF of $x$, $P$ is a conditional PDF of $y$ given $x$ and $V$ is the joint PDF of $x$ and $y$.

2

There are 2 best solutions below

0
On BEST ANSWER

So it seems there is a very easy trick to solve this. We only need to use the Bayes rule to get the marginal distribution $P(y)$: $$ P(y)W(x|y) = P(y|x)W(x)$$

$P(y)$ can then be taken outside the integral and we end up with an expectation over a univariate conditional Gaussian $W(x|y)$ which can be calculated for example as stated in the matrix cookbook.

0
On

So, you have a multinormal vector $(y,x)$ (where both components $y, x$ can be vectors), and want the expectation of the quadratic form $-x^T M X +2x^T m - c$, $M$ an symmetric constant matrix, $m$ a constant vector and $c$ a constant scalar. Since this function (random variable) you want the expectation of, do not depend on $y$, the distribution of $y$ is irrelevant for the expectation. For the first part, you can just apply my answer to sum of squares of dependent gaussian random variables Then $\text{E} x^T m$ is simply $\text{E}(x)^T m$.