570 likes | 717 Views
Chasing peaks. Stochastic motion rate. Ornstein-Uhlenbeck mean. Continuous state. Brian O’Meara http:// www.brianomeara.info. and/or. Directionally: be pulled towards some value. In what ways can a continuous trait change in an instant of time?.
E N D
Chasing peaks Stochastic motion rate Ornstein-Uhlenbeck mean Continuous state Brian O’Meara http://www.brianomeara.info
and/or Directionally: be pulled towards some value In what ways can a continuous trait change in an instant of time? Randomly: increase or decrease slightly by chance
In what ways can a continuous trait change in an instant of time? Randomly: increase or decrease slightly by chance and/or Directionally: be pulled towards some value
In what ways can a continuous trait change in an instant of time? Randomly: increase or decrease slightly by chance and/or Directionally: be pulled towards some value
Rate of wiggle In what ways can a continuous trait change in an instant of time? Randomly: increase or decrease slightly by chance and/or Directionally: be pulled towards some value
In what ways can a continuous trait change in an instant of time? Randomly: increase or decrease slightly by chance and/or Directionally: be pulled towards some value
Adds the entire difference In what ways can a continuous trait change in an instant of time? Randomly: increase or decrease slightly by chance and/or Directionally: be pulled towards some value
In what ways can a continuous trait change in an instant of time? Randomly: increase or decrease slightly by chance and/or Directionally: be pulled towards some value Allows directional change less than 100% (even zero)
In what ways can a continuous trait change in an instant of time? Randomly: increase or decrease slightly by chance and/or Directionally: be pulled towards some value Ornstein-Uhlenbeck process
1 1 1 1 1
2 2 2 2 2
3 3 3 3 3
And so forth... 3 3 3 3 3
Brownian rate OU attraction OU mean
Brownian rate OU attraction OU mean
Brownian rate OU attraction OU mean
Info and collage from Luke Mahler, http://lukemahler.com/. Photos by J. Losos, B.Falk, M. Landestoy, and L. Mahler
∆ Butler & King 2004
Build up (paint regimes one at a time) Merge SURFACE: Ingram & Mahler 2013
Brownian rate OU attraction OU mean
Brownian rate OU attraction OU mean
Brownian motion OU or OU-like Hansen & Martins 1996
What is the probability of heads for our coin? Sim 1: HTTTTHTT Sim 2: TTHHTTTT Sim 3: TTTHHTHH Sim 4: TTTTTTTT Sim 5: TTTTTTTT Sim 6: THHTTTTH Sim 7: TTHHTTHT Sim 8: TTTTHHHT Sim 9: HTTHTTTT Sim 10: THHHTHTT Let p = 0.2
What is the probability of heads for our coin? Sim 1: HTTTTHTT Sim 2: TTHHTTTT Sim 3: TTTHHTHH Sim 4: TTTTTTTT Sim 5: TTTTTTTT Sim 6: THHTTTTH Sim 7: TTHHTTHT Sim 8: TTTTHHHT Sim 9: HTTHTTTT Sim 10: THHHTHTT Let p = 0.2
What is the probability of heads for our coin? Sim 1: HTTTTHTT Sim 2: TTHHTTTT Sim 3: TTTHHTHH Sim 4: TTTTTTTT Sim 5: TTTTTTTT Sim 6: THHTTTTH Sim 7: TTHHTTHT Sim 8: TTTTHHHT Sim 9: HTTHTTTT Sim 10: THHHTHTT Let p = 0.2 P(data) ≈ 1/10
What is the probability of heads for our coin? True Approximation 200 simulations per p
What is the probability of heads for our coin? True Approximation 2,000 simulations per p
What is the probability of heads for our coin? True Approximation 20,000 simulations per p
What is the probability of heads for our coin? True Approximation 200,000 simulations per p
Discrete time taxa[[i]]$nextstates= taxa[[i]]$states + intrinsicFn(…) + extrinsicFn(…)