We consider the probabilistic linear model in which all random variables are normal.
A normal linear model is a probabilistic linear model in which the parameter and noise vectors have normal (Gaussian) densities. The model is also called the Gaussian linear model or the linear model with Gaussian noise.
Let $(x: \Omega \to \R ^d, A \in \R ^{n
\times d}, e: \Omega \to \R ^n)$ be a
probabilistic linear model over the probability
space $(\Omega , \mathcal{A} , \mathbfsf{P} )$ in
which $x$ and $e$ have normal densities.
Recall that a probabilistic linear model has
observation vector $y: \Omega \to \R ^n$ defined
by
\[
y = Ax + e.
\]
Since $x$ and $e$ are normal and independent,
$y$ is normal.1
Moreover, the random vector $(x, y)$ is normal
with covariance
\[
\pmat{
\Sigma _x & \Sigma _{x}A^\top \\
A \Sigma _{x} & A\Sigma _{x}A^\top + \Sigma _{e}
}.
\] \[
\Sigma _{x}A^\top (A\Sigma _{x}A^\top + \Sigma _e)^{-1}\gamma
\] \[
\Sigma _{x} - \Sigma _{x}A^\top (A\Sigma _{x}A^\top +
\Sigma _e)^{-1}A\Sigma _{x}.
\]
We can write the conditional mean as
\[
(\Sigma _{x}^{-1} + A^\top \Sigma _{e}^{-1}A)^{-1} A^\top
\Sigma _{e}^{-1}
\] \[
(\Sigma _x^{-1} + A^\top \Sigma _{e}^{-1} A)^{-1}.
\]
\[ \Sigma _{x} A^\top (A\Sigma _{x}A^\top + \Sigma _{e})^{-1}\gamma . \]