[BioC] help with ath1121501
Martin Olivier
martinol at ensam.inra.fr
Fri Oct 13 14:32:09 CEST 2006
Hi all,
I have (again) some problems to use the package ath1121501.
I don't understand the difference between ath1121501G02PROBE and
ath1121501GO2ALLPROBES...
I made a list listGO2GENE to have for each term GO all the genes that belong
to this term GO. For this, I use ath1121501G02PROBE or
ath1121501GO2ALLPROBES. The problem is that the result depends on the choice
between ath1121501G02PROBE or ath1121501GO2ALLPROBES.
For example, if I use ath1121501G02PROBE, listGO2GENE["GO:0008104"] gives me
three genes in the the GO term GO:0008104
$"GO:0008104"
[1] "AT1G48090" "AT3G50380" "AT5G24740"
and if I use ath1121501GO2ALLPROBES, I obtain more than 570 genes
for this GO term...
So what is the right file I have to use to obtain the list that gives me
all the genes
(in TAIR format) for a particular GO term.
Here under, I give you the code I used to obtain the list
listGO2GENE...may be
the pb is in my code??But it corresponds to a part of the code GOHyperG,
so it
would be right...
goV<-mget(ls(ath1121501GO2PROBE),ath1121501GO2PROBE)
listGO2GENE<- sapply(goV, function(x) {
if (length(x) == 0 || is.na(x))
return(NA)
lls <- unique(unlist(mget(x,
get("ath1121501LOCUSID",mode="environment"))))
lls <- lls[!is.na(lls)]
lls
})
Thanks for your help,
Olivier.
More information about the Bioconductor
mailing list