[BioC] TopGo package
Dan Tenenbaum
dtenenba at fhcrc.org
Wed Mar 19 00:12:50 CET 2014
Hi Qian,
Please keep the Bioconductor list in the loop (hit Reply All) that way we can all teach and learn.
For best results, you need to provide all the code necessary to reproduce the problem. You haven't done that below, because we still don't know what finTab and day7 are.
Also, please run the function
sessionInfo()
...after you get your error, and then paste its output into your email. This tells us about your operating system, R version, and version of packages loaded, which can help us figure out what the problem is. This is what my sessionInfo() looks like:
> sessionInfo()
R Under development (unstable) (2013-10-12 r64048)
Platform: x86_64-apple-darwin10.8.0 (64-bit)
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] parallel stats graphics grDevices utils datasets methods
[8] base
other attached packages:
[1] topGO_2.15.2 SparseM_1.03 GO.db_2.10.1
[4] RSQLite_0.11.4 DBI_0.2-7 AnnotationDbi_1.25.15
[7] GenomeInfoDb_0.99.19 Biobase_2.23.6 BiocGenerics_0.9.3
[10] graph_1.41.3
loaded via a namespace (and not attached):
[1] grid_3.1.0 IRanges_1.21.34 lattice_0.20-27 stats4_3.1.0
[5] tools_3.1.0
Thanks,
Dan
----- Original Message -----
> From: "Qian Wang" <qian.wang725 at gmail.com>
> To: "Dan Tenenbaum" <dtenenba at fhcrc.org>
> Sent: Tuesday, March 18, 2014 4:06:17 PM
> Subject: Re: [BioC] TopGo package
>
>
>
>
> univ <- as.numeric( finTab $ day7 )
>
> univ <- as.factor( univ )
>
> names( univ ) <- finTab $ gene_id
>
> head( univ )
>
>
>
>
> # make topGOdata object
>
> day7GOdata <- new( "topGOdata" ,
>
> description = "Day 7 DE genes: Wt vs. Mut" ,
>
> ontology = "BP" ,
>
> allGenes = univ ,
>
> nodeSize = 10 ,
>
> annot = annFUN.org , mapping = " org.Dm.eg " )
>
> resFisher <- runTest( day7GOdata , alogorithm = "classic" , statistic
> = "fisher" )
>
> allRes <- GenTable( day7GOdata , classicFisher = resFisher , topNodes
> = 20 )
>
> par( fig =c( 0 , 1 , 0 , 1 ), mar =c( 3 , 5 , 2 , 0.5 ))
>
>
>
>
> this is the GO analysis Code. Not sure whether answers your
> questions. I am new to R.
>
>
> Thanks.
> Q
>
>
>
> On Tue, Mar 18, 2014 at 3:30 PM, Dan Tenenbaum < dtenenba at fhcrc.org >
> wrote:
>
>
>
>
>
> ----- Original Message -----
> > From: "Qian Wang" < qian.wang725 at gmail.com >
> > To: bioconductor at r-project.org
> > Sent: Tuesday, March 18, 2014 12:59:53 PM
> > Subject: [BioC] TopGo package
> >
> > Hi,
> >
> > I am using TopGo for my gene ontology analysis of the RNA-seq
> > data.I
> > am able to get the top lists of GO terms. But I also want to print
> > out the significant genes associated with a particular GO term. I
> > am
> > using printGenes.
> >
> >
> >
> > day7infGOdata <- new("topGOdata",
> > description = "Day 7 inf DE genes: DXV vs. PBS",
> > ontology = "BP",
> > allGenes = univ,
> > nodeSize=10,
> > annot = annFUN.org, mapping=" org.Dm.eg ")
>
> What's the value of univ? And what is the output of
>
> sessionInfo()
>
> ?
>
> Dan
>
>
>
> > resFisherD7inf <- runTest(day7infGOdata, alogorithm="classic",
> > statistic="fisher")
> > allResD7inf <- GenTable(day7infGOdata,
> > classicFisher=resFisherD7inf,
> > topNodes=20)
> >
> >
> > > goID <- allResD7inf[1,"GO:0050878"]
> > > gt <- printGenes(day7infGOdata, whichTerms=goID, chip=" org.Dm.eg
> > > ",
> > > numChar=40)
> >
> > The error message is :
> > Error in (function (classes, fdef, mtable) :
> > unable to find an inherited method for function ‘printGenes’ for
> > signature ‘"topGOdata", "NULL", "missing"’
> > >
> >
> > Can anybody help me?
> >
> > Thanks.
> > Q
> > [[alternative HTML version deleted]]
> >
> >
> > _______________________________________________
> > Bioconductor mailing list
> > Bioconductor at r-project.org
> > https://stat.ethz.ch/mailman/listinfo/bioconductor
> > Search the archives:
> > http://news.gmane.org/gmane.science.biology.informatics.conductor
>
>
>
>
> --
> Regards,
>
> Qian Wang
More information about the Bioconductor
mailing list