For a joint Gaussian distribution \(X \sim N(\mu, \Sigma)\), partitioned as
\(X = (X_A, X_B)\), the conditional distribution is:
$$X_B | X_A = x_A \sim N(\mu_{B|A}, \Sigma_{B|A})$$
where:
$$\mu_{B|A} = \mu_B + \Sigma_{BA} \Sigma_{AA}^{-1} (x_A - \mu_A)$$
$$\Sigma_{B|A} = \Sigma_{BB} - \Sigma_{BA} \Sigma_{AA}^{-1} \Sigma_{AB}$$
This is equivalent to the regression formulation used by fippy:
$$\beta = \Sigma_{BA} \Sigma_{AA}^{-1}$$
$$\mu_{B|A} = \mu_B + \beta (x_A - \mu_A)$$
$$\Sigma_{B|A} = \Sigma_{BB} - \beta \Sigma_{AB}$$
Assumptions:
Advantages:
Limitations:
Strong distributional assumption
May produce out-of-range values for bounded features
Cannot handle categorical features
Integer features are treated as continuous and rounded back to integers