ckaksCodon: Compute the conditional kaks(conditional selection pressure)...

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/ckaks_codon.R

Description

Compute the conditional kaks(conditional selection pressure) among codons, the amino mutation in a site will be treated as whole in computation.

Usage

1
ckaksCodon(seq_formated, kaks = TRUE, lod_cut = 2, setPosition = c())

Arguments

seq_formated

Formated alignment sequence. i.e. the result after the treatment of DataFormatCorMut.

kaks

A logical variable to indicate whether kaks is turn on or off, if kaks is TRUE, conditional kaks will be computed only among positive seelction sites, or if kaks is FALSE, conditional kaks will be computed only among all sites of sequences.

lod_cut

The LOD confidence score cutoff, the default value is 2. If lod is larger than 2, it means the positive selection of individual site or the conditional selection pressure among sites are significant.

setPosition

The positions of sequence to compute. setPosition should be a vector of interger type indicating the positions to compute.

Value

A object of ckaks class will be return. ckaks includes two slots of matrix:ckaks and lod, which indicate the ckaks and lod confidence score respectively.

Author(s)

Zhenpeng Li

References

Chen, L., Perlina, A. & Lee, C. J. Positive selection detection in 40,000 human immunodeficiency virus (HIV) type 1 sequences automatically identifies drug resistance and positive fitness mutations in HIV protease and reverse transcriptase. Journal of virology 78,3722-3732 (2004).

See Also

filterSites,plot.ckaks,ckaksAA

Examples

1
2
3
examplefile=system.file("extdata","PI_treatment.aln",package="CorMut")
example=seqFormat(examplefile)
result=ckaksCodon(example)

CorMut documentation built on April 28, 2020, 7:09 p.m.