fiter.WCs: Filter low count WCs.

Description Usage Arguments Value

View source: R/WaveQTL_preprocess_funcs.R

Description

This function filters out some WCs that are computed based on very low counts. This function considers WCs as low count if the total counts used in their computation were less than or equal to ‘meanR.thresh’ per individual on average.

Usage

1
fiter.WCs(Data, meanR.thresh)

Arguments

Data

matrix (or a vector when N = 1) with N (# of samples) by T (# of bps in a region); This matrix contains original data before a wavelet transform; Here, T should be a power of 2.

meanR.thresh

If average reads across individuals <= meanR.thresh, those WCs are filtered out.

Value

filtered.WCs a vector of length T; t-th element indicates whether t-th WC in output from the function FWT is filtered (0) or not (1).


jean997/cfdrSims documentation built on May 18, 2019, 11:43 p.m.