Stochastic Processes MATH5835, P. Del Moral UNSW, School of - - PowerPoint PPT Presentation

stochastic processes
SMART_READER_LITE
LIVE PREVIEW

Stochastic Processes MATH5835, P. Del Moral UNSW, School of - - PowerPoint PPT Presentation

Stochastic Processes MATH5835, P. Del Moral UNSW, School of Mathematics & Statistics Lectures Notes 1 Consultations (RC 5112): Wednesday 3.30 pm 4.30 pm & Thursday 3.30 pm 4.30 pm 1/25 2/25 Stochastic processes Random


slide-1
SLIDE 1

Stochastic Processes

MATH5835, P. Del Moral UNSW, School of Mathematics & Statistics Lectures Notes 1 Consultations (RC 5112): Wednesday 3.30 pm 4.30 pm & Thursday 3.30 pm 4.30 pm

1/25

slide-2
SLIDE 2

2/25

slide-3
SLIDE 3

Stochastic processes

◮ Random dynamical system. ◮ Sequential simulation of random variables.

3/25

slide-4
SLIDE 4

Stochastic processes

◮ Random dynamical system. ◮ Sequential simulation of random variables. ◮ Randomness:

3/25

slide-5
SLIDE 5

Stochastic processes

◮ Random dynamical system. ◮ Sequential simulation of random variables. ◮ Randomness: occurring with undefinite aim/pattern/regularity, odd

and unpredictable, unknown, unidentified, out of place,. . .

◮ Simulation:

3/25

slide-6
SLIDE 6

Stochastic processes

◮ Random dynamical system. ◮ Sequential simulation of random variables. ◮ Randomness: occurring with undefinite aim/pattern/regularity, odd

and unpredictable, unknown, unidentified, out of place,. . .

◮ Simulation: imitation, mimicking, feigning, pretending,

duplicate/replica/clone, counterfeit, fake,. . .

3/25

slide-7
SLIDE 7

Stochastic processes

◮ Random dynamical system. ◮ Sequential simulation of random variables. ◮ Randomness: occurring with undefinite aim/pattern/regularity, odd

and unpredictable, unknown, unidentified, out of place,. . .

◮ Simulation: imitation, mimicking, feigning, pretending,

duplicate/replica/clone, counterfeit, fake,. . . ⇓

The generation of random numbers is too important to be left to chance.

Robert Coveyou [Studies in Applied Mathematics, III (1970), 70-111.]

3/25

slide-8
SLIDE 8

Lost in the Great Sloan Wall

The Tau Zero Foundation interstellar travels in any dimension. Tony Gonzales random travelling plans in the universe lattices at superluminal speeds

◮ dimension 1 and 2 from Reykjavik:

4/25

slide-9
SLIDE 9

Lost in the Great Sloan Wall

The Tau Zero Foundation interstellar travels in any dimension. Tony Gonzales random travelling plans in the universe lattices at superluminal speeds

◮ dimension 1 and 2 from Reykjavik:

Main drawbacks: infinite returns back home....

4/25

slide-10
SLIDE 10

Lost in the Great Sloan Wall

◮ Free trip travel voucher to the 3d- Great Sloan Wall:

5/25

slide-11
SLIDE 11

Lost in the Great Sloan Wall

◮ Free trip travel voucher to the 3d- Great Sloan Wall:

Main advantage: finite mean returns back home.

5/25

slide-12
SLIDE 12

Lost in the Great Sloan Wall

◮ Free trip travel voucher to the 3d- Great Sloan Wall:

Main advantage: finite mean returns back home. Main drawbacks: Wanders off in the infinite universe and never returns back home!

5/25

slide-13
SLIDE 13

Lost in the Great Sloan Wall - Why??

◮ Was it predictable?

6/25

slide-14
SLIDE 14

Lost in the Great Sloan Wall - Why??

◮ Was it predictable? ◮ What is the stochastic model?

6/25

slide-15
SLIDE 15

Lost in the Great Sloan Wall - Why??

◮ Was it predictable? ◮ What is the stochastic model? ◮ How to simulate it?

6/25

slide-16
SLIDE 16

Lost in the Great Sloan Wall - Why??

◮ Was it predictable? ◮ What is the stochastic model? ◮ How to simulate it? ◮ Is there some math. formulae?

6/25

slide-17
SLIDE 17

Lost in the Great Sloan Wall - Why??

