PopGenome (version 2.7.2)

Achaz.stats-methods: Achaz statistic

Description

Achaz statistic

Usage

# S4 method for GENOME
Achaz.stats(object,new.populations=FALSE,new.outgroup=FALSE,subsites=FALSE)

Arguments

object

an object of class "GENOME"

new.populations

list of populations. default:FALSE

new.outgroup

outgroup vector. 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

returned value is a modified object of class "GENOME" -------------------------------------------------------------- The following Slots will be modified in the "GENOME" object --------------------------------------------------------------

Yach

Achaz Y statistic

References

Achaz G.,2008 Testing for neutrality in samples with sequencing errors. Genetics 179: 1409.

Examples

# NOT RUN {
# GENOME.class <- readData("\home\Alignments")
# GENOME.class
# GENOME.class <- Achaz.stats(GENOME.class)
# GENOME.class <- Achaz.stats(GENOME.class,list(1:7,8:12))
# show the result:
# GENOME.class@Yach
 

# }