Easiest way to calculate Ka Ks ratio and divergence time
|
|
The Ka/Ks ratio is used to estimate the nature of evolution among neutral, purifying selection and beneficial mutations acting on a set of homologous protein-coding genes. Although there are several wonderful programs out there to calculate the Ka, Ks and Ka/Ks ratio but kaks function of R package seqinr is the easiest way to do it. The function kaks use the method published by Li et al (J. Mol. Evol., 36:96-99, 1993) to calculate the Ka, Ks and Ka/Ks ratio.
In the following script take aligned the fasta sequences in the clustal format and, finally, R package seqinr will be used to calculate the Ka, Ks and Ka/Ks ratio.
Cheat Sheet to Install and work with R on Ubuntu HERE
Prerequisite
We need the following R libraries to run the script- seqinr
- ape
- phangorn
Ka/Ks_calculator.R
|
Was This Post Useful? Add This To Del.icio.us Share on Facebook StumbleUpon This Add to Technorati Share on Twitter |
Subscribe to:
Post Comments (Atom)
No comments:
Post a Comment
Have Problem ?? Drop a comments here!