◮ Was it predictable? ◮ What is the stochastic model? ◮ How to simulate it? ◮ Is there some math. formulae?

6/25

slide-18
SLIDE 18

Lost in the Great Sloan Wall - The stochastic model

  • Stoch. model =

7/25

slide-19
SLIDE 19

Lost in the Great Sloan Wall - The stochastic model

  • Stoch. model =Simple random walks on Zd:

d = 1 ⇒ Xn = Xn−1 + Un with Un = +1 or − 1 proba 1/2

7/25

slide-20
SLIDE 20

Lost in the Great Sloan Wall - The stochastic model

  • Stoch. model =Simple random walks on Zd:

d = 1 ⇒ Xn = Xn−1 + Un with Un = +1 or − 1 proba 1/2

◮ dimension 2, 3, and any d ?

7/25

slide-21
SLIDE 21

Lost in the Great Sloan Wall - The stochastic model

  • Stoch. model =Simple random walks on Zd:

d = 1 ⇒ Xn = Xn−1 + Un with Un = +1 or − 1 proba 1/2

◮ dimension 2, 3, and any d ? blackboard.

7/25

slide-22
SLIDE 22

Lost in the Great Sloan Wall - The stochastic model

  • Stoch. model =Simple random walks on Zd:

d = 1 ⇒ Xn = Xn−1 + Un with Un = +1 or − 1 proba 1/2

◮ dimension 2, 3, and any d ? blackboard. ◮ simulation ?

7/25

slide-23
SLIDE 23

Lost in the Great Sloan Wall - The stochastic model

  • Stoch. model =Simple random walks on Zd:

d = 1 ⇒ Xn = Xn−1 + Un with Un = +1 or − 1 proba 1/2

◮ dimension 2, 3, and any d ? blackboard. ◮ simulation ? flip coins!

7/25

slide-24
SLIDE 24

Lost in the Great Sloan Wall - The stochastic model

  • Stoch. model =Simple random walks on Zd:

d = 1 ⇒ Xn = Xn−1 + Un with Un = +1 or − 1 proba 1/2

◮ dimension 2, 3, and any d ? blackboard. ◮ simulation ? flip coins! ◮ Math analysis ?

7/25

slide-25
SLIDE 25

Lost in the Great Sloan Wall - The stochastic model

  • Stoch. model =Simple random walks on Zd:

d = 1 ⇒ Xn = Xn−1 + Un with Un = +1 or − 1 proba 1/2

◮ dimension 2, 3, and any d ? blackboard. ◮ simulation ? flip coins! ◮ Math analysis ?

intuition/blackboard.

7/25

slide-26
SLIDE 26

Lost in the Great Sloan Wall - The stochastic model

  • Stoch. model =Simple random walks on Zd:

d = 1 ⇒ Xn = Xn−1 + Un with Un = +1 or − 1 proba 1/2

◮ dimension 2, 3, and any d ? blackboard. ◮ simulation ? flip coins! ◮ Math analysis ?

intuition/blackboard. ⊕ More rigorous analysis for d = 2, 3 = Project No 1

7/25

slide-27
SLIDE 27

Meeting Alice in Wonderland

Alice + White rabbit ∈ Polygonal labyrinth (no communicating edges!):

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 ◮ Chances to meet after 5 moves : > 25% ◮ Chances to meet after 12 moves : > 51% ◮ Chances to meet after 40 moves : > 99%

8/25

slide-28
SLIDE 28

Meeting Alice in Wonderland - Why??

◮ Was it predictable?

9/25

slide-29
SLIDE 29

Meeting Alice in Wonderland - Why??

◮ Was it predictable? ◮ What is the stochastic model?

9/25

slide-30
SLIDE 30

Meeting Alice in Wonderland - Why??

◮ Was it predictable? ◮ What is the stochastic model? ◮ How to simulate it?

9/25

slide-31
SLIDE 31

Meeting Alice in Wonderland - Why??

◮ Was it predictable? ◮ What is the stochastic model? ◮ How to simulate it? ◮ Is there some math. formulae?

9/25

slide-32
SLIDE 32

Meeting Alice in Wonderland - Why??

◮ Was it predictable? ◮ What is the stochastic model? ◮ How to simulate it? ◮ Is there some math. formulae?

9/25

slide-33
SLIDE 33

Meeting Alice in Wonderland- The stochastic model

  • Stoch. model =

10/25

