For a multivariate normal vector $\mathbf{x} = (\mathbf{x}_1, \mathbf{x}_2, \mathbf{x}_3, \mathbf{x}_4)$ with mean 0 and variance $\mathbf{\Sigma} = \begin{bmatrix} \mathbf{\Sigma}_{11} & \mathbf{\Sigma}_{12} & \mathbf{\Sigma}_{13} & \mathbf{\Sigma}_{14} \\ \mathbf{\Sigma}_{21} & \mathbf{\Sigma}_{22} & \mathbf{\Sigma}_{23} & \mathbf{\Sigma}_{24} \\ \mathbf{\Sigma}_{31} & \mathbf{\Sigma}_{32} & \mathbf{\Sigma}_{33} & \mathbf{\Sigma}_{34} \\ \mathbf{\Sigma}_{41} & \mathbf{\Sigma}_{42} & \mathbf{\Sigma}_{43} & \mathbf{\Sigma}_{44} \end{bmatrix}$, is it true that
- $(\mathbf{x}_1, \mathbf{x}_3)$ is normally distributed with mean 0 and variance $\begin{bmatrix} \mathbf{\Sigma}_{11} & \mathbf{\Sigma}_{13} \\ \mathbf{\Sigma}_{31} & \mathbf{\Sigma}_{33} \end{bmatrix}$?
- $\mathbf{x}_3 | \mathbf{x}_1$ is normally distributed with mean 0 and variance $\mathbf{\Sigma}_{33} - \mathbf{\Sigma}_{31} \mathbf{\Sigma}_{11}^{-1} \mathbf{\Sigma}_{13}$?
I believe (1) is true because of the accepted answer at Deriving the conditional distributions of a multivariate normal distribution. I think (2) follows from (1), if (1) is true. However, I wonder why (1) is true - I guess (1) could be verified mathematically by integrating out $\mathbf{x}_2, \mathbf{x}_4$?