View source: R/filter_allele_balance.R
filter_allele_balance | R Documentation |
This function requires a vcfR object as input, and returns a vcfR object filtered to convert heterozygous sites with an allele balance falling outside of the specified ratio to 'NA'. If no ratio is specified, a default .25-.75 limit is imposed. From the dDocent filtering page "Allele balance: a number between 0 and 1 representing the ratio of reads showing the reference allele to all reads, considering only reads from individuals called as heterozygous, we expect that the allele balance in our data (for real loci) should be close to 0.5".
filter_allele_balance(vcfR, min.ratio = NULL, max.ratio = NULL)
vcfR |
a vcfR object |
min.ratio |
minumum allele ratio for a called het |
max.ratio |
maximum allele ratio for a called het |
An identical vcfR object, except that genotypes failing the allele balance filter have been converted to 'NA'.
filter_allele_balance(vcfR = SNPfiltR::vcfR.example)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.