slide-34
SLIDE 34

Meeting Alice in Wonderland- The stochastic model

  • Stoch. model =Random walk on (finite & complete) graphs

G := (V, E) = (Vertices, Edges)

10/25

slide-35
SLIDE 35

Meeting Alice in Wonderland- The stochastic model

  • Stoch. model =Random walk on (finite & complete) graphs

G := (V, E) = (Vertices, Edges) Neighborhood systems: On the set of vertices x ∼ y = ⇒ (x, y) ∈ E = set of edges ⇓ Set of neighbors of x ∈ V := N(x) = {y ∈ V : y ∼ x}

10/25

slide-36
SLIDE 36

Meeting Alice in Wonderland- The stochastic model

  • Stoch. model =Random walk on (finite & complete) graphs

G := (V, E) = (Vertices, Edges) Neighborhood systems: On the set of vertices x ∼ y = ⇒ (x, y) ∈ E = set of edges ⇓ Set of neighbors of x ∈ V := N(x) = {y ∈ V : y ∼ x} Stochastic model: Xn uniform r.v. on the set of neighbors N(Xn−1)

  • P (Xn = y | Xn−1 = x) =

1 #N(Xn−1) 1N (Xn−1)(y)

10/25

slide-37
SLIDE 37

Meeting Alice in Wonderland - Simulation ⊕ Analysis

◮ Simulation ?

11/25

slide-38
SLIDE 38

Meeting Alice in Wonderland - Simulation ⊕ Analysis

◮ Simulation ? flip coins!

11/25

slide-39
SLIDE 39

Meeting Alice in Wonderland - Simulation ⊕ Analysis

◮ Simulation ? flip coins! ◮ Math analysis ?

11/25

slide-40
SLIDE 40

Meeting Alice in Wonderland - Simulation ⊕ Analysis

◮ Simulation ? flip coins! ◮ Math analysis ?

blackboard.

11/25

slide-41
SLIDE 41

Meeting Alice in Wonderland - Simulation ⊕ Analysis

◮ Simulation ? flip coins! ◮ Math analysis ?

blackboard. Other application domains?

11/25

slide-42
SLIDE 42

Meeting Alice in Wonderland - Simulation ⊕ Analysis

◮ Simulation ? flip coins! ◮ Math analysis ?

blackboard. Other application domains? Ex.: web-graph analysis (ranking ⇒ recommendations)

11/25

slide-43
SLIDE 43

The MIT Blackjack team

Mr M. card shuffle tracking and P. Diaconis magic number ≤ 5 shuffles = ⇒ possible predictions ≥ 90%, 6 shuffles = ⇒ possible predictions ≥ 40% ! Some questions?

12/25

slide-44
SLIDE 44

The MIT Blackjack team

Mr M. card shuffle tracking and P. Diaconis magic number ≤ 5 shuffles = ⇒ possible predictions ≥ 90%, 6 shuffles = ⇒ possible predictions ≥ 40% ! Some questions?

◮ Perfect shuffling ?

12/25

slide-45
SLIDE 45

The MIT Blackjack team

Mr M. card shuffle tracking and P. Diaconis magic number ≤ 5 shuffles = ⇒ possible predictions ≥ 90%, 6 shuffles = ⇒ possible predictions ≥ 40% ! Some questions?

◮ Perfect shuffling ? ⇔ fully unpredictable cards.

12/25

slide-46
SLIDE 46

The MIT Blackjack team

Mr M. card shuffle tracking and P. Diaconis magic number ≤ 5 shuffles = ⇒ possible predictions ≥ 90%, 6 shuffles = ⇒ possible predictions ≥ 40% ! Some questions?

◮ Perfect shuffling ? ⇔ fully unpredictable cards. ◮ A shuffle ?

12/25

slide-47
SLIDE 47

The MIT Blackjack team

Mr M. card shuffle tracking and P. Diaconis magic number ≤ 5 shuffles = ⇒ possible predictions ≥ 90%, 6 shuffles = ⇒ possible predictions ≥ 40% ! Some questions?

◮ Perfect shuffling ? ⇔ fully unpredictable cards. ◮ A shuffle ? ⇔ A permutation of 52 or 78 cards.

12/25

slide-48
SLIDE 48

The MIT Blackjack team

Mr M. card shuffle tracking and P. Diaconis magic number ≤ 5 shuffles = ⇒ possible predictions ≥ 90%, 6 shuffles = ⇒ possible predictions ≥ 40% ! Some questions?

