Hi all,

I want to cluster some genes based on their GO term annotation information in Arabidopsis and plot the profile. I tried some packages in Bioconductor including 'goProfiles', 'org.Hs.eg.db', but it says "No ancestors found for this GOTermsList ".

The code I tested is as follows:
----------------------
require(goProfiles)
require(org.At.tair.db)

#first, transform tair ID to entrez ID
x <- org.At.tairENTREZID
mapped_genes <- mappedkeys(x) 
tairENID <- as.list(x[mapped_genes])

#then, plot the profile
arab.MF <- basicProfile(genelist=tairENID, onto = "MF", level = 2, orgPackage = "org.At.tair.db")
---------------------

For the example of 'basicProfile' in goProfiles, it works, but it is for human, the code is:
----------------------
require(goProfiles)
data(prostateIds)
welsh.MF <- basicProfile(welsh01EntrezIDs, onto = "MF", level = 2, orgPackage = "org.Hs.eg.db")
singh.MF <- basicProfile(singh01EntrezIDs, onto = "MF",level = 2, orgPackage = "org.Hs.eg.db")
welsh.singh.MF <- mergeProfilesLists(welsh.MF,singh.MF, profNames = c("Welsh", "Singh"))
printProfiles(welsh.singh.MF, percentage = TRUE)
plotProfiles(welsh.MF, aTitle = "Welsh (2001). Prostate cancer data")
----------------------

Can you see what the problem is? Or is there any other way to cluster the genes and plot the profile? 

Thanking you,

Rae

	[[alternative HTML version deleted]]

