Random Walk with Anti-Correlated Steps

Similar documents
An introduction to the theory of SDDP algorithm

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

BBP-type formulas, in general bases, for arctangents of real numbers

Lecture 2 October ε-approximation of 2-player zero-sum games

RC, RL and RLC circuits

Linear Response Theory: The connection between QFT and experiments

Vehicle Arrival Models : Headway

Operating Systems Exercise 3

The field of mathematics has made tremendous impact on the study of

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

1 Review of Zero-Sum Games

Math 333 Problem Set #2 Solution 14 February 2003

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Electrical and current self-induction

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

A Shooting Method for A Node Generation Algorithm

Some Basic Information about M-S-D Systems

Christos Papadimitriou & Luca Trevisan November 22, 2016

Spring Ammar Abu-Hudrouss Islamic University Gaza

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

20. Applications of the Genetic-Drift Model

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

Stability and Bifurcation in a Neural Network Model with Two Delays

ON THE DEGREES OF RATIONAL KNOTS

Asymptotic Equipartition Property - Seminar 3, part 1

Matrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality

Math 2214 Solution Test 1B Fall 2017

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

Supplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence

EXERCISES FOR SECTION 1.5

6.2 Transforms of Derivatives and Integrals.

Predator - Prey Model Trajectories and the nonlinear conservation law

More Digital Logic. t p output. Low-to-high and high-to-low transitions could have different t p. V in (t)

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

d 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3

STATE-SPACE MODELLING. A mass balance across the tank gives:

Chapter 7: Solving Trig Equations

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

Reading from Young & Freedman: For this topic, read sections 25.4 & 25.5, the introduction to chapter 26 and sections 26.1 to 26.2 & 26.4.

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.

Estimation of Poses with Particle Filters

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

R.#W.#Erickson# Department#of#Electrical,#Computer,#and#Energy#Engineering# University#of#Colorado,#Boulder#

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

Longest Common Prefixes

Lab 10: RC, RL, and RLC Circuits

Solutions for Assignment 2

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

Challenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k

5.1 - Logarithms and Their Properties

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

Math-Net.Ru All Russian mathematical portal

Ensamble methods: Bagging and Boosting

Y 0.4Y 0.45Y Y to a proper ARMA specification.

KINEMATICS IN ONE DIMENSION

Notes for Lecture 17-18

10. State Space Methods

Physics for Scientists & Engineers 2

Dynamic Programming 11/8/2009. Weighted Interval Scheduling. Weighted Interval Scheduling. Unweighted Interval Scheduling: Review

Learning Objectives: Practice designing and simulating digital circuits including flip flops Experience state machine design procedure

5.2. The Natural Logarithm. Solution

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

Chapter 2: Principles of steady-state converter analysis

Suggested Solutions to Assignment 4 (REQUIRED) Submisson Deadline and Location: March 27 in Class

Solutionbank Edexcel AS and A Level Modular Mathematics

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

Particle Swarm Optimization

Ensamble methods: Boosting

Problem Set #3: AK models

1. VELOCITY AND ACCELERATION

A Hop Constrained Min-Sum Arborescence with Outage Costs

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

Expert Advice for Amateurs

Orientation. Connections between network coding and stochastic network theory. Outline. Bruce Hajek. Multicast with lost packets

Approximation Algorithms for Unique Games via Orthogonal Separators

ON THE BEAT PHENOMENON IN COUPLED SYSTEMS

Average Number of Lattice Points in a Disk

Testing for a Single Factor Model in the Multivariate State Space Framework

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

On a Fractional Stochastic Landau-Ginzburg Equation

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Zhihan Xu, Matt Proctor, Ilia Voloh

Deep Learning: Theory, Techniques & Applications - Recurrent Neural Networks -

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

Appendix to Online l 1 -Dictionary Learning with Application to Novel Document Detection

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

