are_rare_vars: Which variants are rare in the genomic databases?

In bahlolab/varpr: Prioritises variants from the merged_collab file

Description

are_rare_vars reads in the specified genomic database columns from the vars file and outputs a logical vector indicating which rows/variants have a MAF less than or equal to the specified maf, or are novel.

Usage

are_rare_vars(vars = NULL, gendb = c("both_kg_esp", "kg", "esp"),
  maf = 0.05)

Arguments

vars	The data frame of vars.
gendb	The genomic database to choose (default: both 1KG and ESP6500).

Value

a logical vector with length equal to nrow(vars).

See Also

Logic.

Examples

## Not run: 
big_logical_vec <- are_rare_vars(vars, gendb = "both_kg_esp", maf = 0.05)

## End(Not run)

bahlolab/varpr documentation built on May 11, 2019, 5:26 p.m.