Learn R Programming

bayesm (version 1.1-2)

ghkvec: Compute GHK approximation to Multivariate Normal Integrals

Description

ghkvec computes the GHK approximation to the integral of a multivariate normal density over a half plane defined by a set of truncation points.

Usage

ghkvec(L, trunpt, above, r)

Arguments

L
lower triangular Cholesky root of Covariance matrix
trunpt
vector of truncation points
above
vector of indicators for truncation above(1) or below(0)
r
number of draws to use in GHK

Value

  • approximation to integral

concept

  • multivariate normal distribution
  • integral

References

For further discussion, see Bayesian Statistics and Marketing by Allenby, McCulloch, and Rossi, Chapter 3. http://gsbwww.uchicago.edu/fac/peter.rossi/research/bsm.html

Examples

Run this code
##

Sigma=matrix(c(1,.5,.5,1),ncol=2)
L=t(chol(Sigma))
trunpt=c(0,0,1,1)
above=c(1,1)
ghkvec(L,trunpt,above,100)

Run the code above in your browser using DataLab