We just finished the Erdős-Stone Theorem, and ex(n, F ) (1 1/(χ(F ) 1)) ( n

Some Ramsey results for the n-cube

Math 315: Linear Algebra Solutions to Assignment 6

15. Vector Valued Functions

i L = VT L (16.34) 918a i D v OUT i L v C V - S 1 FIGURE A switched power supply circuit with diode and a switch.

3.1 More on model selection

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS

1. Kinematics I: Position and Velocity

The motions of the celt on a horizontal plane with viscous friction

Transcription:

Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and experimenal daa. The experimenal analysis includes he compuaion of he expeced values of random walks for seps up o 22. The resul shows he expeced value asympoically converging o. Inroducion Le {µ } be he horizonal random walk 3 ha obeys he following rules:. The walk sars a µ 0 2. An urn conains whie marbles and red marbles. A each sep a marble is seleced from ha urn wihou replacemen. 3. If he marble chosen is red µ increases by one (µ + µ + ) 4. If he marble chosen is whieµ decreases by one (µ + µ -) 5. Le p/3 be he probabiliy o draw a red marble. Accordingly, q2/3 is he probabiliy o draw a whie marble. n is he oal number of marbles in ha urn. Noe ha he absolue maximum posiive is while he walk bounded by in he negaive direcion. Noe also he walk will erminae a - This ype of random walks involve ani-correlaed seps, meaning he more red marbles are drawn (posiive direcion) he more likely are seps ino he negaive direcion (whie marbles are drawn). Ani-correlaed walks are implemened simply be no replacing drawn marbles. We are ineresed in he expeced value E µ ] of {µ }, ha is for he expeced value of a random walk wih ani-correlaed seps. ] proves E µ ] O(). We are going o refine his resul by showing ha he expeced value of random walks wihou replacemen is likely o be exacly, hus E µ ]. Applicaions for random walks wih replacemen can be found in caching algorihms. ] inroduces he RMARK algorihm a random walk which cos is E µ ]. 3] I can be shown ha for random walks wih replacemen E µ ]. Deparmen of Compuer Science, California Sae Universiy of Norhridge, CA 9330. Email: jnoga@ecs.csun.edu 2 Deparmen of Compuer Science, California Sae Universiy of Norhridge, CA 9330. Email: wagiboy@gmail.com 3 Sricly pu we consider random processes because he probabiliy of seps changes wih ime.

2 Conjecure E µ ] We conjecure E µ ] for a random walk wih ani-correlaed sep size. Suppor for his conjecure comes from wo plausibiliy argumen and our experimenal daa. Plausibiliy argumen one: Le } be random walk wih replacemen wih p/3, he probabiliy o draw a red marble. Noe ha p remains consan during he walk. 4] shows he expeced value of random walks wih replacemen is exacly, in shor E ν ]. For random walks wihou replacemen p is changing, in fac p is geing smaller as he walk progresses: A p 0 /3. Afer an equal number of red and whie marbles have been drawn µ reurns o zero. A his poin he probabiliy ha a red marble will be drawn nex is no greaer han /3. Le p he probabiliy o draw a red marble, hen p < p0. The probabiliy of subsequen reurns of µ o zero exhibi an even furher decrease such ha p > p > L. Noe, if whie marbles are drawn subsequen p do increase bu a no 0 p2 effec o E µ ] because µ moves in he negaive direcion. Now when he probabiliy p of drawing reds decreases over he course of a random walk, i is plausible ha he expeced value of random walks wihou replacemen is no greaer han for random walks wih replacemen, or E µ ] E ν ]. max max Plausibiliy argumen wo: As grows large drawing marbles does affec p very lile. In fac, as δ p does no change a all. In oher words he random walk wihou replacemen virually becomes a random walk wih replacemen, hus E µ ] E ν ]. max max 3 Experimenal Mehods In his secion we provide experimenal suppor for our conjecure E max µ ] by compuing E µ ] for random walk of up o 30 seps. We implemened and compared several algorihms o efficienly compue he problem, including exhausive compuing, recursively selecing represenaives, ieraively selecing represenaives and linear opimizaions. As experimenal mehod are implemened and compared on a.6 GHz Penium M machine wih 52 MB of main memory. The final experimenal daa has been compued on a 3.2??? GHz Penium 4 wih hyper hreading and 2 GB??? of main memory. 3. Exhausive Mehod We used he following pseudo code o exhausively compue E µ ]. inpu n 3/2* oal_max_µ0; for each of he n! permuaions µ0; max_µ0; for each marble permuaion if marble is red µ++; if( µ>max_µ ) max_µµ // find maximum µ

