getAdFromVcf: get the AD from vcf

In binmapr: Call Marker from Snp Data using Binmap

Description

get the AD from vcf

Usage

getAdFromVcf(file, keep = NULL, chromosome = NULL)

Arguments

file	vcf file path
keep	a vector store high confidence site, format should be "chr_1"
chromosome	vector, which chromosome to use if it is NULL, all chromosome will be include in analysis

Value

a list with AD and CHROM

Author(s)

Zhougeng Xu

Examples

library(vcfR)
data(vcfR_test)
orig_dir <- getwd()
temp_dir <- tempdir()
setwd( temp_dir )
write.vcf( vcfR_test, file = "test.vcf.gz" )
ad <- getAdFromVcf("test.vcf.gz")
ad
# return is full of NA, because the origin vcf don't have INFO/AD
setwd( orig_dir )

binmapr documentation built on Oct. 30, 2019, 12:14 p.m.