ON THE CHOICE OF ALTERNATIVE MEASURES IN IMPORTANCE SAMPLING WITH MARKOV-CHAINS

被引:14
作者
ANDRADOTTIR, S [1 ]
HEYMAN, DP [1 ]
OTT, TJ [1 ]
机构
[1] BELL COMMUN RES INC,MORRISTOWN,NJ
关键词
D O I
10.1287/opre.43.3.509
中图分类号
C93 [管理学];
学科分类号
12 ; 1201 ; 1202 ; 120202 ;
摘要
In the simulation of Markov chains, importance sampling involves replacing the original transition matrix, say P, with a suitably chosen transition matrix Q that tends to visit the states of interest more frequently. The likelihood ratio of P relative to Q is an important random variable in the importance sampling method. It always has expectation one, and for any interesting pair of transition matrices P and Q, there is a sample path length that causes the likelihood ratio to be close to zero with probability close to one. This may cause the variance of the importance sampling estimator to be larger than the variance of the traditional estimator. We develop sufficient conditions for ensuring the tightness of the distribution of the logarithm of the likelihood ratio for all sample path lengths, and we show that when these conditions are satisfied, the likelihood ratio is approximately lognormally distributed with expected value one. These conditions can be used to eliminate some choices of the alternative transition matrix Q that are likely to result in a variance increase. We also show that if the likelihood ratio is to remain well behaved for all sample path lengths, the alternative transition matrix Q has to approach the original transition matrix P as the sample path length increases. The practical significance of this result is that importance sampling can be difficult to apply successfully in simulations that involve long sample paths.
引用
收藏
页码:509 / 519
页数:11
相关论文
共 11 条
[1]   MONTE-CARLO SIMULATION AND LARGE DEVIATIONS THEORY FOR UNIFORMLY RECURRENT MARKOV-CHAINS [J].
BUCKLEW, JA ;
NEY, P ;
SADOWSKY, JS .
JOURNAL OF APPLIED PROBABILITY, 1990, 27 (01) :44-59
[2]  
GLYYN PW, 1989, MANAGE SCI, V35, P1367
[3]   Statistical Analysis and Simulation Study of Video Teleconference Traffic in ATM Networks [J].
Heyman, Daniel P. ;
Tabatabai, Ali ;
Lakshman, T. V. .
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS FOR VIDEO TECHNOLOGY, 1992, 2 (01) :49-59
[4]   PROCESS WITH CHAIN DEPENDENT GROWTH RATE [J].
KEILSON, J ;
SUBBARAO, S .
JOURNAL OF APPLIED PROBABILITY, 1970, 7 (03) :699-&
[5]  
KEILSON J, 1964, P CAMB PHILOS SOC, V60, P547
[6]  
KLEIJNEN JPC, 1978, MANAGE SCI, V24, P1772
[7]   LIKELIHOOD RATIO SENSITIVITY ANALYSIS FOR MARKOVIAN MODELS OF HIGHLY DEPENDABLE SYSTEMS [J].
NAKAYAMA, MK ;
GOYAL, A ;
GLYNN, PW .
OPERATIONS RESEARCH, 1994, 42 (01) :137-157
[8]  
OTT TJ, 1995, DOUBLE SEQUENCE CENT
[9]   A QUICK SIMULATION METHOD FOR EXCESSIVE BACKLOGS IN NETWORKS OF QUEUES [J].
PAREKH, S ;
WALRAND, J .
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, 1989, 34 (01) :54-66
[10]  
SEIGMUND D, 1976, ANN STAT, V4, P673