Hello, I have a Gaussian mixture model with 3 components. The output of my model, of course: allocations, mixing probabilities, component-specific means and variances. In post-processing, I am trying to compute the highest posterior density for my model. How can we do that?