Personal tools
- Log in
Navigation
- Main Page
- Using the software
Tutorials
Coming soon
- Adaptive Dynamics and Population Structure

Price 3.0

From Evolution and Games

(Difference between revisions)

Revision as of 22:39, 18 September 2010

Go to the front-page of Price equation tutorials

1 A model with frequency dependence
2 An A mutant in a population of B's
3 A B mutant in a population of A's
4 Comparing the two covariances
5 The Price equation for this very simple model
6 What would a proper probability theorist do?
7 What would a proper statistician do?
8 So how does the Price equation help here?

A model with frequency dependence

We go back to the asexual model, but now we will take a model with frequency dependence. This will be a bit more complicated, but we can still use a basic, well known model. Suppose there are $N$ individuals. They are paired randomly, and play a game in those pairs. The game gives them payoffs according to the following payoff matrix

$\begin{array}{ccc} & A & B \\ A & 1 & 0 \\ B & 0 & 1 \end{array}$

If $q i = 1$ , then that means that individual $i$ plays strategy $A$ , and $q i = 0$ means that individual $i$ plays strategy $B$ . Depending on the match, each individual gets a payoff $π i$ . The next generation is drawn, as before, one individual at a time. But now the probability that $i$ is drawn as a parent depends on the payoff, and not just on the own genotype.

$\mathbb{P}\left( i\text{ is chosen}\right) =\frac{\alpha \pi _{i}+\beta }{\sum_{j=1}^{N}\left( \alpha \pi _{j}+\beta \right) }$

We do the same procedure as before, but with a few different starting points.

An A mutant in a population of B's

Suppose $q 1 = 1$ and $q i = 0$ for $i = 2,..., N$ . The probability of any individual $i = 2,..., N$ to be matched with individual $1$ is $1 / (N - 1)$ .

The probability with which an individual $i$ with $q i = 1$ , who is matched to an individual $j$ with $q j = 0$ , is chosen is:

$p_{1,0,1}=\frac{\beta }{\left( N-2\right) \alpha +N\beta }$

The probability with which an individual $i$ with $q i = 0$ , who is matched to an individual $j$ with $q j = 0$ , is chosen is:

$p_{0,0,1}=\frac{\alpha +\beta }{\left( N-2\right) \alpha +N\beta }$

The probability with which an individual $i$ with $q i = 0$ , who is matched to an individual $j$ with $q j = 1$ , is chosen is:

$p_{0,1,1}=\frac{\beta }{\left( N-2\right) \alpha +N\beta }$

We can also compute the covariance of these two actual random variables. Let us start by computing expectations...

$\mathbb{E}\left[ X\right] =\frac{1}{N}$

[show details]

For symmetry reasons...

$\mathbb{E}\left[ Y\right] =1$

$\mathbb{E}\left[ XY\right] =\frac{\beta }{\left( N-2\right) \alpha +N\beta }$

[show details]

The Covariance is thus...

$Cov\left( X,Y\right) =\mathbb{E}\left[ XY\right] -\mathbb{E}\left[ X\right] \mathbb{E}\left[ Y\right] =\frac{\beta }{\left( N-2\right) \alpha +N\beta }-\frac{1}{N}$

A B mutant in a population of A's

Suppose $q 1 = 0$ and $q i = 1$ for $i = 2,..., N$ . The probability of any individual $i = 2,..., N$ to be matched with individual $1$ is $1 / (N - 1)$ .

The probability with which an individual $i$ with $q i = 0$ , who is matched to an individual $j$ with $q j = 1$ ...

$p_{0,1,N-1}=\frac{\beta }{\left( N-2\right) \alpha +N\beta }$

The probability with which an individual $i$ with $q i = 1$ , who is matched to an individual $j$ with $q j = 1$ ...

$p_{1,1,N-1}=\frac{\alpha +\beta }{\left( N-2\right) \alpha +N\beta }$

The probability with which an individual $i$ with $q i = 1$ , who is matched to an individual $j$ with $q j = 0$ ...

$p_{1,0,N-1}=\frac{\beta }{\left( N-2\right) \alpha +N\beta }$

