Generate the MPE of "theta" in following GaussianNIW structure: theta|gamma ~ NIW(gamma) x|theta ~ Gaussian(theta) where theta = (mu,Sigma) is the Gaussian parameter, gamma = (m,k,v,S) is the Normal-Inverse-Wishart(NIW) parameter. The model structure and prior parameters are stored in a "GaussianNIW" object. MPE is theta_MPE = mean_theta(theta|gamma,x)
# S3 method for GaussianNIW
MPE(obj, ...)A "GaussianNIW" object.
Additional arguments to be passed to other inherited types.
A named list, the MPE estimate of "theta".
Murphy, Kevin P. "Conjugate Bayesian analysis of the Gaussian distribution." def 1.22 (2007): 16.
Gelman, Andrew, et al. "Bayesian Data Analysis Chapman & Hall." CRC Texts in Statistical Science (2004).