Computes the conditional distribution of a subset of parameters given
observed values for other parameters. Uses the closed-form Schur complement
for the multivariate normal. Returns a dist object (not an mle),
since conditioning loses the MLE context.
Usage
# S3 method for class 'mle'
conditional(x, P = NULL, ..., given_indices = NULL, given_values = NULL)