Online Supplement to Efficient Computation of Overlapping Variance Estimators for Simulation INFORMS Journal on Computing
|
|
- Sydney Spencer
- 5 years ago
- Views:
Transcription
1 Online Supplement to Efficient Computation of Overlapping Variance Estimators for Simulation INFORMS Journal on Computing Christos Alexopoulos H. Milton Stewart School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA, Nilay Tanık Argon Department of Industrial and Systems Engineering, University of Wisconsin Madison, Madison, Wisconsin , USA, David Goldsman H. Milton Stewart School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA, Natalie M. Steiger Maine Business School, University of Maine, Orono, Maine , USA, Gamze Tokol Decision Analytics, Atlanta, Georgia 30306, USA, James R. Wilson Edward P. Fitts Department of Industrial and Systems Engineering, North Carolina State University, Campus Box 7906, Raleigh, North Carolina , USA, The purpose of this online supplement is to present the following: (i linear time algorithms to obtain overlapping area and Cramér von Mises estimators; (ii experimental results on the bias, variance, meansquared error, and efficiency of the NBM and overlapping variance estimators; and (iii analytical results on the bias and variance of the overlapping estimators when they are applied to i.i.d. normal data.. Algorithms We start with the algorithms for computing overlapping area and CvM estimators with several weighting functions in O(n time. Figures and 2 present the algorithms for obtaining the overlapping area estimator A O (f ; b, m for weighting functions f r (t r p0 c pt p and f cos,j (t 8πj cos(2πjt, respectively. Figure 3 presents an algorithm to obtain the overlapping CvM estimator C O (g; b, m with weighting function g(t r p0 c pt p.
2 Figure : Algorithm for the Overlapping Area Estimator with Polynomial Weighting Functions Step : Initialization V 0, a 0 0, a 0, S (p 0 for p 0,...,r, and k Step 2: Calculate S (p 0,m for p 0,...,r; a 0(f r ; a (f r ; and V 0 (f r S (p S (p + k p Y k for p 0,...,r; V V f r (k/ms (0 Until k m + V V + S (0 a /m a 0 a 0 + f r (k/m, k k + a a + kf r (k/m Step 3: Calculate W i (f r, V i (f r for i,...,n m; accumulate A O (f r ; b, m for all overlapping batches A O V 2 /m 3 and i Until i n m + Return A O A O/ (n m + Calculate S (p i,m backwards in p to save storage space p r S (p p ( px x0 ( p x S (x + m p Y i+m Y i I {p0}, p p Until p W r p0 c p S (p /m p, V V + (Y i+m Y i a /m W + Y i a 0 A O A O + V 2 /m 3, i i + 2. Efficiency of the Variance Estimators In Alexopoulos et al. (2006, we summarize the results of an experimental performance evaluation of the new overlapping variance estimators; and we find that they perform efficaciously compared with their nonoverlapping counterparts at least with respect to the measures of bias and variance. In particular, there is little difference between corresponding nonoverlapping and overlapping estimators when it comes to bias; yet overlapping variance estimators often achieve substantial savings compared with nonoverlapping variance estimators when it comes to variance. We have also shown that all of the overlapping estimators studied herein require computation on the same order of magnitude as that of NBM, the simplest estimator in terms of computation. A natural question to ask is: just how equal are the O(n computational efforts required by the various variance estimators? We might reasonably expect that among these O(n estimators, NBM requires less work than any of the batched STS estimators, which in turn require less work than the corresponding overlapping STS estimators. The good news is that, for all practical purposes, the differences in computational effort 2
3 Figure 2: Algorithm for the Overlapping Area Estimator with Trigonometric Weighting Functions Step : Initialization j user-selected positive integer, V 0, a 0 0, a 0, Q 0, Q 0, S (0 0, and k Step 2: Calculate S (0 0,m, Q m, Q m, a 0(f cos,j, a (f cos,j, and V 0 (f cos,j S (0 S (0 + Y k, V V f cos,j (k/ms (0 Q Q + cos(α j ky k, Q Q + sin(α j ky k a 0 a 0 + f cos,j (k/m, a a + kf cos,j (k/m, k k + Until k m + V V + S (0 a /m Step 3: Calculate Q i, Q i, W i(f cos,j, V i (f cos,j for i,...,n m; accumulate A O (f cos,j ; b, m for all overlapping batches A O V 2 /m 3 and i Q Q + cos(α j (i + my i+m cos(α j iy i Q Q + sin(α j (i + my i+m sin(α j iy i W 8πj (cos(α j iq + sin(α j iq V V + (Y i+m Y i a /m W + Y i a 0 A O A O + V 2 /m 3, i i + Until i n m + Return A O A O /(n m + are meaningless, especially when one regards the observations as expensive compared to the subsequent computational costs once the observations are in hand. Nevertheless, we undertook a great deal of Monte Carlo work to evaluate the actual effort necessary to compute the various variance estimators. In one of our simplest experiments, we generated n 60,000 observations of the waiting time process for the M/M/ queue with 80% traffic intensity; and we organized the resulting time series into batches of size m 4,000, 8,000, and 6,000 so that we could obtain precise estimates of the total computing time C, bias B, variance V, and mean-squared error M of the corresponding variance estimators N (b, m, A O (f 0 ; b, m, A O (f 2 ; b, m, C O (g 0 ; b, m, and C O (g2 ; b, m for b 40, 20, and 0. Table reports the results based on 6,000 i.i.d. replications of each of the variance estimators under study. Each total computing-time entry C in Table is the sum of two components:. On an IBM T43p laptop computer with an Intel Pentium M processor having a speed of 2.3 GHz and 3
4 Figure 3: Algorithm for the Overlapping CvM Estimator with Polynomial Weighting Functions Step : Initialization S 0, M 0, J 0, a p 0 for p {0,, 2}, S (p 0 for p 0,...,r + Z (p 0 for p 0,...,r, l, and k Step 2: Calculate S (p 0,m for p 0,...,r +, U l for l 2,...,m+, M(, and a p (g for p {0,, 2}; initialize the second term of (4 from Alexopoulos et al. (2006 S (p S (p + k p Y k for p 0,...,r +, J J + g(k/mks (0 /m if l U 2 g(/my else Calculate Z (p l,,l backwards in p to save storage space p r Z (p p ( px x0 Z (x + Y l, p p Until p U l+ U l + Y l a 0 r p0 c p Z (p /m p + g(l/ms (0 endif a p a p + k p g(k/m for p {0,, 2}, M M + g(k/ms (0, l l +, k l Until l m + Step 3: Calculate W i (g r+ for i,...,n m, U l for l m + 2,...n, M(l for l 2,...,n m + ; accumulate the third term of (4 from Alexopoulos et al. (2006 and C O (g; b, m K M, G a 0, R Y M, C O a 2 (S (0 /m 2 2S (0 J, i Calculate S (p i,m backwards in p to save storage space p r + S (p p ( px x0 ( p x S (x + m p Y i+m Y i I {p0}, p p Until p W r p0 c p S (p+ /m p+, J J + W Y i a /m K K + r p0 c p S (p /m p a 0 Y l m, M M + K U l m+, R R + Y l m+ M if l n m + Calculate Z (p l,,m backwards in p to save storage space p r Z (p p ( px x0 Z (x + Y l (m + p Y l I {l m+2}, p p Until p U l+ U l + Y l a 0 r p0 c p Z (p /m p else G G g((l n + m 2/m, S S + Y l, p r Calculate Z (p l,l (n m,m backwards in p to save storage space Z (p p ( px x0 Z (x + Y l I {ln m+2} (m + p Y l (m+, p p Until p U l+ U l + Y l G r p0 c p Z (p /m p g((l n + m/ms endif C O C O + a 2 (S (0 /m 2 2S (0 J, l l +, i i + Until l n + l n m + 2 Step 4: Calculate M(l for l n m + 2,...,n; accumulate the third term of (4 and C O (g; b, m M M U l, R R + Y l M, l l + Until l n + Return C O (C O + R/(m 2 (n m + 4
5 Performance Table : Efficiency Analysis of the Variance Estimators Variance Estimator Measure N (b, m A O (f 0 ; b, m A O (f 2 ; b, m C O (g 0 ; b, m C O (g2 ; b, m Results for m 4,000 and b 40 C v.est (sec C (sec B V 430, ,27 332, , ,708 M 43, , ,62 25,74 326,96 Q Q Results for m 8,000 and b 20 C v.est (sec C (sec B V 637, , ,47 33,398 47,670 M 637,87 38,64 426,48 37,036 47,77 Q Q Results for m 6,000 and b 0 C v.est (sec C (sec B V,084, ,573 68,96 430, ,620 M,084, ,092 68,936 43, ,762 Q Q
6 .00 GB of RAM, we computed the average execution time C sim 4.22 seconds per run for an Arena simulation model (Kelton, Sadowski, and Sturrock 2004 of the M/M/ queue in which 60,000 customer waiting times are merely computed and tabulated in an Arena Record module before the associated entity is disposed. 2. On the same laptop computer, we computed the average execution time C v.est (in seconds to compute the variance estimators N (b, m, A O (f 0 ; b, m, A O (f 2 ; b, m, C O (g 0 ; b, m, and C O (g2 ; b, m from a set of n 60,000 queue waiting times for the M/M/ queue that were already resident in computer memory. To compute C v.est for, say, the NBM variance estimator N (b, m, we used the PORTLIB function DTIME of Fortran PowerStation (Microsoft Corp. 995 to obtain the total CPU time required to compute 6,000 replications of N (b, m; and then a similar experiment was performed for the overlapping variance estimators A O (f 0 ; b, m, A O (f 2 ; b, m, C O (g 0 ; b, m, and C O (g2 ; b, m separately. Thus for each variance estimator in Table, we report C sim, C v.est, and C C sim + C v.est. We take the NBM variance estimator N (b, m as the baseline with performance measures C 0, B 0, V 0, and M 0 ; and for any other variance estimator with corresponding performance measures C, B, V, and M, we see that the relative efficiency statistics are given by Q C V /( C0 V 0 and Q C M /( C0 M 0. ( On any reasonably sized practical applications, the estimator computation times are insignificant compared to the times required to obtain the observations, thus rendering NBM noncompetitive in terms of the relative efficiency measures Q and Q. 3. Exact Results for i.i.d. Normal Case The purpose here is simply to see whether the new estimators perform as advertised on an example for which we are able to calculate the expected value and variance of the variance estimators exactly. Suppose that we sample from an i.i.d. standard normal process, {Y i : i,...,n}. After some tedious algebra that is detailed below, we obtain the results displayed in Table 2 describing the bias and variance properties of some of the proposed overlapping variance estimators. In Table 2 all the bias results are exact; on the other hand, in the tabulated expressions for the variance of the overlapping estimators, we omitted terms of the form 6
7 O(/m (for fixed b. We see that the entries in Table 2 indeed match up nicely with those in Table of Alexopoulos et al. (2006 with γ 0. Of course, in this very special case, the natural estimator for σ 2 would have been the usual sample variance S 2, which is unbiased and has variance 2/(n. What follows are the derivations of the entries appearing in Table 2. Table 2: Bias and Variance for Variance Estimators in I.i.d. Normal Case Estimator Bias Variance A O (f 0 ; b, m m 2 24b 3 35(b 2 A O (f 2 ; b, m m 2 2m 4 m 6 354b (b 2 C O (g 0 ; b, m m 2 88b 5 20(b 2 C O (g 2 ; b, m 4 m 2 5 m 4 O(b, m b (b 2 4b 3 b 2 +4b+6 3(b 4 Overlapping Area Estimator We first obtain exact results for the overlapping area estimator. Using the arguments found in, e.g., 4 of Foley and Goldsman (999, one can apply a little algebraic elbow grease to obtain A O i (f ; m m ] 2 h(f ; j, my m 3 i+j, i,...,n m +, where h(f ; j,m m m l m f(l m f( l m, l lj and where we notice that the h(f ; j,m s need only be calculated once beforehand. We now calculate the mean and variance of the overlapping area estimator for the current example. First, 7
8 we give a general expression for the mean. Since the Y i s are i.i.d. with unit variance, we have EA O (f ; b, m] EA O (f ; m] ( m E m 2 ] h(f ; j, my 3 j ( m Var m h(f ; j, my 3 j m h 2 (f ; j, m. (2 m 3 Before moving ahead to the variance of the overlapping area estimator, we state a useful lemma. Lemma (Patel and Read 996 If (Z,Z 2 is bivariate normal with marginal means equal to zero, then Cov(Z 2,Z2 2 2Cov2 (Z,Z 2. In addition, define for k 0,, 2,..., R A k Cov(A O (f ; m, AO +k (f ; m Cov ( m 2 h(f ; i, my m 3 i, ( m 2 h(f ; j, my m 3 k+j 2 ( m m 6 Cov2 h(f ; i, my i, 2 m 6 m 2 m 6 m ik+ m h(f ; j, my k+j m h(f ; i, mh(f ; j,mcov(y i,y k+j ] 2 (by Lemma h(f ; i, mh(f ; i k, m] 2, (3 since the Y i s are i.i.d. Then the definition of A O (f ; b, m and stationarity imply that (n m + Var(A O (f ; b, m n m + k n m R0 A + 2 ( k n m + k R0 A + 2 ( k n m + again since the underlying observations are i.i.d. Cov(A O i (f ; m, AO j (f ; m R A k R A k, (4 8
9 Now we see what happens to the area and variance of the overlapping area estimator when we apply specific weighting functions to the i.i.d. standard normal process. Example For the overlapping area estimator with constant weighting function f 0 (t 2, it is easy to show that h(f 0 ; j,m ( 2 j m +. 2 Plugging this into Equations (2 (4 leads to EA O (f 0 ; b, m] m 2, which matches up with the appropriate entry in Table of Alexopoulos et al. (2006 (taking into account that γ 0, and Var(A O (f 0 ; b, m (m2 m 5 (24b m m 3 b + o(m 3 b ] 35m 5 m(b + ] 2. 24b 3 35(b 2, for large m and b, which is in synch with the appropriate variance entry from Table of Alexopoulos et al. (2006. Example 2 Consider the first-order unbiased overlapping area estimator with quadratic weights f 2 (t 840(3t 2 3t + /2. After some algebra, we have h(f 2 ; i, m 20 4i 3 6i 2 (m + + 2i(m 2 + 3m + m(m + ] 2m 2. Again using (2 (4, we eventually obtain and EA O (f 2 ; b, m] 2m6 + 7m m m 6, Var(A O (f 2 ; b, m 354m + O(/m 4290(bm m + 26m2 + O(. 354b (bm m (b, 2 for large m and b, in agreement with the results found in Table of Alexopoulos et al. (2006. Remark Some additional algebra not reported here allows us to verify the overlapping area estimator results for the AR( process analyzed via Monte Carlo in 4. of Alexopoulos et al. (
10 Overlapping CvM Estimator We next look at the overlapping CvM estimator for an i.i.d. standard normal process. As for the overlapping area estimator, we begin with a general expression for the expected value of the variance estimator. Since the Y i s are i.i.d. with unit variance, we have EC O (g; b, m] EC O (g; m] m g( k m 2 m k2 E(Ȳ,m O Ȳ,k O 2 ] k m g( k m 2 m k2 Var(Ȳ,m Ȳ,k k m g( k m 2 m k2 Var(Ȳ,m + Var(Ȳ,k 2Cov(Ȳ,m, Ȳ,k ] k m ( g( k m 2 m k2 m + k 2 m k m ( g( k m 2 m k2 k. (5 m k Before undertaking our calculation of the overlapping CvM estimator s variance, we start off with a preliminary calculation: For i, j,..., m and k 0,, 2,..., Cov ( Ȳ,m O Ȳ,i O, Ȳ +k,m O Ȳ +k,j O m m Cov(Y m 2 α,y k+β jm im α β i α β m Cov(Y α,y k+β + ij m α β i α β j Cov(Y α,y k+β j Cov(Y α,y k+β max(m k, 0 min(j, max(m k, 0 max(i k, 0 + min(j, max(i k, 0 m2 jm im ij (since the Y i s are i.i.d. with unit variance ( + k m( i m, if k<iand j i k k + ( ( k m 2 i j m k ( i ( m j, if k<iand m k<j m, if k<iand i k<j m k k, if i k<mand j m k m 2 ( ( k m j, if i k<mand m k<j m 0, if m k. (6 0
11 For k 0,, 2,..., let Rk C Cov(C O (g; m, CO +k (g; m Cov m g( i m 2 m i2( Ȳ,m O Ȳ,i O 2, m g( j m 2 m j 2( Ȳ+k,m O Ȳ +k,j O 2 2 g( i m 4 m g( j m i2 j 2 Cov (Ȳ 2,m O Ȳ,i O, Ȳ +k,m O Ȳ +k,j O, the last step following from Lemma. Evocation of Equation (6 then reveals that Rk C 2 k 2 k m k g( i m 4 m g( j m i2 j 2 + ( k 2 k ( g( i m m g( j m i2 j m m 4 + ( + k m + 2 ik+ ji k+ +k 2 ik+ m k ik+ jm k+ i k g( i m g( j m j ( 2 i m jm k+ 2 g( i m g( j m i2 j 2 ( k m 2 + ( k i g( i m g( j m ( i m ( 2 j m 2 ( j m 2, (7 0 for k m. Finally, as in the derivation of Equation (4, we have for k 0,,...,m, and Rk C (n m + Var(C O (g; b, m R0 C + 2 ( k 2 k Rk C n m +. (8 With Equations (5, (7, and (8 in hand, we can calculate the mean and variance of the overlapping CvM estimator when specific weighting functions are applied to the i.i.d. normal process. Example 3 For the overlapping CvM estimator with constant weighting function g 0 (t 6, we obtain EC O (g 0 ; b, m] m 2, which is consistent with Table of Alexopoulos et al. (2006 (with γ 0, and Var(C O (g 0 ; b, m (88b 5m7 + 88m 6 + m 5 (588b O(m 4 b 20m 5 (m(b b 5 20(b 2, for large m and b, which is again consistent with Table of Alexopoulos et al. (2006. Example 4 For the overlapping CvM estimator with quadratic weighting function g2 (t t 50t 2,wehave EC O (g 2 ; b, m] + 4 m 2 5 m 4,
12 which agrees with Table of Alexopoulos et al. (2006 (with γ 0; and Var(C O (g2 ; b, m m0 (bm m + O(m 9 b m 9 (m(b b (b, 2 for large m and b, which agrees with Table of Alexopoulos et al. (2006. Overlapping Batch Means Estimator Finally, we examine the OBM estimator for an i.i.d. standard normal process. The expected value is given by (n m + (n m nm EO(b, m] E(Ȳ O i,m Ȳ n 2 ] Var(Ȳ O i,m Ȳ n ] Var(Ȳi,m O 2Cov(Ȳ i,m O, Ȳ n + Var(Ȳ n (n m + m 2 n + ], n so that EO(b, m], in agreement with Table of Alexopoulos et al. (2006 (with γ 0. As for the variance of the OBM estimator, we have (n m + 2 (n m 2 Var(O(b, m 2n 2 m 2 ( Cov (Ȳi,m O 2 Ȳ n 2,(Ȳj,m O Ȳ n 2 Cov 2 (Ȳ O i,m Ȳ n, Ȳ O j,m Ȳ n (by Lemma Cov(Ȳi,m O, Ȳ j,m O Cov(Ȳ i,m O, Ȳ n Cov(Ȳj,m O, Ȳ n + Var(Ȳ n Cov(Ȳi,m O, Ȳ j,m O n (n m + Var(Ȳ,m n] O ] 2 (n m + m ] 2 m j + 2 (n m + j n m 2 ] 2 (n m + j Cov(Ȳ,m O, Ȳ +j,m O n 2 n ] n 2 jm ] 2 (n m + j,
13 so that, after carrying out some algebra, Var(O(b, m 4b3 b 2 + 4b + 6 2b 2 + 3(b 4 3m(b 2b(b4 3b 3 + 5b 2 0b (n m + (b 4 + b2 (b 2 6b + 8 3(n m + 2 (b 4 4b3 b 2 + 4b + 6 3(b 4 + O(/m, which is in agreement with Table of Alexopoulos et al. (2006. References Alexopoulos, C., N. T. Argon, D. Goldsman, N. M. Steiger, G. Tokol, and J. R. Wilson Efficient computation of overlapping variance estimators for simulation. Technical report. School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, GA. Foley, R. D., and D. Goldsman Confidence intervals using orthonormally weighted standardized times series. ACM Transactions on Modeling and Computer Simulation Kelton, W. D., R. P. Sadowski, and D. T. Sturrock Simulation with Arena, 3rd ed. McGraw-Hill, New York. Microsoft Corp Fortran PowerStation Reference. Microsoft Corp., Redmond, WA. Patel, J. K., and C. B. Read Handbook of the Normal Distribution, 2nd ed. Marcel Dekker, New York. 3
Variance Parameter Estimation Methods with Data Re-use
Variance Parameter Estimation Methods with Data Re-use Christos Alexopoulos, David Goldsman, Melike Meterelliyoz Claudia Antonini James R. Wilson Georgia Institute of Technology, Atlanta, GA, USA Universidad
More informationA SEQUENTIAL PROCEDURE FOR ESTIMATING THE STEADY-STATE MEAN USING STANDARDIZED TIME SERIES. Christos Alexopoulos David Goldsman Peng Tang
Proceedings of the 2013 Winter Simulation Conference R. Pasupathy, S.-H. Kim, A. Tolk, R. Hill, and M. E. Kuhl, eds. A SEQUENTIAL PROCEDURE FOR ESTIMATING THE STEADY-STATE MEAN USING STANDARDIZED TIME
More informationFolded standardized time series area variance estimators for simulation
Submitted to IIE Transactions. Folded standardized time series area variance estimators for simulation CLAUDIA ANTONINI 1 Departamento de Matemáticas Puras y Aplicadas, Universidad Simón Bolivar, Sartenejas,
More informationArea variance estimators for simulation using folded standardized time series
To appear, IIE Transactions. Area variance estimators for simulation using folded standardized time series CLAUDIA ANTONINI 1 Departamento de Matemáticas Puras y Aplicadas, Universidad Simón Bolívar, Sartenejas,
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 informationARD: AN AUTOMATED REPLICATION-DELETION METHOD FOR SIMULATION ANALYSIS
Proceedings of the 2013 Winter Simulation Conference R. Pasupathy, S.-H. Kim, A. Tolk, R. Hill, and M. E. Kuhl, eds. ARD: AN AUTOMATED REPLICATION-DELETION METHOD FOR SIMULATION ANALYSIS Emily K. Lada
More informationSIMULATION OUTPUT ANALYSIS
1 / 64 SIMULATION OUTPUT ANALYSIS Dave Goldsman School of ISyE Georgia Tech Atlanta, Georgia, USA sman@gatech.edu www.isye.gatech.edu/ sman April 15, 2016 Outline 1 Introduction 2 Finite-Horizon Simulation
More informationOptimal Linear Combinations of Overlapping Variance Estimators for Steady-State Simulation
Optimal Linear Combinations of Overlapping Variance Estimators for Steady-State Simulation Tûba Aktaran-Kalaycı, Christos Alexopoulos, David Goldsman, and James R. Wilson Abstract To estimate the variance
More informationProceedings of the 2008 Winter Simulation Conference S. J. Mason, R. R. Hill, L. Mönch, O. Rose, T. Jefferson, J. W. Fowler eds.
Proceedings of the 2008 Winter Simulation Conference S. J. Mason, R. R. Hill, L. Mönch, O. Rose, T. Jefferson, J. W. Fowler eds. A DISTRIBUTION-FREE TABULAR CUSUM CHART FOR CORRELATED DATA WITH AUTOMATED
More informationSIMULATION OUTPUT ANALYSIS
1 / 71 SIMULATION OUTPUT ANALYSIS Dave Goldsman School of ISyE Georgia Tech Atlanta, Georgia, USA sman@gatech.edu www.isye.gatech.edu/ sman 10/19/17 2 / 71 Outline 1 Introduction 2 A Mathematical Interlude
More informationARD: AN AUTOMATED REPLICATION-DELETION METHOD FOR SIMULATION ANALYSIS
ARD: AN AUTOMATED REPLICATION-DELETION METHOD FOR SIMULATION ANALYSIS Emily K. Lada Anup C. Mokashi SAS Institute Inc. 100 SAS Campus Drive, R5413 Cary, NC 27513-8617, USA James R. Wilson North Carolina
More informationMonitoring Autocorrelated Processes Using A Distribution-Free Tabular CUSUM Chart With Automated Variance Estimation
To appear, IIE Transactions. Monitoring Autocorrelated Processes Using A Distribution-Free Tabular CUSUM Chart With Automated Variance Estimation JOONGSUP (JAY) LEE 1 CHRISTOS ALEXOPOULOS 2 DAVID GOLDSMAN
More informationAsummary and an analysis are given for an experimental performance evaluation of WASSP, an automated
INFORMS Journal on Computing Vol. 19, No. 2, Spring 2007, pp. 150 160 issn 1091-9856 eissn 1526-5528 07 1902 0150 informs doi 10.1287/ijoc.1050.0161 2007 INFORMS Performance of a Wavelet-Based Spectral
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. SEQUEST: A SEQUENTIAL PROCEDURE FOR ESTIMATING STEADY-STATE QUANTILES
More informationA New Model-Free CuSum Procedure for Autocorrelated Processes
A New Model-Free CuSum Procedure for Autocorrelated Processes Seong-Hee Kim, Christos Alexopoulos, David Goldsman, and Kwok-Leung Tsui School of Industrial and Systems Engineering Georgia Institute of
More informationISyE 6644 Fall 2014 Test 3 Solutions
1 NAME ISyE 6644 Fall 14 Test 3 Solutions revised 8/4/18 You have 1 minutes for this test. You are allowed three cheat sheets. Circle all final answers. Good luck! 1. [4 points] Suppose that the joint
More informationMonitoring autocorrelated processes using a distribution-free tabular CUSUM chart with automated variance estimation
IIE Transactions (2009) 41, 979 994 Copyright C IIE ISSN: 0740-817X print / 1545-8830 online DOI: 10.1080/07408170902906035 Monitoring autocorrelated processes using a distribution-free tabular CUSUM chart
More informationESTIMATION AND OUTPUT ANALYSIS (L&K Chapters 9, 10) Review performance measures (means, probabilities and quantiles).
ESTIMATION AND OUTPUT ANALYSIS (L&K Chapters 9, 10) Set up standard example and notation. Review performance measures (means, probabilities and quantiles). A framework for conducting simulation experiments
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 information7 Variance Reduction Techniques
7 Variance Reduction Techniques In a simulation study, we are interested in one or more performance measures for some stochastic model. For example, we want to determine the long-run average waiting time,
More informationA SPACED BATCH MEANS PROCEDURE FOR SIMULATION ANALYSIS
Proceedings of the 2007 Winter Simulation Conference S. G. Henderson, B. Biller, M.-H. Hsieh, J. Shortle, J. D. Tew, and R. R. Barton, eds. SBatch: A SPACED BATCH MEANS PROCEDURE FOR SIMULATION ANALYSIS
More informationPerformance of Skart: A Skewness- and Autoregression-Adjusted Batch-Means Procedure for Simulation Analysis
INFORMS Journal on Computing, to appear. Performance of Skart: A Skewness- and Autoregression-Adjusted Batch-Means Procedure for Simulation Analysis Ali Tafazzoli, James R. Wilson Edward P. Fitts Department
More informationEconometrics Summary Algebraic and Statistical Preliminaries
Econometrics Summary Algebraic and Statistical Preliminaries Elasticity: The point elasticity of Y with respect to L is given by α = ( Y/ L)/(Y/L). The arc elasticity is given by ( Y/ L)/(Y/L), when L
More informationLecture 6: Linear models and Gauss-Markov theorem
Lecture 6: Linear models and Gauss-Markov theorem Linear model setting Results in simple linear regression can be extended to the following general linear model with independently observed response variables
More informationA distribution-free tabular CUSUM chart for autocorrelated data
IIE Transactions (007) 39, 317 330 Copyright C IIE ISSN: 0740-817X print / 1545-8830 online DOI: 10.1080/07408170600743946 A distribution-free tabular CUSUM chart for autocorrelated data SEONG-HEE KIM
More informationPERFORMANCE COMPARISON OF MSER-5 AND N-SKART ON THE SIMULATION START-UP PROBLEM
Proceedings of the 2010 Winter Simulation Conference B. Johansson, S. Jain, J. Montoya-Torres, J. Hugan, and E. Yücesan, eds. PERFORMACE COMPARISO OF MSER-5 AD -SKART O THE SIMULATIO START-UP PROBLEM Anup
More informationN-SKART: A NONSEQUENTIAL SKEWNESS- AND AUTOREGRESSION-ADJUSTED BATCH-MEANS PROCEDURE FOR SIMULATION ANALYSIS
Proceedings of the 2009 Winter Simulation Conference M. D. Rossetti, R. R. Hill, B. Johansson, A. Dunkin, and R. G. Ingalls, eds. N-SKART: A NONSEQUENTIAL SKEWNESS- AND AUTOREGRESSION-ADJUSTED BATCH-MEANS
More informationMore on Input Distributions
More on Input Distributions Importance of Using the Correct Distribution Replacing a distribution with its mean Arrivals Waiting line Processing order System Service mean interarrival time = 1 minute mean
More informationA Randomized Algorithm for the Approximation of Matrices
A Randomized Algorithm for the Approximation of Matrices Per-Gunnar Martinsson, Vladimir Rokhlin, and Mark Tygert Technical Report YALEU/DCS/TR-36 June 29, 2006 Abstract Given an m n matrix A and a positive
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 informationSolution manual of Simulation with Arena 6th Edition by W. David Kelton, Randall Sadowski and Nancy Zupick
Solution manual of Simulation with Arena 6th Edition by W. David Kelton, Randall Sadowski and Nancy Zupick Link download full: http://testbankair.com/download/simulation-with-arena-6th-edition-solutions-by-kelton/
More informationCourse information: Instructor: Tim Hanson, Leconte 219C, phone Office hours: Tuesday/Thursday 11-12, Wednesday 10-12, and by appointment.
Course information: Instructor: Tim Hanson, Leconte 219C, phone 777-3859. Office hours: Tuesday/Thursday 11-12, Wednesday 10-12, and by appointment. Text: Applied Linear Statistical Models (5th Edition),
More informationSimple Linear Regression
Simple Linear Regression In simple linear regression we are concerned about the relationship between two variables, X and Y. There are two components to such a relationship. 1. The strength of the relationship.
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationQueueing Review. Christos Alexopoulos and Dave Goldsman 10/25/17. (mostly from BCNN) Georgia Institute of Technology, Atlanta, GA, USA
1 / 26 Queueing Review (mostly from BCNN) Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/25/17 2 / 26 Outline 1 Introduction 2 Queueing Notation 3 Transient
More informationQueueing Review. Christos Alexopoulos and Dave Goldsman 10/6/16. (mostly from BCNN) Georgia Institute of Technology, Atlanta, GA, USA
1 / 24 Queueing Review (mostly from BCNN) Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/6/16 2 / 24 Outline 1 Introduction 2 Queueing Notation 3 Transient
More informationQUANTILE AND HISTOGRAM ESTIMATION
Proceedings of the 2001 Winter Simulation Conference B.A.Peters,J.S.Smith,D.J.Medeiros,andM.W.Rohrer,eds. QUANTILE AND HISTOGRAM ESTIMATION E. Jack Chen Mathematical Modeling Group BASF Corporation 3000
More informationHand Simulations. Christos Alexopoulos and Dave Goldsman 9/2/14. Georgia Institute of Technology, Atlanta, GA, USA 1 / 26
1 / 26 Hand Simulations Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 9/2/14 2 / 26 Outline 1 Monte Carlo Integration 2 Making Some π 3 Single-Server Queue 4
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 Siulation Conference L. Yilaz, W. K. V. Chan, I. Moon, T. M. K. Roeder, C. Macal, and M. D. Rossetti, eds. JACKKNIFED VARIANCE ESTIMATORS FOR SIMULATION OUTPUT ANALYSIS Keal
More informationTHE ROYAL STATISTICAL SOCIETY 2009 EXAMINATIONS SOLUTIONS GRADUATE DIPLOMA MODULAR FORMAT MODULE 3 STOCHASTIC PROCESSES AND TIME SERIES
THE ROYAL STATISTICAL SOCIETY 9 EXAMINATIONS SOLUTIONS GRADUATE DIPLOMA MODULAR FORMAT MODULE 3 STOCHASTIC PROCESSES AND TIME SERIES The Society provides these solutions to assist candidates preparing
More informationSPC Monitoring and Variance Estimation
SPC Monitoring and Variance Estimation C. Alexopoulos D. Goldsman K.-L. Tsui School of Industrial and Systems Engineering Georgia Institute of Technology Atlanta, GA 30332-0205 W. Jiang INSIGHT, AT&T Labs
More informationCOMPARING SYSTEMS. Dave Goldsman. November 15, School of ISyE Georgia Institute of Technology Atlanta, GA, USA
1 / 103 COMPARING SYSTEMS Dave Goldsman School of ISyE Georgia Institute of Technology Atlanta, GA, USA sman@gatech.edu November 15, 2017 Outline 1 Introduction and Review of Classical Confidence Intervals
More informationConjugate Gradient (CG) Method
Conjugate Gradient (CG) Method by K. Ozawa 1 Introduction In the series of this lecture, I will introduce the conjugate gradient method, which solves efficiently large scale sparse linear simultaneous
More informationDIMACS Workshop on Parallelism: A 2020 Vision Lattice Basis Reduction and Multi-Core
DIMACS Workshop on Parallelism: A 2020 Vision Lattice Basis Reduction and Multi-Core Werner Backes and Susanne Wetzel Stevens Institute of Technology 29th March 2011 Work supported through NSF Grant DUE
More informationInitial Transient Phase of Steady state Simulation: Methods of Its Length Detection and Their Evaluation in Akaroa2
Initial Transient Phase of Steady state Simulation: Methods of Its Length Detection and Their Evaluation in Akaroa2 A Thesis Submitted in Partial Fulfilment of the Requirement for the Degree of Master
More informationWeb-based Supplementary Materials for A Robust Method for Estimating. Optimal Treatment Regimes
Biometrics 000, 000 000 DOI: 000 000 0000 Web-based Supplementary Materials for A Robust Method for Estimating Optimal Treatment Regimes Baqun Zhang, Anastasios A. Tsiatis, Eric B. Laber, and Marie Davidian
More informationSTATISTICAL SCREENING, SELECTION, AND MULTIPLE COMPARISON PROCEDURES IN COMPUTER SIMULATION
Proceedings of the 1998 Winter Simulation Conference D.J. Medeiros E.F. Watson J.S. Carson and M.S. Manivannan eds. STATISTICAL SCREENING SELECTION AND MULTIPLE COMPARISON PROCEDURES IN COMPUTER SIMULATION
More informationA Multivariate Two-Sample Mean Test for Small Sample Size and Missing Data
A Multivariate Two-Sample Mean Test for Small Sample Size and Missing Data Yujun Wu, Marc G. Genton, 1 and Leonard A. Stefanski 2 Department of Biostatistics, School of Public Health, University of Medicine
More informationSTAT 520: Forecasting and Time Series. David B. Hitchcock University of South Carolina Department of Statistics
David B. University of South Carolina Department of Statistics What are Time Series Data? Time series data are collected sequentially over time. Some common examples include: 1. Meteorological data (temperatures,
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 informationA QUICKER ASSESSMENT OF INPUT UNCERTAINTY. Eunhye Song Barry L. Nelson
Proceedings of the 2013 Winter Simulation Conference R. Pasupathy, S.-H. Kim, A. Tolk, R. Hill, and M. E. Kuhl, eds. A QUICKER ASSESSMENT OF INPUT UNCERTAINTY Eunhye Song Barry. Nelson Department of Industrial
More informationSTAT 100C: Linear models
STAT 100C: Linear models Arash A. Amini June 9, 2018 1 / 56 Table of Contents Multiple linear regression Linear model setup Estimation of β Geometric interpretation Estimation of σ 2 Hat matrix Gram matrix
More informationMarkov Chain Model for ALOHA protocol
Markov Chain Model for ALOHA protocol Laila Daniel and Krishnan Narayanan April 22, 2012 Outline of the talk A Markov chain (MC) model for Slotted ALOHA Basic properties of Discrete-time Markov Chain Stability
More informationGeneralized Network Flow Techniques for Dynamic Voltage Scaling in Hard Real-Time Systems
Generalized Network Flow Techniques for Dynamic Voltage Scaling in Hard Real-Time Systems Vishnu Swaminathan and Krishnendu Chakrabarty Department of Electrical & Computer Engineering Duke University Durham,
More informationPart 4: Multi-parameter and normal models
Part 4: Multi-parameter and normal models 1 The normal model Perhaps the most useful (or utilized) probability model for data analysis is the normal distribution There are several reasons for this, e.g.,
More informationMA 575 Linear Models: Cedric E. Ginestet, Boston University Revision: Probability and Linear Algebra Week 1, Lecture 2
MA 575 Linear Models: Cedric E Ginestet, Boston University Revision: Probability and Linear Algebra Week 1, Lecture 2 1 Revision: Probability Theory 11 Random Variables A real-valued random variable is
More informationOn Expected Gaussian Random Determinants
On Expected Gaussian Random Determinants Moo K. Chung 1 Department of Statistics University of Wisconsin-Madison 1210 West Dayton St. Madison, WI 53706 Abstract The expectation of random determinants whose
More informationDiscrete Mathematics and Probability Theory Spring 2016 Rao and Walrand Note 26. Estimation: Regression and Least Squares
CS 70 Discrete Mathematics and Probability Theory Spring 2016 Rao and Walrand Note 26 Estimation: Regression and Least Squares This note explains how to use observations to estimate unobserved random variables.
More informationExpectations and Variance
4. Model parameters and their estimates 4.1 Expected Value and Conditional Expected Value 4. The Variance 4.3 Population vs Sample Quantities 4.4 Mean and Variance of a Linear Combination 4.5 The Covariance
More informationSlides 12: Output Analysis for a Single Model
Slides 12: Output Analysis for a Single Model Objective: Estimate system performance via simulation. If θ is the system performance, the precision of the estimator ˆθ can be measured by: The standard error
More informationEE/PEP 345. Modeling and Simulation. Spring Class 11
EE/PEP 345 Modeling and Simulation Class 11 11-1 Output Analysis for a Single Model Performance measures System being simulated Output Output analysis Stochastic character Types of simulations Output analysis
More information[y i α βx i ] 2 (2) Q = i=1
Least squares fits This section has no probability in it. There are no random variables. We are given n points (x i, y i ) and want to find the equation of the line that best fits them. We take the equation
More informationCPSC 531: System Modeling and Simulation. Carey Williamson Department of Computer Science University of Calgary Fall 2017
CPSC 531: System Modeling and Simulation Carey Williamson Department of Computer Science University of Calgary Fall 2017 Quote of the Day A person with one watch knows what time it is. A person with two
More informationClosed-Form Estimators for the Gamma Distribution Derived from Likelihood Equations
The American Statistician ISSN: 3-135 (Print) 1537-2731 (Online) Journal homepage: http://www.tandfonline.com/loi/utas2 Closed-Form Estimators for the Gamma Distribution Derived from Likelihood Equations
More informationProperties of Summation Operator
Econ 325 Section 003/004 Notes on Variance, Covariance, and Summation Operator By Hiro Kasahara Properties of Summation Operator For a sequence of the values {x 1, x 2,..., x n, we write the sum of x 1,
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 informationEconometrics of financial markets, -solutions to seminar 1. Problem 1
Econometrics of financial markets, -solutions to seminar 1. Problem 1 a) Estimate with OLS. For any regression y i α + βx i + u i for OLS to be unbiased we need cov (u i,x j )0 i, j. For the autoregressive
More informationGenerating Uniform Random Numbers
1 / 44 Generating Uniform Random Numbers Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/29/17 2 / 44 Outline 1 Introduction 2 Some Lousy Generators We Won t
More informationCHAPTER 2 EXTRACTION OF THE QUADRATICS FROM REAL ALGEBRAIC POLYNOMIAL
24 CHAPTER 2 EXTRACTION OF THE QUADRATICS FROM REAL ALGEBRAIC POLYNOMIAL 2.1 INTRODUCTION Polynomial factorization is a mathematical problem, which is often encountered in applied sciences and many of
More information2.1 Linear regression with matrices
21 Linear regression with matrices The values of the independent variables are united into the matrix X (design matrix), the values of the outcome and the coefficient are represented by the vectors Y and
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 informationTopics in Probability and Statistics
Topics in Probability and tatistics A Fundamental Construction uppose {, P } is a sample space (with probability P), and suppose X : R is a random variable. The distribution of X is the probability P X
More informationOnline publication date: 22 December 2010 PLEASE SCROLL DOWN FOR ARTICLE
This article was downloaded by: Wilson, James R.] On: 22 December 2010 Access details: Access Details: subscription number 931526073] Publisher Taylor & Francis Informa Ltd Registered in England and Wales
More informationReminders. Thought questions should be submitted on eclass. Please list the section related to the thought question
Linear regression Reminders Thought questions should be submitted on eclass Please list the section related to the thought question If it is a more general, open-ended question not exactly related to a
More informationProblem Set 1 Solution Sketches Time Series Analysis Spring 2010
Problem Set 1 Solution Sketches Time Series Analysis Spring 2010 1. Construct a martingale difference process that is not weakly stationary. Simplest e.g.: Let Y t be a sequence of independent, non-identically
More informationST 371 (IX): Theories of Sampling Distributions
ST 371 (IX): Theories of Sampling Distributions 1 Sample, Population, Parameter and Statistic The major use of inferential statistics is to use information from a sample to infer characteristics about
More informationADJUSTED POWER ESTIMATES IN. Ji Zhang. Biostatistics and Research Data Systems. Merck Research Laboratories. Rahway, NJ
ADJUSTED POWER ESTIMATES IN MONTE CARLO EXPERIMENTS Ji Zhang Biostatistics and Research Data Systems Merck Research Laboratories Rahway, NJ 07065-0914 and Dennis D. Boos Department of Statistics, North
More informationX i. X(n) = 1 n. (X i X(n)) 2. S(n) n
Confidence intervals Let X 1, X 2,..., X n be independent realizations of a random variable X with unknown mean µ and unknown variance σ 2. Sample mean Sample variance X(n) = 1 n S 2 (n) = 1 n 1 n i=1
More informationMA 575 Linear Models: Cedric E. Ginestet, Boston University Non-parametric Inference, Polynomial Regression Week 9, Lecture 2
MA 575 Linear Models: Cedric E. Ginestet, Boston University Non-parametric Inference, Polynomial Regression Week 9, Lecture 2 1 Bootstrapped Bias and CIs Given a multiple regression model with mean and
More informationAbsolute Value Programming
O. L. Mangasarian Absolute Value Programming Abstract. We investigate equations, inequalities and mathematical programs involving absolute values of variables such as the equation Ax + B x = b, where A
More informationE-Companion to Fully Sequential Procedures for Large-Scale Ranking-and-Selection Problems in Parallel Computing Environments
E-Companion to Fully Sequential Procedures for Large-Scale Ranking-and-Selection Problems in Parallel Computing Environments Jun Luo Antai College of Economics and Management Shanghai Jiao Tong University
More informationProceedings of the 2016 Winter Simulation Conference T. M. K. Roeder, P. I. Frazier, R. Szechtman, E. Zhou, T. Huschka, and S. E. Chick, eds.
Proceedings of the 016 Winter Simulation Conference T. M. K. Roeder, P. I. Frazier, R. Szechtman, E. Zhou, T. Huschka, and S. E. Chick, eds. SPEEDING UP PAIRWISE COMPARISONS FOR LARGE SCALE RANKING AND
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 informationNew Negative Latin Square Type Partial Difference Sets in Nonelementary Abelian 2-groups and 3-groups
New Negative Latin Square Type Partial Difference Sets in Nonelementary Abelian 2-groups and 3-groups John Polhill Department of Mathematics, Computer Science, and Statistics Bloomsburg University Bloomsburg,
More informationNonresponse weighting adjustment using estimated response probability
Nonresponse weighting adjustment using estimated response probability Jae-kwang Kim Yonsei University, Seoul, Korea December 26, 2006 Introduction Nonresponse Unit nonresponse Item nonresponse Basic strategy
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 informationNikolai K. Krivulin. St.Petersburg State University
A MAX-ALGEBRA APPROACH TO MODELING AND SIMULATION OF TANDEM QUEUEING SYSTEMS Nikolai K Krivulin StPetersburg State University Abstract Max-algebra models of tandem single-server queueing systems with both
More informationExpectations and Variance
4. Model parameters and their estimates 4.1 Expected Value and Conditional Expected Value 4. The Variance 4.3 Population vs Sample Quantities 4.4 Mean and Variance of a Linear Combination 4.5 The Covariance
More informationCosmology & CMB. Set5: Data Analysis. Davide Maino
Cosmology & CMB Set5: Data Analysis Davide Maino Gaussian Statistics Statistical isotropy states only two-point correlation function is needed and it is related to power spectrum Θ(ˆn) = lm Θ lm Y lm (ˆn)
More informationGenerating Uniform Random Numbers
1 / 41 Generating Uniform Random Numbers Christos Alexopoulos and Dave Goldsman Georgia Institute of Technology, Atlanta, GA, USA 10/13/16 2 / 41 Outline 1 Introduction 2 Some Lousy Generators We Won t
More informationENEE 621 SPRING 2016 DETECTION AND ESTIMATION THEORY THE PARAMETER ESTIMATION PROBLEM
c 2007-2016 by Armand M. Makowski 1 ENEE 621 SPRING 2016 DETECTION AND ESTIMATION THEORY THE PARAMETER ESTIMATION PROBLEM 1 The basic setting Throughout, p, q and k are positive integers. The setup With
More informationA Single-Pass Algorithm for Spectrum Estimation With Fast Convergence Han Xiao and Wei Biao Wu
4720 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL 57, NO 7, JULY 2011 A Single-Pass Algorithm for Spectrum Estimation With Fast Convergence Han Xiao Wei Biao Wu Abstract We propose a single-pass algorithm
More informationUC Berkeley Department of Electrical Engineering and Computer Sciences. EECS 126: Probability and Random Processes
UC Berkeley Department of Electrical Engineering and Computer Sciences EECS 6: Probability and Random Processes Problem Set 3 Spring 9 Self-Graded Scores Due: February 8, 9 Submit your self-graded scores
More information7.1 Coupling from the Past
Georgia Tech Fall 2006 Markov Chain Monte Carlo Methods Lecture 7: September 12, 2006 Coupling from the Past Eric Vigoda 7.1 Coupling from the Past 7.1.1 Introduction We saw in the last lecture how Markov
More informationReview: Second Half of Course Stat 704: Data Analysis I, Fall 2014
Review: Second Half of Course Stat 704: Data Analysis I, Fall 2014 Tim Hanson, Ph.D. University of South Carolina T. Hanson (USC) Stat 704: Data Analysis I, Fall 2014 1 / 13 Chapter 8: Polynomials & Interactions
More informationProblem Y is an exponential random variable with parameter λ = 0.2. Given the event A = {Y < 2},
ECE32 Spring 25 HW Solutions April 6, 25 Solutions to HW Note: Most of these solutions were generated by R. D. Yates and D. J. Goodman, the authors of our textbook. I have added comments in italics where
More informationGradient-based Adaptive Stochastic Search
1 / 41 Gradient-based Adaptive Stochastic Search Enlu Zhou H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of Technology November 5, 2014 Outline 2 / 41 1 Introduction
More informationDISPLAYING STATISTICAL POINT ESTIMATORS: THE LEADING-DIGIT PROCEDURE
Proceedings of the 2008 Winter Simulation Conference S. J. Mason, R. R. Hill, L. Mönch, O. Rose, T. Jefferson, J. W. Fowler eds. DISPLAYING STATISTICAL POINT ESTIMATORS: THE LEADING-DIGIT PROCEDURE Wheyming
More information14 Random Variables and Simulation
14 Random Variables and Simulation In this lecture note we consider the relationship between random variables and simulation models. Random variables play two important roles in simulation models. We assume
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 information