◮ Perfect shuffling ? ⇔ fully unpredictable cards. ◮ A shuffle ? ⇔ A permutation of 52 or 78 cards. ◮ Unpredictable decks ?

12/25

slide-49
SLIDE 49

The MIT Blackjack team

Mr M. card shuffle tracking and P. Diaconis magic number ≤ 5 shuffles = ⇒ possible predictions ≥ 90%, 6 shuffles = ⇒ possible predictions ≥ 40% ! Some questions?

◮ Perfect shuffling ? ⇔ fully unpredictable cards. ◮ A shuffle ? ⇔ A permutation of 52 or 78 cards. ◮ Unpredictable decks ? ⇔ Uniform distribution on Gd (d = 52, 78).

12/25

slide-50
SLIDE 50

The MIT Blackjack team

Mr M. card shuffle tracking and P. Diaconis magic number ≤ 5 shuffles = ⇒ possible predictions ≥ 90%, 6 shuffles = ⇒ possible predictions ≥ 40% ! Some questions?

◮ Perfect shuffling ? ⇔ fully unpredictable cards. ◮ A shuffle ? ⇔ A permutation of 52 or 78 cards. ◮ Unpredictable decks ? ⇔ Uniform distribution on Gd (d = 52, 78). ◮ How to sample the uniform distribution P (σ) = 1/d! ?

12/25

slide-51
SLIDE 51

The MIT Blackjack team

Mr M. card shuffle tracking and P. Diaconis magic number ≤ 5 shuffles = ⇒ possible predictions ≥ 90%, 6 shuffles = ⇒ possible predictions ≥ 40% ! Some questions?

◮ Perfect shuffling ? ⇔ fully unpredictable cards. ◮ A shuffle ? ⇔ A permutation of 52 or 78 cards. ◮ Unpredictable decks ? ⇔ Uniform distribution on Gd (d = 52, 78). ◮ How to sample the uniform distribution P (σ) = 1/d! ? ◮ How many shuffles?

12/25

slide-52
SLIDE 52

The MIT Blackjack team

Mr M. card shuffle tracking and P. Diaconis magic number ≤ 5 shuffles = ⇒ possible predictions ≥ 90%, 6 shuffles = ⇒ possible predictions ≥ 40% ! Some questions?

◮ Perfect shuffling ? ⇔ fully unpredictable cards. ◮ A shuffle ? ⇔ A permutation of 52 or 78 cards. ◮ Unpredictable decks ? ⇔ Uniform distribution on Gd (d = 52, 78). ◮ How to sample the uniform distribution P (σ) = 1/d! ? ◮ How many shuffles?

52! ≃ 8.053 1067 78! ≃ 1.13 10115 >> number 1080of particles in the universe

12/25

slide-53
SLIDE 53

Why shuffling cards?

Application domains?

13/25

slide-54
SLIDE 54

Why shuffling cards?

Application domains?

◮ Software security-design:

Online gambling, iPod songs shuffles.

13/25

slide-55
SLIDE 55

Why shuffling cards?

Application domains?

◮ Software security-design:

Online gambling, iPod songs shuffles.

◮ Cryptography:

Encrypted codes, (pseudo)-random key generators,

13/25

slide-56
SLIDE 56

Why shuffling cards?

Application domains?

◮ Software security-design:

Online gambling, iPod songs shuffles.

◮ Cryptography:

Encrypted codes, (pseudo)-random key generators,

◮ Random search algo:

Simulated annealing (traveling salesman).

13/25

slide-57
SLIDE 57

Why shuffling cards?

Application domains?

◮ Software security-design:

Online gambling, iPod songs shuffles.

◮ Cryptography:

Encrypted codes, (pseudo)-random key generators,

◮ Random search algo:

Simulated annealing (traveling salesman).

◮ Computer sciences:

Reallocations/balancing techniques.

13/25

slide-58
SLIDE 58

Why shuffling cards?

Application domains?

◮ Software security-design:

Online gambling, iPod songs shuffles.

◮ Cryptography:

Encrypted codes, (pseudo)-random key generators,

◮ Random search algo:

Simulated annealing (traveling salesman).

◮ Computer sciences:

Reallocations/balancing techniques.

13/25

slide-59
SLIDE 59

Shuffling cards - The stochastic model?

Random walks on the symmetric group Gd with d = 52, 78, . . .

