(Difference between revisions)
 Revision as of 00:48, 17 October 2013 (view source) (→Simulating evolution and then pooled sequencing)← Previous diff Revision as of 00:50, 17 October 2013 (view source) (→Simulating evolution and then pooled sequencing)Next diff → Line 16: Line 16: * '''Simulate pooled sequencing''' * '''Simulate pooled sequencing''' - (1) draw a coverage from $Pois(\lambda)$ where $\lambda \approx 40-50$ + (1) draw a coverage (n) from $Pois(\lambda)$ where $\lambda \approx 40-50$ - + (2) $f_i$ (frequency of the ith SNP) $= B(n, f_i^{true})$

## Revision as of 00:50, 17 October 2013

Analyzing pooled sequenced data with selection Main project page
Previous entry

## Simulating evolution and then pooled sequencing

• We start with a population of N individuals with a certain SNP distribution P(M) ~ N(μ,Σ).

To do this we can use Dick Hudson's MS to lay down neutral variants (but the distribution really does not matter).

• We then evolve the populations for g generations with recombination

To do this we can use Kessner's Forqs. Forqs allows the user to specify the recombination rate file which might be useful in the future.

• Simulate pooled sequencing

(1) draw a coverage (n) from Pois(λ) where $\lambda \approx 40-50$ (2) fi (frequency of the ith SNP) $= B(n, f_i^{true})$