filter.noncovered.bins: Filter non-covered bins

Description Usage Arguments Details Value Author(s)

View source: R/filter.noncovered.bins.R

Description

Find bins that are not covered and remove them from the input 'bins.df'. This function is useful for targeted sequencing when the binning was unaware of the targeted regions. In practice, the bin counts across a few samples are used to assess if a bin is covered or not.

Usage

1
2
filter.noncovered.bins(bins.df, files.df, nb.samples = 10,
  bc.med.min = NULL, plot = TRUE)

Arguments

bins.df

the data.frame with the bins information that will be subsetted.

files.df

the data.frame with the file information, e.g. location of the bin count files.

nb.samples

the number of samples to use. Default is 10.

bc.med.min

If non-NULL, the minimum median coverage to use for the filtering. If NULL (default), this value will be estimated.

plot

Display a graph ? Default is TRUE.

Details

If 'bc.med.min' is NULL, the function will quickly model the coverage due to off-target reads and set the minimum threshold as (median+3sd) of the off-target coverage. The off-tagert coverage is approximated by the maximum non-null coverage in bins with median coverage of 0. If a bin has 0 read in more than half of the samples, it's likely a non-covered bins. However sometimes a few samples have some reads there. The distribution of these noise-coverage is used to fix the minimum threshold in median coverage.

Value

a subset of the input 'bins.df' data.frame

Author(s)

Jean Monlong


jmonlong/PopSV documentation built on Dec. 9, 2018, 10:10 a.m.