Learn R Programming

bbricks (version 0.1.1)

MPE.GaussianNIG: MPE estimate of a "GaussianNIG" object

Description

Generate the MPE estimate of "theta" in following Gaussian-NIG structure: beta,sigma^2|gamma ~ NIG(gamma) x|beta,sigma^2,X ~ Gaussian(X where gamma = (m,V,a,b) is the Normal-Inverse-Gamma(NIG) parameter, "m" is a numeric "location" parameter; "V" is a symmetric positive definite matrix representing the "scale" parameters; "a" and "b" are the "shape" and "rate" parameter of the Inverse Gamma distribution. The model structure and prior parameters are stored in a "GaussianNIG" object. MPE is beta,sigma^2_MPE = E(beta,sigma^2|gamma,x,X), E() is the expectation function.

Usage

# S3 method for GaussianNIG
MPE(obj, ...)

Arguments

obj

A "GaussianNIG" object.

...

Additional arguments to be passed to other inherited types.

Value

A named list, the MPE estimate of beta and sigma^2.

References

Banerjee, Sudipto. "Bayesian Linear Model: Gory Details." Dowloaded from http://www. biostat. umn. edu/~ ph7440 (2008).

See Also

GaussianNIG