Simulation of Markovian models using Bootstrap method Summer Computer Simulation Conference (SCSC) 2010
Ricardo M. Czekster, Paulo Fernandes, Afonso Sales, Dione Taschetto and Thais Webber Pontifícia Universidade Católica do Rio Grande do Sul (PUCRS) PaleoProspec Project - PUCRS/Petrobras Funded also by CAPES and CNPq - Brazil
Outline
Introduction Markovian model simulation Bootstrap simulation Numerical results Conclusion
2 / 47
Approach modeling formalism
System
solve
Performance indices
Model
improve
Well-known modeling formalism Markov Chains
3 / 47
Approach modeling formalism
System
solve
Performance indices
Model
improve
Well-known modeling formalism Markov Chains
4 / 47
Markov Chains formalism Simple primitives: states and transitions ◮
Biology Physics Social sciences Business Computer science and telecommunication
5 / 47
Markov Chains Advantage: intuitive and easy to model “small” systems Limitation: representing huge models by state/transition form becomes intractable (state space explosion problem)
Structured formalisms Model description by components
Structured mathematical description
Examples Queueing Networks (QN) [Little61, Basket et al. 75, Reiser et al. 80] Stochastic Petri Nets (SPN) [Florin et al. 85] Performance Evaluation Process Algebra (PEPA) [Hillston95] Stochastic Automata Networks (SAN) [Plateau84] 6 / 47
Model solution It plays equally an important role Numerical solution (iterative methods): ◮ ◮ ◮
Power method [Stewart94] Arnoldi [Arnoldi51] GMRES [Saad and Schultz86]
Simulation (model’s events firing transitions according to pseudorandom number generator): ◮ ◮ ◮
Traditional [Ross96] Monte Carlo [Häggström02] Backward [Propp and Wilson 96]
7 / 47
Simulation methods They suffer from a lack of precision in the results due to the very nature of simulation ◮ ◮
Based on a pseudorandom number generator When compared to the analytical solution
Bootstrap method It is a well known statistical method applied to many fields to improve accuracy when performing sample estimations for complex distributions
Our idea Apply the Bootstrap method in the simulation context in order to reduce the “noise” that were observed from the samples
8 / 47
Bootstrap method ˜ Λ
Λ
¯xΛ
K1
K2
¯xK1
¯xK2
...
Kz
¯xKz
¯xK
˜ infinite sized set Λ: ˜ Λ: subset (sampling) of Λ K: bootstrap (with the same size of Λ)
These set of results, in average, are better than the average of a single sampling 9 / 47
Bootstrap method ˜ Λ
Λ
¯xΛ
K1
K2
¯xK1
¯xK2
...
Kz
¯xKz
¯xK
˜ infinite sized set Λ: ˜ Λ: subset (sampling) of Λ K: bootstrap (with the same size of Λ)
These set of results, in average, are better than the average of a single sampling 10 / 47
each visited state = sample mean permanence ′ probability π = πn
28 / 47
Bootstrap simulation States
Transition Matrix
0 Initial state
0 Time 0
1
2
0 0.10 0.65 0.25 1 0.25 0.55 0.20 2 0.30 0.25 0.45 n = trajectory length K: bootstrap z: number of bootstraps mean permanence Pz ¯ xi probability π = i=1 z
n = trajectory length K: bootstrap z: number of bootstraps mean permanence Pz ¯ xi probability π = i=1 z
38 / 47
Bootstrap simulation States
2
Transition Matrix 1 ...
1 Initial state
0
0 0
1
K1
0 1 2
0
K2
0 1 2
2
3
¯x1
Kz
0 ... 1 2
Time 5 ... n
4
normalize
0 1 2
¯x2
0 1 2
¯xz
0 ... 1 2
0 1 2 0 0.10 0.65 0.25 1 0.25 0.55 0.20 2 0.30 0.25 0.45 n = trajectory length K: bootstrap z: number of bootstraps mean permanence Pz ¯ xi probability π = i=1 z
39 / 47
Bootstrap simulation States
2
Transition Matrix 1 ...
1 Initial state
0
0 0
1
K1
0 1 2
0 2
K2
0 1 2
3
¯x1
Kz
0 ... 1 2
Time 5 ... n
4
normalize
0 1 2
¯x2
π ¯x1[0] + ¯x2[0] + ··· + ¯xz[0] z ¯x1[1] + ¯x2[1] + ··· + ¯xz[1] z
= π0
¯x1[2] + ¯x2[2] + ··· + ¯xz[2] z
= π2
= π1
¯xz
0 ... 1 2
0 1 2
0 1 2 0 0.10 0.65 0.25 1 0.25 0.55 0.20 2 0.30 0.25 0.45 n = trajectory length K: bootstrap z: number of bootstraps mean permanence Pz ¯ xi probability π = i=1 z
40 / 47
Examples Alternate Service Patterns (ASP): describes an Open Queueing Network with servers that map P different service patterns. First Available Server (FAS): indicates the availability of N servers, where tasks are firstly assigned to the first server, case it is available. If the server is busy, the task must be sent to the second server and so on. Resource Sharing (RS): maps R shared resources to P processes.
Mean relative error It is computed comparing the simulation results with the numerical results obtained by the Power method
41 / 47
Results: ASP model 1
0.0082 0.0069
0.001
104
105
106
107
108
0.0009
0.0028 0.0061
0.0902 0.0285
0.01
0.0272 0.0110
0.3740
0.2608 0.1017
Mean relative error
0.1
Traditional Bootstrap
109
Trajectory length (n)
42 / 47
Results: FAS model
0.01
104
105
106
107
108
0.0146
0.0485 0.0383
0.1429 0.0582
0.1
0.4680 0.1552
0.8672 0.4214
1.2876
Mean relative error
1
1.5182 0.9986
Traditional Bootstrap
109
Trajectory length (n)
43 / 47
Results: RS model 1
104
105
106
107
108
0.0011
0.001
0.0034 0.0010
0.0108 0.0033
0.01
0.0307 0.0110
0.1071 0.0325
0.3487 0.1063
0.3151
Mean relative error
0.1
Traditional Bootstrap
109
Trajectory length (n)
44 / 47
Conclusion Simulation is an alternative solution to solve Markovian models when analytical solution becomes intractable Due to the nature of simulation (pseudorandom number generator) there is a lack of precision in the results A large quantity of samples is necessary in order to achieve better precision results Bootstrap simulation can produce a more precise result using a smaller amount of samples than Traditional simulation
Future works Model classification: determination of classes of models more suitable for the Bootstrap simulation Parallel solution: distribute the bootstraps for different processors/computers in a grid Compare with other methods: Monte Carlo, Backward, ... 45 / 47
Conclusion Simulation is an alternative solution to solve Markovian models when analytical solution becomes intractable Due to the nature of simulation (pseudorandom number generator) there is a lack of precision in the results A large quantity of samples is necessary in order to achieve better precision results Bootstrap simulation can produce a more precise result using a smaller amount of samples than Traditional simulation
Future works Model classification: determination of classes of models more suitable for the Bootstrap simulation Parallel solution: distribute the bootstraps for different processors/computers in a grid Compare with other methods: Monte Carlo, Backward, ... 46 / 47
equally important role when examining complex systems. There are ... Although numerical solution produces reliable results, it is .... transition matrix storage [9].
packet network traffic. Recent studies had led to the conclusion that the Ethernet, ATM traffic, telnet,. FTP and variable-bit-rate (VBR) video traffic can be more ...
The paper presents a simple mathematical model for water pollution. We consider the advection diffusion equation as an Initial Boundary Value Problem (IBVP) for the estimation of water pollution. For the numerical solution of the IBVP, the derivation
IJRIT International Journal of Research in Information Technology, Volume 2, .... discrete approximation to this functions) which satisfies a given relationship ...
IJRIT International Journal of Research in Information Technology, Volume 2, .... discrete approximation to this functions) which satisfies a given relationship ...
take place freely in the open area. Indeed ... Software Organization & MonIToring of Ambient Systems Work- shop. ..... SUMO is an open-source application im-.
of mobile nodes tracking a particular target, which may be either static or move ... it uses real maps obtained from the TIGER/Lines database. [11]. For each route ...
Jan 10, 2012 - Contact email: [email protected]. The authors ..... before, we add shadow shocks to each of these equations (giving a total of ..... Papers. CEPREMAP, April. http://ideas.repec.org/p/cpm/dynare/001.html.
Keywords: Partial identification, incomplete specification test, duality, dilation ...... 15 minutes for n = 1000 on an AMD Opteron (tm) Processor 250 with 4G of ...
Corresponding to a bootstrap resample Ïâ is a bootstrap replication ...... univariate bootstrap estimation of bias and variance for an arbitrary statistic, theta. It.
Jan 27, 2005 - parametric neural netWork classi?ers that assume little a priori information. See K. ..... in the spectral mea surement. 2=PT(ttT)'1P+diag(o2). (9).
Jan 27, 2005 - tion sequence organizes spectral data into clusters having a high degree ..... HoWever, cluster analysis does not utilize a priori informa tion and ...
These DNS are usually implemented via spectral (pseu- dospectral ... In this domain, the flow is naturally periodic along azimuthal (θ) direction and taken to be ...
LMP calculated from the ACOPF algorithm and outperforms the conventional lossless DCOPF algorithm. This is reasonable since the FND model considers the ...
Try one of the apps below to open or edit this item. pdf-1464\steady-state-simulation-of-an-oil-refinery-using-commercial-software-by-gerald-l-kaes.pdf.
However, even quadratic speedup is considerable when N is large. Like all quantum computer algorithms, Grover's algorithm is probabilistic, in the sense that it.
Jul 20, 2010 - models from medical imaging data have improved substantially in the last ..... real-time visualization of segmentation results. ...... 14th Annual Conference on Computer Graphics and Interactive Techniques, 1987; 163â169.
Modern RF chips for mobile devices, for instance, typically .... Using this simple and numerically cheap method, one can do many kinds of analysis ... Floor plan with relocation option that was considered after nonlinear phase noise analysis ...
CFD platform, in: accepted the 12th International Topical Meeting on. Nuclear Reactor Thermal Hydraulics (NURETH-12), Pittsburgh, Pennsyl- vania, USA ...