(Difference between revisions)
 Revision as of 00:19, 17 October 2013 (view source) (→Simulating evolution and then pooled sequencing)← Previous diff Revision as of 00:48, 17 October 2013 (view source) (→Simulating evolution and then pooled sequencing)Next diff → Line 12: Line 12: * '''We then evolve the populations for g generations with recombination''' * '''We then evolve the populations for g generations with recombination''' - To do this we can use Kessner's [http://arxiv.org/abs/1310.3234 Forqs]. Forqs allows the user to specify the recombination rate file. Which might be useful in the future. + To do this we can use Kessner's [http://arxiv.org/abs/1310.3234 Forqs]. Forqs allows the user to specify the recombination rate file which might be useful in the future. * '''Simulate pooled sequencing''' * '''Simulate pooled sequencing''' + + (1) draw a coverage from $Pois(\lambda)$ where $\lambda \approx 40-50$

## Revision as of 00:48, 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 from Pois(λ) where $\lambda \approx 40-50$