Skip to content

Derivation of the continuous-time Price equation

June 26, 2012

The standard discrete-time Price equation is really useful, but its derivation is much less beautiful than the continuous-time version in my opinion. The beautiful thing about the continuous-time version is that it can be used to illustrate the clear relationship between the Price equation and the chain rule.

Model and notation

Consider a set of groups of individual organisms and label the groups with i = 1, 2, .... Groups are characterized by a particular trait; every individual belonging to group i is characterized by the same trait-value, x_i. The number, n_i, of individuals in group i changes with time according to the equation,

\frac{\mathrm{d} n_i}{\mathrm{d} t} = r_i n_i

where r_i is the per capita growth rate of the ith group, which may itself change over time.

Objective

Given this model, our goal is to derive an equation for the time dynamics of the weighted average of the trait, which is defined as,

\mathrm{E}(x) = \sum_i p_i x_i

where,

p_i = \frac{n_i}{\sum_j n_j}

is the relative abundance of the ith group.

The chain rule

As I mentioned earlier, the continuous-time Price equation is essentially a special case of the chain rule. In particular, by the chain rule we can write the time-derivative of the weighted average trait value as,

\frac{\mathrm{dE}(x)}{\mathrm{d}t} = \sum_i \frac{\partial \mathrm{E}(x)}{\partial p_i} \frac{\partial p_i}{\partial t} +  \sum_i \frac{\partial \mathrm{E}(x)}{\partial x_i} \frac{\partial x_i}{\partial t}

Substituting the partial derivatives of the weighted average trait value,

\frac{\mathrm{dE}(x)}{\mathrm{d}t} = \sum_i x_i \frac{\partial p_i}{\partial t} +  \sum_i p_i \frac{\partial x_i}{\partial t}

or with dots,

\frac{\mathrm{dE}(x)}{\mathrm{d}t} = \sum_i x_i \dot{p}_i +  \sum_i p_i \dot{x}_i

This is the most important equation in the post. Essentially its the continuous-time Price equation, but in a non-standard form. The first summation gives the selection effect: it will be large if groups with high trait values tend to be increasing in relative abundance. The second summation gives the trait dynamics effect: it will be large if abundant groups tend to be increasing in their trait values. What I like about this form of the Price equation is that it demonstrates a deep connection between evolutionary biology and many other general equations in science that are also based on the chain rule (e.g. Hamilton’s equations). The rest of the post will be about how to get from this ‘chain-rule form’ to the standard ‘moment form’.

Aside: relative abundance dynamics

Before getting to the standard form, we need another result. The chain-rule shows up again,

\frac{\partial p_i}{\partial t} = \sum_j \frac{\partial p_i}{\partial n_j} \frac{\partial n_j}{\partial t}

The first partial in each term comes from differentiating the definition for the relative abundances (we need the quotient rule here),

\frac{\partial p_i}{\partial n_j} = \begin{cases} \frac{-p_i}{N}, & i \neq j \\ \frac{1-p_i}{N}, & i = j \end{cases}

