[BioC] outlier removal from gene chip
Weiwei Shi
helprhelp at gmail.com
Tue Sep 19 19:02:44 CEST 2006
dear listers:
I have a question on whether bioconductor has some tool-kit to detect
outliers and remove them.
my original dataset looks like this:
V1 V51 V53 V55 V57
1 -493249600 1.459459 -3.069444 -1.300000 1.935484
2 -1613096495 -1.139269 -5.525281 -16.592593 -1.831978
3 1626196571 -3.500000 -1.011662 2.223881 3.921053
4 -1397009217 -3.571429 1.685714 -1.180297 -6.807692
5 1428659728 -1.405405 -1.469004 -4.779754 -1.033708
6 459853658 -2.158879 -7.510823 -1.085581 -9.382979
7 530182506 -1.431677 -1.336343 -3.126437 4.878788
8 1173842263 1.215385 1.856410 -2.059794 -6.020833
9 28847 2.407895 -2.048889 -1.730337 -1.178947
10 -1961875610 2.864159 -2.301234 -4.733264 -1.172058
V1: internal probe id
the rests are different samples. the cells are fold-change of disease/normal.
summary of the sample columns( V51, ... V57) gives the following:
V51 V53 V55 V57
Min. :-482.000 Min. : -55.7342 Min. :-122.074 Min. :-14086.750
1st Qu.: -2.159 1st Qu.: -1.7312 1st Qu.: -2.125 1st Qu.: -1.831
Median : -1.199 Median : -1.0416 Median : -1.200 Median : -1.080
Mean : -0.918 Mean : 0.1662 Mean : -1.027 Mean : -1.874
3rd Qu.: 1.441 3rd Qu.: 1.5721 3rd Qu.: 1.419 3rd Qu.: 1.521
Max. : 198.434 Max. :1478.1639 Max. : 95.768 Max. : 683.519
My question is, is there any package which can detect those outliers
(like -14086.750)and remove them and get an "average" for each gene
(instead of each probe)?
Thank you.
Weiwei
--
Weiwei Shi, Ph.D
Research Scientist
GeneGO, Inc.
"Did you always know?"
"No, I did not. But I believed..."
---Matrix III
More information about the Bioconductor
mailing list