# Difference between revisions of "Drummond:PopGen"

Dadrummond (talk | contribs) (→Evolution is linear on a log-odds scale) |
Dadrummond (talk | contribs) (→Diffusion approximation: Added statistical analysis section.) |
||

Line 139: | Line 139: | ||

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

Insert math here. | Insert math here. | ||

+ | |||

+ | ==Statistical analysis of relative growth rates== | ||

+ | We have three strains, <math>i</math>, <math>j</math> and <math>r</math>, where <math>r</math> is a reference strain. | ||

+ | Strains <math>i</math> and <math>j</math> have fitness <math>w_i = e^{r_i}</math> and <math>w_j=e^{r_j}</math>. Define the selection coefficient <math>s_{ij} = \ln \frac{w_i}{w_j} = r_i - r_j</math> as usual. | ||

+ | We have data consisting of triples (number of generations, number of cells of type <math>i</math>, number of cells of type <math>r</math>). | ||

+ | |||

+ | What is the best estimate, and error, on <math>s_{ij}</math>? | ||

+ | |||

+ | ===Model=== | ||

+ | Given | ||

+ | Let <math>\Pr(s_{ij}=t) = \mathcal{N}(t;\mu_{ij}, \sigma^2_{ij})</math>. | ||

+ | |||

+ | ===Maximum-likelihood approach=== | ||

+ | |||

+ | |||

+ | ===Bayesian approach=== |

## Revision as of 07:43, 22 April 2010

## 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 be the number of organisms of type at time , and let be the *per-capita reproductive rate per generation*. If counts generations, then

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

where the last simplification follows from L'Hôpital's rule. Explicitly, let . Then

The solution to the equation

*instantaneous growth rate*for convenience.

## Continuous rate of change

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

Let and be the instantaneous rates of increase of type 1 and type 2, respectively. Then

This result says that the proportion of type 1 changes most rapidly when and most slowly when is very close to 0 or 1.

## Evolution is linear on a log-odds scale

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 ,

This differential equation has the solution

showing that the log-odds of finding type 1 changes linearly in time, increasing if and decreasing if .

## Diffusion approximation

Insert math here.

## Statistical analysis of relative growth rates

We have three strains, , and , where is a reference strain. Strains and have fitness and . Define the selection coefficient as usual. We have data consisting of triples (number of generations, number of cells of type , number of cells of type ).

What is the best estimate, and error, on ?

### Model

Given Let .