Communications for Statistical Applications and Methods

The Korean Statistical Society (한국통계학회)
권호 표지
DOI QR코드

DOI QR Code

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

  • Received : 2021.08.02
  • Accepted : 2021.10.18
  • Published : 2021.11.30
Download PDF
⟨ Previous Next ⟩

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

Acknowledgement

This work was supported by Hanshin University research grant.

References

  1. Burkardt J (2014). The truncated normal distribution, Online document, Available from: http://people.sc.fsu.edu/jburkardt/presentations/truncatednormal.pdf
  2. Eddelbuettel D and Sanderson C (2014). RcppArmadillo: Accelerating R with high-performance C++ linear algebra, Computational Statistics & Data Analysis, 71, 1054-1063. https://doi.org/10.1016/j.csda.2013.02.005
  3. 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
  4. Genz A (1992). Numerical computation of multivariate normal probabilities, Journal of Computational and Graphical Statistics, 1, 141-150. https://doi.org/10.2307/1390838
  5. Genz A (1993). Comparison of methods for the computation of multivariate normal probabilities, Computing Science and Statistics, 25, 400-405.
  6. 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
  7. Jackel P (2005). A note on multivariate Gauss-Hermite quadrature, London: ABN-Amro, Available from: http://www.pjaeckel.webspace.virginmedia.com/ANoteOnMultivariate GaussHermiteQuadrature.pdf
  8. Kan R and Robotti C (2017). On moments of folded and truncated multivariate normal distributions, Journal of Computational and Graphical Statistics, 26, 930-934. https://doi.org/10.1080/10618600.2017.1322092
  9. Lee SC (2021). moments calculation for truncated multivariate normal in nonlinear generalized mixed models, Communications for Statistical Applications and Methods, 27, 377-383. https://doi.org/10.29220/CSAM.2020.27.3.377
  10. Manjunath BG and Wilhelm S (2012). Moments calculation for the doubly truncated multivariate normal density, SSRN Electronic Journal, arXiv preprint arXiv:1206.5387.
  11. 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/.
  12. Wilhelm S and Manjunath BG (2010). tmvtnorm: Truncated multivariate normal and Student t distribution, The R Journal, 2.