Using properties of the multivariate Normal distribution, the conditional expectation and covariance matrix of zˆN given zˆC are readily available E[zˆN|zˆC]=RNCRCC−1zˆC V[zˆN|zˆC]=RNN−RNCRCC−1RCN and do not depend on Σγ. We rewrite the marginal likelihood p(y|γ,X)=N(zˆ|0,R+RΣγR) in terms of the marginal distribution of zˆC and conditional distribution of zˆN given zˆC to obtain the following expression