Entering edit mode
Xiaohui Wu
▴
280
@xiaohui-wu-4141
Last seen 10.3 years ago
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]]