Likelihood Ratio Gradient Estimation for Stochastic Systems
P. W. Glynn
Communications of the ACM, Vol. 33, 75-84 (1990)
Consider a computer system having a CPU that feeds jobs to two input/output (I/O) devices having different speeds. Let θ be the fraction of jobs routed to the first I/O device, so that 1 - θ is the fraction routed to the second. Suppose that α = α(θ) is the steady-sate amount of time that a job spends in the system. Given that θ is a decision variable, a designer might wish to minimize α(θ) over θ. Since α(·) is typically difficult to evaluate analytically, Monte Carlo optimization is an attractive methodology. By analogy with deterministic mathematical programming, efficient Monte Carlo gradient estimation is an important ingredient of simulation-based optimization algorithms. As a consequence, gradient estimation has recently attracted considerable attention in the simulation community. It is our goal, in this article, to describe one efficient method for estimating gradients in the Monte Carlo setting, namely the likelihood ratio method (also known as the efficient score method). This technique has been previously described (in less general settings than those developed in this article) in [6, 16, 18, 21]. An alternative gradient estimation procedure is infinitesimal perturbation analysis; see [11, 12] for an introduction. While it is typically more difficult to apply to a given application than the likelihood ratio technique of interest here, it often turns out to be statistically more accurate. In this article, we first describe two important problems which motivate our study of efficient gradient estimation algorithms. Next, we will present the likelihood ratio gradient estimator in a general setting in which the essential idea is most transparent. The section that follows then specializes the estimator to discrete-time stochastic processes. We derive likelihood-ratio-gradient estimators for both time-homogeneous and non-time homogeneous discrete-time Markov chains. Later, we discuss likelihood ratio gradient estimation in continuous time. As examples of our analysis, we present the gradient estimators for time-homogeneous continuous-time Markov chains; non-time homogeneous continuous-time Markov chains; semi-Markov processes; and generalized semi-Markov processes. (The analysis throughout these sections assumes the performance measure that defines α(θ) corresponds to a terminating simulation.) Finally, we conclude the article with a brief discussion of the basic issues that arise in extending the likelihood ratio gradient estimator to steady-state performance measures.
 Glynn, P. W. Stochastic approximation for Monte Carlo optimization. In Proceedings of the 1986 Winter Simulation Conference, (1986), pp. 356-365.
 Ho, Y. C. Performance evaluation and perturbation analysis of discrete event dynamic systems. IEEE Tram. on Auto. Control AC-32, (1987), 563-572.
 Ho, Y. C. and Cassandras, C. A new approach to the analysis of discrete event dynamic systems. Automatica 29, (1983), 149-167.
 Reiman, M. I. and Weiss, A. Sensitivity analysis for simulations via likelihood ratios. Oper. Res. 37, (1989), 830-844.
 Rubinstein, R. Y. Sensitivity analysis and performance extrapolation for computer simulation models. Oper. Rex 37, (1989), 72-81.