User:Hussein Alasadi/Notebook/stephens/2013/10/16: Difference between revisions
From OpenWetWare
(fix raw html notebook nav) |
|||
| (4 intermediate revisions by one other user not shown) | |||
| Line 2: | Line 2: | ||
|- | |- | ||
|style="background-color: #EEE"|[[Image:owwnotebook_icon.png|128px]]<span style="font-size:22px;"> Analyzing pooled sequenced data with selection</span> | |style="background-color: #EEE"|[[Image:owwnotebook_icon.png|128px]]<span style="font-size:22px;"> Analyzing pooled sequenced data with selection</span> | ||
|style="background-color: #F2F2F2" align="center"| | |style="background-color: #F2F2F2" align="center"|[[File:Report.png|frameless|link={{#sub:{{FULLPAGENAME}}|0|-11}}]][[{{#sub:{{FULLPAGENAME}}|0|-11}}|Main project page]]<br />{{#if:{{#lnpreventry:{{FULLPAGENAME}}}}|[[File:Resultset_previous.png|frameless|link={{#lnpreventry:{{FULLPAGENAME}}}}]][[{{#lnpreventry:{{FULLPAGENAME}}}}{{!}}Previous entry]] }}{{#if:{{#lnnextentry:{{FULLPAGENAME}}}}|[[{{#lnnextentry:{{FULLPAGENAME}}}}{{!}}Next entry]][[File:Resultset_next.png|frameless|link={{#lnnextentry:{{FULLPAGENAME}}}}]]}} | ||
|- | |- | ||
| colspan="2"| | | colspan="2"| | ||
| Line 8: | Line 8: | ||
==Simulating evolution and then pooled sequencing== | ==Simulating evolution and then pooled sequencing== | ||
* '''We start with a population of N individuals with a | * '''We start with a population of N individuals with a SNP distribution <math> P(M)</math> ~ <math> N(\mu, \Sigma) </math>.''' | ||
We can use Dick Hudson's MsSel (MS with selection) to generate a SNP distribution. | |||
* '''We then evolve the populations for g generations with selection for SNP j''' | |||
In MsSel, we can specify an allele frequency trajectory for a site. We can generate a trajectory based on a simple W-F simulation with selection. | |||
Or, we may have to use [http://arxiv.org/abs/1310.3234 forqs]. | |||
* '''Simulate pooled sequencing''' | * '''Simulate pooled sequencing''' | ||
| Line 19: | Line 24: | ||
(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> | ||
Latest revision as of 23:29, 26 September 2017
 Analyzing pooled sequenced data with selection
|
 Main project page Previous entry
|
Simulating evolution and then pooled sequencing
We can use Dick Hudson's MsSel (MS with selection) to generate a SNP distribution.
In MsSel, we can specify an allele frequency trajectory for a site. We can generate a trajectory based on a simple W-F simulation with selection.
(1) draw a coverage (n) from [math]\displaystyle{ Pois(\lambda) }[/math] where [math]\displaystyle{ \lambda \approx 40-50 }[/math] (2) [math]\displaystyle{ f_i }[/math] (frequency of the ith SNP) ~ [math]\displaystyle{ B(n, f_i^{true}) }[/math]
| |
