= Solution
Let $K\in\mathbb R^{n\times n}$ have entries $K_{ij}=\kappa(x_i,x_j)$, let $k_*\in\mathbb R^n$ have entries $(k_*)_i=\kappa(x_i,x_*)$, and let $k_{**}=\kappa(x_*,x_*)$. The <Gaussian process> prior and independent <Gaussian noise> imply
$$
\begin{pmatrix}f(x_*)\\y\end{pmatrix}
\sim N\!\left[
\begin{pmatrix}0\\0_n\end{pmatrix},
\begin{pmatrix}
k_{**}&k_*^T\\
k_*&K+\sigma^2I_n
\end{pmatrix}
\right].
$$
Applying the <conditional multivariate normal distribution> gives the <Gaussian process regression posterior>
$$
f(x_*)\mid y,X,x_*
\sim N(m_*,v_*),
$$
where
$$
m_*=k_*^T(K+\sigma^2I_n)^{-1}y,
\qquad
v_*=k_{**}-k_*^T(K+\sigma^2I_n)^{-1}k_*.
$$
Back to article page
