transversions_only: Remove transversions (C->T and G->A substitutions)

In admixr: An Interface for Running 'ADMIXTOOLS' Analyses

Description

Remove substitutions that are more likely to be a result of ancient DNA damage (C->T and G->A substitutions).

Usage

transversions_only(data, outfile = tempfile(fileext = ".snp"))

Arguments

data	EIGENSTRAT data object.
outfile	Path to an output snp file with coordinates of excluded sites.

Value

Updated S3 EIGENSTRAT data object with an additional 'exclude' slot specifying the path to the set of SNPs to be removed from a downstream analysis.

Examples

## Not run: # download an example genomic data set and prepare it for analysis
snps <- eigenstrat(download_data(dirname = tempdir()))

# perform the calculation only on transversions
snps_tv <- transversions_only(snps)
results_d <- d(W = "French", X = "Dinka", Y = "Altai", Z = "Chimp", data = snps_tv)

## End(Not run)

admixr documentation built on July 8, 2020, 6:19 p.m.