A COMPARISON OF OUTPUT-ANALYSIS METHODS FOR SIMULATIONS OF PROCESSES WITH MULTIPLE REGENERATION SEQUENCES. James M. Calvin Marvin K.
|
|
- Geraldine Brooke Nichols
- 5 years ago
- Views:
Transcription
1 Proceedings of the 00 Winter Simulation Conference E. Yücesan, C.-H. Chen, J. L. Snowdon, and J. M. Charnes, eds. A COMPARISON OF OUTPUT-ANALYSIS METHODS FOR SIMULATIONS OF PROCESSES WITH MULTIPLE REGENERATION SEQUENCES James M. Calvin Marvin K. Nakayama Department of Computer Science New Jersey Institute of Technology Newark, NJ 070, U.S.A. ABSTRACT We compare several simulation estimators for a performance measure of a process having multiple regeneration sequences. We examine the setting of two regeneration sequences. We compare two existing estimators, the permuted estimator and the semi-regenerative estimator, and two new estimators, a type of U-statistic estimator and a type of V -statistic estimator. The last two estimators are obtained by resampling trajectories without and with replacement, respectively. The permuted estimator and the U-statistic estimator turn out to be equivalent, but the others are in general different. We show that when estimating the second moment of a cumulative cycle reward, the semi-regenerative and V -statistic estimators have non-negative bias, with the semi-regenerative bias being larger. The permuted estimator was previously shown to be unbiased. Although some of the estimators have different small-sample properties, they all satisfy central limit theorems with the same asymptotic variance constant. INTRODUCTION A regenerative process is a process that has an infinite sequence of stopping times, known as regeneration points, at which times the process probabilistically restarts. The evolution of the process between two successive regeneration points is called a regenerative cycle, and cycles are i.i.d. e.g., see Shedler 99). For example, for an irreducible, positiverecurrent Markov chain, the successive hitting times to a fixed state constitute a regeneration sequence. The standard regenerative method Crane and Iglehart 975) exploits the i.i.d. cycle structure in regenerative processes to construct asymptotically valid confidence intervals for a large class of performance measures. Many regenerative processes possess several different regeneration sequences. For example, for a Markov chain, each fixed state yields a regeneration sequence, but there are many choices for the state to use. The standard regenerative method only uses one of these sequences, and in this paper, we examine several estimators that exploit more than one sequence. Throughout the paper, we consider estimating the second moment of the cumulative cycle reward of a discretetime Markov chain on a discrete state space, but the ideas can also be applied to other performance measures and more general processes having multiple regeneration sequences. Also, we consider only two regeneration sequences. Several methods to exploit multiple regeneration sequences have been proposed. 998, 000) developed permuted regenerative estimators, which can be described in the context of regeneration sequences as follows. Run a simulation of a fixed number of cycles from one sequence, and based on the resulting sample path, calculate an estimate of the performance measure. Construct a new sample path from the original one by permuting the cycles from both sequences, and calculate an estimate based on the new path. Averaging over all permutations yields the permuted estimator. Calculation of this estimator does not require constructing all permutations since a simple closedform representation for the estimator can be derived. Calvin and Nakayama 000) show that the permuted estimator always has the same expectation as the standard regenerative estimator, so if the standard estimator is unbiased, so is the permuted estimator. Moreover, compared to the standard estimator, the permuted estimator always has no larger, and typically strictly smaller, variance. Calvin, Glynn, and Nakayama 00) developed another approach, the semi-regenerative method, which gets its name from its close connection to semi-regenerative processes see Section 0.6 of Çinlar 975). The idea is to fix a set A of states, and to break up the sample path into trajectories determined by successive entrances into A. Using the semi-regenerative structure, one derives an alternative representation of the performance measure in terms of expectations of random quantities defined over trajectories, and then constructs an estimator by replacing each
2 expectation by its natural estimator. Calvin, Glynn, and Nakayama 00 also present a stratified semi-regenerative estimator in which trajectories are sampled in an i.i.d. fashion, but we will not consider that approach in this paper.) Other methods for simulating processes with multiple regeneration sequences include the almost regenerative method Gunther and Wolff 980) and A-segments Zhang and Ho 99). In this paper, we present two new ways to construct estimators. Both approaches start with a sample path of a fixed number of cycles with respect to one regeneration sequence. The path is then broken into trajectories determined by visits into a set A, and the new estimators are obtained by resampling the trajectories and averaging over all possible resamples. One estimator considers resampling without replacement, which is a type of U-statistic, and the other resamples with replacement, which is a type of V -statistic. See Serfling 980, Chapter 5, for details on U- and V -statistics.) We compare four estimators: the permuted, semiregenerative, U-statistic, and V -statistic estimators. We show that the U-statistic and permuted estimators are equivalent. The other two estimators are different, and we derive the exact differences of all the estimators based on a sample path of a fixed number m of cycles from one regeneration sequence. It turns out that the semi-regenerative and V - statistic estimators are biased high for any m, with larger bias for the semi-regenerative estimator. Asymptotically, all four estimators are equivalent and satisfy central limit theorems with the same asymptotic variance constant. The rest of the paper has the following organization. In Section we describe the mathematical model and the performance measure considered. We derive the semiregenerative estimator in Section, and we present the permuted estimator from 000) in Section 4. Section 5 contains the resampled estimators, and we compare the four estimators in Section 6. We state some conclusions in Section 7. All proofs of theorems stated in this paper are given in 00). FRAMEWORK For simplicity, we specialize the development to discretetime Markov chains, but the results hold for more general processes with multiple regeneration sequences. Let X = X j : j = 0,,,...) be an irreducible positive-recurrent discrete-time Markov chain on a discrete state space S = {,,,...}. For any state x S, we can define a sequence of regeneration points Tx) = T k x) : k = 0,,,...), where T 0 x) = inf{j 0 : X j = x} and T k x) = inf{j > T k x) : X j = x}, k. We call the sample-path segment X j : T k x) j<t k x)) the kth x-cycle. Fix a state w S. Forx S, define αx) = τ fx k ), where τ = inf{k : X k = w} and f : S Ris a reward function. Our goal is to estimate αw). We will compare the estimator of αw) using the semi-regenerative method and the permuted estimator of 000), and we will present two other estimators of αw) based on resampling. Fix a set A S with w A. Define T 0 = inf{j : X j A} and T k = inf{j T k : X j A} for k. We will carry out our analysis with A ={, } and w =. Also, we will assume that we simulate one sample path of a fixed number m of -cycles, and let the sample path be X m = X j : j = 0,,...,T m )), where X 0 =. Note that τ = T ). Also, let W = W k : k = 0,,,...) with W k = X Tk for k 0, which is the embedded chain of X on visits to the set A. Let M = m {T k ), k 0 : T k ) <T m )}, which is the number of returns to the set A up to time T m ), sot M = T m ). The standard estimator of α) based on the sample path X m is where for k =,,...,m. t STD X m ) = m Y k = T k ) j=t k ) m Yk, ) k= fx j ) SEMI-REGENERATIVE ESTIMATOR We now derive the semi-regenerative estimator for α). The key to applying the semi-regenerative method is deriving an alternative representation for α) based on returns to the set A. Let T = inf{k : X k A}. Forx A, note that [ T αx) = fx k ) τ Iτ > T,X T = y) fx k ) y A
3 Expanding this formula for αx), we see that it is composed of three terms: αx) = T fx k ) τ T Iτ > T,X T = y) fx k ) fx j ) y A j=0 τ Iτ > T,X T = y) fx k ). ) y A We now derive new expressions for the second and third terms in ). Each summand in the second term in ) satisfies τ Iτ > T,X T = y) = Fx,y)ey), T fx k ) j=0 fx j ) where ] T Fx,y) = [Iτ > T,X T = y) fx k ), Also so letting [ τ ] ey) = E y fx k ). [ T ] ex) = fx k ) y A [Iτ > T,X T = y)] ey), [ T ] qx) = fx k ), and Gx, y) = [Iτ > T,X T = y)],e= ex) : x A), q = qx) : x A), and G = Gx, y) : x,y A), we see that e = q Ge, ore = I G) q. Thus, each summand in the second term in ) satisfies τ Iτ > T,X T = y) = Fx,y)I G) q)y). T fx k ) j=0 X j For the summand in the third term in ), we make the observation that τ Iτ > T,X T = y) fx k ) Now define = Gx, y) αy). cx) = T fx k ). Let α = αx) : x A), c = cx) : x A), and F = F x, y) : x,y A). Then, putting this all together, we get α = c FI G) q Gα, so α = I G) c FI G) q). Now suppose A = {, } and w =. Note that Gx, ) = Fx,) = 0 for x A. Let J = J x, y) : x,y A) with J = I G), and note that Thus, J = = ) G, ) 0 G, ) G, )/ G, )) 0 / G, )) )[ ) J, ) c) α = 0 J, ) c) ) 0 F, ) J, ) 0 F, ) 0 J, ) ) c) J, )c) = J, )c) so ). ) q) q) F, ) J, )F, ))J, )q) J, F, )q) )] ), α) = c) J, )c) [F, ) J, )F, )]J, )q). ) To derive an estimator based on ) from the sample path X m of m -cycles, we need estimates of the quantities on the right-hand side of ). Recall that M was defined such that T M = T m ). Forx,y A, we call a trajectory that begins in state x and ends in state y an x, y)-trajectory, and define hx, y) = M IW k = x,w k = y) as the number of such trajectories along the sample path X m. Let Hx) =
4 M IW k = x), which is the number of trajectories in X m that begin in state x. Note that hx, )hx, ) = Hx) for x A. Also, h, ) = h, ) since X m is a path of a fixed number of -cycles, and H) = m. For x,y A, let T x, y) = inf{t k : k 0, W k = x, W k = y} and T i x, y) = inf{t k T i x, y) : k 0, W k = x, W k = y} for i. Also, let T x, y) = inf{t k : k 0, W k = x, W k = y} and T i x, y) = inf{t k : k 0, T k T i x, y), W k = x, W k = y} for i. Thus, for k =,,...,hx,y), the kth trajectory starting in state x A and ending in state y A begins at time T k x, y) and finishes at time T k x, y), and we define Y k x, y) as the sum of the rewards along that trajectory; i.e., Y k x, y) = T k x,y) j=t k x,y) fx j ). Define S j x, y) = hx,y) k= Y k x, y) j. Now for x A, observe that S x, )S x, ))/H x) and S x, ) S x, ))/H x) are the natural estimators of cx) and qx), respectively. Also, we use hx, )/H x) and S x, )/H x) as estimators for Gx, ) and Fx,), respectively. Our estimator of J, ) = G, )/ G, )) is then ˆ J, ) = = ) h, ) h, ) H) H) ) h, ) H) = H) H) h, ) H) since h, ) h, ) = H) and h, ) = h, ). Similarly, our estimator of J, ) = / G, )) is J, ˆ ) = H)/h, ). Thus, the semi-regenerative estimator of α) based on the sample path X m is t SR X m ) = H) S, )S, )) ) H) S, ) S, ) H) H) [ ) ] S, ) H) S, ) H) H) H) ) H) S, ) S, ) h, ) H) = S, ) S, ) S, ) S, ) m h, ) [S, )S, ) S, )S, ) S, )S, ) S, ]) 4) since H) = m. To simplify the comparison with other estimators, we define so Q = m S, ) S, ) S, ) S, ) [ S, )S, ) S, )S, ) h, ) S, )S, ) ]), 5) t SR X m ) = Q 4 PERMUTED ESTIMATOR mh, ) S,. 6) The permuted estimator of α) based on the sample path X m is obtained as follows. Recall we defined the standard estimator of α) as t STD X m ) in ). The permuted estimator of α) is obtained by permuting -cycles and -cycles to obtain a new sample path X m, and averaging t STD X m ) over all possible permuted paths X m. Calvin and Nakayama 000) showed that the permuted estimator is t P X m ) m S, ) S, ) S, ) [ S, )S, ) h, ) ] S, )S, ) S, )S, ) h, ) h, ) S, ) h, ) S, = Q mh, ) ) where Q is defined in 5). 5 RESAMPLED ESTIMATORS S, S, ) ) ), 7) We now present two estimators obtained via resampling the x, y)-trajectories for x,y A. One estimator is based on resampling the x, y)-trajectories without replacement i.e., a permutation of trajectories), and the other resamples with replacement. Fix a sample path X m. For x,y A, define x, y) as the set of permutations of,,...,hx,y)), and let = x,y A x, y). Let x, y) ={i,i,...,i hx,y) ) : i j hx, y) for j =,,...,hx,y)}, which is the set of all hx, y)-dimensional vectors in which the components are selected with replacement from
5 {,,...,hx,y)}. Let = x,y A x, y). Define = {i,i,...,i h,) ) : i =, i h,) h, ), i j i j for j =,,...,h, ) }. Let Ɣ = and Ɣ =, and let K, D) denote a generic element of Ɣ, with K = Kx, y) : x,y A), Kx,y) = K i x, y) : i =,,...,hx,y)) x, y), and D = D,D,...,D h,) ). Also, define D h,) = h, ) and the function g : R as gk, D) = [ h,) Y Ki,), m i= Y Ki,), ) Y Ki,), ) h,) i= D i ] Y Kj,), ) j=d i = [ h,) Y Ki,), m i= Y Ki,), Y Ki,), h,) i= D i Y Kj,), ) j=d i Y Ki,), )Y Ki,), ) { Y Ki,), ) Y Ki,), ) } D i )] Y Kj,), ). 8) j=d i Observe that the function g is conditional on the original sample path X m. Note that x, y) resp., x, y)) is the set of all possible orderings of hx, y) x, y)-trajectories, where the x, y)- trajectories are chosen without resp., with) replacement, and in 8), Kx,y) x, y) is one such ordering of hx, y) x, y)-trajectories. Also, note that is the set of vectors specifying which, )-trajectories fall in between each pairing of, )-trajectory and, )-trajectory determined by K, ) and K, ); i.e., for D,...,D h,) ), slots D j to D j are the locations of the, )-trajectories in the ordered list K, ), ) of, )-trajectories that fall between the, )-trajectory indexed by K j, ) and the, )-trajectory indexed by K j, ). Observe that if D j = D j, then there are no, )-trajectories that fall between the, )-trajectory indexed by K j, ) and the, )-trajectory indexed by K j, ). Suppose we set K x, y) =,,...,hx,y)) for all x,y A, and let K = K x, y) : x,y A). Note that K x, y) x, y) for all x,y A, and K. Define D = D,D,...,D h,) ) with D i = min{j : T j, ) >T i, )} for i. Note that D. Then gk,d ) = t STD X m ), the standard estimator of α). We now define the resampled estimators and t U X m ) = Ɣ t V X m ) = Ɣ K,D) Ɣ K,D) Ɣ gk, D) 9) gk, D), 0) where t U X m ) is the estimator based on resampling without replacement, and t V X m ) is the estimator with replacement. In addition to resampling trajectories, both estimators also average over all possible allocations of h, ), )- trajectories to the -cycles containing a hit to state, which is accomplished in 9) and 0) by averaging over D. Note that t U X m ) is a type of U-statistic and t V X m ) is a type of V -statistic. Figure presents an example of a sample path X m top) and a path middle) obtained from X m by resampling without replacement the x, y)-trajectories, for x, y A. State corresponds to the horizontal axis, and state corresponds to the dashed horizontal line. The original path X m, which is shown on top, has m = 5 -cycles. Also, the original path has h, ) =, )-trajectories, h, ) =, )- trajectories, h, ) =, )-trajectories, and h, ) =, )-trajectories. Moreover, K, ) =, K, ) = ; K, ) =, K, ) =, K, ) = ; K, ) =, K, ) =, K, ) = ; K, ) =, K, ) =, K, ) = ; D =, D =, D = 4, D 4 = 4. The resampled path shown on the bottom has K, ) =, K, ) = ; K, ) =, K, ) =, K, ) = ; K, ) =, K, ) =, K, ) = ; K, ) =, K, ) =, K, ) = ; D =, D =, D =, D 4 = 4. The bottom path in Figure is one possible path obtained from X m by resampling with replacement the trajectories. The resampled path has K, ) =, K, ) = ; K, ) =, K, ) =, K, ) = ; K, ) =, K, ) =, K, ) = ; K, ) =, K, ) =, K, ) = ; D =, D =, D =, D 4 = 4. In the resampled path, Y, ) appears twice, and Y, ) does not appear. Also, for the, )-trajectories, Y, ) appears once, Y, ) does not appear, and Y, ) appears twice. For the, )-trajectories, each Y k, ) appears once. For
6 Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Y,) Figure : A Sample Path top) and Paths Obtained by Resampling Trajectories Without Replacement middle) and With Replacement bottom)
7 the, )-trajectories, Y, ) does not appear, Y, ) appears once, and Y, ) appears twice. Theorem For all X m and m, t U X m ) = t P X m ) and t V X m ) = Q h, ) ) mh, ) h, ) ) S,, ) where Q is defined in 5). Thus, since t U X m ) = t P X m ), we see that permuting trajectories and reallocating, )-trajectories to the -cycles leads to the same estimator as permuting cycles. 6 COMPARING THE ESTIMATORS We start by comparing the semi-regenerative estimator t SR X m ) in 6) and the permuted estimator t P X m ) in 7). In the second term in the right-hand side of 7), the sum of the squares of the Y k, ) i.e., S, )) is subtracted from the square of the sum i.e., S, ), whereas the second term in the right-hand side of 6) only has the square of the sum. A similar situation occurs when comparing the V -statistic estimator t V X m ) in ) and the permuted estimator t P X m ). The estimators of α) satisfy a complete ordering. Theorem For all X m and m, where t SR X m ) t V X m ) t U X m ) = t P X m ), ) the first inequality in ) is strict if and only if S, ) = 0; the second inequality in ) is strict if and only if not all the Y k, ), k =,,...,h, ), are the same i.e., there exists i, j {,,..., h, )} such that Y i, ) = Y j, )); t SR X m )>t P X m ) if and only if Y k, ) = 0 for some k =,,...,h, ). Moreover, where E[t SR X m )] E[t V X m )] E[t U X m )] ) = E[t P X m )] = α), the first inequality in ) is strict if and only if P [S, ) = 0] < ; the second inequality in ) is strict if and only if Var[Y, )] > 0; E[t SR X m )] > E[t P X m )] if and only if P [Y, ) = 0] <. We now compare some asymptotic properties of our D estimators. Let denote convergence in distribution and Na,b) denote a normal distribution with mean a and variance b. Theorem Assume E[ T ) j=0 fx j ) ) 4 ] <. Then for t X m ) defined as either t P X m ), t SR X m ), or t V X m ), we have that t X m ) converges with probability one to α) while a suitably standardized version of t X m ) converges in distribution to a normal distribution with mean zero and variance σα, i.e., and t X m ) α), a.s., m [ t X m ) α)] D N0,σ α ) as m, for some constant σα > 0, where σ α is the same for all three estimators and is given in Calvin and Nakayama 000). For all of the estimators of α, one can estimate the asymptotic variance constant σα from one simulation run. Thus, one does not need to run multiple replications to obtain estimates of the asymptotic variance. See Calvin and Nakayama 000) for details. 7 CONCLUSIONS In this paper we showed that a wide spectrum of estimators can be obtained in simulations of processes having multiple regeneration sequences. We compared some of these estimators, two of which were proposed previously the permuted and semi-regenerative estimators) and two are new the U- and V -statistic estimators). We showed that the U-statistic and permuted estimators are equivalent. All of the other estimators are in general different, but they are asymptotically equivalent and they all satisfy the same central limit theorem. When estimating the second moment of a cycle reward, we showed that the estimators satisfy a complete ordering, with the semi-regenerative estimator being the largest, the V -statistic estimator being next largest, and the U-statistic and permuted estimators being the smallest. More work is needed to analyze the behavior of these estimators for other performance measures and to determine which is the best perhaps in terms of smallest mean squared error) for a specific measure in the small-sample context. ACKNOWLEDGMENTS This work was supported by the National Science Foundation under grants DMI and DMI Thus, t SR X m ) and t V X m ) are biased high.
8 REFERENCES Calvin, J. M. and M. K. Nakayama Using permutations in regenerative simulations to reduce variance. ACM Transactions on Modeling and Computer Simulation 8: 5 9. Calvin, J. M. and M. K. Nakayama Simulation of processes with multiple regeneration sequences. Probability in the Engineering and Informational Sciences, 4: Calvin, J. M. and M. K. Nakayama. 00. Resampled regenerative estimators. Forthcoming. Calvin, J. M., P. W. Glynn and M. K. Nakayama. 00. The semi-regenerative method of simulation output analysis. Forthcoming. Çinlar, E Introduction to Stochastic Processes. Englewood Cliffs, N.J.: Prentice-Hall. Crane, M. and D. L. Iglehart Simulating stable stochastic systems, III: Regenerative processes and discrete-event simulations. Operations Research : 45. Gunther, F. L. and R. W. Wolff The almost regenerative method for stochastic system simulations. Operations Research 8: Serfling, R.J Approximation Theorems of Mathematical Statistics. New York: Wiley. Shedler, G. S. 99. Regenerative Stochastic Simulation. San Diego: Academic Press. Zhang, B. and Y.-C. Ho. 99. Improvements in the likelihood ratio method for steady-state sensitivity analysis and simulation. Performance Analysis 5: interests include applied probability, statistics, simulation and modeling. AUTHOR BIOGRAPHIES JAMES M. CALVIN is an associate professor in the Department of Computer Science at the New Jersey Institute of Technology. He received a Ph.D. in operations research from Stanford University and is an associate editor for ACM Transactions on Modeling and Computer Simulation. Besides simulation output analysis, his research interests include global optimization and probabilistic analysis of algorithms. MARVIN K. NAKAYAMA is an associate professor in the Department of Computer Science at the New Jersey Institute of Technology. He received a Ph.D. in operations research from Stanford University. He won second prize in the 99 George E. Nicholson Student Paper Competition sponsored by INFORMS and is a recipient of a CAREER Award from the National Science Foundation. He is the area editor for the Stochastic Modeling Area of ACM Transactions on Modeling and Computer Simulation and an associate editor for Informs Journal on Computing. His research
The Semi-Regenerative Method of Simulation Output Analysis
The Semi-Regenerative Method of Simulation Output Analysis JAMES M. CALVIN New Jersey Institute of Technology PETER W. GLYNN Stanford University and MARVIN K. NAKAYAMA New Jersey Institute of Technology
More informationIMPORTANCE SAMPLING USING THE SEMI-REGENERATIVE METHOD. Marvin K. Nakayama
Proceedings of the 2001 Winter Simulation Conference B.A.Peters,J.S.Smith,D.J.Medeiros,andM.W.Rohrer,eds. IMPORTANCE SAMPLING USING THE SEMI-REGENERATIVE METHOD James M. Calvin Department of Computer Science
More informationA GENERAL FRAMEWORK FOR THE ASYMPTOTIC VALIDITY OF TWO-STAGE PROCEDURES FOR SELECTION AND MULTIPLE COMPARISONS WITH CONSISTENT VARIANCE ESTIMATORS
Proceedings of the 2009 Winter Simulation Conference M. D. Rossetti, R. R. Hill, B. Johansson, A. Dunkin, and R. G. Ingalls, eds. A GENERAL FRAMEWORK FOR THE ASYMPTOTIC VALIDITY OF TWO-STAGE PROCEDURES
More informationProceedings of the 2014 Winter Simulation Conference A. Tolk, S. Y. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds.
Proceedings of the 2014 Winter Simulation Conference A. Tolk, S. Y. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds. CONSTRUCTING CONFIDENCE INTERVALS FOR A QUANTILE USING BATCHING AND
More informationOn Derivative Estimation of the Mean Time to Failure in Simulations of Highly Reliable Markovian Systems
On Derivative Estimation of the Mean Time to Failure in Simulations of Highly Reliable Markovian Systems Marvin K. Nakayama Department of Computer and Information Science New Jersey Institute of Technology
More informationOn Finding Optimal Policies for Markovian Decision Processes Using Simulation
On Finding Optimal Policies for Markovian Decision Processes Using Simulation Apostolos N. Burnetas Case Western Reserve University Michael N. Katehakis Rutgers University February 1995 Abstract A simulation
More informationProceedings of the 2015 Winter Simulation Conference L. Yilmaz, W. K. V. Chan, I. Moon, T. M. K. Roeder, C. Macal, and M. D. Rossetti, eds.
Proceedings of the 2015 Winter Simulation Conference L. Yilmaz, W. K. V. Chan, I. Moon, T. M. K. Roeder, C. Macal, and M. D. Rossetti, eds. UNBIASED MONTE CARLO COMPUTATION OF SMOOTH FUNCTIONS OF EXPECTATIONS
More informationREFLECTED VARIANCE ESTIMATORS FOR SIMULATION. Melike Meterelliyoz Christos Alexopoulos David Goldsman
Proceedings of the 21 Winter Simulation Conference B. Johansson, S. Jain, J. Montoya-Torres, J. Hugan, and E. Yücesan, eds. REFLECTE VARIANCE ESTIMATORS FOR SIMULATION Melike Meterelliyoz Christos Alexopoulos
More informationRegenerative Processes. Maria Vlasiou. June 25, 2018
Regenerative Processes Maria Vlasiou June 25, 218 arxiv:144.563v1 [math.pr] 22 Apr 214 Abstract We review the theory of regenerative processes, which are processes that can be intuitively seen as comprising
More informationTHE BALANCED LIKELIHOOD RATIO METHOD FOR ESTIMATING PERFORMANCE MEASURES OF HIGHLY RELIABLE SYSTEMS
Proceedings of the 1998 Winter Simulation Conference D.J. Medeiros, E.F. Watson, J.S. Carson and M.S. Manivannan, eds. THE BALANCED LIKELIHOOD RATIO METHOD FOR ESTIMATING PERFORMANCE MEASURES OF HIGHLY
More informationCONFIDENCE INTERVALS FOR QUANTILES WITH STANDARDIZED TIME SERIES. James M. Calvin Marvin K. Nakayama
Proceedings of the 2013 Winter Simulation Conference R. Pasupathy, S.-H. Kim, A. Tolk, R. Hill, and M. E. Kuhl, eds. CONFIDENCE INTERVALS FOR QUANTILES WITH STANDARDIZED TIME SERIES James M. Calvin Marvin
More informationCentral limit theorems for ergodic continuous-time Markov chains with applications to single birth processes
Front. Math. China 215, 1(4): 933 947 DOI 1.17/s11464-15-488-5 Central limit theorems for ergodic continuous-time Markov chains with applications to single birth processes Yuanyuan LIU 1, Yuhui ZHANG 2
More informationLecture 2: Convex Sets and Functions
Lecture 2: Convex Sets and Functions Hyang-Won Lee Dept. of Internet & Multimedia Eng. Konkuk University Lecture 2 Network Optimization, Fall 2015 1 / 22 Optimization Problems Optimization problems are
More informationSEQUENTIAL ESTIMATION OF THE STEADY-STATE VARIANCE IN DISCRETE EVENT SIMULATION
SEQUENTIAL ESTIMATION OF THE STEADY-STATE VARIANCE IN DISCRETE EVENT SIMULATION Adriaan Schmidt Institute for Theoretical Information Technology RWTH Aachen University D-5056 Aachen, Germany Email: Adriaan.Schmidt@rwth-aachen.de
More informationNEW ESTIMATORS FOR PARALLEL STEADY-STATE SIMULATIONS
roceedings of the 2009 Winter Simulation Conference M. D. Rossetti, R. R. Hill, B. Johansson, A. Dunkin, and R. G. Ingalls, eds. NEW ESTIMATORS FOR ARALLEL STEADY-STATE SIMULATIONS Ming-hua Hsieh Department
More informationOn the estimation of the mean time to failure by simulation
On the estimation of the mean time to failure by simulation Peter Glynn, Marvin Nakayama, Bruno Tuffin To cite this version: Peter Glynn, Marvin Nakayama, Bruno Tuffin. On the estimation of the mean time
More informationGeneral Glivenko-Cantelli theorems
The ISI s Journal for the Rapid Dissemination of Statistics Research (wileyonlinelibrary.com) DOI: 10.100X/sta.0000......................................................................................................
More information6 Markov Chain Monte Carlo (MCMC)
6 Markov Chain Monte Carlo (MCMC) The underlying idea in MCMC is to replace the iid samples of basic MC methods, with dependent samples from an ergodic Markov chain, whose limiting (stationary) distribution
More informationP i [B k ] = lim. n=1 p(n) ii <. n=1. V i :=
2.7. Recurrence and transience Consider a Markov chain {X n : n N 0 } on state space E with transition matrix P. Definition 2.7.1. A state i E is called recurrent if P i [X n = i for infinitely many n]
More informationSAMPLE CHAPTERS UNESCO EOLSS STATIONARY PROCESSES. E ξ t (1) K.Grill Institut für Statistik und Wahrscheinlichkeitstheorie, TU Wien, Austria
STATIONARY PROCESSES K.Grill Institut für Statistik und Wahrscheinlichkeitstheorie, TU Wien, Austria Keywords: Stationary process, weak sense stationary process, strict sense stationary process, correlation
More informationSection Properties of Rational Expressions
88 Section. - Properties of Rational Expressions Recall that a rational number is any number that can be written as the ratio of two integers where the integer in the denominator cannot be. Rational Numbers:
More informationHITTING TIME IN AN ERLANG LOSS SYSTEM
Probability in the Engineering and Informational Sciences, 16, 2002, 167 184+ Printed in the U+S+A+ HITTING TIME IN AN ERLANG LOSS SYSTEM SHELDON M. ROSS Department of Industrial Engineering and Operations
More informationSelected Exercises on Expectations and Some Probability Inequalities
Selected Exercises on Expectations and Some Probability Inequalities # If E(X 2 ) = and E X a > 0, then P( X λa) ( λ) 2 a 2 for 0 < λ
More informationTheory and Applications of Stochastic Systems Lecture Exponential Martingale for Random Walk
Instructor: Victor F. Araman December 4, 2003 Theory and Applications of Stochastic Systems Lecture 0 B60.432.0 Exponential Martingale for Random Walk Let (S n : n 0) be a random walk with i.i.d. increments
More informationStochastic Realization of Binary Exchangeable Processes
Stochastic Realization of Binary Exchangeable Processes Lorenzo Finesso and Cecilia Prosdocimi Abstract A discrete time stochastic process is called exchangeable if its n-dimensional distributions are,
More informationAverage Reward Parameters
Simulation-Based Optimization of Markov Reward Processes: Implementation Issues Peter Marbach 2 John N. Tsitsiklis 3 Abstract We consider discrete time, nite state space Markov reward processes which depend
More informationIBM Almaden Research Center, San Jose, California, USA
This article was downloaded by: [Stanford University] On: 20 July 2010 Access details: Access Details: [subscription number 731837804] Publisher Taylor & Francis Informa Ltd Registered in England and Wales
More informationLet (Ω, F) be a measureable space. A filtration in discrete time is a sequence of. F s F t
2.2 Filtrations Let (Ω, F) be a measureable space. A filtration in discrete time is a sequence of σ algebras {F t } such that F t F and F t F t+1 for all t = 0, 1,.... In continuous time, the second condition
More informationAMCS243/CS243/EE243 Probability and Statistics. Fall Final Exam: Sunday Dec. 8, 3:00pm- 5:50pm VERSION A
AMCS243/CS243/EE243 Probability and Statistics Fall 2013 Final Exam: Sunday Dec. 8, 3:00pm- 5:50pm VERSION A *********************************************************** ID: ***********************************************************
More information798 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 44, NO. 10, OCTOBER 1997
798 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL 44, NO 10, OCTOBER 1997 Stochastic Analysis of the Modulator Differential Pulse Code Modulator Rajesh Sharma,
More informationTRANSIENT SIMULATION VIA EMPIRICALLY BASED COUPLING
Probability in the Engineering and Informational Sciences, 13, 1999, 147 167+ Printed in the U+S+A+ TRANSIENT SIMULATION VIA EMPIRICALLY BASED COUPLING EUGENE W. WONG, PETER W. GLYNN, AND DONALD L. IGLEHART*
More informationLecture 3: Policy Evaluation Without Knowing How the World Works / Model Free Policy Evaluation
Lecture 3: Policy Evaluation Without Knowing How the World Works / Model Free Policy Evaluation CS234: RL Emma Brunskill Winter 2018 Material builds on structure from David SIlver s Lecture 4: Model-Free
More informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationAnswers and expectations
Answers and expectations For a function f(x) and distribution P(x), the expectation of f with respect to P is The expectation is the average of f, when x is drawn from the probability distribution P E
More informationMS&E 321 Spring Stochastic Systems June 1, 2013 Prof. Peter W. Glynn Page 1 of 10
MS&E 321 Spring 12-13 Stochastic Systems June 1, 2013 Prof. Peter W. Glynn Page 1 of 10 Section 3: Regenerative Processes Contents 3.1 Regeneration: The Basic Idea............................... 1 3.2
More informationAdvanced Computer Networks Lecture 2. Markov Processes
Advanced Computer Networks Lecture 2. Markov Processes Husheng Li Min Kao Department of Electrical Engineering and Computer Science University of Tennessee, Knoxville Spring, 2016 1/28 Outline 2/28 1 Definition
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables THE UNIVERSITY OF MANCHESTER. 21 June :45 11:45
Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS 21 June 2010 9:45 11:45 Answer any FOUR of the questions. University-approved
More informationChapter 2: Markov Chains and Queues in Discrete Time
Chapter 2: Markov Chains and Queues in Discrete Time L. Breuer University of Kent 1 Definition Let X n with n N 0 denote random variables on a discrete space E. The sequence X = (X n : n N 0 ) is called
More informationThis paper discusses the application of the likelihood ratio gradient estimator to simulations of highly dependable systems. This paper makes the foll
Likelihood Ratio Sensitivity Analysis for Markovian Models of Highly Dependable Systems Marvin K. Nakayama IBM Thomas J. Watson Research Center P.O. Box 704 Yorktown Heights, New York 10598 Ambuj Goyal
More informationEstimation theory. Parametric estimation. Properties of estimators. Minimum variance estimator. Cramer-Rao bound. Maximum likelihood estimators
Estimation theory Parametric estimation Properties of estimators Minimum variance estimator Cramer-Rao bound Maximum likelihood estimators Confidence intervals Bayesian estimation 1 Random Variables Let
More informationOverall Plan of Simulation and Modeling I. Chapters
Overall Plan of Simulation and Modeling I Chapters Introduction to Simulation Discrete Simulation Analytical Modeling Modeling Paradigms Input Modeling Random Number Generation Output Analysis Continuous
More informationTo Be (a Circle) or Not to Be?
To Be (a Circle) or Not to Be? Hassan Boualem and Robert Brouzet Hassan Boualem (hassan.boualem@univ-montp.fr) received his Ph.D. in mathematics from the University of Montpellier in 99, where he studied
More informationSTOCHASTIC PETRI NETS FOR DISCRETE-EVENT SIMULATION
STOCHASTIC PETRI NETS FOR DISCRETE-EVENT SIMULATION Peter J. Haas IBM Almaden Research Center San Jose, CA Part I Introduction 2 Peter J. Haas My Background Mid 1980 s: PhD student studying discrete-event
More informationExact Simulation of the Stationary Distribution of M/G/c Queues
1/36 Exact Simulation of the Stationary Distribution of M/G/c Queues Professor Karl Sigman Columbia University New York City USA Conference in Honor of Søren Asmussen Monday, August 1, 2011 Sandbjerg Estate
More information(b) What is the variance of the time until the second customer arrives, starting empty, assuming that we measure time in minutes?
IEOR 3106: Introduction to Operations Research: Stochastic Models Fall 2006, Professor Whitt SOLUTIONS to Final Exam Chapters 4-7 and 10 in Ross, Tuesday, December 19, 4:10pm-7:00pm Open Book: but only
More informationBivariate Paired Numerical Data
Bivariate Paired Numerical Data Pearson s correlation, Spearman s ρ and Kendall s τ, tests of independence University of California, San Diego Instructor: Ery Arias-Castro http://math.ucsd.edu/~eariasca/teaching.html
More informationFINDING THE BEST IN THE PRESENCE OF A STOCHASTIC CONSTRAINT. Sigrún Andradóttir David Goldsman Seong-Hee Kim
Proceedings of the 2005 Winter Simulation Conference M. E. Kuhl, N. M. Steiger, F. B. Armstrong, and J. A. Joines, eds. FINDING THE BEST IN THE PRESENCE OF A STOCHASTIC CONSTRAINT Sigrún Andradóttir David
More informationMS&E 226: Small Data. Lecture 6: Bias and variance (v2) Ramesh Johari
MS&E 226: Small Data Lecture 6: Bias and variance (v2) Ramesh Johari ramesh.johari@stanford.edu 1 / 47 Our plan today We saw in last lecture that model scoring methods seem to be trading o two di erent
More informationNon-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings. Krzysztof Burdzy University of Washington
Non-Essential Uses of Probability in Analysis Part IV Efficient Markovian Couplings Krzysztof Burdzy University of Washington 1 Review See B and Kendall (2000) for more details. See also the unpublished
More informationMonte Carlo Method for Finding the Solution of Dirichlet Partial Differential Equations
Applied Mathematical Sciences, Vol. 1, 2007, no. 10, 453-462 Monte Carlo Method for Finding the Solution of Dirichlet Partial Differential Equations Behrouz Fathi Vajargah Department of Mathematics Guilan
More informationAsymptotic properties of the maximum likelihood estimator for a ballistic random walk in a random environment
Asymptotic properties of the maximum likelihood estimator for a ballistic random walk in a random environment Catherine Matias Joint works with F. Comets, M. Falconnet, D.& O. Loukianov Currently: Laboratoire
More informationApproximating Martingales in Continuous and Discrete Time Markov Processes
Approximating Martingales in Continuous and Discrete Time Markov Processes Rohan Shiloh Shah May 6, 25 Contents 1 Introduction 2 1.1 Approximating Martingale Process Method........... 2 2 Two State Continuous-Time
More informationINTRODUCTION TO MARKOV CHAIN MONTE CARLO
INTRODUCTION TO MARKOV CHAIN MONTE CARLO 1. Introduction: MCMC In its simplest incarnation, the Monte Carlo method is nothing more than a computerbased exploitation of the Law of Large Numbers to estimate
More informationA Proof of Aldous Spectral Gap Conjecture. Thomas M. Liggett, UCLA. Stirring Processes
A Proof of Aldous Spectral Gap Conjecture Thomas M. Liggett, UCLA Stirring Processes Given a graph G = (V, E), associate with each edge e E a Poisson process Π e with rate c e 0. Labels are put on the
More informationStatistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation
Statistics - Lecture One Charlotte Wickham wickham@stat.berkeley.edu http://www.stat.berkeley.edu/~wickham/ Outline 1. Basic ideas about estimation 2. Method of Moments 3. Maximum Likelihood 4. Confidence
More informationOn Differentiability of Average Cost in Parameterized Markov Chains
On Differentiability of Average Cost in Parameterized Markov Chains Vijay Konda John N. Tsitsiklis August 30, 2002 1 Overview The purpose of this appendix is to prove Theorem 4.6 in 5 and establish various
More informationMS&E 321 Spring Stochastic Systems June 1, 2013 Prof. Peter W. Glynn Page 1 of 7
MS&E 321 Spring 12-13 Stochastic Systems June 1, 213 Prof. Peter W. Glynn Page 1 of 7 Section 9: Renewal Theory Contents 9.1 Renewal Equations..................................... 1 9.2 Solving the Renewal
More informationAlgebra II Vocabulary Word Wall Cards
Algebra II Vocabulary Word Wall Cards Mathematics vocabulary word wall cards provide a display of mathematics content words and associated visual cues to assist in vocabulary development. The cards should
More informationCourse: ESO-209 Home Work: 1 Instructor: Debasis Kundu
Home Work: 1 1. Describe the sample space when a coin is tossed (a) once, (b) three times, (c) n times, (d) an infinite number of times. 2. A coin is tossed until for the first time the same result appear
More informationPOLYNOMIAL ACCELERATION OF MONTE-CARLO GLOBAL SEARCH. James M. Calvin
roceedings of the 999 Winter Simulation Conference. A. Farrington, H. B. Nembhard, D. T. Sturrock, and G. W. Evans, eds. OLYNOMIAL ACCELERATION OF MONTE-CARLO GLOBAL SEARCH James M. Calvin Department of
More informationA Few Notes on Fisher Information (WIP)
A Few Notes on Fisher Information (WIP) David Meyer dmm@{-4-5.net,uoregon.edu} Last update: April 30, 208 Definitions There are so many interesting things about Fisher Information and its theoretical properties
More informationStochastic Process (ENPC) Monday, 22nd of January 2018 (2h30)
Stochastic Process (NPC) Monday, 22nd of January 208 (2h30) Vocabulary (english/français) : distribution distribution, loi ; positive strictement positif ; 0,) 0,. We write N Z,+ and N N {0}. We use the
More informationThe Optimal Stopping of Markov Chain and Recursive Solution of Poisson and Bellman Equations
The Optimal Stopping of Markov Chain and Recursive Solution of Poisson and Bellman Equations Isaac Sonin Dept. of Mathematics, Univ. of North Carolina at Charlotte, Charlotte, NC, 2822, USA imsonin@email.uncc.edu
More information8. Statistical Equilibrium and Classification of States: Discrete Time Markov Chains
8. Statistical Equilibrium and Classification of States: Discrete Time Markov Chains 8.1 Review 8.2 Statistical Equilibrium 8.3 Two-State Markov Chain 8.4 Existence of P ( ) 8.5 Classification of States
More informationWhen is a Markov chain regenerative?
When is a Markov chain regenerative? Krishna B. Athreya and Vivekananda Roy Iowa tate University Ames, Iowa, 50011, UA Abstract A sequence of random variables {X n } n 0 is called regenerative if it can
More informationImproved Methods for Monte Carlo Estimation of the Fisher Information Matrix
008 American Control Conference Westin Seattle otel, Seattle, Washington, USA June -3, 008 ha8. Improved Methods for Monte Carlo stimation of the isher Information Matrix James C. Spall (james.spall@jhuapl.edu)
More informationThe Two Headed Disk: Stochastic Dominance of the Greedy Policy
The Two Headed Disk: Stochastic Dominance of the Greedy Policy Sridhar Seshadri Leonard N. Stern School of Business New York University Doron Rotem Information & Computing Sciences Division Lawrence Berkeley
More informationFormulas for probability theory and linear models SF2941
Formulas for probability theory and linear models SF2941 These pages + Appendix 2 of Gut) are permitted as assistance at the exam. 11 maj 2008 Selected formulae of probability Bivariate probability Transforms
More informationStochastic Processes (Week 6)
Stochastic Processes (Week 6) October 30th, 2014 1 Discrete-time Finite Markov Chains 2 Countable Markov Chains 3 Continuous-Time Markov Chains 3.1 Poisson Process 3.2 Finite State Space 3.2.1 Kolmogrov
More informationarxiv: v1 [math.st] 12 Aug 2017
Existence of infinite Viterbi path for pairwise Markov models Jüri Lember and Joonas Sova University of Tartu, Liivi 2 50409, Tartu, Estonia. Email: jyril@ut.ee; joonas.sova@ut.ee arxiv:1708.03799v1 [math.st]
More informationNotes, March 4, 2013, R. Dudley Maximum likelihood estimation: actual or supposed
18.466 Notes, March 4, 2013, R. Dudley Maximum likelihood estimation: actual or supposed 1. MLEs in exponential families Let f(x,θ) for x X and θ Θ be a likelihood function, that is, for present purposes,
More informationLarge Sample Properties of Estimators in the Classical Linear Regression Model
Large Sample Properties of Estimators in the Classical Linear Regression Model 7 October 004 A. Statement of the classical linear regression model The classical linear regression model can be written in
More informationReview Quiz. 1. Prove that in a one-dimensional canonical exponential family, the complete and sufficient statistic achieves the
Review Quiz 1. Prove that in a one-dimensional canonical exponential family, the complete and sufficient statistic achieves the Cramér Rao lower bound (CRLB). That is, if where { } and are scalars, then
More informationLessons in Estimation Theory for Signal Processing, Communications, and Control
Lessons in Estimation Theory for Signal Processing, Communications, and Control Jerry M. Mendel Department of Electrical Engineering University of Southern California Los Angeles, California PRENTICE HALL
More informationProbability Models for Bayesian Recognition
Intelligent Systems: Reasoning and Recognition James L. Crowley ENSIAG / osig Second Semester 06/07 Lesson 9 0 arch 07 Probability odels for Bayesian Recognition Notation... Supervised Learning for Bayesian
More informationA QUICK ASSESSMENT OF INPUT UNCERTAINTY. Bruce E. Ankenman Barry L. Nelson
Proceedings of the 2012 Winter Simulation Conference C. Laroque, J. Himmelspach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds. A QUICK ASSESSMENT OF INPUT UNCERTAINTY Bruce E. Ankenman Barry L. Nelson
More informationSTA 2201/442 Assignment 2
STA 2201/442 Assignment 2 1. This is about how to simulate from a continuous univariate distribution. Let the random variable X have a continuous distribution with density f X (x) and cumulative distribution
More informationIEOR 6711, HMWK 5, Professor Sigman
IEOR 6711, HMWK 5, Professor Sigman 1. Semi-Markov processes: Consider an irreducible positive recurrent discrete-time Markov chain {X n } with transition matrix P (P i,j ), i, j S, and finite state space.
More informationThe Nonparametric Bootstrap
The Nonparametric Bootstrap The nonparametric bootstrap may involve inferences about a parameter, but we use a nonparametric procedure in approximating the parametric distribution using the ECDF. We use
More informationMarkov chain Monte Carlo
Markov chain Monte Carlo Karl Oskar Ekvall Galin L. Jones University of Minnesota March 12, 2019 Abstract Practically relevant statistical models often give rise to probability distributions that are analytically
More informationON-LINE ERROR BOUNDS FOR STEADY-STATE APPROXIMATIONS: A POTENTIAL SOLUTION TO THE INITIALIZATION BIAS PROBLEM
Proceedings of the 2001 Winter Simulation Conference B. A. Peters, J. S. Smith, D. J. Medeiros, and M. W. Rohrer, eds. ON-LINE ERROR BOUNDS FOR STEADY-STATE APPROXIMATIONS: A POTENTIAL SOLUTION TO THE
More informationA Tour of Reinforcement Learning The View from Continuous Control. Benjamin Recht University of California, Berkeley
A Tour of Reinforcement Learning The View from Continuous Control Benjamin Recht University of California, Berkeley trustable, scalable, predictable Control Theory! Reinforcement Learning is the study
More informationProceedings of the 2014 Winter Simulation Conference A. Tolk, S. Y. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds.
Proceedings of the 2014 Winter Simulation Conference A. Tolk, S. Y. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds. BOOTSTRAP RANKING & SELECTION REVISITED Soonhui Lee School of Business
More informationAsymptotic distribution of the sample average value-at-risk
Asymptotic distribution of the sample average value-at-risk Stoyan V. Stoyanov Svetlozar T. Rachev September 3, 7 Abstract In this paper, we prove a result for the asymptotic distribution of the sample
More informationFast Estimation of the Statistics of Excessive Backlogs in Tandem Networks of Queues.
Fast Estimation of the Statistics of Excessive Backlogs in Tandem Networks of Queues. M.R. FRATER and B.D.O. ANDERSON Department of Systems Engineering Australian National University The estimation of
More informationACCOUNTING FOR INPUT-MODEL AND INPUT-PARAMETER UNCERTAINTIES IN SIMULATION. <www.ie.ncsu.edu/jwilson> May 22, 2006
ACCOUNTING FOR INPUT-MODEL AND INPUT-PARAMETER UNCERTAINTIES IN SIMULATION Slide 1 Faker Zouaoui Sabre Holdings James R. Wilson NC State University May, 006 Slide From American
More informationAsymptotics and Simulation of Heavy-Tailed Processes
Asymptotics and Simulation of Heavy-Tailed Processes Department of Mathematics Stockholm, Sweden Workshop on Heavy-tailed Distributions and Extreme Value Theory ISI Kolkata January 14-17, 2013 Outline
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.262 Discrete Stochastic Processes Midterm Quiz April 6, 2010 There are 5 questions, each with several parts.
More information2. Variance and Covariance: We will now derive some classic properties of variance and covariance. Assume real-valued random variables X and Y.
CS450 Final Review Problems Fall 08 Solutions or worked answers provided Problems -6 are based on the midterm review Identical problems are marked recap] Please consult previous recitations and textbook
More informationStatistics 150: Spring 2007
Statistics 150: Spring 2007 April 23, 2008 0-1 1 Limiting Probabilities If the discrete-time Markov chain with transition probabilities p ij is irreducible and positive recurrent; then the limiting probabilities
More informationN.G.Bean, D.A.Green and P.G.Taylor. University of Adelaide. Adelaide. Abstract. process of an MMPP/M/1 queue is not a MAP unless the queue is a
WHEN IS A MAP POISSON N.G.Bean, D.A.Green and P.G.Taylor Department of Applied Mathematics University of Adelaide Adelaide 55 Abstract In a recent paper, Olivier and Walrand (994) claimed that the departure
More informationMGR-815. Notes for the MGR-815 course. 12 June School of Superior Technology. Professor Zbigniew Dziong
Modeling, Estimation and Control, for Telecommunication Networks Notes for the MGR-815 course 12 June 2010 School of Superior Technology Professor Zbigniew Dziong 1 Table of Contents Preface 5 1. Example
More informationEfficient Rare-event Simulation for Perpetuities
Efficient Rare-event Simulation for Perpetuities Blanchet, J., Lam, H., and Zwart, B. We consider perpetuities of the form Abstract D = B 1 exp Y 1 ) + B 2 exp Y 1 + Y 2 ) +..., where the Y j s and B j
More informationProbability & Statistics - FALL 2008 FINAL EXAM
550.3 Probability & Statistics - FALL 008 FINAL EXAM NAME. An urn contains white marbles and 8 red marbles. A marble is drawn at random from the urn 00 times with replacement. Which of the following is
More informationIIT JAM : MATHEMATICAL STATISTICS (MS) 2013
IIT JAM : MATHEMATICAL STATISTICS (MS 2013 Question Paper with Answer Keys Ctanujit Classes Of Mathematics, Statistics & Economics Visit our website for more: www.ctanujit.in IMPORTANT NOTE FOR CANDIDATES
More informationSolution: (Course X071570: Stochastic Processes)
Solution I (Course X071570: Stochastic Processes) October 24, 2013 Exercise 1.1: Find all functions f from the integers to the real numbers satisfying f(n) = 1 2 f(n + 1) + 1 f(n 1) 1. 2 A secial solution
More informationTest Code: STA/STB (Short Answer Type) 2013 Junior Research Fellowship for Research Course in Statistics
Test Code: STA/STB (Short Answer Type) 2013 Junior Research Fellowship for Research Course in Statistics The candidates for the research course in Statistics will have to take two shortanswer type tests
More informationBandit Algorithms. Zhifeng Wang ... Department of Statistics Florida State University
Bandit Algorithms Zhifeng Wang Department of Statistics Florida State University Outline Multi-Armed Bandits (MAB) Exploration-First Epsilon-Greedy Softmax UCB Thompson Sampling Adversarial Bandits Exp3
More informationECE 3800 Probabilistic Methods of Signal and System Analysis
ECE 3800 Probabilistic Methods of Signal and System Analysis Dr. Bradley J. Bazuin Western Michigan University College of Engineering and Applied Sciences Department of Electrical and Computer Engineering
More informationLinear Regression. Machine Learning CSE546 Kevin Jamieson University of Washington. Oct 5, Kevin Jamieson 1
Linear Regression Machine Learning CSE546 Kevin Jamieson University of Washington Oct 5, 2017 1 The regression problem Given past sales data on zillow.com, predict: y = House sale price from x = {# sq.
More information