# Drummond:PopGen

### From OpenWetWare

(→Per-generation and instantaneous growth rates) |
(→Continuous rate of change) |
||

Line 113: | Line 113: | ||

|} | |} | ||

- | The logit function <math>\mathrm{logit} (p) = \ln {p \over 1-p}</math>, which takes <math>p \in [0,1] \to \mathbb{R}</math>, induces a more natural space for considering changes in frequencies. In logit terms, | + | The logit function <math>\mathrm{logit} (p) = \ln {p \over 1-p}</math>, which takes <math>p \in [0,1] \to \mathbb{R}</math>, induces a more natural space for considering changes in frequencies. Rather than tracking the proportion of type 1 or 2, we instead track their log odds. In logit terms, with <math>L_p(t) \equiv \mathrm{logit} (p(t))</math>, |

+ | |||

+ | <math>{\partial L_p(t) \over \partial t} = {\partial \over \partial t}\left(\ln {p(t) \over 1-p(t)}\right)</math> | ||

+ | |||

+ | <math>{\partial L_p(t) \over \partial t} = s t</math>. | ||

+ | |||

+ | That is, <math>L_p(t)</math>, the log-odds of finding type 1 in a random draw from the population, changes linearly in time with slope <math>s</math>. | ||

==Diffusion approximation== | ==Diffusion approximation== | ||

Insert math here. | Insert math here. |

## Revision as of 09:29, 2 April 2009

## Introduction

Here I will treat some basic questions in population genetics. For personal reasons, I tend to include all the algebra.

## Per-generation and instantaneous growth rates

What is the relationship between per-generation growth rates and the Malthusian parameter, the instantaneous rate of growth?

Let *n*_{i}(*t*) be the number of organisms of type *i* at time *t*, and let *R* be the *per-capita reproductive rate per generation*. If *t* counts generations, then

Now we wish to move to the case where *t* is continuous and real-valued.
As before,

where the last simplification follows from L'Hôpital's rule. Explicitly, let ε = Δ*t*. Then

The solution to the equation

*t*. We can define the

*instantaneous growth rate*

*r*= ln

*R*for convenience.

## Continuous rate of change

If two organisms grow at different rates, how do their proportions in the population change over time?

Let *r*_{1} and *r*_{2} be the instantaneous rates of increase of type 1 and type 2, respectively. Then

*p*(

*t*).

The logit function , which takes , induces a more natural space for considering changes in frequencies. Rather than tracking the proportion of type 1 or 2, we instead track their log odds. In logit terms, with ,

.

That is, *L*_{p}(*t*), the log-odds of finding type 1 in a random draw from the population, changes linearly in time with slope *s*.