BioNet package network analysis important questions for tutorial
Entering edit mode
svlachavas ▴ 780
Last seen 5 days ago
Germany/Heidelberg/German Cancer Resear…

Dear All,

i started to read and use the tutorial for the BioNet package for the functional analysis of biological networks(, in order to implemented afterwards for finding interesting functional modules in my colon cancer microarray data analysis DE genes. More specifically, in page 15 of the vignette i used instead of the command  ” -e ALL edges 001.txt -n ALL nodes 001.txt -N True -E False” the runHeinz command :

runHeinz(heinz.folder="heinz_1.67", heinz.e.file="ALL_edges_001", heinz.n.file="ALL_nodes_001", N=TRUE, E=FALSE, diff=-1, n=1)

where heinz_1.67 is the name of the heinz folder, but i got the following warning:

Warning message:
running command 'cd heinz_1.67; -e ALL_edges_001 -n ALL_nodes_001 -N True -E False -d -1 -s 1' had status 127 

Afterwards the tutorial continued without error, but 

datadir <- file.path(.path.package("BioNet"), "extdata") 

Error: '.path.package' is defunct.
Use 'path.package' instead.
See help("Defunct")

And when i used datadir <- file.path(path.package("BioNet"), "extdata") 

it didnt showed an error. These are some bugs with Bioconductor ? i have also downloaded CPLEX and installed it.

BioNet bioconductor affymetrixchip network analysis • 947 views
Entering edit mode

I am not the maintainer of BioNet (but I will alert them to this thread) but yes, .path.package() should not be used. I have updated it in release (version 1.26.1) and devel (version 1.27.1) so these will be available late tomorrow morning (Seattle time) via biocLite() and will clear up the error you're seeing (if you are using the release (3.0) or devel (3.1) versions of Bioconductor). I am not sure if this will clear up the other warning you saw, but hopefully the maintainer can speak to that.

Entering edit mode

i understand. Thus far,  i didnt noticed a problem in the output, also in the file for the cytoscape--just some small number differences in the number of nodes-but i guess this is due to possible updates of limma and other packages. Moreover, i would like also to ask you one important question regarding the function nodeData in page 16:

 nodeData(module, n = nodes(module), attr = "diff") <- fit2$coefficients[nodes(module), + 1] 

Because i also implemented limma in my DE analysis prior to BioNet :

conditions <- data.imputed.eset$condition
condition <- factor(conditions, levels(condition)[c(2,1)])
pairs <- factor(rep(1:13, each = 2))

design <- model.matrix(~condition+pairs)
fit <- lmFit(data.imputed.eset, design)
fit2 <- eBayes(fit, trend=TRUE)

symbols <- unlist(mget(featureNames(data.imputed.eset), env=hgu133aSYMBOL))
top2 <- topTable(fit2, coef="conditionCancer", number=nrow(fit2), genelist=symbols, adjust.method="fdr","none")

So, how i would transform the above function as in contrast with the vignette i have more than 1 coefficients and i want the fit2$coefficient= conditionCancer ?


Entering edit mode

Also i would like to ask one more important question regarding my dataset:

instead of importing the final image in  cytoscape, i used the function:

plotModule(module, diff.expr = diff, scores = scores)--but the nodes arent coloured. Should i import the xgmml image in cytoscape ?


Login before adding your answer.

Traffic: 258 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6