E ] This mehod enabled us o comforably compue max µ for he values A he program runime on he (slower).6 GHz Penium M machine amouns o 3 ½ hours. As exhausive compuing represens only a crude firs ake on he problem no furher commens need o be made. 3.2 Recursive Selecion Represenaive Mehod Compuing E µ ] involves generaing permuaions of {µ }, for example µ r, w, w 2, r 2, r 3, w 3, w 4, w 5. These permuaions can be grouped by repeaing subsequences of w s and r s, for example µ r w 2 r 2 w 3. Each group represens permuaions of which each generaes he same value µ max. In oher words he selecing one represenaive from group speeds up he compuaion of E µ ] by a facor of We void developed Selec_k_ouof_n( a recursive algorihm in level, o generae in he loopsar selecion ) of elemens ou of n { E ] which is shown here: Implemenaion of his algorihm enables us o compue max µ for( in iloopsar; i<m_n-m_k+2+level; i++ ) for up { o he wihin a few hours on a.6 GHz Penium machine. m_selecionlevel]i; 3.3 Ieraive Selecion Represenaive Mehod SimulaeRandomWalk( m_selecion ); While he recursive else selecion mehod is shor and concise, i suffers from he overhead incurred by frequen Selec_k_ouof_n( he recursive funcion level+, calls and i+ exiss. ); For his reason we convered } our recursive selecion mehod ino an ieraive. } if( level m_finallevel ) Uilizing an ieraive mehod we here able o push he envelope o for compuaions lasing less han one hour. 3.4 Linear Opimizaions End Walk a There is no need o le he algorihm simulae all drawings of µ. A 3/2 of he walk has been compleed. A his poin here is no chance ha E µ ] can grow, simply because here are no enough red marbles lef o advance µ pas 0 ino he posiive range. This implemenaion opimizaion speeds up he compuaion linearly by /3. Skip E max µ ] 0 Selecions The selecion generaing algorihm produces is sequences in lexicographic order. This fac enables us rim down he selecions from he poin forward when sufficienly whie marbles have been drawn prevening µ o become posiive. More precisely, he firs red marble is drawn a /2 followed by a whie marble. This implemenaion opimizaion speeds up he compuaion by???.

4 Experimenal Resuls We have compued E µ ] as a funcion of he up o 30. As he graph in figure shows E µ ] appears o be converging o he asympoic boundary. Figure also shows he graph of E ν ]. Noe ha for he 2 graphs are geing closer and closer, appearing boh o be converging o. This converging behavior experimenally suppors plausibiliy argumen 2 ha in he limi E µ ] E ν ] max max Expeced Maximum Value of Mu.2 0.8 Emax mu] 0.6 0.4 Emax mu] Emax nu] 0.2 0 2 4 6 8 0 2 4 6 8 20 22 dela Table conains he numerical daa used o creae he graph of E µ ] in figure. I also conains he columns for he number of red and whie marbles. The resuls of E ν ] are augmened wih exac raios are heir respecive run imes. #whie #red Run ime E µ ] E µ ] 2 2 0.30 secs 5/27 0.333333 4 4 2 0.30 secs 524/729 0.466666 6 6 3 0.30 secs 607/9683 0.54769 8 8 4 0.30 secs 0.604040 0 0 5 0.30 secs 0.646354 2 2 6.26 secs 0.679595

4 4 7 3.64 secs 0.706570 6 6 8 2.66 secs 0.728995 8 8 9 69 secs 0.747990 (3 min) 20 20 0 342 secs 0.764323 (22 min) 22 22 930 secs (2 hrs 35 m) 0.77854 Table 2 presen daa for E ν ] #whie #red Run ime E ν ] E ν ] 2 2 0.30 secs /3 0.555556 4 4 2 0.30 secs 7/5 0.78793 6 6 3 0.30 secs 46/84 0.83748 8 8 4 0.30 secs 0.86869 0 0 5 0.30 secs 0.906453 2 2 6.26 secs 0.9339 4 4 7 3.64 secs 0.9499 6 6 8 2.66 secs 0.9667 8 8 9 69 secs (3 min) 0.9707 20 20 0 342 secs (22 min) References ] Rajeev Mowani, Randomized Algorihms, Cambridge Universiy, 995 2] Marek Chrobak, Elias Kousoupias, John Noga. More on Randomized On-line Algorihms for Caching. Theoreical Compuer Science, Denmark, volume 3, pages 997-2008, 2003 3] Yair Baral, Marek Chrobak, John Noga, Prabhakar Raghavan. More on Random Walks, Elecrical Neworks and Harmonic k-server Algorihms. Informaion Processing Leers, volume 84, pages 27-276, 2002 4] Feller. An Inroducion o Probabiliy Theroy and Is Applicaions. Wiley & Sons, Inc. New York, volume, pages 282 o 288, 950