kaks {seqinr} | R Documentation |
Ks and Ka are respectively the number of substitutions per synonymous site and per nonsynonymous site between two protein-coding genes. The ratio of nonsynonymous (Ka) to synonymous (Ks) nucleotide substitution rates is an indicator of selective pressures on genes. A ratio significantly greater than 1 indicates positive selective pressure. A ratio around 1 indicates either neutral evolution at the protein level or an averaging of sites under positive and negative selective pressures. A ratio less than 1 indicates pressures to conserve protein sequence (i.e. purifying selection). This function estimates the Ka and Ks values for a set of aligned sequences using the method published by Li (1993) and gives the associated variance matrix.
kaks(x)
x |
An object of class alignment |
ks |
matrix of ks values |
ka |
matrix of ka values |
vks |
variance matrix of ks |
vka |
variance matrix of ka |
When the alignment does not contain enough information (i.e we approach saturation), the Ka and Ks values take the value 10. Negative values indicate that Ka and Ks can not be computed.
D. Charif
Li WH. (1993) Unbiased estimation of the rates of synonymous and nonsynonymous substitution.
J Mol Evol., Jan;36(1):96-9.
Hurst LD. (2002) The Ka/Ks ratio: diagnosing the form of sequence evolution.
Trends Genet., Sept;18(9):486.
The C programm implementing this method was provided by M. Gouy.
To have an overview of the seqinR's functionnality, please consult this vignette:
Charif, D., Lobry, J.R. (2005) SeqinR: a contributed package to the R project for statistical
computing devoted to biological sequences retrieval and analysis. Springer Verlag, Biological and Medical Physics/Biomedical Series, in preparation.
s = read.alignment(File=system.file("sequences/test.phylip", package = "seqinr"),format="phylip") kaks(s)