DOI QR코드

DOI QR Code

An importance sampling for a function of a multivariate random variable

  • Jae-Yeol Park (Department of Statistics, University of Seoul) ;
  • Hee-Geon Kang (Department of Statistics, University of Seoul) ;
  • Sunggon Kim (Department of Statistics, University of Seoul)
  • Received : 2023.09.11
  • Accepted : 2023.10.24
  • Published : 2024.01.31

Abstract

The tail probability of a function of a multivariate random variable is not easy to estimate by the crude Monte Carlo simulation. When the occurrence of the function value over a threshold is rare, the accurate estimation of the corresponding probability requires a huge number of samples. When the explicit form of the cumulative distribution function of each component of the variable is known, the inverse transform likelihood ratio method is directly applicable scheme to estimate the tail probability efficiently. The method is a type of the importance sampling and its efficiency depends on the selection of the importance sampling distribution. When the cumulative distribution of the multivariate random variable is represented by a copula and its marginal distributions, we develop an iterative algorithm to find the optimal importance sampling distribution, and show the convergence of the algorithm. The performance of the proposed scheme is compared with the crude Monte Carlo simulation numerically.

Keywords

Acknowledgement

The authors would like to thank the anonymous reviewers for their comments and suggestions on the first draft of this paper. Their suggestions have greatly improve the quality of the paper. This work was supported by the 2021 sabbatical year research grant of the University of Seoul.

References

  1. Asmussen S, Binswanger K, and Hojgaard B (2000). Rare events simulation for heavy-tailed distributions, Bernoulli, 6, 303-322.  https://doi.org/10.2307/3318578
  2. Asmussen S and Kroese DP (2006). Improved algorithms for rare event simulation with heavy tails, Advances in Applied Probability, 38, 545-558.  https://doi.org/10.1239/aap/1151337084
  3. Asmussen S, Blanchet J, Juneja S, and Rojas-Nandayapa L (2011). Efficient simulation of tail probabilities of sums of correlated lognormals, Annals of Operations Research, 189, 5-23.  https://doi.org/10.1007/s10479-009-0658-5
  4. Blanchet J and Rojas-Nandayapa L (2011). Efficient simulation of tail probabilities of sums of dependent random variables, Journal of Applied Probability, 48, 147-164.  https://doi.org/10.1239/jap/1318940462
  5. De Boer P-T, Kroese DP, Mannor S, and Rubinstein RY (2005). A tutorial on the cross-entropy method, Annals of Operations Research, 134, 19-67.  https://doi.org/10.1007/s10479-005-5724-z
  6. Dupuis P, Leder K, and Wang H (2007). Importance sampling for sums of random variables with regularly varying tails, ACM Transactions on Modeling and Computer Simulation (TOMACS), 17, 14-es. 
  7. Glynn PW and Iglehart DL (1989). Importance sampling for stochastic simulations, Management Science, 35, 1367-1392.  https://doi.org/10.1287/mnsc.35.11.1367
  8. Glynn PW and Whitt W (1992). The asymptotic efficiency of simulation estimators, Operations Research, 40, 505-520.  https://doi.org/10.1287/opre.40.3.505
  9. Homem-de-Mello T and Rubinstein RY (2002). Estimation of rare event probabilities using cross-entropy. In Proceedings of the Winter Simulation Conference, San Diego, CA, USA, 310-319. 
  10. Juneja S and Shahabuddin P (2002). Simulating heavy tailed processes using delayed hazard rate twisting, ACM Transactions on Modeling and Computer Simulation (TOMACS), 12, 94-118.  https://doi.org/10.1145/566392.566394
  11. Juneja S and Shahabuddin P (2006). Rare-event simulation techniques: An introduction and recent advances, Handbooks in Operations Research and Management Science, 13, 291-350.  https://doi.org/10.1016/S0927-0507(06)13011-X
  12. Kroese DP and Rubinstein RY (2004). The transform likelihood ratio method for rare event simulation with heavy tails, Queueing Systems, 46, 317-351.  https://doi.org/10.1023/B:QUES.0000027989.97672.be
  13. Lieber D (1998). Rare-event estimation via cross-entropy and importance sampling (Doctoral dissertation), Technion-Israel Institute of Technology, Haifa.
  14. Nelsen RB (2006). An Introduction to Copulas, Springer, New York. 
  15. Rached NB, Benkhelifa F, Alouini M-S, and Tempone R (2015). A fast simulation method for the log-normal sum distribution using a hazard rate twisting technique. In Proceedings of 2015 IEEE International Conference on Communications (ICC), London, 4259-4264. 
  16. Rached NB, Benkhelifa F, Kammoun A, Alouini M-S, and Tempone R (2018). On the generalization of the hazard rate twisting-based simulation approach, Statistics and Computing, 28, 61-75.  https://doi.org/10.1007/s11222-016-9716-4
  17. Rubinstein RY (1999). The cross-entropy method for combinatorial and continuous optimization, Methodology and Computing in Applied Probability, 1, 127-190.  https://doi.org/10.1023/A:1010091220143
  18. Rubinstein RY (2002). Cross-entropy and rare events for maximal cut and partition problems, ACM Transactions on Modeling and Computer Simulation (TOMACS), 12, 27-53.  https://doi.org/10.1145/511442.511444
  19. Rubinstein RY and Kroese DP (2016). Simulation and the Monte Carlo Method, John Wiley & Sons, Hoboken, New Jersey. 
  20. Rubinstein RY and Shapiro A (1993). Discrete Event Systems: Sensitivity Analysis and Stochastic Optimization by the Score Function Method, Wiley, Chichester, England. 
  21. Sadowsky JS (1993). On the optimality and stability of exponential twisting in Monte Carlo estimation, IEEE Transactions on Information Theory, 39, 119-128.  https://doi.org/10.1109/18.179349
  22. Sadowsky JS and Bucklew JA (1990). On large deviations theory and asymptotically efficient Monte Carlo estimation, IEEE Transactions on Information Theory, 36, 579-588.  https://doi.org/10.1109/18.54903
  23. Sak H and Hormann W (2012). Fast simulations in credit risk, Quantitative Finance, 12, 1557-1569.  https://doi.org/10.1080/14697688.2011.564199