Learn R Programming

bayesm (version 2.0-9)

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
  • GHK method
  • integral

References

For further discussion, see Bayesian Statistics and Marketing by Rossi,Allenby and McCulloch, Chapter 2. http://faculty.chicagogsb.edu/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