14/25

slide-60
SLIDE 60

Shuffling cards - The stochastic model?

Random walks on the symmetric group Gd with d = 52, 78, . . . Xn(i) = Value of the i-th card at time n Yn(i) = X −1

n (i) = Position of the card with label i at time n

14/25

slide-61
SLIDE 61

Shuffling cards - The stochastic model?

Random walks on the symmetric group Gd with d = 52, 78, . . . Xn(i) = Value of the i-th card at time n Yn(i) = X −1

n (i) = Position of the card with label i at time n

Change of order σn+1 or change of values i Yn(i) = position

  • σn+1 (Yn(i)) = Yn+1(i) = new position

i Xn(i) = value

  • Xn
  • σ−1

n+1(i)

  • = Xn+1(i) = new value

14/25

slide-62
SLIDE 62

Shuffling cards - The stochastic model?

Random walks on the symmetric group Gd with d = 52, 78, . . . Xn(i) = Value of the i-th card at time n Yn(i) = X −1

n (i) = Position of the card with label i at time n

Change of order σn+1 or change of values i Yn(i) = position

  • σn+1 (Yn(i)) = Yn+1(i) = new position

i Xn(i) = value

  • Xn
  • σ−1

n+1(i)

  • = Xn+1(i) = new value

⇓ Xn = Xn−1 ◦ σ−1

n

  • r

Yn = σn ◦ Yn−1 with some i.i.d. r.v. σn in some class : transpositions, top-in, riffles,. . .

14/25

slide-63
SLIDE 63

Shuffling cards – Simulation+Analysis

◮ Simulation ?

15/25

slide-64
SLIDE 64

Shuffling cards – Simulation+Analysis

◮ Simulation ? flip coins! ∼ permutations/riffles/top-in

(blackboard).

15/25

slide-65
SLIDE 65

Shuffling cards – Simulation+Analysis

◮ Simulation ? flip coins! ∼ permutations/riffles/top-in

(blackboard).

◮ Math analysis ?

15/25

slide-66
SLIDE 66

Shuffling cards – Simulation+Analysis

◮ Simulation ? flip coins! ∼ permutations/riffles/top-in

(blackboard).

◮ Math analysis ?

convergence to stationarity. Xn and Yn − →n→∞ limiting r.v. X∞ and Y∞ with for any uniform r.v. σ on Gd X∞

in law

= X∞ ◦ σ−1 and Y∞

in law

= σ ◦ Y∞

15/25

slide-67
SLIDE 67

Shuffling cards – Simulation+Analysis

◮ Simulation ? flip coins! ∼ permutations/riffles/top-in

(blackboard).

◮ Math analysis ?

convergence to stationarity. Xn and Yn − →n→∞ limiting r.v. X∞ and Y∞ with for any uniform r.v. σ on Gd X∞

in law

= X∞ ◦ σ−1 and Y∞

in law

= σ ◦ Y∞ ⊕ Top-in & Transpositions = Project No 2

15/25

slide-68
SLIDE 68

Kruskal count (with the classroom!)

16/25

slide-69
SLIDE 69

Kruskal count - Stochastic model ⊕ Simulation ⊕ Analysis

Martin David Kruskal (1925-2006)

◮ Stochastic model ?

  • 17/25
slide-70
SLIDE 70

Kruskal count - Stochastic model ⊕ Simulation ⊕ Analysis

Martin David Kruskal (1925-2006)

◮ Stochastic model ?

Deterministic walker on a random environment.

◮ Simulation ?

17/25

slide-71
SLIDE 71

Kruskal count - Stochastic model ⊕ Simulation ⊕ Analysis

Martin David Kruskal (1925-2006)

◮ Stochastic model ?

Deterministic walker on a random environment.

◮ Simulation ? Permutation sampling!

17/25

slide-72
SLIDE 72

Kruskal count - Stochastic model ⊕ Simulation ⊕ Analysis

Martin David Kruskal (1925-2006)

◮ Stochastic model ?

Deterministic walker on a random environment.

◮ Simulation ? Permutation sampling! ◮ Math analysis ?

17/25

slide-73
SLIDE 73

Kruskal count - Stochastic model ⊕ Simulation ⊕ Analysis

Martin David Kruskal (1925-2006)

◮ Stochastic model ?

Deterministic walker on a random environment.

◮ Simulation ? Permutation sampling! ◮ Math analysis ?

