Price equation examples

In the theory of evolution and natural selection, the Price equation expresses certain relationships between a number of statistical measures on parent and child populations. Without training in the understanding of these measures, the meaning of Price equation is rather opaque. For the less experienced person, simple and particular examples are vital to gaining an intuitive understanding of these statistical measures as they apply to populations, and the relationship between them as expressed in the Price equation.

Evolution of sight

As an example of the simple Price equation, consider a model for the evolution of sight. Suppose ''z''_''i'' is a real number measuring the visual acuity of an organism. An organism with a higher ''z''_''i'' will have better sight than one with a lower value of ''z''_''i''. Let us say that the fitness of such an organism is ''W''_''i''=''z''_''i'' which means the more sighted it is, the fitter it is, that is, the more children it will produce. Beginning with the following description of a parent population composed of 3 types: (''i'' = 0,1,2) with sightedness values ''z''_''i'' = 3,2,1:
:

Using Equation  (4), the child population $n_i'$ (assuming the character ''z''_''i'' doesn't change; that is, ''z''_''i'' = ''z''_''i''')
:

We would like to know how much average visual acuity has increased or decreased in the population. From Equation  (3), the average sightedness of the parent population is ''z'' = 5/3. The average sightedness of the child population is ''z = 2, so that the change in average sightedness is:

:$\Delta z \;\stackrel\; z' - z = \frac$

which indicates that the trait of sightedness is increasing in the population. (Note that the covariance formula used below is not the standard covariance formula commonly used in mainstream textbooks. Refer to Equation (2) for Price's definition of covariance in this context.) Applying the Price equation we have (since ''z''′_''i''= ''z''_''i''):

:$\Delta z = \frac\operatorname\left(w_i, z_i\right) = \frac$

Evolution of sickle cell anemia

As an example of dynamical sufficiency, consider the case of sickle cell anemia. Each person has two sets of genes, one set inherited from the father, one from the mother. Sickle cell anemia is a blood disorder which occurs when a particular pair of genes both carry the 'sickle-cell trait'. The reason that the sickle-cell gene has not been eliminated from the human population by selection is because when there is only one of the pair of genes carrying the sickle-cell trait, that individual (a "carrier") is highly resistant to malaria, while a person who has neither gene carrying the sickle-cell trait will be susceptible to malaria. Let's see what the Price equation has to say about this.

Let ''z''_''i''=''i'' be the number of sickle-cell genes that organisms of type ''i'' have so that ''z''_''i'' = 0 or 1 or 2. Assume the population sexually reproduces and matings are random between type 0 and 1, so that the number of 0–1 matings is ''n''₀''n''₁/(''n''₀+''n''₁) and the number of ''i''–''i'' matings is ''n''²_''i''/ (''n''₀+''n''₁)where ''i'' = 0 or 1. Suppose also that each gene has a 1/2 chance of being passed to any given child and that the initial population is ''n''_''i''= 'n''₀,''n''₁,''n''₂ If ''b'' is the birth rate, then after reproduction there will be

:$b\left(\frac\right)$ type 0 children (unaffected)
:$b\left(\frac\right)$ type 1 children (carriers)
:$b\left(\frac\right)$ type 2 children (affected)

Suppose a fraction ''a'' of type 0 reproduce, the loss being due to malaria. Suppose all of type 1 reproduce, since they are resistant to malaria, while none of the type 2 reproduce, since they all have sickle-cell anemia. The fitness coefficients are then:

:$\begin w_0 &= ab\left(\frac\right) \\ w_1 &= b\left(\frac \right) \\ w_2 &= 0 \end$

To find the concentration ''n''₁ of carriers in the population at equilibrium, with the equilibrium condition of Δ ''z''=0, the simple Price equation is used:

:$0 = \operatorname\left(\frac, z_i\right) = \frac$

where ''f''=''n''₁/''n''₀. At equilibrium then, assuming ''f'' is
not zero:

:$f = \frac = \frac$

In other words, the ratio of carriers to non-carriers will be equal to the above constant non-zero value. In the absence of malaria, ''a''=1 and ''f''=0 so that the sickle-cell gene is eliminated from the gene pool. For any presence of malaria, ''a'' will be smaller than unity and the sickle-cell gene will persist.

The situation has been effectively determined for the infinite (equilibrium) generation. This means that there is dynamical sufficiency with respect to the Price equation, and that there is an equation relating higher moments to lower moments. For example, for the second moments:

:$\begin \langle w_i^2z_i \rangle &= \frac \\ \langle w_i^2 \rangle &= \frac \end$

Sex ratios

In a 2-sex species or deme with sexes 1 and 2 where $z_1= 1$, $z_2 = 0$, $z$ is the relative frequency of sex 1. Since all individuals have one parent of each sex, the fitness of each sex is proportional to the number of the other sex. Consider proportionality constants $a$ and $b$ such that $w_1 = a(1 - z)$and $w_2 = bz$. Under this scenario, ''a'' is the number of children a male would have if there he were the only male and unlimited number of females, while ''b'' is the number of children a female would have if she were the only female and unlimited number of males. This gives $w = (a + b)z(1 - z)$ and $cov(w_i, z_i) + wz = E(w_i z_i) = az(1 - z)$, so cov(w_i, z_i) = az(1 - z) - (a + b)z^2(1 - z) = z - za - (a + b)z]. Hence, $\Delta z = \frac - z$ so that $z' = \frac$.

Under another scenario, every woman has a maximum number of children ($w_0=b$) so that $b\, n_f$ children are created per generation, and every male is responsible for an equal number of children so that $w_1=b\, n_f/n_m$ where $n_f$ and $n_m$ are the total number of females and males respectively. In this case, the sex ratio stabilizes at $z=1/2$.

Evolution of mutability

Suppose there is an environment containing two kinds of food. Let α be the amount of the first kind of food and β be the amount of the second kind. Suppose an organism has a single allele which allows it to utilize a particular food. The allele has four gene forms: ''A''₀, ''A''_''m'', ''B''₀, and ''B''_''m''. If an organism's single food gene is of the ''A'' type, then the organism can utilize ''A''-food only, and its survival is proportional to α. Likewise, if an organism's single food gene is of the ''B'' type, then the organism can utilize ''B''-food only, and its survival is proportional to β. ''A''₀ and ''A''_''m'' are both ''A''-alleles, but organisms with the ''A''₀ gene produce offspring with ''A''₀-genes
only, while organisms with the ''A''_''m'' gene produce, among their ''n'' offspring, (''n''−3''m'') offspring with the ''A''_''m'' gene, and ''m'' organisms of the remaining three gene types. Likewise, ''B''₀ and ''B''_''m'' are both ''B''-alleles, but organisms with the ''B''₀ gene produce offspring with ''B''₀-genes only, while
organisms with the ''B''_''m'' gene produce (''n''−3''m'') offspring with the ''B''_''m'' gene, and ''m'' organisms of the remaining three gene types.

Let ''i''=0,1,2,3 be the indices associated with the ''A''₀, ''A''_''m'', ''B''₀, and ''B''_''m'' genes respectively. Let ''w''_''ij'' be the number of viable type-''j'' organisms produced per type-''i'' organism. The ''w''_''ij'' matrix is: (with ''i'' denoting rows and ''j'' denoting columns)

:$\mathbf = \begin \alpha & 0 & 0 & 0 \\ m\alpha & (n - 3m)\alpha & m\beta & m\beta \\ 0 & 0 & \beta & 0 \\ m\alpha & m\alpha & m\beta & (n - 3m)\beta \end$

Mutators are at a disadvantage when the food supplies α and β are constant. They lose every generation compared to the non-mutating genes. But when the food supply varies, even though the mutators lose relative to an ''A'' or ''B'' non-mutator, they may lose less than them over the long run because, for example, an ''A'' type loses a lot when α is low. In this way, "purposeful" mutation may be selected for. This may explain the redundancy in the genetic code, in which some amino acids are encoded by more than one codon in the DNA. Although the codons produce the same amino acids, they have an effect on the mutability of the DNA, which may be selected for or against under certain conditions.

With the introduction of mutability, the question of identity versus lineage arises. Is fitness measured by the number of children an individual has, regardless of the children's genetic makeup, or is fitness the child/parent ratio of a particular genotype?. Fitness is itself a characteristic, and as a result, the Price equation will handle both.

Suppose we want to examine the evolution of mutator genes. Define the ''z''-score as:

:z_i = \left , 1, 0, 1\right/math>

in other words, 0 for non-mutator genes, 1 for mutator genes. There are
two cases:

Evolution of altruism

To study the evolution of a genetic predisposition to altruism, altruism will be defined as the genetic predisposition to behavior which decreases individual fitness while increasing the average fitness of the group to which the individual belongs. First specifying a simple model, which will only require the simple Price equation. Specify a fitness ''w''_''i'' by a model equation:

:$w_i = \frac = k - a z_i + b z$

where ''z''_''i'' is a measure of altruism, the ''az''_''i'' term is the decrease in fitness of an individual due to altruism towards the group and ''bz'' is the increase in fitness of an individual due to the altruism of the group towards an individual.
Assume that ''a'' and ''b'' are both greater than zero. From the Price equation:

:$w\Delta z = -a\operatorname\left(z_i\right)$

where var(''z''_''i'') is the variance of ''z''_''i'' which is just the covariance of ''z''_''i'' with itself:

:$\operatorname(z_i) \;\stackrel\; \operatorname(z_i^2) - \operatorname(z_i)^2$

It can be seen that, by this model, in order for altruism to persist it must be uniform throughout the group. If there are two altruist types the average altruism of the group will decrease, and the more altruistic will lose out to the less altruistic.

Now assuming a hierarchy of groups which will require the full Price equation. The population will be divided into groups, labelled with index ''i'' and then each group will have a set of subgroups labelled by index ''j''. Individuals will thus be identified by two indices,
''i'' and ''j'', specifying which group and subgroup they belong to. ''n''_''ij'' will specify the number of individuals of type ''ij''. Let ''z''_''ij'' be the degree of altruism expressed by individual ''j'' of group ''i'' towards the members of group ''i''. Let's specify the fitness ''w''_''ij'' by a model equation:

:$w_ = \frac = k - a z_ + b z_i$

The ''a z''_''ij'' term is the fitness the organism loses by being altruistic and is
proportional to the degree of altruism ''z''_''ij'' that it expresses towards members
of its own group. The ''b z''_''i'' term is the fitness that the organism gains from the altruism of the members of its group, and is proportional to the average altruism ''z''_''i'' expressed by the group towards its members. Again, in studying altruistic (rather than spiteful) behavior, it is expected that ''a'' and ''b'' are positive numbers. Note that the above behavior is altruistic only when ''az''_''ij'' >''bz''_''i''. Defining the group averages:

:\begin
n_i &= \sum_j n_ \\
z_i &= \frac\sum_j z_n_ \\
w_i &= \frac\sum_j w_n_ = k + (b - a)z_i \\
n_i'&= \sum_j n_' = n_i + (b - a)z_i\\
z_i'&= \frac\sum_j z_n_'
\end

and global averages:

:$\begin n &= \sum_ n_ = \sum_i n_i \\ z &= \frac\sum_ z_n_ = \frac\sum_i z_in_i \\ w &= \frac\sum_ w_n_ = \frac\sum_i w_in_i \\ n'&= \sum_ n_' = \sum_i n_i' \\ z'&= \frac\sum_ z_n_' = \frac\sum_i z_i'n_i' \end$

It can be seen that since the ''z''_''i'' and ''z''_''i'' are now averages over a particular group, and since these groups are subject to selection, the value of Δ''z''_''i'' = ''z''′_''i''−''z''_''i'' will not necessarily be zero, and the full Price equation will be needed.

:$\Delta z = \operatorname\left(\frac, z_i\right) + \operatorname\left(w_i\,\Delta \frac\right)$

In this case, the first term isolates the advantage to each group conferred by having altruistic members. The second term isolates the loss of altruistic members from their group due to their altruistic behavior. The second term will be negative. In other words, there will be an average loss of altruism due to the in-group loss of altruists, assuming that the altruism is not uniform across the group. The first term is:

:$\operatorname\left(\frac, z_i\right) = \left(b - a\right)\operatorname{var}(z_i)$

In other words, for ''b''>''a'' there may be a positive contribution to the average altruism as a result of a group growing due to its high number of altruists and this growth can offset in-group losses, especially if the variance of the in-group altruism is low. In order for this effect to be significant, there must be a spread in the average altruism of the groups.

Equations
Evolutionary dynamics
Evolutionary biology
Population genetics