So...

$\mathbb{P}\left( Y=k|X=0\right) =\binom{N}{k}\left( p_{0,1,N-1}\right)^{k}\left( 1-p_{0,1,N-1}\right) ^{N-k}=$

$\mathbb{P}\left( Y=k|X=1\right) =\frac{1}{N-1}\binom{N}{k}\left(p_{1,0,N-1}\right) ^{k}\left( 1-p_{1,0,N-1}\right) ^{N-k}+\frac{N-2}{N-1}\binom{N}{k}\left( p_{1,1,N-1}\right) ^{k}\left( 1-p_{1,1,N-1}\right) ^{N-k}$

This is pretty well defined, so we can also compute the covariance of these two actual random variables. Let us start with the expectations... $\mathbb{E}\left[ X\right] =\frac{N-1}{N}$

[show details]

For symmetry reasons...

$\mathbb{E}\left[ Y\right] =1$

$\mathbb{E}\left[ XY\right] = \frac{\left( N-2\right) \alpha +\left( N-1\right) \beta }{\left(N-2\right) \alpha +N\beta }$

[show details]

Thus....

$Cov\left( X,Y\right) =\mathbb{E}\left[ XY\right] -\mathbb{E}\left[ X\right] \mathbb{E}\left[ Y\right] =\frac{\left( N-2\right) \alpha +\left( N-1\right)\beta }{\left( N-2\right) \alpha +N\beta }-\frac{N-1}{N}$

Comparing the two covariances

If we compare these two covariances, for instance for $β = 0$ and $α = 1$ , then in the first case $C o v (X, Y) = - 1 / N$ , while in the second case $C o v (X, Y) = + 1 / N$ . This should not feel strange; it is a coordination game, and if everyone plays $B$ , then playing $A$ is a bad idea, but if everyone else plays $A$ too, then playing $A$ is a good idea. The model covariance is therefore not constant, and changes with the frequency of strategy $A$ .

The Price equation for this very simple model

Here we have a nested sampling procedure; you can choose a starting point (by clicking on each individual in the parent population), keep it fixed, and draw a few transitions. Then change your starting point, and draw a few transitions from there. You will notice that the Price equation varies with every transition, along with the likelyhood of the different transitions that depend on the initial population.

What would a proper probability theorist do?

Try to figure out the transition probabilities for the suggested model as a function of the current state. It seems that Nowak, Liebermann, Ohtsuki and Tarnita are a proper probability theorists, because that is what they did already.

What would a proper statistician do?

Sample often for different starting points. Fit data to the model, estimate $α$ and $β$ . Test if the payoff matrix is as assumed.

So how does the Price equation help here?

Not. Can't think of anything remotely useful about the Price equation. The extra problem here is the frequency dependence. Sometimes this is referred to as the problem that the Price Equation is not dynamically sufficient. That is a step in the right direction, but still inaccurate, because it only describes one of the symptoms, and not the cause. It is the model we have here that gives us, after a bit of computing, the transition probabilities for all frequencies of $A$ , and thereby all properties of the dynamical system. The models in versions 1.0 and 2.0 just concerned transitions leaving from one given starting point. A model can therefore be described as dynamically sufficient or insufficient. Because the Price equation can only deal with one step at at time, it tends to inspire people to think of implicit assumptions that pertain to dynamically insufficient models. That, however, does not mean that the Price Equation cannot be used in combination with a dynamically sufficient model. It can; as shown in the applet above. But the question is: used for what? And the answer is: for nothing.

Retrieved from "http://evolutionandgames.nl/wiki/index.php?title=Price_3.0"

Category: Price equation

@@ Line 1: / Line 1: @@
 [[Category:Price equation]]
+[[Price | Go to the front-page of Price equation tutorials ]]
 == A model with frequency dependence ==

Price 3.0

From Evolution and Games

Revision as of 22:39, 18 September 2010

Contents

A model with frequency dependence

An A mutant in a population of B's

A B mutant in a population of A's

Comparing the two covariances

The Price equation for this very simple model

What would a proper probability theorist do?

What would a proper statistician do?

So how does the Price equation help here?