[Bioc-sig-seq] extract non-zero rows
Estefania Mancini
estefania.mancini at indear.com
Thu Aug 25 16:39:03 CEST 2011
Dear all
I have loaded and analyzed properly 4 454 dataset, corresponding to control and stress samples with their biological replicates.
I would like to know if is possible to filter, in my DGEList object
-which tags dont have zero in any column,
-which of these tags could be consider "housekeeping" (at least with logFC near 0)
The object DGEList looks like this:
>dup.data
An object of class "DGEList"
$samples
group lib.size norm.factors
A8_control control 77953 1
A8_stress stress 176860 1
mq_control control 98109 1
mq_stress stress 145839 1
pi_control control 132479 1
pi_stress stress 142484 1
tj_control control 65827 1
tj_stress stress 144278 1
I have tried to filter using the suggested function:
>dup.de.filter <- dup.data[rowSums(dup.data$counts) >= 0, ]
or with
>dup.de.filter <- dup.data[rowSums(dup.data$counts) >= 1, ]
but have no changes at all. I have many rows which 0 and 1 read in some column which should be excluded.
Also:
dup.de.com
An object of class "DGEExact"
$table
logConc logFC p.value
Glyma13g11940.8 -2.588833 0.26176050 0.7348221
Glyma13g11900.1 -2.875548 0.03020441 0.9688072
Glyma09g24780.1 -3.501041 -0.12108619 0.8754371
Glyma13g12050.1 -3.224648 0.03036675 0.9691009
Glyma13g12070.1 -3.743064 0.14416487 0.8521188
19860 more rows ...
$comparison
[1] "control" "stress"
$genes
NULL
Thanks in advance,
Estefania
More information about the Bioc-sig-sequencing
mailing list