This generic function calculates the Four-Gamete test (Hudson 1985).
Usage
## S4 method for signature 'GENOME'
recomb.stats(object,new.populations=FALSE,subsites=FALSE)
## S4 method for signature 'GENOME'
get.recomb(object)
Arguments
object
an object of class "GENOME"
new.populations
list of populations. default:FALSE
subsites
"transitions": SNPs that are transitions. "transversions": SNPs that are transversions. "syn": synonymous sites. "nonsyn": nonsynonymous sites. "exon": SNPs in exon regions. "intron": SNPs in intron regions. "coding": SNPs in coding regions (CDS). "utr": SNPs in UTR regions. "gene": SNPs in genes.
default:FALSE
Value
The return value is a modified object of class "GENOME"
————————————————————–
The following slots will be modified in the "GENOME" object
————————————————————–
Hudson.RM
Four-gamete test
References
Hudson, R. K. (1985). Statistical Properties of the Number of Recombination Events in the History of a Sample of DNA Sequences
Genetics 111 (1): 147-164.
Examples
# GENOME.class <- readData("homeAlignments")
# GENOME.class
# GENOME.class <- recomb.stats(GENOME.class)
# GENOME.class <- recomb.stats(GENOME.class,list(1:7,8:12))
# show the result:
# recomb.values <- get.recomb(GENOME.class)
# recomb.values[[1]] # first population !
# GENOME.class@region.stats@Hudson.RM