After integrating out the multivariate normal distribution , the marginal distribution is
Up to terms independent of , the log-likelihood is
Differentiating and setting the result to zero gives the maximum marginal likelihood estimator
This is an Empirical Bayes method because the estimated hyperparameter is then inserted into the prior and posterior distributions.
By the stationary distribution property, and have the same marginal distribution . Expanding the square gives
This is the probabilistic representation of the Dirichlet form of a Markov chain.