Simplified model on blackboard.

17/25

slide-74
SLIDE 74

Kruskal count - Stochastic model ⊕ Simulation ⊕ Analysis

Martin David Kruskal (1925-2006)

◮ Stochastic model ?

Deterministic walker on a random environment.

◮ Simulation ? Permutation sampling! ◮ Math analysis ?

Simplified model on blackboard.

17/25

slide-75
SLIDE 75

The magic fern from Daisetsuzan

18/25

slide-76
SLIDE 76

The magic fern from Daisetsuzan

◮ What is the stochastic model?

18/25

slide-77
SLIDE 77

The magic fern from Daisetsuzan

◮ What is the stochastic model? ◮ How to simulate it?

18/25

slide-78
SLIDE 78

The magic fern from Daisetsuzan

◮ What is the stochastic model? ◮ How to simulate it? ◮ Is there some math. formulae?

18/25

slide-79
SLIDE 79

The magic fern from Daisetsuzan

◮ What is the stochastic model? ◮ How to simulate it? ◮ Is there some math. formulae?

18/25

slide-80
SLIDE 80

The magic fern - The stochastic model

  • Stoch. model = Iterated random functions (IRF):

◮ Some functions x ∈ R2 → fi(x) = Ai x + bi ∈ R2, with

i ∈ I = {1, . . . , d}.

◮ Some i.i.d. r.v. ǫn with some law µ on {1, . . . , d}

19/25

slide-81
SLIDE 81

The magic fern - The stochastic model

  • Stoch. model = Iterated random functions (IRF):

◮ Some functions x ∈ R2 → fi(x) = Ai x + bi ∈ R2, with

i ∈ I = {1, . . . , d}.

◮ Some i.i.d. r.v. ǫn with some law µ on {1, . . . , d}

⇓ Xn = fǫn (Xn−1)

19/25

slide-82
SLIDE 82

The magic fern - The stochastic model

  • Stoch. model = Iterated random functions (IRF):

◮ Some functions x ∈ R2 → fi(x) = Ai x + bi ∈ R2, with

i ∈ I = {1, . . . , d}.

◮ Some i.i.d. r.v. ǫn with some law µ on {1, . . . , d}

⇓ Xn = fǫn (Xn−1) =

  • fǫn ◦ fǫn−1 ◦ . . . ◦ fǫ1
  • (x0)

19/25

slide-83
SLIDE 83

The magic fern - The stochastic model

  • Stoch. model = Iterated random functions (IRF):

◮ Some functions x ∈ R2 → fi(x) = Ai x + bi ∈ R2, with

i ∈ I = {1, . . . , d}.

◮ Some i.i.d. r.v. ǫn with some law µ on {1, . . . , d}

⇓ Xn = fǫn (Xn−1) =

  • fǫn ◦ fǫn−1 ◦ . . . ◦ fǫ1
  • (x0)

Do you believe this? ⊕ Fractal & IRF = Project No 3

19/25

slide-84
SLIDE 84

The magic fern from Daisetsuzan - An illustration

Scilab program: fractal.tree.sce A1 = c

  • b1 =

1/2

  • A2 =

r cos(ϕ) −r sin(ϕ) r sin(ϕ) r cos(ϕ)

  • b2 =
  • 1

2 − r 2 cos(ϕ)

c − r

2 sin(ϕ)

  • et

A3 = q cos(ψ) −r sin(ψ) q sin(ψ) r cos(ψ)

  • b3 =
  • 1

2 − q 2 cos(ψ) 3c 5 − q 2 sin(ψ)

  • with ǫn i.i.d. uniform on {1, 2, 3}, and with the parameters

c = 0.255, r = 0.75, q = 0.625 ϕ = −π 8 , ψ = π 5 , |X0| ≤ 1.

20/25

slide-85
SLIDE 85

The Kepler-22b Eve

Migration in 2457 of selected 1000 humans: Reproduction rate (20 years)

◮ 5.5 thousands years ⇒ 25% population ∈ same family. ◮ 6.5 thousands years ⇒ 68% population ∈ same family. ◮ 10 thousands years ⇒ more than 99% population ∈ same family.

21/25

slide-86
SLIDE 86

The Kepler-22b Eve

Migration in 2457 of selected 1000 humans: Reproduction rate (20 years)

