prep_vcf_for_est_m_rd {whoa}R Documentation

convert a VCF into an 012,-1 matrix and read_depth bin matrix for estimation

Description

convert a VCF into an 012,-1 matrix and read_depth bin matrix for estimation

Usage

prep_vcf_for_est_m_rd(v, DF = "DP", minBin)

Arguments

v

a vcfR object into which a VCF file has been read

DF

Field to use for obtaining total read depth. Currently all that is available is DP.

minBin

minimum number of observations for each read depth bin

Value

This sends back a list with mat012: the 012,-1 matrix of genotypes, and dp_bins_list: the list returned by bin_depths().

Examples

pv <- prep_vcf_for_est_m_rd(lobster_buz_2000, minBin = 1000)

[Package whoa version 0.0.2 Index]