User:Hussein Alasadi/Notebook/stephens/2013/10/16

From OpenWetWare

(Difference between revisions)
Jump to: navigation, search
(Simulating evolution and then pooled sequencing)
(Simulating evolution and then pooled sequencing)
Line 18: Line 18:
(1) draw a coverage (n) from <math> Pois(\lambda) </math> where <math> \lambda \approx 40-50 </math>
(1) draw a coverage (n) from <math> Pois(\lambda) </math> where <math> \lambda \approx 40-50 </math>
-
(2) <math> f_i </math> (frequency of the ith SNP) <math> = B(n, f_i^{true}) </math>
+
(2) <math> f_i </math> (frequency of the ith SNP) ~ <math> B(n, f_i^{true}) </math>

Revision as of 11:52, 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})





Personal tools