filter.allele_ages: Filter allele_ages by quality of tree

In leospeidel/relater: What the Package Does (One Line, Title Case)

Description

This function filters out selection pvalues on bad trees using summary statistics on tree quality stored in qual.

Usage

filter.allele_ages(allele_ages, qual)

Arguments

allele_ages	data.table. Obtained from get.allele_ages()
qual	data.frame. qual file

Value

Returns a data table.

Examples

# get allele ages
mut         <- read.mut(system.file("extdata/example.mut.gz", package = "relater"), CHR = 1)
sele        <- read.sele(system.file("extdata/example.sele.gz", package = "relater"), CHR = 1)
freq        <- read.freq(system.file("extdata/example.freq.gz", package = "relater"), CHR = 1)
allele_ages <- get.allele_ages(mut, freq, sele)

# read quality file
qual        <- read.qual(system.file("extdata/example.qual.gz", package = "relater"), CHR = 1)
filter.allele_ages(allele_ages, qual)

leospeidel/relater documentation built on Feb. 3, 2022, 8:22 p.m.