I have a posterior distribution over a set of parameters denoted by $\theta = \{w, \phi, \lambda\}$ and the posterior is the joint distribution given the observed data i.e. $P(\theta | D)$. Now, this posterior is a multivariate normal distribution (due to the approximate way it is estimated).
Now I want to find out the marginal distributions from this joint. So, for example, $P(\lambda)$ could be found my marginalizing over $w$ and $\phi$. Is there an easy way to do this for such MVN distributions? Also, can I find the mode for this marginal distribution easily?