filter.allele_ages: Filter allele_ages by quality of tree

Description Usage Arguments Value Examples

View source: R/manipulate_sele.R

Description

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

Usage

1
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

1
2
3
4
5
6
7
8
9
# 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 Jan. 6, 2020, 7:11 a.m.