Communications for Statistical Applications and Methods

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

DOI QR Code

Moments calculation for truncated multivariate normal in nonlinear generalized mixed models

  • Received : 2020.02.10
  • Accepted : 2020.03.07
  • Published : 2020.05.31
Download PDF
⟨ Previous Next ⟩

Abstract

The likelihood-based inference in a nonlinear generalized mixed model often requires computing moments of truncated multivariate normal random variables. Many methods have been proposed for the computation using a recurrence relation or the moment generating function; however, these methods rely on high dimensional numerical integrations. The numerical method is known to be inefficient for high dimensional integral in accuracy. Besides the accuracy, the methods demand too much computing time to use them in practical analyses. In this note, a moment calculation method is proposed under an assumption of a certain covariance structure that occurred mostly in generalized mixed models. The method needs only low dimensional numerical integrations.

Keywords

References

  1. Amemiya T (1974). Multivariate regression and simultaneous equation models when the dependent variables are truncated normal, Econometrica, 42, 999-1012. https://doi.org/10.2307/1914214
  2. Arismendi JC (2013). Multivariate truncated moments, Journal of Multivariate Analysis, 117, 41-75. https://doi.org/10.1016/j.jmva.2013.01.007
  3. Begier MH and Hamdan MA (1971). Correlation in a bivariate normal distribution with truncation in both variables, Australian Journal of Statistics, 13, 77-82. https://doi.org/10.1111/j.1467-842X.1971.tb01245.x
  4. Burkardt J (2014). The truncated normal distribution, Online Document. Available from: http://people.sc.fsu.edu/jburkardt/presentations/truncated normal.pdf
  5. Croissant Y (2017). Package pglm: Panel generalized linear models, R package version 0.2-1. Available from: https://cran.r-project.org/web/packages/pglm/pglm.pdf
  6. Eddelbuettel D, Francois R, Allaire J, Ushey K, Kou Q, Russell N, Bates D, and Chambers J (2018). Seamless R and C++ Integration. Available from: http://www.rcpp.org, http://dirk.eddelbuettel.com/code/rcpp.html
  7. Genz A (1992). Numerical computation of multivariate normal probabilities, Journal of Computational and Graphical Statistics, 1, 141-150. https://doi.org/10.2307/1390838
  8. Genz A (1993). Comparison of methods for the computation of multivariate normal probabilities, Computing Science and Statistics, 25, 400-405.
  9. Genz A, Bretz F, Miwa T, Mi X, Leisch F, Scheipl F, Bornkamp B, Maechler M, and Hothorn T (2020). mvtvnorm: Multivariate Normal and t Distributions, R package version 1.0-12. Available from: https://cran.r-project.org/web/packages/mvtnorm/mvtnorm.pdf
  10. 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
  11. Manjunath BG and Wilhelm S (2012). Moments calculation for the doubly truncated multivariate normal density, arXiv preprint arXiv:1206.5387.
  12. Lee LF (1979). On the first and second moments of the truncated multi-normal distribution and a simple estimator, Economics Letters, 3, 165-169. https://doi.org/10.1016/0165-1765(79)90111-3
  13. Lee M and Lee SC (2019). Restricted maximum likelihood estimation of a censored random effects panel regression model, Communications for Statistical Applications and Methods, 26, 371-383. https://doi.org/10.29220/CSAM.2019.26.4.371
  14. McCulloch CE (1994). Maximum likelihood variance components estimation for binary data, Journal of the American Statistical Association, 89, 330-335. https://doi.org/10.1080/01621459.1994.10476474
  15. Muthen B (1990). Moments of the censored and truncated bivariate normal distribution, British Jour-nal of Mathematical and Statistical Psychology, 43, 131-143. https://doi.org/10.1111/j.2044-8317.1990.tb00930.x
  16. 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/.
  17. Rosenbaum S (1961). Moments of a truncated bivariate normal distribution, Journal of the Royal Statistical Society: Series B (Methodological), 23, 405-408. https://doi.org/10.1111/j.2517-6161.1961.tb00422.x
  18. Sarrias M (2016). Discrete choice models with random parameters in R: The Rchoice Package, Journal of Statistical Software, 74, 1-31. https://doi.org/10.18637/jss.v074.i10
  19. Shah SM and Parikh NT (1964). Moments of singly and doubly truncated standard bivariate normal distribution, Vidya, 7, 81-91.
  20. SAS (2015). SAS/STAT 14.1 Users Guide: The NLMIXED Procedure.
  21. Tallis GM (1961). The moment generating function of the truncated multi-normal distribution, Journal of the Royal Statistical Society: Series B (Methodological), 23, 223-229. https://doi.org/10.1111/j.2517-6161.1961.tb00408.x
  22. Wilhelm S and Manjunath BG (2010). tmvtnorm: Truncated multivariate normal and Student t distribution, The R Journal, 2.