mnormalexpgamma {lestat}R Documentation

A Multivariate Normal-ExpGamma Distribution

Description

Creates an object representing a multivariate Normal-ExpGamma distribution. If (x,y)(x,y) has a multivariate Normal-ExpGamma distribution, then the marginal distribution of yy is an ExpGamma distribution, and the conditional distribution of xx given yy is multivariate normal.

Usage

mnormalexpgamma(mu=c(0,0), P, alpha, beta)

Arguments

mu

The mu parameter. It must be a vector of length at least 2. The default value is (0,0).

P

The P parameter.

alpha

The alpha parameter.

beta

The beta parameter.

Details

If (x,y)(x,y) has a multivariate Normal-ExpGamma distribution with parameters μ\mu, PP, α\alpha, and β\beta, then the marginal distribution of yy has an ExpGamma distribution with parameters α\alpha, β\beta, and -2, and conditionally on yy, xx has a multivariate normal distribution with expectation μ\mu and precision matrix e2yPe^{-2y}P. The probability density is proportional to

f(x,y)=exp((2α+k)ye2y(β+(xμ)tP(xμ)/2)) f(x,y)=\exp(-(2\alpha + k)y - e^{-2y}(\beta + (x-\mu)^tP(x-\mu)/2))

where kk is the dimension.

Value

A multivariate Normal-ExpGamma probability distribution.

Author(s)

Petter Mostad <mostad@chalmers.se>

See Also

gamma,normal,expgamma, normalgamma,normalexpgamma mnormal,mnormalgamma

Examples

plot(mnormalexpgamma(alpha=3, beta=3))

[Package lestat version 1.9 Index]