◮ 5.5 thousands years ⇒ 25% population ∈ same family. ◮ 6.5 thousands years ⇒ 68% population ∈ same family. ◮ 10 thousands years ⇒ more than 99% population ∈ same family.

The Seven Daughters of Eve of Bryan Sykes:

21/25

slide-87
SLIDE 87

The Kepler-22b Eve

Migration in 2457 of selected 1000 humans: Reproduction rate (20 years)

◮ 5.5 thousands years ⇒ 25% population ∈ same family. ◮ 6.5 thousands years ⇒ 68% population ∈ same family. ◮ 10 thousands years ⇒ more than 99% population ∈ same family.

The Seven Daughters of Eve of Bryan Sykes: Mitochondrial Eve ≃ 140 − 200 thousands years (missprint in manuscript)

21/25

slide-88
SLIDE 88

The Kepler-22b Eve - Stochastic model

Random walks on functions a : i ∈ {1, . . . , d} → {1, . . . , d}

22/25

slide-89
SLIDE 89

The Kepler-22b Eve - Stochastic model

Random walks on functions a : i ∈ {1, . . . , d} → {1, . . . , d} Birth and death = selection of the parents/ancestors/types

A0

  • A1
  • A2
  • A3
  • 1

1 1 1 1 2 2 2 2 2 3 3 3 3 3 4 4 4 4 4 5 5 5 5 5 6 6 6 6 6

22/25

slide-90
SLIDE 90

The Kepler-22b Eve - Simulation ⊕ Analysis

◮ Simulation?

23/25

slide-91
SLIDE 91

The Kepler-22b Eve - Simulation ⊕ Analysis

◮ Simulation?

Neutral genetic model Each A(i) uniform on {1, . . . , d}

◮ Analysis?

23/25

slide-92
SLIDE 92

The Kepler-22b Eve - Simulation ⊕ Analysis

◮ Simulation?

Neutral genetic model Each A(i) uniform on {1, . . . , d}

◮ Analysis? Iterated random mappings.

The label of the ”surviving” ancestors of the population at time n: Xn = Xn−1 ◦ An = A0 ◦ . . . ◦ An : i ∈ {1, . . . , d} → {1, . . . , d} Intuitively: Xn − →n→∞ X∞ = Constant random function

23/25

slide-93
SLIDE 93

The Kepler-22b Eve - Simulation ⊕ Analysis

◮ Simulation?

Neutral genetic model Each A(i) uniform on {1, . . . , d}

◮ Analysis? Iterated random mappings.

The label of the ”surviving” ancestors of the population at time n: Xn = Xn−1 ◦ An = A0 ◦ . . . ◦ An : i ∈ {1, . . . , d} → {1, . . . , d} Intuitively: Xn − →n→∞ X∞ = Constant random function ⊕ More rigorous analysis = Project No 4

23/25

slide-94
SLIDE 94

Poisson typos (reminder of stats course)

◮ Arrival times Tn of events :

misprints in a text, trades counts, bus arrivals, machine failures, catastrophies, number of tries in rugby games, . . .

◮ Stat. model on the time axis = i.i.d. exponential inter-arrivals

P ((Tn+1 − Tn) ∈ dt | Tn) = e− t 1[0,∞[(t) dt

24/25

slide-95
SLIDE 95

◮ Consequences (admitted, cf. pb.1 p. 71)

∀n ≥ 0 P (#{Ti ∈ [0, λ]} = n) = e−λ λn n! = Poisson r.v. and T1 Tn+1 , T2 Tn+1 , . . . , Tn Tn+1

  • = ordered unif. stats on [0, 1] ⊥ Tn+1

◮ Extension to any state space (⊃ the book of D. Poisson)?

25/25

slide-96
SLIDE 96

◮ Consequences (admitted, cf. pb.1 p. 71)

∀n ≥ 0 P (#{Ti ∈ [0, λ]} = n) = e−λ λn n! = Poisson r.v. and T1 Tn+1 , T2 Tn+1 , . . . , Tn Tn+1

  • = ordered unif. stats on [0, 1] ⊥ Tn+1

◮ Extension to any state space (⊃ the book of D. Poisson)?

Counting processes with intensity x ∈ S → f (x) ∈ [0, ∞[ X =

  • 1≤i≤N

δX i with

◮ N Poisson r.v. with λ =

  • f (x) dx.

◮ Given N = n, X 1, . . . , X n i.i.d. with density g(x) = f (x)/λ. 25/25