https://github.com/brentp/bio-playground
miscellaneous scripts for bioinformatics/genomics that dont merit their own repo.
https://github.com/brentp/bio-playground
Last synced: 3 months ago
JSON representation
miscellaneous scripts for bioinformatics/genomics that dont merit their own repo.
- Host: GitHub
- URL: https://github.com/brentp/bio-playground
- Owner: brentp
- License: mit
- Created: 2010-03-02T21:41:07.000Z (over 15 years ago)
- Default Branch: master
- Last Pushed: 2018-11-26T20:36:32.000Z (over 6 years ago)
- Last Synced: 2025-03-24T08:25:01.642Z (3 months ago)
- Language: C
- Homepage:
- Size: 979 KB
- Stars: 182
- Watchers: 15
- Forks: 104
- Open Issues: 1
-
Metadata Files:
- Readme: README.rst
- License: LICENSE
Awesome Lists containing this project
README
.. contents::
Miscellaneous scripts for bioinformatics that dont merit their own repo.
All under MIT License unless otherwise specified.
Ks Calc
--------
Abnormal nucleotide frequency tends to throw off normal procedures for estimating `evolutionary models `_. A practical situation is when calculating the Ks values for the grass genes where a significant portion of them are high-GC genes (see details `here `_). In the case of high GC genes, most of the substitutions will be either G or C, therefore the Jukes-Cantor correction will under-estimate the Ks values. The codon models in PAML, on the contrary, tend to over-estimate Ks values. The Ks calculator we want to implement here, ignores the inference of models (where it is difficult anyway, since you have very few sites to estimate the parameters in the model). Instead, we ask this: **given biased substitutions, lengths, run simulations and try to fit an evolutionary model based on the simulations.**
.. image:: http://chart.apis.google.com/chart?cht=lc&chls=8|8&chd=t2:65,65,65|75,75,75|40,50,80&chs=300x200&chm=V,FFFFFF,0,,25|@tObserved+alignment,,0,.05:.87,10|@twith+difference+D,,0,.05:.8,10|@tSimulate+alignments,,0,.55:.87,10|@twith+various+Ks,,0,.55:.8,10|@tProb(D)=0.3,,0,.3:.45,10|@tProb(D)=0.6,ff0000,0,.3:.37,10|@tProb(D)=0.4,,0,.3:.3,10|@tKs=0.1,,0,.15:.45,10|@tKs=0.2,ff0000,0,.15:.37,10|@tKs=0.3,,0,.15:.3,10|@tKs=...,808080,0,.15:.23,10|@tMaximum+Likelihood+Estimate,ff0000,0,.5:.37,10|a,990066,2,1,9.0&chma=0,0,30,0
:alt: method