Question: simLL method in GOstats to measure distances between yeast genes
0
gravatar for Maria Persico
14.4 years ago by
Maria Persico100
Maria Persico100 wrote:
Dear All, I did some modifications in function simLL to allow this function access to ygs98 annotation: this is the code: simLL2<-function (ll1, ll2, Ontology = "MF", measure = "LP", dropCodes = NULL) { require(ygs98) #or library(ygs98) wh = match.arg(Ontology, c("MF", "BP", "CC")) ll1GO = .getWHEC(ll1, wh, dropCodes) ll2GO = .getWHEC(ll2, wh, dropCodes) dataenv = get(paste("GO", wh, "PARENTS", sep = ""), mode = "environment") g1 = GOGraph(ll1GO, dataenv) g2 = GOGraph(ll2GO, dataenv) ..... ....etc etc Then, I run the function with these 2 yeast geneID: aa = simLL2("855134", "851236", "BP") and I obtain > aa [1] NA If I run the function with 2 human llID(aa = simLL2("9184", "3547", "BP")), the function give me the right result.This means that the function reads the hgu95av2 library. How can I adapt the simmLL function to work with yeast annotation? Thanks, Maria Maria Persico MINT database, Cesareni Group Universita' di Tor Vergata, via della Ricerca Scientifica 00133 Roma, Italy Tel: +39 0672594315 FAX: +39 0672594766 e-mail: maria@cbm.bio.uniroma2.it
yeast hgu95av2 ygs98 • 659 views
ADD COMMENTlink modified 14.4 years ago by rgentleman5.5k • written 14.4 years ago by Maria Persico100
Answer: simLL method in GOstats to measure distances between yeast genes
0
gravatar for rgentleman
14.4 years ago by
rgentleman5.5k
United States
rgentleman5.5k wrote:
Hi Maria, The simLL function relies on have Entrez gene (used to be LocusLink IDs). Which seem to be the type of identifiers you are providing - so you have no need of ygs98 - that I can see. Now the first one of your identifiers: "855134" if you check does not have any GO annotation - so you won't get a GO graph for it. And the human one works because both genes are annotated for the GO category you are interested in. The result does not mean that the function reads any other library - it does just what it says it does. Locus/Entrez to GO - but if there is no GO anntation then you will get NA. Best wishes, Robert Maria Persico wrote: > Dear All, > I did some modifications in function simLL to allow this function access > to ygs98 annotation: > > this is the code: > > simLL2<-function (ll1, ll2, Ontology = "MF", measure = "LP", dropCodes = > NULL) > { > require(ygs98) > #or library(ygs98) > wh = match.arg(Ontology, c("MF", "BP", "CC")) > ll1GO = .getWHEC(ll1, wh, dropCodes) > ll2GO = .getWHEC(ll2, wh, dropCodes) > dataenv = get(paste("GO", wh, "PARENTS", sep = ""), mode = > "environment") > g1 = GOGraph(ll1GO, dataenv) > g2 = GOGraph(ll2GO, dataenv) > ..... > ....etc etc > > Then, I run the function with these 2 yeast geneID: > > aa = simLL2("855134", "851236", "BP") > > and I obtain > >>aa > > [1] NA > > If I run the function with 2 human llID(aa = simLL2("9184", "3547", > "BP")), the function give me the right result.This means that the function > reads the hgu95av2 library. > > How can I adapt the simmLL function to work with yeast annotation? > > Thanks, > > Maria > > > > > Maria Persico > MINT database, Cesareni Group > Universita' di Tor Vergata, via della Ricerca Scientifica > 00133 Roma, Italy > Tel: +39 0672594315 > FAX: +39 0672594766 > e-mail: maria@cbm.bio.uniroma2.it > > _______________________________________________ > Bioconductor mailing list > Bioconductor@stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor >
ADD COMMENTlink written 14.4 years ago by rgentleman5.5k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 333 users visited in the last hour