where N = \sum_i n_i is the total abundance. After (as people tend to say) `some’ algebra,

\frac{\partial p_i}{\partial t} = p_i (r_i - \sum_j p_j r_j)

which implies, given our definition of weighted averaging (i.e. mathematical expectation),

\frac{\partial p_i}{\partial t} = p_i (r_i - \mathrm{E}(r))

This says that groups with per capita growth rates greater than the average will increase in relative abundance, which makes sense.

Selection and covariance

This is the title of Price’s original article. It refers to the mathematical equivalence between the concept of selection from evolutionary theory and the concept of covariance from probability theory. The original work by Price was in discrete-time. Here I’m interested in continuous-time and the chain rule. Using the time-derivative of the relative abundance derived above, we find that the first term of the Price equation can be expressed as,

\sum_i x_i \frac{\partial p_i}{\partial t} = \sum_i x_i p_i (r_i - \mathrm{E}(r))

Again using weighted averaging,

\sum_i x_i \frac{\partial p_i}{\partial t} = \mathrm{E} (x_i (r_i - \mathrm{E}(r)))

From the definition of covariance,

\sum_i x_i \frac{\partial p_i}{\partial t} = \mathrm{Cov}(x, r)

Average trait-value dynamics

Now for the second term of the continuous-time Price equation. Basically, its just the weighted average trait dynamics,

\sum_i p_i \frac{\partial x_i}{\partial t} = \mathrm{E}(\dot{x})

where \dot{x_i} = \frac{\partial x_i}{\partial t}.

Standard form of the continuous-time Price equation

Putting all this together gives,

\frac{\mathrm{dE}(x)}{\mathrm{d}t} = \mathrm{Cov}(x, r) + \mathrm{E}(\dot{x})

So…the rate of trait evolution depends on (1) the covariance between the trait and per capita growth rate and (2) the expected rate of change of the trait values. This is a beautiful equation but the chain-rule form is also beautiful, and should get more air time. Here it is again:

\frac{\mathrm{dE}(x)}{\mathrm{d}t} = \sum_i x_i \dot{p}_i +  \sum_i p_i \dot{x}_i

Hope I didn’t mess anything up! Further reading here (scroll down to the Day and Gandon 2005 book chapter). Also, Jeremy Fox of the Oikos blog has some great Price-equation-esque work in ecology (e.g.).

UPDATE: Check out the great references below that commenters have pointed out.

Advertisements
15 Comments leave one →
  1. Jeremy Fox permalink
    August 2, 2012 10:10 pm

    This is very useful, even for a supposed Price equation “expert” like me. Because I’ve only ever had occasion to use to the discrete time Price equation in my own work, I’ve never really sat down and worked my way through the continuous time version. Steve Ellner’s recently done some nice stuff with the continuous time version in Ecology Letters, using it as a way to partition ‘evolutionary’ and ‘ecological’ changes in mean phenotype over time.

    • August 2, 2012 10:28 pm

      I had been hoping you’d find this! Glad you found it useful. I had been playing around with some simple models of functional diversity dynamics, and found that I was just re-discovering the continuous-time Price equation. One of these days I’ll catch up to the state of the art! Thanks for the Ellner reference too. I think I saw him presenting that stuff at ESA one year but never got around to looking at the paper.

      • Jeremy Fox permalink
        August 2, 2012 10:48 pm

        If memory serves, Norberg et al. (PNAS 2000 or 2001?) is another application of the continuous time Price equation (specifically, the special case of Fisher’s Fundamental Theorem), in an ecological context. Norberg et al. were interested in how a diverse set of species can better “track” environmental fluctuations than any single species and so maintain higher total biomass on average, basically by harboring variation on which selection can act.

      • August 4, 2014 7:32 pm

        Another application of the continuous-time version of the Price Equation: I used it to track symbionts as characters that can be inherited vertically from one’s parents or horizontally from other individuals in a population (Am Nat 2007 170:542).

  2. August 12, 2012 12:28 pm

    Another good reference is ‘Applying population-genetic models in theoretical evolutionary epidemiology’ by Troy Day and Sylvain Gandon. It’s also from Ecology Letters.

    • August 12, 2012 3:34 pm

      Thanks for this great reference…I’m going to give it a close reading as soon as I have some time.

    • Jeremy Fox permalink
      August 13, 2012 6:15 pm

      As long as we’re forcing you to read everything cool that everyone’s ever done with the Price equation…

      Kerr and Godfray-Smith 2009 Evolution. Extending the (discrete) Price equation to handle immigration into the descendant population.

      Steven Frank’s recent review in the Journal of Evolutionary Biology. Sincerely flattered to have had some of my stuff cited in that one.

  3. Michael permalink
    November 11, 2014 5:43 pm

    Can you not just use the product rule from the off and skip the chain rule to derive the main equation?

    • November 11, 2014 5:53 pm

      Yes. Well … the sum rule followed by the product rule, which are both special cases of the chain rule. But yes, I was probably being overly complex.

      • Michael permalink
        November 12, 2014 1:49 pm

        Thanks Steven, just thought I would check I wasn’t missing the point. If you could help me understand a couple of other points I would really appreciate it! I am getting a bit lost at the start of the section ‘Aside: relative abundance dynamics’.

        1) In the chain rule for dpi/dt, where does the sum_j come from?

        2) I always end up with an N^2 in the denominator following the quotient rule. I have a derivative of N in the numerator, but I don’t see why it cancels.

        3) I’m also unsure how you then used the ‘some algebra’ following the expression for dpi/dnj.

        Thanks!

        Michael

      • November 12, 2014 5:57 pm

        1) This is a ‘total derivative’:

        2&3) Now that I think about, its probably easier to use the product rule on:

        n_i = p_i N
        \frac{dn_i}{dt} = \frac{dp_i}{dt} N + p_i\frac{dN}{dt}
        \frac{dp_i}{dt} =  \frac{1}{N}\left[\frac{dn_i}{dt} - p_i\frac{dN}{dt}\right]
        \frac{dp_i}{dt} =  \frac{1}{N}\left[r_in_i - p_i\sum_j r_jn_j\right]
        \frac{dp_i}{dt} =  \left[r_ip_i - p_i\sum_j r_jp_j\right]
        \frac{dp_i}{dt} =  \left[r_ip_i - p_iE(r)\right]
        \frac{dp_i}{dt} =  p_i\left[r_i - E(r)\right]

  4. Michael permalink
    November 13, 2014 9:25 am

    Thanks Steven!

  5. February 23, 2017 11:02 pm

    Very nice derivation. There is a tiny typo: I think that in the two equations below “chain rule”, it should be dp/dt, dx/dt full derivatives, not partial. (Though it seems wikipdeia on chain rule also has the same mistake). See, for example http://mathworld.wolfram.com/Derivative.html eq 45. If the derivative on the left hand side is partial, then you end with partials on the right, and if the derivative is full on left, you end with full on right. But this is really nitpicking.

Trackbacks

  1. Still lovin’ the chain rule… « ecology & stats
  2. Derivation of the continuous time Price equation | Dynamic Ecology

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: