Browse > Article
http://dx.doi.org/10.29220/CSAM.2021.28.6.673

trunmnt: An R package for calculating moments in a truncated multivariate normal distribution  

Lee, Seung-Chun (Department of Applied Statistics, Hanshin University)
Publication Information
Communications for Statistical Applications and Methods / v.28, no.6, 2021 , pp. 673-679 More about this Journal
Abstract
The moment calculation in a truncated multivariate normal distribution is a long-standing problem in statistical computation. Recently, Kan and Robotti (2017) developed an algorithm able to calculate all orders of moment under different types of truncation. This result was implemented in an R package MomTrunc by Galarza et al. (2021); however, it is difficult to use the package in practical statistical problems because the computational burden increases exponentially as the order of the moment or the dimension of the random vector increases. Meanwhile, Lee (2021) presented an efficient numerical method in both accuracy and computational burden using Gauss-Hermit quadrature. This article introduces trunmnt implementation of Lee's work as an R package. The Package is believed to be useful for moment calculations in most practical statistical problems.
Keywords
moment of truncated multivariate normal; MomTrunc; Gauss-Hermite quadrature;
Citations & Related Records
연도 인용수 순위
  • Reference
1 Genz A (1993). Comparison of methods for the computation of multivariate normal probabilities, Computing Science and Statistics, 25, 400-405.
2 Genz et al. (2021). mvtvnorm: Multivariate Normal and t Distributions, R package version 1.1, Available from: https://cran.r-project.org/web/packages/mvtnorm/mvtnorm.pdf
3 Kan R and Robotti C (2017). On moments of folded and truncated multivariate normal distributions, Journal of Computational and Graphical Statistics, 26, 930-934.   DOI
4 Lee SC (2021). moments calculation for truncated multivariate normal in nonlinear generalized mixed models, Communications for Statistical Applications and Methods, 27, 377-383.   DOI
5 Manjunath BG and Wilhelm S (2012). Moments calculation for the doubly truncated multivariate normal density, SSRN Electronic Journal, arXiv preprint arXiv:1206.5387.
6 R Core Team (2019). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/.
7 Wilhelm S and Manjunath BG (2010). tmvtnorm: Truncated multivariate normal and Student t distribution, The R Journal, 2.
8 Burkardt J (2014). The truncated normal distribution, Online document, Available from: http://people.sc.fsu.edu/jburkardt/presentations/truncatednormal.pdf
9 Jackel P (2005). A note on multivariate Gauss-Hermite quadrature, London: ABN-Amro, Available from: http://www.pjaeckel.webspace.virginmedia.com/ANoteOnMultivariate GaussHermiteQuadrature.pdf
10 Genz A (1992). Numerical computation of multivariate normal probabilities, Journal of Computational and Graphical Statistics, 1, 141-150.   DOI
11 Galarza CE, Kan R, and Lachos VH (2021). MomTrunc: Moments of folded and doubly truncated multivariate distributions, Online document, Available from: https://cran.r-project.org/web/packages/MomTrunc/MomTrunc.pdf
12 Eddelbuettel D and Sanderson C (2014). RcppArmadillo: Accelerating R with high-performance C++ linear algebra, Computational Statistics & Data Analysis, 71, 1054-1063.   DOI