Proceedings of the 2014 Winter Simulation Conference A. Tolk, S. Y. Diallo, I. O. Ryzhov, L. Yilmaz, S. Buckley, and J. A. Miller, eds.
|
|
- Oswin Benson
- 5 years ago
- Views:
Transcription
1 Proceedings of the 2014 Winter Simuation Conference A. Tok, S. Y. Diao, I. O. Ryzhov, L. Yimaz, S. Buckey, and J. A. Mier, eds. STATISTICAL UNCERTAINTY ANALYSIS FOR STOCHASTIC SIMULATION WITH DEPENDENT INPUT MODELS Wei Xie Rensseaer Poytechnic Institute Troy, NY , USA Barry L. Neson Northwestern University Evanston, IL 60208, USA Russe R. Barton The Pennsyvania State University University Park, PA 16802, USA ABSTRACT When we use simuation to estimate the performance of a stochastic system, ack of fideity in the random input modes can ead to poor system performance estimates. Since the components of many compex systems coud be dependent, we want to buid input modes that faithfuy capture such key properties. In this paper, we use the fexibe NORma To Anything (NORTA) representation for dependent inputs. However, to use the NORTA representation we need to estimate the margina distribution parameters and a correation matrix from rea-word data, introducing input uncertainty. To quantify this uncertainty, we empoy the bootstrap to capture the parameter estimation error and an equation-based stochastic kriging metamode to propagate the input uncertainty to the output mean. Asymptotic anaysis provides theoretica support for our approach, whie an empirica study demonstrates that it has good finite-sampe performance. 1 INTRODUCTION Stochastic simuation is used to characterize the behavior of compex systems that are driven by random input modes. The choice of input modes directy impacts the system performance estimates. It is common practice to treat the input modes as a coection of independent univariate distributions. However, these simpe input modes do not aways faithfuy represent the physica input processes. For exampe, in a manufacturing system the processing times for a singe workpiece at a series of machining stations coud be dependent due to characteristics of that particuar workpiece; and in a suppy chain system the customer demands over different products from muti-product warehouses coud be reated. Ignoring such dependence can ead to poor estimates of system performance measures (Livny et a. 1993). Thus, it is desirabe to buid input modes that faithfuy capture the dependence. Considering the amount of information needed to construct fu joint distributions, in this paper we focus on input modes characterized by margina distributions and correation matrix. The margina distributions have known parametric famiies with unknown parameter vaues. The dependence between different components of the input modes can be measured by various criteria (Bier and Ghosh 2006); we focus on the Spearman rank correation. Even though product-moment correation is widey used in engineering appications, it cannot capture compex noninear dependence and its definition needs the variances of the components to be finite. The Spearman rank correation can avoid these probems. Since the input modes are estimated from finite sampes of rea-word data, a compete statistica characterization of stochastic system performance requires quantifying both simuation and input estimation error. We buid on Xie et a. (2014a), which used a metamode-assisted bootstrapping approach to form /14/$ IEEE 674
2 Xie, Neson, and Barton a confidence interva (CI) accounting for the impact of input and simuation uncertainty on the system s mean performance estimates. However, their study assumes that the input distributions are univariate and mutuay independent. To efficienty and correcty account for the dependence between various components of input modes, we introduce a more genera metamode-assisted bootstrapping framework that can quantify the impact of dependent input modes and simuation estimation error on system performance estimates, whie simutaneousy reducing the infuence of simuation estimation error compared with direct simuation methods. Specificay, we estimate the key properties of a fexibe input mode with rea-word data. The bootstrap is then used to quantify the estimation error of input margina distributions and dependence measures. An equation-based stochastic kriging metamode propagates the input uncertainty to the output mean. From this, we can derive a CI that accounts for both simuation and input uncertainty. Notice that we are interested in the dependence between different components of the input distributions instead the dependence among the estimated input-mode parameter vaues. There are two centra contributions of this paper: First, we generaize the metamode-assisted bootstrap framework of Xie et a. (2014a) to stochastic simuations with dependent input modes. Second, we propose a rigorous anaysis for cases where the dependence is measured by Spearman rank correation. The next section reviews some previous work for input uncertainty in simuation. This is foowed by a forma description of the probem of interest in Section 3. In Section 4, we propose a generaized metamodeassisted bootstrapping framework and provide an agorithm to buid a CI accounting for both input and simuation estimation error on the system mean performance estimates. Our approach is supported with asymptotic anaysis. We then report resuts of finite-sampe behavior from an empirica study in Section 5 and concude the paper in Section 6. 2 BACKGROUND Various approaches to account for input uncertainty in stochastic simuations have been proposed, see Barton (2012) for a review. The metamode-assisted bootstrapping approach was introduced in Barton et a. (2014). In that paper the input modes are a coection of independent univariate parametric distributions with known parametric famiies and unknown parameters vaues. The input parameters are estimated by rea-word data and the bootstrap is used to quantify the input-parameter uncertainty. Based on the simuation outputs at a few design points, a fexibe stochastic kriging (SK) metamode is buit to propagate the input uncertainty to the output mean. A CI is derived to quantify the impact of input uncertainty. Compared with the direct simuation approach, the metamode can reduce the impact of simuation estimation error. However, Barton et a. (2014) assumed that the simuation budget is not tight and the metamode uncertainty can be ignored. If the true mean response surface is compex, especiay for high-dimensiona probems with many input distributions, and the computationa budget is tight, then the impact of metamode uncertainty can no onger be ignored. The metamode-assisted bootstrapping approach was improved in Xie et a. (2014a) to buid a CI accounting for the impact from both input and metamode uncertainty on the system mean estimates. Further, a variance decomposition was proposed to estimate the reative contribution of input to overa uncertainty, which is very usefu for decision makers to determine where to put more effort to reduce the system uncertainty. The metamode-assisted bootstrapping approach demonstrated robust performance even when there is a tight computationa budget and simuation estimation error is arge. However, Xie et a. (2014a) aso assumed that the input modes are mutuay independent univariate distributions. To faithfuy capture the dependence between different components of input modes, Cario and Neson (1997) proposed a fexibe NORma To Anything (NORTA) distribution to represent and generate random vectors with amost arbitrary margina distributions and correation matrix. Since NORTA representations are estimated from finite sampes of rea-word data, Bier and Coru (2011) proposed a Bayesian approach to account for parameter uncertainty. Specificay, the uncertainty about the NORTA distribution parameter estimates is quantified by posterior distributions. For compex stochastic systems with a arge number of correated inputs, a fast agorithm was proposed to draw sampes from these posterior distributions 675
3 Xie, Neson, and Barton to quantify the input uncertainty. Then, the direct simuation method was used to propagate the input uncertainty to the output mean. However, when the simuated system is compex and the computationa budget is tight, the direct simuation method can not efficienty use the computationa budget and incurs substantia simuation estimation error. The good performance of metamode-assisted bootstrapping for stochastic simuations with independent univariate input distributions in Barton et a. (2014) and Xie et a. (2014a) motivates us to extend it to more compex cases with dependence in the input modes. Therefore, in this paper, we use fexibe NORTA representations for unknown dependent inputs. Since the NORTA method is based on estimating margina distribution parameters and a correation matrix from rea-word data, metamode-assisted bootstrapping is generaized to quantify the impact of both NORTA parameters and simuation estimation error whie simutaneousy reducing the infuence of simuation estimation error due to output variabiity. 3 PROBLEM STATEMENT The stochastic simuation output is a function of random numbers and the input modes denoted by F, where F is composed of a coection of input distributions used to drive the simuation. For notation simpification, we do not expicity incude the random numbers. The output from the jth repication of a simuation with input modes F can be written as Y j (F) = µ(f) + ε j (F) where µ(f) = E[Y j (F)] denotes the unknown output mean and ε j (F) represents the simuation error with mean zero. Notice that the simuation output depends on the choice of input modes. Let F c denote the true correct input modes, which are unknown and estimated from finite sampes of rea-word data. Our goa is to quantify the impact of the statistica error on system mean performance estimates by finding a (1 α)100% CI [Q L,Q U ] such that Pr{µ(F c ) [Q L,Q U ]} = 1 α. Let F = {F 1,F 2,...,F L } with F 1,F 2,...,F L mutuay independent; F coud be composed of univariate and mutivariate joint distributions. In this paper, we do not consider time-series input processes. Let F be a d dimensiona distribution having margina distributions denoted by {F,1,F,2,...,F,d } with d 1. For the th distribution F with d > 1, et X F be a d 1 random vector having a d d Spearman rank correation matrix denoted by R X with eements R X (i, j) = cor(f,i (X,i ),F, j (X, j )) = E[F,i(X,i )F, j (X, j )] E[F,i (X,i )]E[F, j (X, j )] Var(F,i (X,i ))Var(F, j (X, j )) with i, j = 1,2,...,d. Since the correation matrix is symmetric and diagona terms are 1, we can view a d d correation matrix as an eement of a d d (d 1)/2 dimensiona vector space. Then, the Spearman rank correation matrix can be uniquey specified by a d 1 vector denoted by V X. For F with d = 1, V X is empty and d = 0. For F with = 1,2,...,L, we assume that the famiies of margina distributions {F,1,F,2,...,F,d } are known, but not their parameter vaues. Let an h,i 1 vector θ,i denote the unknown parameters for the ith margina distribution F,i. By stacking θ,i with i = 1,2,...,d together, we have a d 1 dimensiona parameter vector θ (θ,1,θ,2,...,θ,d ) with d d i=1 h,i. Suppose the input mode is characterized by margina distributions and correation matrix, which are specified by ϑ {(θ ;V X ), = 1,2,...,L} that incudes d L =1 d + L =1 d eements. We ca ϑ the input mode parameters. By abusing notation, we can rewrite µ(f) as µ(ϑ). The true input parameters ϑ c are unknown and estimated from finite sampes of rea-word data. Thus, our goa can be restated as finding a (1 α)100% CI [Q L,Q U ] such that Pr{µ(ϑ c ) [Q L,Q U ]} = 1 α. (1) 676
4 Xie, Neson, and Barton { Let m denote the } number of i.i.d. rea-word observations avaiabe from the th input process X,m X (1),X (2),...,X (m ) with d 1 random vector X (i) i.i.d F c for i = 1,2,...,m. Let X m = {X,m, = 1,2,...,L} be the coection of sampes from a L input distributions in F c, where m = (m 1,m 2,...,m L ). The estimators for the input mode parameters are denoted by ϑ m {( θ,m ; V X,m ), = 1,2,...,L} and are a function of X m. The rea-word data are a particuar reaization of X m, say x (0) m. Input uncertainty is quantified by the samping distribution of µ( ϑ m ). Further, since the underying response surface µ( ) is unknown, for any ϑ et µ(ϑ) denote the corresponding mean response estimator. Thus, there are both input and simuation estimation errors in the system mean performance estimates. This paper focuses on estimating the performance of compex stochastic systems with input modes incuding mutivariate distributions. Our objective is to propose an approach to quantify the overa impact of both input and simuation estimation error on system mean performance estimates and then buid a CI satisfying Equation (1). Further, since each simuation run coud be computationay expensive and we have computationa budget denoted by N, we want to reduce the infuence of simuation estimation error due to output variabiity. 4 METAMODEL-ASSISTED BOOTSTRAPPING FRAMEWORK For simuations with parametric input distributions that are univariate and mutuay independent, a metamode-assisted bootstrapping framework was used to account for the impact of both input and simuation estimation errors in Xie et a. (2014a). In this section, we generaize the metamode-assisted bootstrapping approach to aow NORTA input modes. A procedure to buid a CI for the system mean performance that accounts for overa uncertainty is described and supported by asymptotic anaysis. 4.1 Bootstrap for Input Uncertainty The way we choose to represent input modes pays an important roe in the impementation of metamodeassisted bootstrapping. For the th input distribution F, instead of using the natura parameters θ to characterize the margina distributions, we can use moments; see Barton et a. (2014) for an expanation. Suppose that the parametric margina distribution F,i can be uniquey characterized by its first h,i finite moments denoted by the h,i 1 vector ψ,i for i = 1,2,...,d. By stacking ψ,i with i = 1,2,...,d together, we have a d 1 dimensiona vector of margina moments ψ (ψ,1,ψ,2,...,ψ,d ). Therefore, the input modes can be characterized by the moments M {(ψ ;V X ), = 1,2,...,L}. Abusing notation, we rewrite µ(ϑ) as µ(m ). The true moments of dependent input modes, denoted by M c, are unknown and estimated using X m. We use moment estimators for margina distributions, denoted by ψ,m. The estimation error of input modes can be quantified by the samping distribution of Mm {( ψ,m ; V X,m ), = 1,2,...,L}, denoted by FM c m. Therefore, the input uncertainty can be measured by the samping distribution of µ( M m ) with M m FM c m. Since it is hard to derive the samping distribution FM c m, we use bootstrap resamping to approximate it (Shao and Tu 1995). Since F 1,F 2,...,F L are mutuay independent, we can do bootstrapping for each distribution separatey. For distribution F, et A {1,2,...,m } be the index set of the rea-word observations. The procedure to quantify the input uncertainty by the bootstrap is as foows: 1. For the th distribution F with = 1,2,...,L, draw m sampes with repacement from set A and obtain indexes {i 1,i 2,...,i m }; choose corresponding sampes from rea-word data x (0) m and get x (1),m {x (i 1),x (i 2),...,x (i m ) }. Denote the coection of bootstrap sampes by x (1) m = 677
5 { x (1) Xie, Neson, and Barton,m, = 1,2,...,L} and use it to cacuate the bootstrapped moment estimates, denoted by { } M m (1) ( ψ (1),m ;Ṽ (1) X,m ), = 1,2,...,L. 2. Repeat the previous step B times to generate M m (b) for b = 1,2,...,B. The bootstrap resamped moments are drawn from the bootstrap distribution denoted by F Mm ( x (0) m ) with M m (b) F Mm ( x (0) m ). For estimation of a CI, B is recommended to be a few thousand. In this paper, a denotes a quantity estimated from rea-word data, whie a denotes a quantity estimated from bootstrapped data. 4.2 NORTA Representation In this paper, mutivariate input modes are characterized by their margina distributions and correation matrix. Since this partia characterization does not uniquey determine the joint distributions except in specia cases, we use the NORTA representation introduced by Cario and Neson (1997); NORTA vectors are competey specified by their margina distributions and correation matrix. In this section we describe how to appy NORTA in our metamode-assisted bootstrapping approach. Since F = {F 1,F 2,...,F L } with F 1,F 2,...,F L mutuay independent, we ony need to appy NORTA separatey for the distributions F with d > 1. Specificay, we represent X as a transformation of a d dimensiona standard mutivariate norma (MVN) vector Z = (Z,1,Z,2,...,Z,d ) with product-moment correation matrix ρ Z, ( ) X = F,1 1 [Φ(Z,1);θ,1 ],F,2 1 [Φ(Z,2);θ,2 ],...,F,d 1 [Φ(Z,d );θ,d ]. (2) If the margina distribution famiies are given, as we assume here, then the NORTA representation for F is (θ,ρ Z ). Let R Z denote a d d Spearman rank correation matrix for Z. Notice that there is a cosed form reationship between product-moment and Spearman rank correations for the standard mutivariate norma distribution (Cemen and Reiy 1999): For any i, j = 1,2,...,d, we have R Z (i, j) = 6 π sin 1 ( ρz (i, j) 2 ). (3) In this paper, we focus on continuous margina distributions with stricty increasing cdfs {F,1,F,2,...,F,d } for = 1,2,...,L. Since the Spearman rank correation is invariant under the monotone one-to-one transformation F,i 1 [Φ( )], we have R X = R Z. Given ϑ = {(θ ;V X ), = 1,2,...,L}, the procedure to generate the NORTA random variates is as foows: 1. From V X, get the Spearman rank correation matrix R Z = R X. 2. Cacuate ρ Z (i, j) = 2sin(πR Z (i, j)/6) for i, j = 1,2,...,d. 3. i.i.d. Generate Z MVN(0,ρ Z ) and obtain X by using Equation (2). 4. Repeat Steps 1 3 for a F with d > 1. Notice that when we use Spearman rank correation to measure the input mode dependence, the choice of margina distributions and correation is separabe. Thus, for any feasibe ϑ, we can find the corresponding NORTA representations. 4.3 Stochastic Kriging Metamode Instead of running simuations for systems with input modes represented by bootstrap moment sampes { M m (1), M (2) m,..., M (B) m }, we can choose a sma number of design points and run simuations there to buid a metamode; see Xie et a. (2014a). We then propagate the input uncertainty to the output mean 678
6 Xie, Neson, and Barton by using the metamode. In this paper, we use a fexibe stochastic kriging (SK) metamode introduced by Ankenman et a. (2010). We briefy review it in this section. For more detaied information, see Ankenman et a. (2010). Suppose that the underying true (but unknown) response surface is a reaization of a stationary Gaussian Process (GP). We mode the simuation output Y by Y j (x) = β 0 +W(x) + ε j (x). (4) This mode incudes two sources of uncertainty: the simuation output uncertainty ε j (x) and the mean response uncertainty W(x). Since stochastic systems with dependent input modes having simiar key properties tend to have mean responses cose to each other, SK uses a mean-zero, second-order stationary GP W( ) to account for this spatia dependence of the response surface. The uncertainty about the unknown true response surface µ(x) is represented by a GP M(x) β 0 +W(x) (note that β 0 can be repaced by a more genera trend term f(x) β). For many, but not a, simuation settings the output is an average of a arge number of more basic outputs, so a norma approximation can be appied: ε(x) N(0,σε 2 (x)). In SK, the covariance between W(x) and W(x ) quantifies how knowedge of the surface at one ocation affects the prediction at another ocation. In this paper, a parametric form of the covariance is used to capture this spatia dependence, Σ(x,x ) = Cov[W(x),W(x )] = τ 2 r(x x ), where τ 2 denotes the variance and r( ) is a correation function that depends ony on the distance x x. Using prior information about the smoothness of µ( ), we can choose the form of correation function. Based on previous study (Xie et a. 2010), we use the product-form Gaussian correation function ( r(x x d ) ) = exp φ j (x j x j) 2 (5) j=1 for the empirica evauation in Section 5. Let φ = (φ 1,φ 2,...,φ d ) represent the correation function parameters. So in summary, before having any simuation resuts the uncertainty about µ(x) is represented by a Gaussian process M(x) GP(β 0,τ 2 r(x x )). To reduce the uncertainty about µ(x) we choose an experiment design consisting of pairs D {(x i,n i ),i = 1,2,...,k} at which to run simuations, where (x i,n i ) denotes the ocation and the number of repications, respectivey, at the ith design point. The simuation outputs at D are Y D {(Y 1 (x i ),Y 2 (x i ),...,Y ni (x i ));i = 1,2,...,k} and the sampe mean at design point x i is Ȳ (x i ) = n i j=1 Y j(x i )/n i. Let the sampe means at a k design points be Ȳ D = (Ȳ (x 1 ),Ȳ (x 2 ),...,Ȳ (x k )) T. Since the use of common random numbers is detrimenta to prediction (Chen et a. 2012), the simuations at different design points are independent and the variance of Ȳ D is represented by a k k diagona matrix C = diag { σε 2 (x 1 )/n 1,σε 2 (x 2 )/n 2,...,σε 2 } (x k )/n k. Let Σ be the k k spatia covariance matrix of the design points and et Σ(x, ) be the k 1 spatia covariance vector between each design point and a fixed prediction point x. If the parameters (τ 2,φ,C) are known, then the metamode uncertainty can be characterized by a refined GP M p (x) that denotes the conditiona distribution of M(x) given a simuation outputs, M p (x) GP(m p (x),σ 2 p(x)) (6) where m p ( ) is the minimum mean squared error (MSE) inear unbiased predictor and the corresponding variance is m p (x) = β 0 + Σ(x, ) (Σ +C) 1 (Ȳ D β 0 1 k 1 ), (7) σ 2 p(x) = τ 2 Σ(x, ) (Σ +C) 1 Σ(x, ) + η [1 k 1(Σ +C) 1 1 k 1 ] 1 η (8) 679
7 Xie, Neson, and Barton where β 0 = [1 k 1 (Σ +C) 1 1 k 1 ] 1 1 k 1 (Σ +C) 1 Ȳ D and η = 1 1 k 1 (Σ +C) 1 Σ(x, ) (Ankenman et a. 2010). Since in reaity the spatia correation parameters τ 2 and φ are unknown, maximum ikeihood estimates are typicay used for prediction, and the sampe variance is used as an estimate for the simuation variance at design points C (Ankenman et a. 2010). By pugging these into Equations (7) and (8) we can obtain the estimated mean m p (x) and variance σ 2 p(x). Thus, the metamode we use is µ(x) = m p (x) with variance estimated by σ 2 p(x). Notice that since accounting for the parameter estimation error is intractabe, this pug-in estimator is commony used in the kriging iterature. In the metamode-assisted bootstrapping approach, the dependent input mode characterized by moments M can be interpreted as a ocation x in a d dimensiona space. 4.4 Procedure to Buid CI Since there are both input and simuation estimation errors in the system mean performance estimates, in this section, we propose a procedure to buid a CI quantifying the overa uncertainty by using our generaized metamode-assisted bootstrapping approach. Asymptotic anaysis shows that as m, B, the CI is consistent. Based on a hierarchica approach, we propose the foowing procedure to buid a (1 α)100% bootstrap percentie CI: 1. Identify a design space E for input mode moments M over which to fit the metamode. Since the metamode is used to propagate the input uncertainty measured by the bootstrapped moments Mm to the output mean, the design space is chosen to be the smaest eipsoid covering most ikey bootstrapped moments. See Barton et a. (2014) for more detaied information. 2. To obtain an experiment design D = {(M i,n i ),i = 1,2,...,k}, use a Latin hypercube sampe to embed k design points into the design space E and assign equa repications to these points to exhaust N. 3. For i = 1 to k (oop through k design points) (a) Use moment matching to cacuate the margina parameters θ (i) for = 1,2,...,L. (b) By foowing the description in Section 4.2, find the NORTA representation with parameters ( θ (i) ),ρ(i) Z for F with = 1,2,...,L and d > 1. Next i 4. At a design points, generate sampes of X by using NORTA representations for F with d > 1 and using standard approaches (Neson 2013) for F with d = 1, = 1,2,...,L. Use these sampes to drive simuations (see Section 4.2) and obtain outputs y D. Compute the sampe average ȳ(m i ) and sampe variance s 2 (M i ) of the simuation outputs, i = 1,2,...,k. Fit a SK metamode to obtain m p ( ) and σ p( ) 2 using ( ȳ(m i ),s 2 ) (M i ),M i, i = 1,2,...,k. 5. For b = 1 to B (a) Generate bootstrap moments M (b) ( m by foowing the procedure in Section 4.1. (b) Draw M b N m p ( M m (b) ), σ p( 2 M ) m (b) ). Next b [ ] 6. Report CI: M ( B α 2 ), M ( B(1 α 2 ) ) where, M (1) M (2) M (B) are the sorted vaues. If SK parameters (τ 2,φ,C) are known and we repace M b in Step 5(b) of the CI procedure with M b N ( m p ( M m (b) ),σp( 2 M m (b) ) ), then we can show that the CI obtained [M ( B α 2 ),M ( B(1 α 2 ) )] is asymptoticay consistent by Theorem 1. Further, this CI characterizes the impact from both input and metamode uncertainty on the system performance estimate. Theorem 1 Suppose that the foowing assumptions hod. 680
8 Xie, Neson, and Barton Figure 1: A stochastic activity network. 1. We have i.i.d. observations X ( j) i.i.d F c for j = 1,2,...,m and = 1,2,...,L. 2. The ε j (x) i.i.d. N(0,σε 2 (x)) for any x, and M(x) is a stationary, separabe GP with a continuous correation function satisfying 1 r(x x c ) og( x x 2 ) 1+γ for a x x 2 δ for some c > 0, γ > 0 and δ < 1, where x x 2 = d j=1 (x j x j )2. 3. The input processes, simuation noise ε j (x) and GP M(x) are mutuay independent and the bootstrap process is independent of a of them. As m, we have m /m c with = 1,2,...,L for a constant c > 0. Then the interva [M ( B α 2 ),M ( B(1 α 2 ) )] is asymptoticay consistent, meaning the iterated imit im m im B Pr{M ( Bα/2 ) M p (M c ) M ( B(1 α/2) ) } = 1 α. (9) Proof: By Theorem 1 in Xie et a. (2014b), we have M a.s. m M c as m. Then, foowing simiar steps as in proving Theorem 1 in Xie et a. (2014a), we can show Equation (9). 5 EMPIRICAL STUDY Since F 1,F 2,...,F L are mutuay independent, without oss of generaity we consider an exampe with L = 1 and suppress the subscript for the th input distribution. We use a stochastic activity network as shown in Figure 1 to examine the finite-sampe performance of our metamode-assisted bootstrapping approach. Suppose that the time required to compete task (arc) i is denoted by X i for i = 1,2,...,5 and X = (X 1,X 2,...,X 5 ). We wish to compute the time to compete the project, which is the ongest path through the network, Y = max{x 1 + X 2 + X 5,X 1 + X 4,X 3 + X 5 }. We are interested in the mean response E[Y ]. We assume that F c is NORTA. The margina distributions are X i exp(θi c ) for i = 1,2,...,5 with means θ c = (10,5,12,11,5). The Spearman rank correation matrix is R c X =
9 Xie, Neson, and Barton Therefore, the number of parameters characterizing the input mode F is d = d +d = 5+(5 4)/2 = 15. Since the true mean response µ(ϑ c ) is unknown, we run 10 7 repications and obtain the estimated true mean response with standard error To evauate our metamode-assisted bootstrapping approach, we pretend that the input-mode parameters (θ c,r c X ) are unknown and they are estimated by m i.i.d. observations from Fc ; this represents obtaining rea-word data. The goa is to buid a CI quantifying the impact of both input and simuation estimation error on the system mean response estimate. We compare metamode-assisted bootstrapping to the conditiona CI and direct bootstrapping. For the conditiona CI, we fit the input distribution to the rea-word data by moment matching and aocate a computationa budget N repications to simuating the resuting system. In direct bootstrapping, we run N/B repications of the simuation at each bootstrap moment Ȳ b = Ȳ ( M (b) ), and report the percentie CI m [Ȳ( B α 2 ),Ȳ ( B(1 α 2 ) ) ] m. M (b), record the average simuation output Tabe 1 shows the statistica performance of conditiona and direct bootstrapping CIs and metamodeassisted bootstrapping with k = 80 design points, m = 100, 500, 1000 rea-word observations, and computationa budget of N = 10 3, 10 4, and 10 5 repications. We ran 1000 macro-repications of the entire experiment. In each macro-repication, we first generate m mutivariate observations by using NORTA with parameters (θ c,r c X ). Then, for the conditiona CI, we run N repications at the estimated parameters ( θ m, R X,m ) and buid CIs with nomina 95% coverage of the response mean. For direct bootstrapping and metamode-assisted bootstrapping, we use bootstrapping to generate B = 1000 sampes moments to quantify the input uncertainty. Since µ( ) is unknown, we use the fixed computationa budget N to propagate the input uncertainty either via direct simuation or the SK metamode to buid percentie CIs with nomina 95% coverage. From Tabe 1 we observe that under the same computationa budget N, the conditiona CIs that ony account for the simuation uncertainty tend to have undercoverage. The CIs obtained by direct simuation are much wider and they typicay have obvious over-coverage. The CIs obtained by metamode-assisted bootstrapping have coverage much coser to the nomina eve of 95%. As N increases and simuation estimation error decreases, the undercoverage probem for the conditiona CI becomes worse. Since direct bootstrapping and metamode-assisted bootstrapping use the same set of bootstrapped sampes to quantify input uncertainty, the overcoverage for the direct bootstrap represents the additiona uncertainty introduced whie propagating the input uncertainty to the output mean. Tabe 1 shows that the metamode can effectivey use the computationa budget and reduce the impact from simuation estimation error. Further, as the computationa budget increases, the difference between the CIs obtained by the two methods diminishes. Figure 2 shows a scatter pot of conditiona CIs and CIs obtained by direct simuation and metamodeassisted bootstrapping when m = 500, k = 80 and N = It incudes resuts from 100 macro-repications. The horizonta axis represents (Q L +Q U )/2 that gives a point estimate of system mean performance, where Q L and Q U are the ower and upper bounds of the CIs. The vertica axis is (Q U Q L )/2 that gives haf width of the CIs. Region 1 contains points that correspond to CIs having underestimation and Region 3 contains points corresponding to overestimation, whie Region 2 contains CIs that cover µ(f c ) (Kang and Schmeiser 1990). Conditiona CIs have short width and their centers have arge variance. Therefore, they have serious undercoverage. For resuts from metamode-assisted bootstrapping, the proportion of CIs in Region 2 is cose to 95% and CIs outside tend to have underestimation. 6 CONCLUSIONS In this paper, we extended the metamode-assisted bootstrapping approach to cases with dependence in the input modes. The input modes are characterized by their margina distributions and Spearman rank correation, which are estimated from rea-word data. Metamode-assisted bootstrapping uses the bootstrap to quantify the input uncertainty and propagates it to the output mean by using a SK metamode. This 682
10 Xie, Neson, and Barton Tabe 1: Resuts for CIs when m = 100,500,1000. m = 100 N = 10 3 N = 10 4 N = 10 5 conditiona CI coverage 43% 14.1% 4.4% CI width (mean) CI width (SD) direct simuation coverage 100% 100% 99.4% CI width (mean) CI width (SD) metamode-assisted coverage 97.6% 96.1% 94.8% CI width (mean) CI width (SD) m = 500 N = 10 3 N = 10 4 N = 10 5 conditiona CI coverage 76.9% 32.1% 12% CI width (mean) CI width (SD) direct simuation coverage 100% 100% 100% CI width (mean) CI width (SD) metamode-assisted coverage 96.3% 98.7% 97.2% CI width (mean) CI width (SD) m = 1000 N = 10 3 N = 10 4 N = 10 5 conditiona CI coverage 83.4% 42.2% 16.7% CI width (mean) CI width (SD) direct simuation coverage 100% 100% 100% CI width (mean) CI width (SD) metamode-assisted coverage 97.1% 94% 96.6% CI width (mean) CI width (SD)
11 Xie, Neson, and Barton Figure 2: Scatter pot of conditiona CIs and CIs obtained by direct simuation and metamode-assisted bootstrapping when m = 500, k = 80 and N = approach deivers a CI quantifying the overa uncertainty of the system performance estimate. Compared with the direct bootstrap, the metamode can make effective use of the simuation budget. An empirica study demonstrated that our metamode-assisted bootstrap has good finite-sampe performance under various quantities of rea-word data and simuation budget. ACKNOWLEDGMENTS This paper is based upon work supported by the Nationa Science Foundation under Grant No. CMMI REFERENCES Ankenman, B. E., B. L. Neson, and J. Staum Stochastic Kriging for Simuation Metamodeing. Operations Research 58: Barton, R. R Tutoria: Input Uncertainty in Output Anaysis. In Proceedings of the 2012 Winter Simuation Conference, edited by C. Laroque, J. Himmespach, R. Pasupathy, O. Rose, and A.M. Uhrmacher, 67 78: IEEE. Barton, R. R., B. L. Neson, and W. Xie Quantifying Input Uncertainty via Simuation Confidence Intervas. Informs Journa on Computing 26: Bier, B., and C. G. Coru Accounting for Parameter Uncertainty in Large-Scae Stochastic Simuations with Correated Inputs. Operations Research 59: Bier, B., and S. Ghosh Mutivariate Input Processes. In Handbooks in Operations Research and Management Science: Simuation, edited by S. Henderson and B. L. Neson, Chapter 5. Esevier. 684
12 Xie, Neson, and Barton Cario, M. C., and B. L. Neson Modeing and Generating Random Vectors with Arbitrary Margina Distributions and Correation Matrix. Technica report, Department of Industria Engineering and Management Sciences, Northwestern University. Chen, X., B. E. Ankenman, and B. L. Neson The Effect of Common Random Numbers on Stochastic Kriging Metamodes. ACM Transactions on Modeing and Computer Simuation 22:7. Cemen, R. T., and T. Reiy Correations and Copuas for Decision and Risk Anaysis. Management Science 45: Kang, K., and B. Schmeiser Graphica Methods for Evauation and Comparing Confidence-Interva Procedures. Operations Research 38 (3): Livny, M., B. Meamed, and A. K. Tsiois The Impact of Autocorreation on Queueing Systems. Management Science 39: Neson, B. L Foundations and Methods of Stochastic Simuation: A First Course. Springer-Verag. Shao, J., and D. Tu The Jackknife and Bootstrap. Springer-Verag. Xie, W., B. L. Neson, and R. R. Barton. 2014a. Statistica Uncertainty Anaysis for Stochastic Simuation. Technica report, Department of Industria Engineering and Management Sciences, Northwestern University. Xie, W., B. L. Neson, and R. R. Barton. 2014b. Statistica Uncertainty Anaysis for Stochastic Simuation with Dependent Input Modes. Technica report, Department of Industria Engineering and Management Sciences, Northwestern University. Xie, W., B. L. Neson, and J. Staum The Infuence of Correation Functions on Stochastic Kriging Metamodes. In Proceedings of the 2010 Winter Simuation Conference, edited by B. Johansson, S. Jain, J. Montoya-Torres, J. Hugan, and E. Yucesan, : IEEE. AUTHOR BIOGRAPHIES WEI XIE is an assistant professor in the Department of Industria and Systems Engineering at Rensseaer Poytechnic Institute. Her research interests are in computer simuation, appied statistics and data anaytics. Her emai address is weixie2013@u.northwestern.edu. BARRY L. NELSON is the Water P. Murphy Professor and Chair of the Department of Industria Engineering and Management Sciences at Northwestern University. He is a Feow of INFORMS and IIE. His research centers on the design and anaysis of computer simuation experiments on modes of stochastic systems, and he is the author of Foundations and Methods of Stochastic Simuation: A First Course, from Springer. His e-mai address is nesonb@northwestern.edu. RUSSELL R. BARTON is Professor of Suppy Chain and Information Systems and Professor of Industria Engineering at the Pennsyvania State University. He currenty serves as senior associate dean for research and facuty in the Smea Coege of Business. He received a B.S. degree in Eectrica Engineering from Princeton University and M.S. and Ph.D. degrees in Operations Research from Corne University. Before entering academia, he spent tweve years in industry. He is a past president of the INFORMS Simuation Society, and currenty serves as the I-Sim representative on the INFORMS Subdivisions Counci and chairs the QSR Advisory Board. He is a senior member of IIE and IEEE and a Certified Anaytics Professiona. His research interests incude appications of statistica and simuation methods to system design and to product design, manufacturing and deivery. His emai address is rbarton@psu.edu. 685
A Multivariate Input Uncertainty in Output Analysis for Stochastic Simulation
A Mutivariate Input Uncertainty in Output Anaysis for Stochastic Simuation WEI XIE, Rensseaer Poytechnic Institute BARRY L. NELSON, Northwestern University RUSSELL R. BARTON, The Pennsyvania State University
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 Simuation Conference A. Tok, S. Y. Diao, I. O. Ryzhov, L. Yimaz, S. Buckey, and J. A. Mier, eds. ADVANCED TUTORIAL: INPUT UNCERTAINTY QUANTIFICATION Eunhye Song Barry L.
More informationA. Distribution of the test statistic
A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationSydU STAT3014 (2015) Second semester Dr. J. Chan 18
STAT3014/3914 Appied Stat.-Samping C-Stratified rand. sampe Stratified Random Samping.1 Introduction Description The popuation of size N is divided into mutuay excusive and exhaustive subpopuations caed
More informationStochastic Variational Inference with Gradient Linearization
Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,
More informationSTA 216 Project: Spline Approach to Discrete Survival Analysis
: Spine Approach to Discrete Surviva Anaysis November 4, 005 1 Introduction Athough continuous surviva anaysis differs much from the discrete surviva anaysis, there is certain ink between the two modeing
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 informationFitting Algorithms for MMPP ATM Traffic Models
Fitting Agorithms for PP AT Traffic odes A. Nogueira, P. Savador, R. Vaadas University of Aveiro / Institute of Teecommunications, 38-93 Aveiro, Portuga; e-mai: (nogueira, savador, rv)@av.it.pt ABSTRACT
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna
More informationDo Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix
VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix
More informationEfficiently Generating Random Bits from Finite State Markov Chains
1 Efficienty Generating Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown
More informationNEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION
NEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION Hsiao-Chang Chen Dept. of Systems Engineering University of Pennsyvania Phiadephia, PA 904-635, U.S.A. Chun-Hung Chen
More informationAppendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model
Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this
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 informationTwo-Stage Least Squares as Minimum Distance
Two-Stage Least Squares as Minimum Distance Frank Windmeijer Discussion Paper 17 / 683 7 June 2017 Department of Economics University of Bristo Priory Road Compex Bristo BS8 1TU United Kingdom Two-Stage
More informationIIE Transactions Publication details, including instructions for authors and subscription information:
This artice was downoaded by: [Professor Barry Neson] On: 24 March 205, At: 09:50 Pubisher: Tayor & Francis Informa Ltd Registered in Engand and Waes Registered Number: 072954 Registered office: Mortimer
More informationASummaryofGaussianProcesses Coryn A.L. Bailer-Jones
ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe
More informationBayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with?
Bayesian Learning A powerfu and growing approach in machine earning We use it in our own decision making a the time You hear a which which coud equay be Thanks or Tanks, which woud you go with? Combine
More informationMustafa H. Tongarlak Bruce E. Ankenman Barry L. Nelson
Proceedings of the 0 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelspach, K. P. White, and M. Fu, eds. RELATIVE ERROR STOCHASTIC KRIGING Mustafa H. Tongarlak Bruce E. Ankenman Barry L.
More informationA Comparison Study of the Test for Right Censored and Grouped Data
Communications for Statistica Appications and Methods 2015, Vo. 22, No. 4, 313 320 DOI: http://dx.doi.org/10.5351/csam.2015.22.4.313 Print ISSN 2287-7843 / Onine ISSN 2383-4757 A Comparison Study of the
More information6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7
6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17 Soution 7 Probem 1: Generating Random Variabes Each part of this probem requires impementation in MATLAB. For the
More informationStatistical Learning Theory: A Primer
Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More informationAlberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain
CORRECTIONS TO CLASSICAL PROCEDURES FOR ESTIMATING THURSTONE S CASE V MODEL FOR RANKING DATA Aberto Maydeu Oivares Instituto de Empresa Marketing Dept. C/Maria de Moina -5 28006 Madrid Spain Aberto.Maydeu@ie.edu
More informationAST 418/518 Instrumentation and Statistics
AST 418/518 Instrumentation and Statistics Cass Website: http://ircamera.as.arizona.edu/astr_518 Cass Texts: Practica Statistics for Astronomers, J.V. Wa, and C.R. Jenkins, Second Edition. Measuring the
More informationProcess Capability Proposal. with Polynomial Profile
Contemporary Engineering Sciences, Vo. 11, 2018, no. 85, 4227-4236 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ces.2018.88467 Process Capabiity Proposa with Poynomia Profie Roberto José Herrera
More informationA proposed nonparametric mixture density estimation using B-spline functions
A proposed nonparametric mixture density estimation using B-spine functions Atizez Hadrich a,b, Mourad Zribi a, Afif Masmoudi b a Laboratoire d Informatique Signa et Image de a Côte d Opae (LISIC-EA 4491),
More informationA Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)
A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,
More informationSome Measures for Asymmetry of Distributions
Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester
More informationExpectation-Maximization for Estimating Parameters for a Mixture of Poissons
Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Brandon Maone Department of Computer Science University of Hesini February 18, 2014 Abstract This document derives, in excrutiating
More informationFORECASTING TELECOMMUNICATIONS DATA WITH AUTOREGRESSIVE INTEGRATED MOVING AVERAGE MODELS
FORECASTING TEECOMMUNICATIONS DATA WITH AUTOREGRESSIVE INTEGRATED MOVING AVERAGE MODES Niesh Subhash naawade a, Mrs. Meenakshi Pawar b a SVERI's Coege of Engineering, Pandharpur. nieshsubhash15@gmai.com
More informationStatistics for Applications. Chapter 7: Regression 1/43
Statistics for Appications Chapter 7: Regression 1/43 Heuristics of the inear regression (1) Consider a coud of i.i.d. random points (X i,y i ),i =1,...,n : 2/43 Heuristics of the inear regression (2)
More informationOptimality of Inference in Hierarchical Coding for Distributed Object-Based Representations
Optimaity of Inference in Hierarchica Coding for Distributed Object-Based Representations Simon Brodeur, Jean Rouat NECOTIS, Département génie éectrique et génie informatique, Université de Sherbrooke,
More informationExplicit overall risk minimization transductive bound
1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,
More informationAutomobile Prices in Market Equilibrium. Berry, Pakes and Levinsohn
Automobie Prices in Market Equiibrium Berry, Pakes and Levinsohn Empirica Anaysis of demand and suppy in a differentiated products market: equiibrium in the U.S. automobie market. Oigopoistic Differentiated
More informationSUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS
ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia
More informationTwo-sample inference for normal mean vectors based on monotone missing data
Journa of Mutivariate Anaysis 97 (006 6 76 wwweseviercom/ocate/jmva Two-sampe inference for norma mean vectors based on monotone missing data Jianqi Yu a, K Krishnamoorthy a,, Maruthy K Pannaa b a Department
More informationMARKOV CHAINS AND MARKOV DECISION THEORY. Contents
MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After
More informationFRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)
1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using
More informationAvailable online at ScienceDirect. Procedia Computer Science 96 (2016 )
Avaiabe onine at www.sciencedirect.com ScienceDirect Procedia Computer Science 96 (206 92 99 20th Internationa Conference on Knowedge Based and Inteigent Information and Engineering Systems Connected categorica
More informationEfficient Generation of Random Bits from Finite State Markov Chains
Efficient Generation of Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown
More informationAn Algorithm for Pruning Redundant Modules in Min-Max Modular Network
An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai
More informationA Solution to the 4-bit Parity Problem with a Single Quaternary Neuron
Neura Information Processing - Letters and Reviews Vo. 5, No. 2, November 2004 LETTER A Soution to the 4-bit Parity Probem with a Singe Quaternary Neuron Tohru Nitta Nationa Institute of Advanced Industria
More informationAsymptotic Properties of a Generalized Cross Entropy Optimization Algorithm
1 Asymptotic Properties of a Generaized Cross Entropy Optimization Agorithm Zijun Wu, Michae Koonko, Institute for Appied Stochastics and Operations Research, Caustha Technica University Abstract The discrete
More informationAALBORG UNIVERSITY. The distribution of communication cost for a mobile service scenario. Jesper Møller and Man Lung Yiu. R June 2009
AALBORG UNIVERSITY The distribution of communication cost for a mobie service scenario by Jesper Møer and Man Lung Yiu R-29-11 June 29 Department of Mathematica Sciences Aaborg University Fredrik Bajers
More informationTechnical Appendix for Voting, Speechmaking, and the Dimensions of Conflict in the US Senate
Technica Appendix for Voting, Speechmaking, and the Dimensions of Confict in the US Senate In Song Kim John Londregan Marc Ratkovic Juy 6, 205 Abstract We incude here severa technica appendices. First,
More informationImproving the Accuracy of Boolean Tomography by Exploiting Path Congestion Degrees
Improving the Accuracy of Booean Tomography by Expoiting Path Congestion Degrees Zhiyong Zhang, Gaoei Fei, Fucai Yu, Guangmin Hu Schoo of Communication and Information Engineering, University of Eectronic
More informationEfficient Bayesian Multifidelity Approach in Metamodeling of an Electromagnetic Simulator Response
Progress In Eectromagnetics Research M, Vo. 54, 47 54, 2017 Efficient Bayesian Mutifideity Approach in Metamodeing of an Eectromagnetic Simuator Response Tarek Bdour 1, 2, *, Christophe Guiffaut 1, and
More informationConsistent linguistic fuzzy preference relation with multi-granular uncertain linguistic information for solving decision making problems
Consistent inguistic fuzzy preference reation with muti-granuar uncertain inguistic information for soving decision making probems Siti mnah Binti Mohd Ridzuan, and Daud Mohamad Citation: IP Conference
More informationarxiv: v2 [stat.ap] 28 Mar 2016
Bayesian modeing of networks in compex business inteigence probems Daniee Durante, Say Paganin, Bruno Scarpa and David B. Dunson arxiv:0.00v [stat.ap] Mar 0 University of Padua, Department of Statistica
More information8 Digifl'.11 Cth:uits and devices
8 Digif'. Cth:uits and devices 8. Introduction In anaog eectronics, votage is a continuous variabe. This is usefu because most physica quantities we encounter are continuous: sound eves, ight intensity,
More informationAdjustment of automatic control systems of production facilities at coal processing plants using multivariant physico- mathematical models
IO Conference Series: Earth and Environmenta Science AER OEN ACCESS Adjustment of automatic contro systems of production faciities at coa processing pants using mutivariant physico- mathematica modes To
More informationStatistical Learning Theory: a Primer
??,??, 1 6 (??) c?? Kuwer Academic Pubishers, Boston. Manufactured in The Netherands. Statistica Learning Theory: a Primer THEODOROS EVGENIOU AND MASSIMILIANO PONTIL Center for Bioogica and Computationa
More informationAsynchronous Control for Coupled Markov Decision Systems
INFORMATION THEORY WORKSHOP (ITW) 22 Asynchronous Contro for Couped Marov Decision Systems Michae J. Neey University of Southern Caifornia Abstract This paper considers optima contro for a coection of
More informationIn-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017
In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative
More informationII. PROBLEM. A. Description. For the space of audio signals
CS229 - Fina Report Speech Recording based Language Recognition (Natura Language) Leopod Cambier - cambier; Matan Leibovich - matane; Cindy Orozco Bohorquez - orozcocc ABSTRACT We construct a rea time
More informationBayesian Nonparametric Covariance Regression
Journa of Machine Learning Research 16 (215) 251-2542 Submitted 2/15; Revised 7/15; Pubished 12/15 Bayesian Nonparametric Covariance Regression Emiy B. Fox Department of Statistics University of Washington
More informationRadar/ESM Tracking of Constant Velocity Target : Comparison of Batch (MLE) and EKF Performance
adar/ racing of Constant Veocity arget : Comparison of Batch (LE) and EKF Performance I. Leibowicz homson-csf Deteis/IISA La cef de Saint-Pierre 1 Bd Jean ouin 7885 Eancourt Cede France Isabee.Leibowicz
More informationDISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE
DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,
More informationSymbolic models for nonlinear control systems using approximate bisimulation
Symboic modes for noninear contro systems using approximate bisimuation Giordano Poa, Antoine Girard and Pauo Tabuada Abstract Contro systems are usuay modeed by differentia equations describing how physica
More informationMoreau-Yosida Regularization for Grouped Tree Structure Learning
Moreau-Yosida Reguarization for Grouped Tree Structure Learning Jun Liu Computer Science and Engineering Arizona State University J.Liu@asu.edu Jieping Ye Computer Science and Engineering Arizona State
More informationUniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete
Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity
More informationAkaike Information Criterion for ANOVA Model with a Simple Order Restriction
Akaike Information Criterion for ANOVA Mode with a Simpe Order Restriction Yu Inatsu * Department of Mathematics, Graduate Schoo of Science, Hiroshima University ABSTRACT In this paper, we consider Akaike
More informationProceedings of the 2018 Winter Simulation Conference M. Rabe, A. A. Juan, N. Mustafee, A. Skoogh, S. Jain, and B. Johansson, eds.
Proceedings of the 2018 Winter Simulation Conference M. Rabe, A. A. Juan, N. Mustafee, A. Skoogh, S. Jain, and B. Johansson, eds. METAMODEL-ASSISTED RISK ANALYSIS FOR STOCHASTIC SIMULATION WITH INPUT UNCERTAINTY
More informationA Sparse Covariance Function for Exact Gaussian Process Inference in Large Datasets
A Covariance Function for Exact Gaussian Process Inference in Large Datasets Arman ekumyan Austraian Centre for Fied Robotics The University of Sydney NSW 26, Austraia a.mekumyan@acfr.usyd.edu.au Fabio
More informationCollective organization in an adaptative mixture of experts
Coective organization in an adaptative mixture of experts Vincent Vigneron, Christine Fuchen, Jean-Marc Martinez To cite this version: Vincent Vigneron, Christine Fuchen, Jean-Marc Martinez. Coective organization
More informationAnalysis of Emerson s Multiple Model Interpolation Estimation Algorithms: The MIMO Case
Technica Report PC-04-00 Anaysis of Emerson s Mutipe Mode Interpoation Estimation Agorithms: The MIMO Case João P. Hespanha Dae E. Seborg University of Caifornia, Santa Barbara February 0, 004 Anaysis
More informationVALIDATED CONTINUATION FOR EQUILIBRIA OF PDES
VALIDATED CONTINUATION FOR EQUILIBRIA OF PDES SARAH DAY, JEAN-PHILIPPE LESSARD, AND KONSTANTIN MISCHAIKOW Abstract. One of the most efficient methods for determining the equiibria of a continuous parameterized
More informationHigh-order approximations to the Mie series for electromagnetic scattering in three dimensions
Proceedings of the 9th WSEAS Internationa Conference on Appied Mathematics Istanbu Turkey May 27-29 2006 (pp199-204) High-order approximations to the Mie series for eectromagnetic scattering in three dimensions
More informationarxiv: v1 [math.ca] 6 Mar 2017
Indefinite Integras of Spherica Besse Functions MIT-CTP/487 arxiv:703.0648v [math.ca] 6 Mar 07 Joyon K. Boomfied,, Stephen H. P. Face,, and Zander Moss, Center for Theoretica Physics, Laboratory for Nucear
More informationApproximate Bandwidth Allocation for Fixed-Priority-Scheduled Periodic Resources (WSU-CS Technical Report Version)
Approximate Bandwidth Aocation for Fixed-Priority-Schedued Periodic Resources WSU-CS Technica Report Version) Farhana Dewan Nathan Fisher Abstract Recent research in compositiona rea-time systems has focused
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationMinimizing Total Weighted Completion Time on Uniform Machines with Unbounded Batch
The Eighth Internationa Symposium on Operations Research and Its Appications (ISORA 09) Zhangiaie, China, September 20 22, 2009 Copyright 2009 ORSC & APORC, pp. 402 408 Minimizing Tota Weighted Competion
More informationGeneral Certificate of Education Advanced Level Examination June 2010
Genera Certificate of Education Advanced Leve Examination June 2010 Human Bioogy HBI6T/Q10/task Unit 6T A2 Investigative Skis Assignment Task Sheet The effect of using one or two eyes on the perception
More informationhttps://doi.org/ /epjconf/
HOW TO APPLY THE OPTIMAL ESTIMATION METHOD TO YOUR LIDAR MEASUREMENTS FOR IMPROVED RETRIEVALS OF TEMPERATURE AND COMPOSITION R. J. Sica 1,2,*, A. Haefee 2,1, A. Jaai 1, S. Gamage 1 and G. Farhani 1 1 Department
More informationCryptanalysis of PKP: A New Approach
Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in
More informationTRAVEL TIME ESTIMATION FOR URBAN ROAD NETWORKS USING LOW FREQUENCY PROBE VEHICLE DATA
TRAVEL TIME ESTIMATIO FOR URBA ROAD ETWORKS USIG LOW FREQUECY PROBE VEHICLE DATA Erik Jeneius Corresponding author KTH Roya Institute of Technoogy Department of Transport Science Emai: erik.jeneius@abe.kth.se
More informationGlobal sensitivity analysis using low-rank tensor approximations
Goba sensitivity anaysis using ow-rank tensor approximations K. Konaki 1 and B. Sudret 1 1 Chair of Risk, Safety and Uncertainty Quantification, arxiv:1605.09009v1 [stat.co] 29 May 2016 ETH Zurich, Stefano-Franscini-Patz
More informationBayesian Unscented Kalman Filter for State Estimation of Nonlinear and Non-Gaussian Systems
Bayesian Unscented Kaman Fiter for State Estimation of Noninear and Non-aussian Systems Zhong Liu, Shing-Chow Chan, Ho-Chun Wu and iafei Wu Department of Eectrica and Eectronic Engineering, he University
More informationSupport Vector Machine and Its Application to Regression and Classification
BearWorks Institutiona Repository MSU Graduate Theses Spring 2017 Support Vector Machine and Its Appication to Regression and Cassification Xiaotong Hu As with any inteectua project, the content and views
More informationSubstructuring Preconditioners for the Bidomain Extracellular Potential Problem
Substructuring Preconditioners for the Bidomain Extraceuar Potentia Probem Mico Pennacchio 1 and Vaeria Simoncini 2,1 1 IMATI - CNR, via Ferrata, 1, 27100 Pavia, Itay mico@imaticnrit 2 Dipartimento di
More informationTELECOMMUNICATION DATA FORECASTING BASED ON ARIMA MODEL
TEECOMMUNICATION DATA FORECASTING BASED ON ARIMA MODE Anjuman Akbar Muani 1, Prof. Sachin Muraraka 2, Prof. K. Sujatha 3 1Student ME E&TC,Shree Ramchandra Coege of Engineering, onikand,pune,maharashtra
More informationTesting for the Existence of Clusters
Testing for the Existence of Custers Caudio Fuentes and George Casea University of Forida November 13, 2008 Abstract The detection and determination of custers has been of specia interest, among researchers
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationNonlinear Gaussian Filtering via Radial Basis Function Approximation
51st IEEE Conference on Decision and Contro December 10-13 01 Maui Hawaii USA Noninear Gaussian Fitering via Radia Basis Function Approximation Huazhen Fang Jia Wang and Raymond A de Caafon Abstract This
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationPREDICTION OF DEFORMED AND ANNEALED MICROSTRUCTURES USING BAYESIAN NEURAL NETWORKS AND GAUSSIAN PROCESSES
PREDICTION OF DEFORMED AND ANNEALED MICROSTRUCTURES USING BAYESIAN NEURAL NETWORKS AND GAUSSIAN PROCESSES C.A.L. Baier-Jones, T.J. Sabin, D.J.C. MacKay, P.J. Withers Department of Materias Science and
More informationBP neural network-based sports performance prediction model applied research
Avaiabe onine www.jocpr.com Journa of Chemica and Pharmaceutica Research, 204, 6(7:93-936 Research Artice ISSN : 0975-7384 CODEN(USA : JCPRC5 BP neura networ-based sports performance prediction mode appied
More informationResearch Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation
Appied Mathematics and Stochastic Anaysis Voume 007, Artice ID 74191, 8 pages doi:10.1155/007/74191 Research Artice On the Lower Bound for the Number of Rea Roots of a Random Agebraic Equation Takashi
More informationarxiv: v1 [cs.lg] 31 Oct 2017
ACCELERATED SPARSE SUBSPACE CLUSTERING Abofaz Hashemi and Haris Vikao Department of Eectrica and Computer Engineering, University of Texas at Austin, Austin, TX, USA arxiv:7.26v [cs.lg] 3 Oct 27 ABSTRACT
More informationTurbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University
Turbo Codes Coding and Communication Laboratory Dept. of Eectrica Engineering, Nationa Chung Hsing University Turbo codes 1 Chapter 12: Turbo Codes 1. Introduction 2. Turbo code encoder 3. Design of intereaver
More informationIncremental Reformulated Automatic Relevance Determination
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 60, NO. 9, SEPTEMBER 22 4977 Incrementa Reformuated Automatic Reevance Determination Dmitriy Shutin, Sanjeev R. Kukarni, and H. Vincent Poor Abstract In this
More informationarxiv: v1 [cs.db] 1 Aug 2012
Functiona Mechanism: Regression Anaysis under Differentia Privacy arxiv:208.029v [cs.db] Aug 202 Jun Zhang Zhenjie Zhang 2 Xiaokui Xiao Yin Yang 2 Marianne Winsett 2,3 ABSTRACT Schoo of Computer Engineering
More informationCategories and Subject Descriptors B.7.2 [Integrated Circuits]: Design Aids Verification. General Terms Algorithms
5. oward Eicient arge-scae Perormance odeing o Integrated Circuits via uti-ode/uti-corner Sparse Regression Wangyang Zhang entor Graphics Corporation Ridder Park Drive San Jose, CA 953 wangyan@ece.cmu.edu
More informationStatistical Astronomy
Lectures for the 7 th IAU ISYA Ifrane, nd 3 rd Juy 4 p ( x y, I) p( y x, I) p( x, I) p( y, I) Statistica Astronomy Martin Hendry, Dept of Physics and Astronomy University of Gasgow, UK http://www.astro.ga.ac.uk/users/martin/isya/
More informationThe EM Algorithm applied to determining new limit points of Mahler measures
Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationSmoothness equivalence properties of univariate subdivision schemes and their projection analogues
Numerische Mathematik manuscript No. (wi be inserted by the editor) Smoothness equivaence properties of univariate subdivision schemes and their projection anaogues Phiipp Grohs TU Graz Institute of Geometry
More informationNumerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process
ISIJ Internationa Vo 54 (2014) No 2 pp 254 258 Numerica Simuation for Optimizing Temperature Gradients during Singe Crysta Casting Process Aeksandr Aeksandrovich INOZEMTSEV 1) Aeksandra Sergeevna DUBROVSKAYA
More information