summary.DNAmat_class: Summary statistics

Description Usage Arguments Value See Also Examples

View source: R/summary_DNAmat_class.R

Description

Calculates summary statistics useful for analyzing DNA.

Usage

1
2
## S3 method for class 'DNAmat_class'
summary(DNAmat)

Arguments

DNAmat

matrix; either a single-nucleotide polymorphism matrix or a segregating sites matrix of class 'DNAmat_class'.

Value

List containing the following components:

SFS

The site frequency spectrum

Watterson

Watterson's estimator.

pairwDiff

The pairwise difference estimator.

TajimaD

Tajima's D

See Also

SFS, mutRate, TajimaD

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
# Using the function simDNAseq:
summary(simDNAseq(n = 12, seqLen = 20, mutRate = 5,
               popType = "varPop", expRate = 1.5))

# Creating a segregating sites matrix by hand:
DNAmat <- matrix(c(0,1,0,0,
                   0,1,0,1,
                   0,0,0,0,
                   0,0,0,0), 4,4, byrow=TRUE)
class(DNAmat) <- "DNAmat_class"
summary(DNAmat)

aumath-advancedr2019/simDNA documentation built on Nov. 27, 2019, 11 a.m.