[BioC] Bug in getOntoly-function (package annotate)?
Jukka Hiissa
jahiis at utu.fi
Thu Oct 2 09:00:48 CEST 2008
Dear BioC users,
I am having problems with the getOntology-function (annotate-package v.
1.18 ) The command
getOntology(getGO( "7650504", "illuminaHumanBCv2")[[1]], "BP" )
[1] "GO:0032515" "GO:0032516"
returns two GO identifiers, but they are not correct since their
ontologies are in fact NA
$`GO:0032515`
$`GO:0032515`$GOID
[1] "GO:0032515"
$`GO:0032515`$Evidence
[1] "IDA"
$`GO:0032515`$Ontology
[1] NA
$`GO:0032516`
$`GO:0032516`$GOID
[1] "GO:0032516"
$`GO:0032516`$Evidence
[1] "IDA"
$`GO:0032516`$Ontology
[1] NA
I looked at the getOntology source code and I think that the problem is
that the NAs are removed from the index vector but not from the list
containing the GO-elements. The last two lines in the getOntology
function are
onts = onts[!is.na(onts)]
unique(names(inlist[onts %in% which]))
where the onts contains the ontology labels, e.g.
GO:0000159 GO:0005524 GO:0005625 GO:0008160 GO:0008601 GO:0008601
NA "MF" "CC" "MF" "MF" "MF"
GO:0016887 GO:0019211 GO:0032515 GO:0032516 GO:0035307 GO:0035308
"MF" "MF" NA NA "BP" "BP"
GO:0042803 GO:0046982 GO:0051721
"MF" "MF" "MF"
After the NAs are removed from the onts vector the indexes of the two
BP-ontology cases are 8 and 9, but these are not the correct indexes for
the inlist list, which still contains the three missing ontology cases.
names( inlist )
[1] "GO:0000159" "GO:0005524" "GO:0005625" "GO:0008160" "GO:0008601"
[6] "GO:0008601" "GO:0016887" "GO:0019211" "GO:0032515" "GO:0032516"
[11] "GO:0035307" "GO:0035308" "GO:0042803" "GO:0046982" "GO:0051721
Best regards,
Jukka Hiissa
More information about the Bioconductor
mailing list