iqrxFilter: iqrxFilter Function to filter expression data within a...

Description Usage Arguments Details Value Author(s) References Examples

View source: R/iqrxFilter.R

Description

This function filters the expression of the supplied countDat object. iqrxFilter is a filtering function used to remove rows (genes) of various expression data.

Usage

1
iqrxFilter(cD, iqr_multi = 1.5, MEDIAN = FALSE, na.rm = TRUE)

Arguments

cD

A countDat object.

iqr_multi

Numeric multiplier; removes any outliers that are iqr_multi times the mid-50 percentile distance greater or less than the 25th and 75th percentiles, by default

MEDIAN

Boolean, Calculate RowMeans or RowMedians.

na.rm

Boolean, NA removal.

Details

This function filters the expression of the supplied cD object, based on a selected percentage cutoff and selected interquartile range multiplier. The function iqrxFilter will: 1) log-base two transform all RPKM values (obligatory); (2) remove any outliers that were 1.5 times the mid-50 percentile distance greater or less than the 75th and 25th percentiles (by default), respectively; and (3) uses mean values and instead of median values (by default).

Value

Returns a filtered countDat object.

Author(s)

AJ Vaestermark, JR Walters.

References

Jue et al. BMC Genomics 2013 14:150

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
data(hmel.data.doser)
reps <- c("Male", "Male", "Male", "Female", "Female", "Female")
annotxn <- data.frame("Chromosome" = factor(hmel.dat$chromosome,
levels = 1:21))
hm.tr<-hmel.dat$trxLength
hm<-new("countDat",data=hmel.dat$readcounts,seglens=hm.tr,
annotation=annotxn)
replicates(hm) <- reps
libsizes(hm) <- getLibsizes2(hm, estimationType = "total")
rpkm(hm) <- make_RPKM(hm)
f_hm <- iqrxFilter(hm)

avastermark19/doseR_0.99.4 documentation built on May 6, 2019, 12:08 p.m.