Additional high-scoring subnetworks in BioNet
1
0
Entering edit mode
enricoferrero ▴ 660
@enricoferrero-6037
Last seen 3.0 years ago
Switzerland

I'd like to use BioNet to identify high-scoring subnetworks in my data.

The tutorial clearly explains how to identify the top scoring subnetwork (I don't have the CPLEX library so I'm using runFastHeinz()):

library(BioNet)
library(DLBCL)
data(dataLym)
data(interactome)

pvals <- cbind(t=dataLym$t.pval, s=dataLym$s.pval)
rownames(pvals) <- dataLym$label
pval <- aggrPvals(pvals, order=2, plot=FALSE)

subnet <- subNetwork(dataLym$label, interactome)
subnet <- rmSelfLoops(subnet)
subnet
A graphNEL graph with undirected edges
Number of Nodes = 2559
Number of Edges = 7788

fb <- fitBumModel(pval, plot=FALSE)
scores <- scoreNodes(subnet, fb, fdr=0.001)

module <- runFastHeinz(subnet, scores)
module
A graphNEL graph with undirected edges
Number of Nodes = 37
Number of Edges = 44

But what about other high-scoring subnetworks? It is likely that in a network with 2259 nodes there will be other interesting active subnetworks to explore other than the first one. How can I retrieve these?

Bonus question: are there any parameters that I can set to increase or decrease the size of the resulting subnetworks?

Thank you.

bionet network • 1.5k views
ADD COMMENT
1
Entering edit mode
@branislav-misovic-4248
Last seen 5.6 years ago
Netherlands/Amsterdam

hi  Enrico 
 Last year ,during the Networks course,  I  asked same question Gunnar W. Klau  (one of  the authors who worked with Heinz &  link to Cytoscape Examine plugin )   ,   he said :  "you can remove your 1st Module nodes from you PPI DB  and rerun the Heinz to get 2nd Module ... in the past Heinz  gave more than 1  Module" .
  I do not know the reason why they stopped providing that output .

ADD COMMENT
0
Entering edit mode

Thanks Branislav, that sounds like a reasonable (although not ideal) workaround.

Hopefully someone from the dev team can step in and comment.

ADD REPLY

Login before adding your answer.

Traffic: 527 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6