Price equation 1076205 224517986 2008-07-09T05:44:30Z Lombar2 3540431 The '''Price equation''' (also known as '''Price's equation''') is a covariance equation which is a mathematical description of [[evolution]] and [[natural selection]]. The Price equation was derived by [[George R. Price]], working in London to rederive [[W.D. Hamilton]]'s work on [[kin selection]]. Suppose there is a population of <math>n</math> individuals over which the amount of a particular characteristic varies. Those <math>n</math> individuals can be grouped by the amount of the characteristic that each displays. In this case, at most there will be <math>n</math> groups of <math>n</math> distinct values of the characteristic, and at least there will be 1 group of a single shared value of the characteristic. Index each group with <math>i</math> so that the number of members in the group is <math>n_i</math> and the value of the characteristic shared among all members of the group is <math>z_i</math>. Now assume that having <math>z_i</math> of the characteristic is associated with having a [[Fitness (biology)|fitness]] <math>w_i</math> where the product <math>w_i n_i</math> represents the number of offspring in the next generation. Denote this number of offspring from group <math>i</math> by <math>n_i'</math> so that <math>w_i=n_i'/n_i</math>. Let <math>z_i'</math> be the amount of the characteristic displayed by the offspring from group <math>i</math>. Denote the amount of change in characteristic in group <math>i</math> by <math>\Delta z_i</math> defined by :<math> \Delta{z_i} \ \stackrel{\mathrm{def}}{=}\ z_i' - z_i </math> Now take <math>z</math> to be the [[average]] characteristic value in this population and <math>z'</math> to be the average characteristic value in the next generation. Define the change in average characteristic by <math>\Delta{z}</math>. That is, :<math> \Delta{z} \ \stackrel{\mathrm{def}}{=}\ z' - z </math> Note that this is ''not'' the average value of <math>\Delta{z_i}</math>. Also take <math>w</math> to be the average fitness of this population. The Price equation states: :<math> w\,\Delta{z}=\operatorname{cov}(w_i,z_i)+\operatorname{E}(w_i\,\Delta z_i) , \,\! </math> where the functions <math>\operatorname{E}</math> and <math>\operatorname{cov}</math> are respectively defined in Equations (1) and (2) below but are ''like'' the sample versions of the [[expected value]] and [[covariance]] operators from [[probability]] (see [[Sample mean and covariance]]). Note that this is really a [[difference equation]] relating the average value of a characteristic in one generation to the average value of the characteristic in the very next generation. In fact, assuming that <math>w</math> is not zero, it is often useful to write it as :<math> \Delta{z}=\frac{\operatorname{cov}(w_i,z_i)}{w}+\frac{\operatorname{E}(w_i\,\Delta z_i)}{w}\, </math> In the specific case that characteristic <math>z_i = w_i</math> (i.e., the fitness is the characteristic of interest), then Price's equation reformulates [[Fisher's fundamental theorem of natural selection]]. Price's equation is, importantly, a [[Tautology (logic)|tautology]]. It is a statement of mathematical fact between certain variables, and its value lies in the insight gained by assigning certain values encountered in evolutionary genetics to the variables. For example, the statement "if every pair of birds has two offspring, then among ten pairs of birds there will be twenty offspring" is a tautology. It doesn't really impart any new information about birds so much as it organizes our concepts about birds and their offspring. The Price equation is much more sophisticated than the above statement, but at its core, it too is a mathematically provable tautology. The Price equation also has applications in [[economics]]. ==Proof of the Price equation== To prove the Price equation, the following definitions are needed. If <math>n_i</math> is the number of occurrences of a pair of real numbers <math>x_i</math> and <math>y_i</math>, then: * The average or [[expected value]] of the <math>x_i</math> values is: {| border="0" cellpadding="0" cellspacing="0" width="100%" |- | align="left" | :<math> \operatorname{E}(x_i)\ \stackrel{\mathrm{def}}{=}\ \frac{\sum_i x_i n_i}{\sum_i n_i} </math> | align="right" | <math>(1)\,</math> |} * The [[covariance]] between the <math>x_i</math> and <math>y_i</math> values is: {| border="0" cellpadding="0" cellspacing="0" width="100%" |- | align="left" | :<math> \operatorname{cov}(x_i,y_i) \ \stackrel{\mathrm{def}}{=}\ \frac{\sum_i n_i~[x_i-\operatorname{E}(x_i)]~[y_i-\operatorname{E}(y_i)]}{\sum_i n_i} = \operatorname{E}(x_iy_i)-\operatorname{E}(x_i)\operatorname{E}(y_i) </math> | align="right" | <math>(2)\,</math> |} The notation <math>\langle x_i \rangle = \operatorname{E}( x_i )</math> will also be used when convenient. Suppose there is a population of organisms all of which have a genetic characteristic described by some real number. For example, high values of the number represent an increased visual acuity over some other organism with a lower value of the characteristic. Groups can be defined in the population which are characterized by having the same value of the characteristic. Let subscript <math>i</math> identify the group with characteristic <math>z_i</math> and let <math>n_i</math> be the number of organisms in that group. The total number of organisms is then <math>n</math> where: :<math> n = \sum_i n_i\, </math> The average value of the characteristic is <math>z</math> defined as: {| border="0" cellpadding="0" cellspacing="0" width="100%" |- | align="left" | :<math> z \ \stackrel{\mathrm{def}}{=}\ \operatorname{E}(z_i) = \frac{\sum_i z_i n_i}{n} </math> | align="right" | <math>(3)\,</math> |} Now suppose that the population reproduces, all parents are eliminated, and then there is a selection process on the children, by which less fit children are removed from the reproducing population. After reproduction and selection, the population numbers for the child groups will change to ''n''&prime;<sub>''i''</sub>. Primes will be used to denote child parameters, unprimed variables denote parent parameters. The total number of children is ''n' '' where: :<math> n' = \sum_i n'_i\, </math> The fitness of group ''i'' will be defined to be the ratio of children to parents: {| border="0" cellpadding="0" cellspacing="0" width="100%" |- | align="left" | :<math> w_i = \frac{n_i'}{n_i} </math> | align="right" | <math>(4)\,</math> |} with average fitness of the population being {| border="0" cellpadding="0" cellspacing="0" width="100%" |- | align="left" | :<math> w \ \stackrel{\mathrm{def}}{=}\ \operatorname{E}(w_i) = \frac{\sum_i w_i n_i}{n} = \frac{\sum_i \frac{n_i'}{n_i} n_i}{n} = \frac{\sum_i n_i'}{n} = \frac{n'}{n} </math> | align="right" | <math>(5)\,</math> |} The average value of the child characteristic will be ''z' '' where: {| border="0" cellpadding="0" cellspacing="0" width="100%" |- | align="left" | :<math> z' = \frac{\sum_i z'_i n_i'}{n'} </math> | align="right" | <math>(6)\,</math> |} where ''z''&prime;<sub>''i''</sub> are the (possibly new) values of the characteristic in the child population. Equation (2) shows that: {| border="0" cellpadding="0" cellspacing="0" width="100%" |- | align="left" | :<math> \operatorname{cov}(w_iz_i)=\operatorname{E}(w_iz_i)-wz </math> | align="right" | <math>(7)\,</math> |} Call the change in characteristic value from parent to child populations <math>\Delta z_i</math> so that <math>\Delta z_i = z'_i - z_i</math>. As seen in Equation (1), the expected value operator <math>\operatorname{E}</math> is [[linear]], so {| border="0" cellpadding="0" cellspacing="0" width="100%" |- | align="left" | :<math> \operatorname{E}(w_i\,\Delta z_i)=\operatorname{E}(w_iz'_i)-\operatorname{E}(w_iz_i) </math> | align="right" | <math>(8)\,</math> |} Combining Equations (7) and (8) leads to {| border="0" cellpadding="0" cellspacing="0" width="100%" |- | align="left" | :<math>\operatorname{cov}(w_iz_i)+\operatorname{E}(w_i\,\Delta z_i) = \bigl(\operatorname{E}(w_iz_i)-wz\bigr) + \bigl(\operatorname{E}(w_iz'_i)-\operatorname{E}(w_iz_i)\bigr) = \operatorname{E}(w_iz'_i)-wz</math> | align="right" | <math>(9)\,</math> |} but from Equation (1) gives: :<math>\operatorname{E}(w_iz'_i)=\frac{\sum_i w_iz'_in_i}{n}</math> and from Equation (4) gives: {| border="0" cellpadding="0" cellspacing="0" width="100%" |- | align="left" | :<math>\operatorname{E}(w_iz'_i)=\frac{\sum_i\frac{n'_i}{n_i}z'_in_i}{n}=\frac{\sum_in'_iz'_i}{n}=\frac{n'}{n}~\frac{\sum_i z'_in'_i}{n'}</math> | align="right" | <math>(10)\,</math> |} Applying Equations (5) and (6) to Equation (10) and then applying the result to Equation (9) gives the Price Equation: :<math>\operatorname{cov}(w_iz_i)+\operatorname{E}(w_i\,\Delta z_i)=wz'-wz=w\,\Delta z\,</math> ==Simple Price equation== When the characters ''z''<sub>''i''</sub> do not change from the parent to the child generation, the second term in the Price equation becomes zero resulting in a simplified version of the Price equation: :<math> w\,\Delta z = \operatorname{cov}\left(w_i,z_i\right) </math> which can be restated as: :<math> \Delta z = \operatorname{cov}\left(v_i,z_i\right) </math> where ''v''<sub>''i''</sub> is the fractional fitness: ''v''<sub>''i''</sub>= ''w''<sub>''i''</sub>/''w''. This simple Price equation can be proven using the definition in Equation (2) above. It makes this fundamental and tautological statement about evolution: "If a certain inheritable characteristic is correlated with an increase in fractional fitness, the average value of that characteristic in the child population will be increased over that in the parent population." ===Example: Evolution of sight=== As an example of the simple Price equation, consider a model for the evolution of sight. Suppose ''z''<sub>''i''</sub> is a real number measuring the visual acuity of an organism. An organism with a higher ''z''<sub>''i''</sub> will have better sight than one with a lower value of ''z''<sub>''i''</sub>. Let us say that the fitness of such an organism is ''w''<sub>''i''</sub>=''z''<sub>''i''</sub> 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''<sub>''i''</sub> = 3,2,1: :{| class="wikitable" |- |align="center"|'''''i''''' |align="center"|'''0''' |align="center"|'''1''' |align="center"|'''2''' |- |align="center"|'''''n''<sub>''i''</sub>''' |align="center"|10 |align="center"|20 |align="center"|30 |- |align="center"|'''''z''<sub>''i''</sub>''' |align="center"|3 |align="center"|2 |align="center"|1 |} Using Equation (4), the child population (assuming the character ''z''<sub>''i''</sub> doesn't change) :{| class="wikitable" |- |align="center"|'''''i''''' |align="center"|'''0''' |align="center"|'''1''' |align="center"|'''2''' |- |align="center"|'''''n''<sub>''i''</sub>''' |align="center"|30 |align="center"|40 |align="center"|30 |- |align="center"|'''''z''<sub>''i''</sub>''' |align="center"|3 |align="center"|2 |align="center"|1 |} 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: :<math> \Delta z \ \stackrel{\mathrm{def}}{=}\ z'-z = 1/3 \,\! </math> which indicates that the trait of sightedness is increasing in the population. Applying the Price equation we have (since ''z''&prime;<sub>''i''</sub>= ''z''<sub>''i''</sub>): :<math> \Delta z = \operatorname{cov}\left(w_i,z_i\right)/w = 1/3 </math> ===Dynamical sufficiency and the simple Price equation=== Sometimes the genetic model being used encodes enough information into the parameters used by the Price equation to allow the calculation of the parameters for all subsequent generations. This property is referred to as dynamical sufficiency. For simplicity, the following looks at dynamical sufficiency for the simple Price equation, but is also valid for the full Price equation. Referring to the definition in Equation (2), the simple Price equation for the character ''z'' can be written: :<math> w(z'-z)=\langle w_i z_i\rangle - wz </math> For the second generation: :<math> w'(z''-z')=\langle w'_i z'_i\rangle - w'z' </math> The simple Price equation for ''z'' only gives us the value of ''z'' '&nbsp; for the first generation, but does not give us the value of ''w'' '&nbsp; and 〈''w'' '<sub>i</sub> ''z'' '<sub>i</sub>&nbsp;〉 which needs to be calculated ''z'' ' '&nbsp; for the second generation. ''w'' '&nbsp; and 〈''w'' '<sub>i</sub> ''z'' '<sub>i</sub>&nbsp;〉 can both be thought of characteristics of the first generation, so the Price equation can be used to calculate them as well: :<math>w(w'-w)=\langle w_i^2\rangle - w^2</math> :<math> w(\langle w'_i z'_i\rangle-\langle w_i z_i\rangle)=\langle w_i ^2 z_i\rangle - w\langle w_i z_i\rangle </math> The five 0-generation variables ''w'', ''z'', 〈''w''<sub>''i''</sub>&nbsp; ''z''<sub>''i''</sub>&nbsp;〉, 〈''w''<sup>2</sup><sub>''i''</sub>&nbsp;〉, and 〈''w''<sup>2</sup><sub>''i''</sub>&nbsp;''z''<sub>''i''</sub>&nbsp;〉, which must be known before proceeding to calculate the three first generation variables ''w'' ', ''z'' ', and 〈''w'' '<sub>''i''</sub>&nbsp;''z'' '<sub>''i''</sub>&nbsp;〉, which are needed to calculate ''z'' '&nbsp; for the second generation. It can be seen that in general the Price equation cannot be used to propagate forward in time unless there is a way of calculating the higher moments (〈''w''<sup>''n''</sup><sub>''i''</sub>&nbsp;〉 and 〈''w''<sup>''n''</sup><sub>''i''</sub> ''z''<sub>''i''</sub>&nbsp;〉) from the lower moments in a way that is independent of the generation. Dynamical sufficiency means that such equations can be found in the genetic model, allowing the Price equation to be used alone as a propagator of the dynamics of the model forward in time. ===Example: Evolution of sickle cell anemia=== [[Image:autorecessive.svg|300px|thumb|An example of autosomal recessive inheritance. In the sickle cell case, the two parents are "carriers" who are resistant to malaria. Their children have one chance in four of inheriting both sickle cell genes and suffering sickle cell anemia, two chances in four of being a carrier themselves, and being resistant to malaria like their parents, and one chance in four of not inheriting the gene from either parent, and being susceptible to malaria.]] 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''<sub>''i''</sub>=''i'' be the number of sickle-cell genes that organisms of type ''i'' have so that ''z''<sub>''i''</sub>&nbsp;= 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''<sub>0</sub>''n''<sub>1</sub>/(''n''<sub>0</sub>+''n''<sub>1</sub>) and the number of ''i''–''i'' matings is ''n''<sup>2</sup><sub>''i''</sub>/[2(''n''<sub>0</sub>+''n''<sub>1</sub>)] where ''i''&nbsp;= 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''<sub>''i''</sub>=[''n''<sub>0</sub>,''n''<sub>1</sub>,''n''<sub>2</sub>]. If ''b'' is the birth rate, then after reproduction there will be :<math>b\left(\frac{n_0^2/2+n_0n_1/2+n_1^2/8}{n_0+n_1}\right)</math> type 0 children (unaffected) :<math>b\left(\frac{n_0n_1/2+n_1^2/4}{n_0+n_1} \right)</math> type 1 children (carriers) :<math>b\left(\frac{n_1^2/8}{n_0+n_1} \right)</math> 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: :<math>w_0=ab\left(\frac{n_0^2/2+n_0n_1/2+n_1^2/8}{n_0(n_0+n_1)}\right)</math> :<math>w_1=b\left(\frac{n_0n_1/2+n_1^2/4}{n_1(n_0+n_1)} \right)</math> :<math>w_2=0\,</math> To find the concentration ''n''<sub>1</sub> of carriers in the population at equilibrium, with the equilibrium condition of Δ ''z''=0, the simple Price equation is used: :<math> 0=\operatorname{cov}(w_i/w,z_i) = \frac{f(2-2a-af)}{(1+f)(2a+2f+af)} </math> where ''f''=''n''<sub>1</sub>/''n''<sub>0</sub>. At equilibrium then, assuming ''f'' is not zero: :<math> f=\frac{n_1}{n_0}=\frac{2(1-a)}{a} </math> 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: :<math> \langle w_i^2z_i \rangle = \frac{\langle w_i z_i \rangle}{z}</math> :<math> \langle w_i^2 \rangle = \frac{-b z^2 w^2 + 2 b z^2 w \langle w_i z_i \rangle + b z \langle w_i z_i \rangle^2 - b z^2 \langle w_i z_i \rangle^2 - 4 \langle w_i z_i \rangle^3} { b z^2 - 4 z \langle w_i z_i \rangle} </math> ===Example: sex ratios=== In a 2-sex species or deme with sexes 1 and 2 where <math> z_1= 1</math>, <math> z_2 = 0</math>, <math>z</math> 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 other sex's size. Consider proportionality constants <math> a</math> and <math> b</math> such that <math> w_1 = a(1 - z)</math> and <math> w_2 = bz</math>. This gives <math> w = (a + b)z(1 - z)</math> and <math> cov(w_i, z_i) + wz = E(w_i,z_i) = az(1 - z)</math>, so <math> cov(w_i, z_i) = az(1 - z) - (a + b)z^2(1 - z) = z(1 - z)(a - (a + b)z) </math>. Hence, <math> \Delta z = a/(a + b) - z </math> so that <math> z' = a/(a + b) </math>. ==Full Price equation== The simple Price equation was based on the assumption that the characters ''z''<sub>''i''</sub> do not change over one generation. If it is assumed that they do change, with ''z''<sub>''i''</sub> being the value of the character in the child population, then the full Price equation must be used. A change in character can come about in a number of ways. The following two examples illustrate two such possibilities, each of which introduces new insight into the Price equation. ===Examples=== ====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''<sub>''i''</sub> by a model equation: :<math> w_i = \frac{n'_i}{n_i} = k - a z_i + b z </math> where ''z''<sub>''i''</sub> is a measure of altruism, the ''az''<sub>''i''</sub> 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: :<math> w\,\Delta z = -a~\operatorname{var}\left(z_i\right) </math> where var(''z''<sub>''i''</sub>) is the [[variance]] of ''z''<sub>''i''</sub> which is just the covariance of ''z''<sub>''i''</sub> with itself: :<math> \operatorname{var}(z_i) \ \stackrel{\mathrm{def}}{=}\ \operatorname{E}(z_i^2)-\operatorname{E}(z_i)^2 </math> 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, 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''<sub>''ij''</sub> will specify the number of individuals of type ''ij''. Let ''z''<sub>''ij''</sub> be the degree of altruism expressed by individual ''j'' of group ''i'' towards the members of group ''i''. Let's specify the fitness ''w''<sub>''ij''</sub> by a model equation: :<math> w_{ij} = \frac{n'_{ij}}{n_{ij}} = k - a z_{ij} + b z_i </math> The ''a z''<sub>''ij''</sub> term is the fitness the organism loses by being altruistic and is proportional to the degree of altruism ''z''<sub>''ij''</sub> that it expresses towards members of its own group. The ''b z''<sub>''i''</sub> 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''<sub>''i''</sub> expressed by the group towards its members. Again, in studying study 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''<sub>''ij''</sub> >''bz''<sub>''i''</sub>. Defining the group averages: :<math>n_i = \sum_j n_{ij}\,</math> :<math>z_i = \frac{\sum_j z_{ij}n_{ij}}{n_i}</math> :<math>w_i = \frac{\sum_j w_{ij}n_{ij}}{n_i}=k+(b-a)z_i</math> :<math>n_i'= \sum_j n_{ij}'=n_i(k+(b-a)z_i)\,</math> :<math>z_i'= \frac{\sum_j z_{ij}n_{ij}'}{n_i'}</math> and global averages: :<math>n = \sum_{ij} n_{ij} = \sum_i n_i\,</math> :<math>z = \frac{\sum_{ij} z_{ij}n_{ij}}{n} = \frac{\sum_i z_in_i}{n}</math> :<math>w = \frac{\sum_{ij} w_{ij}n_{ij}}{n} = \frac{\sum_i w_in_i}{n}</math> :<math>n'= \sum_{ij} n_{ij}' = \sum_i n_i'\,</math> :<math>z'= \frac{\sum_{ij} z_{ij}n_{ij}'}{n'} = \frac{\sum_i z_i'n_i'}{n'}</math> It can be seen that since the ''z''<sub>''i''</sub> and ''z''<sub>''i''</sub> are now averages over a particular group, and since these groups are subject to selection, the value of Δ''z''<sub>''i''</sub>&nbsp;= ''z''&prime;<sub>''i''</sub>−''z''<sub>''i''</sub> will not necessarily be zero, and the full Price equation will be needed. :<math> \Delta z = \operatorname{cov}(w_i/w,z_i)+\operatorname{E}(w_i\,\Delta z_i/w)\, </math> 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: :<math> \operatorname{cov}(w_i/w,z_i)=\left(b-a\right)\operatorname{var}(z_i) </math> 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. ====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''<sub>0</sub>, ''A''<sub>''m''</sub>, ''B''<sub>0</sub>, and ''B''<sub>''m''</sub>. 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''<sub>0</sub> and ''A''<sub>''m''</sub> are both ''A''-alleles, but organisms with the ''A''<sub>0</sub> gene produce offspring with ''A''<sub>0</sub>-genes only, while organisms with the ''A''<sub>''m''</sub> gene produce (1−3''m'') offspring with the ''A''<sub>''m''</sub> gene, and ''m'' organisms of the remaining three gene types. Likewise, ''B''<sub>0</sub> and ''B''<sub>''m''</sub> are both ''B''-alleles, but organisms with the ''B''<sub>0</sub> gene produce offspring with ''B''<sub>0</sub>-genes only, while organisms with the ''B''<sub>''m''</sub> gene produce (1−3''m'') offspring with the ''B''<sub>''m''</sub> gene, and ''m'' organisms of the remaining three gene types. Let ''i''=0,1,2,3 be the indices associated with the ''A''<sub>0</sub>, ''A''<sub>''m''</sub>, ''B''<sub>0</sub>, and ''B''<sub>''m''</sub> genes respectively. Let ''w''<sub>''ij''</sub> be the number of viable type-''j'' organisms produced per type-''i'' organism. The ''w''<sub>''ij''</sub> matrix is: (with ''i'' denoting rows and ''j'' denoting columns) :{| class="wikitable" |- |align="center"|α |align="center"|0 |align="center"|0 |align="center"|0 |- |align="center"|mα |align="center"|(1−3''m'')α |align="center"|''m''β |align="center"|''m''β |- |align="center"|0 |align="center"|0 |align="center"|β |align="center"|0 |- |align="center"|''m''α |align="center"|''m''α |align="center"|''m''β |align="center"|(1−3''m'')β |} :<math>\mathbf{w}= \begin{bmatrix} \alpha & 0 &0 & 0 \\ m\alpha & (1-3m)\alpha & m\beta & m\beta \\ 0 & 0 & \beta & 0 \\ m\alpha & m\alpha & m\beta & (1-3m)\beta \end{bmatrix} </math> 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 acid]]s 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: :<math> z_i = \left[0,1,0,1\right] </math> in other words, 0 for non-mutator genes, 1 for mutator genes. There are two cases: =====Genotype fitness===== Lets focus on the idea of the fitness of the genotype. The index ''i'' indicates the genotype and the number of type ''i'' genotypes in the child population is: :<math> n'_i = \sum_j w_{ji}n_j\, </math> which gives fitness: :<math> w_i=\frac{n'_i}{n_i} </math> Since the individual mutability ''z''<sub>''i''</sub> does not change, the average mutabilities will be: :<math>z = \frac{\sum_i z_i n_i}{n}</math> :<math>z' = \frac{\sum_i z_i n'_i}{n'}</math> with these definitions, the simple Price equation now applies. =====Lineage fitness===== In this case we want to look at the idea that fitness is measured by the number of children an organism has, regardless of their genotype. Note that we now have two methods of grouping, by lineage, and by genotype. It is this complication that will introduce the need for the full Price equation. The number of children an ''i''-type organism has is: :<math> n'_i = n_i\sum_j w_{ij}\, </math> which gives fitness: :<math> w_i=\frac{n'_i}{n_i} = \sum_j w_{ij} </math> We now have characters in the child population which are the average character of the ''i''-th parent. :<math> z'_j = \frac{\sum_i n_i z_i w_{ij} }{\sum_i n_i w_{ij}} </math> with global characters: :<math>z = \frac{\sum_i z_i n_i}{n}</math> :<math>z' = \frac{\sum_i z_i n'_i}{n'}</math> with these definitions, the full Price equation now applies. ==Cultural References== Price's equation features in the plot and title of the 2008 thriller film ''[[WΔZ]]''. ==References== * {{cite journal | last=Frank | first=S.A. | title=George Price's contributions to Evolutionary Genetics | journal=Journal of Theoretical Biology | year=1995 | volume=175 | pages=373–388 | url=http://stevefrank.org/reprints-pdf/95JTB-Price.pdf | doi=10.1006/jtbi.1995.0148}} * {{cite journal | last=Frank | first=S.A. | title=The Price Equation, Fisher's Fundamental Theorem, Kin Selection, and Causal Analysis | journal=Evolution | year=1997 | volume=51 | issue=6 | pages=1712–1729 | url=http://stevefrank.org/reprints-pdf/97Evol-Causal.pdf | doi=10.2307/2410995}} * {{cite journal | last=Gardner | first=A. | title=The Price equation | journal=Curr. Biol. | year=2008 | volume=18 | pages=R198–R202 | url=http://www.biology.ed.ac.uk/research/groups/gardner/publications/Gardner_2008.pdf | doi=10.1016/j.cub.2008.01.005}} * {{cite journal | last=Grafen | first=A. | title=Developments of the Price equation and natural selection under uncertainty | journal=Proc. R. Soc. London B | year=2000 | volume=267 | pages=1223–1227 | url=http://users.ox.ac.uk/~grafen/cv/PriceUnc.pdf | doi=10.1098/rspb.2000.1131}} * {{cite journal | author=Langdon, W. B. | title=Evolution of GP Populations: Price's Selection and Covariance Theorem | journal=Genetic Programming and Data Structures | year=1998 | chapter=8.1 | pages=167–208 }} * {{cite journal | author=Price, G.R. | authorlink = George R. Price | title=Selection and covariance | journal=[[Nature (journal)|Nature]] | year=1970 | volume=227 | pages=520–521 | doi=10.1038/227520a0}} * {{cite journal | author=Price, G.R. | authorlink = George R. Price | title=Fisher's "fundamental theorem" made clear | journal=Annals of Human Genetics | year=1972 | volume=36 | pages=129–140 | doi=10.1111/j.1469-1809.1972.tb00764.x }} * {{cite journal | author=Van Veelen, M. | title=On the use of the Price equation | journal=Journal of Theoretical Biology | year=2005 | volume=237 | pages=412–426 | doi=10.1016/j.jtbi.2005.04.026}} * {{cite journal | last=Day | first=T. | title=Insights from Price's equation into evolutionnary epidemiology | journal=DIMACS Series in Discrete Mathematics and Theoretical Computer Science | year=2006 | volume=71 | pages=23–43}} {{popgen}} [[Category:Equations]] [[Category:Evolutionary dynamics]] [[Category:Evolutionary biology]] [[Category:Population genetics]] [[de:Price-Gleichung]]