[BioC] issues with xps filter
Larry N. Singh
larryns at pcbi.upenn.edu
Wed Jun 30 19:01:29 CEST 2010
Hi everyone,
I'm having some issues with the madFilter in xps. Please see below:
# Please note I've already loaded the RMA data in data.rma
preflr <- Prefilter(mad = c(0.25))
rma.pfr <- prefilter(data.rma, "tmp_prefilter", getwd(),
filter=preflr, minfilters=1, logbase="log2")
tmp <- validData(rma.pfr)
d <- dim(tmp(tmp[, "FLAG"] == 1,])
> d
208 2
#### BUT
x <- validData(data.rma)
x <- log(x)/log(2)
mads <- apply(x, 1, mad)
> sum(mads > 0.25)
[1] 3641
I don't get the same numbers for MAD filtering.
Can someone please tell me what I'm missing here?
Thank you kindly for any help.
-Larry.
More information about the Bioconductor
mailing list