[BioC] DEGraph graph format?
Hamid Bolouri
hbolouri at fhcrc.org
Fri Jun 8 02:43:12 CEST 2012
hello;
Can anyone tell me how to use DEGraph with the pathways in NCIGraphData?
The DEGraph Demo:
>data("Loi2008_DEGraphVignette", package="DEGraph")
>classData <- classLoi2008
>exprData <- exprLoi2008
>annData <- annLoi2008
>grList <- grListKEGG
>res <- testOneGraph(grList[[1]],exprData,classData,verbose=T,prop=0.2)
works fine for me. But replacing grList with NCI.cyList from NCIGraph:
>library(NCIgraphData)
>data("NCI-cyList")
> NCI.cyList[[1]]
A graphNEL graph with directed edges
Number of Nodes = 35
Number of Edges = 40
I get this error:
>res <- testOneGraph(NCI.cyList[[1]],exprData,classData,verbose=T,prop=0.2)
Keeping genes in the graph *and* the expression data set...
35 genes of the graph were not found in the expression data set:
chr [1:35] "6749854621221256793-pid_m_25632-674985462-829166685-pid_m_100726" ...
227 genes of the expression data set are absent from the graph:
chr [1:227] "31" "32" "207" "208" "355" "356" "369" "572" ...
Error: all.equal(dataGN, graphGN) is not TRUE
Keeping genes in the graph *and* the expression data set...done
I get the same error with 'reactome.cyList' graphs and with graphs generated by 'parseNCInetwork'.
Thanks
Hamid Bolouri
> sessionInfo()
R version 2.15.0 (2012-03-30)
Platform: i386-pc-mingw32/i386 (32-bit)
locale:
[1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252
[3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C
[5] LC_TIME=English_United States.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] NCIgraphData_0.99.4 DEGraph_1.8.0 R.utils_1.12.1
[4] R.oo_1.9.3 R.methodsS3_1.2.2
loaded via a namespace (and not attached):
[1] BiocGenerics_0.2.0 graph_1.34.0 grid_2.15.0 KEGGgraph_1.12.0
[5] lattice_0.20-6 mvtnorm_0.9-9992 NCIgraph_1.4.0 RBGL_1.32.0
[9] RCurl_1.91-1.1 RCytoscape_1.6.3 Rgraphviz_1.34.1 rrcov_1.3-01
[13] stats4_2.15.0 tools_2.15.0 XML_3.9-4.1 XMLRPC_0.2-4
More information about the Bioconductor
mailing list