### R code from vignette source 'vignettes/BioNet/inst/doc/Tutorial.Rnw' ################################################### ### code chunk number 1: Tutorial.Rnw:47-49 ################################################### options(width=60) ps.options(family="sans") ################################################### ### code chunk number 2: Tutorial.Rnw:52-56 ################################################### library(BioNet) library(DLBCL) data(dataLym) data(interactome) ################################################### ### code chunk number 3: Tutorial.Rnw:60-63 ################################################### pvals <- cbind(t=dataLym$t.pval, s=dataLym$s.pval) rownames(pvals) <- dataLym$label pval <- aggrPvals(pvals, order=2, plot=FALSE) ################################################### ### code chunk number 4: Tutorial.Rnw:67-70 ################################################### subnet <- subNetwork(dataLym$label, interactome) subnet <- rmSelfLoops(subnet) subnet ################################################### ### code chunk number 5: Tutorial.Rnw:76-78 ################################################### fb <- fitBumModel(pval, plot=FALSE) scores <- scoreNodes(subnet, fb, fdr=0.001) ################################################### ### code chunk number 6: Tutorial.Rnw:85-88 ################################################### module <- runFastHeinz(subnet, scores) logFC <- dataLym$diff names(logFC) <- dataLym$label ################################################### ### code chunk number 7: Tutorial.Rnw:95-96 ################################################### plotModule(module, scores=scores, diff.expr=logFC) ################################################### ### code chunk number 8: Tutorial.Rnw:151-155 ################################################### library(BioNet) library(DLBCL) data(exprLym) data(interactome) ################################################### ### code chunk number 9: Tutorial.Rnw:163-164 ################################################### exprLym ################################################### ### code chunk number 10: Tutorial.Rnw:170-171 ################################################### interactome ################################################### ### code chunk number 11: Tutorial.Rnw:178-180 ################################################### network <- subNetwork(featureNames(exprLym), interactome) network ################################################### ### code chunk number 12: Tutorial.Rnw:185-187 ################################################### network <- largestComp(network) network ################################################### ### code chunk number 13: Tutorial.Rnw:197-201 ################################################### library(genefilter) library(impute) expressions <- impute.knn(exprs(exprLym))$data t.test <- rowttests(expressions, fac=exprLym$Subgroup) ################################################### ### code chunk number 14: Tutorial.Rnw:204-205 ################################################### t.test[1:10, ] ################################################### ### code chunk number 15: Tutorial.Rnw:212-215 ################################################### library(xtable) top.table <- xtable(t.test[1:10,], display=c("s", "f", "f", "f")) print(top.table, floating=FALSE) ################################################### ### code chunk number 16: Tutorial.Rnw:225-230 ################################################### data(dataLym) ttest.pval <- t.test[, "p.value"] surv.pval <- dataLym$s.pval names(surv.pval) <- dataLym$label pvals <- cbind(ttest.pval, surv.pval) ################################################### ### code chunk number 17: Tutorial.Rnw:239-240 ################################################### pval <- aggrPvals(pvals, order=2, plot=FALSE) ################################################### ### code chunk number 18: Tutorial.Rnw:247-249 ################################################### fb <- fitBumModel(pval, plot=FALSE) fb ################################################### ### code chunk number 19: Tutorial.Rnw:252-258 ################################################### dev.new(width=13, height=7) par(mfrow=c(1,2)) hist(fb) plot(fb) dev.off() ################################################### ### code chunk number 20: Tutorial.Rnw:273-274 ################################################### plotLLSurface(pval, fb) ################################################### ### code chunk number 21: Tutorial.Rnw:283-284 ################################################### scores <- scoreNodes(network=network, fb=fb, fdr=0.001) ################################################### ### code chunk number 22: Tutorial.Rnw:291-294 ################################################### network <- rmSelfLoops(network) writeHeinzEdges(network=network, file="lymphoma_edges_001", use.score=FALSE) writeHeinzNodes(network=network, file="lymphoma_nodes_001", node.scores = scores) ################################################### ### code chunk number 23: Tutorial.Rnw:311-313 ################################################### datadir <- file.path(.path.package("BioNet"), "extdata") dir(datadir) ################################################### ### code chunk number 24: Tutorial.Rnw:319-322 ################################################### module <- readHeinzGraph(node.file=file.path(datadir, "lymphoma_nodes_001.txt.0.hnz"), network=network) diff <- t.test[, "dm"] names(diff) <- rownames(t.test) ################################################### ### code chunk number 25: Tutorial.Rnw:325-326 ################################################### plotModule(module, diff.expr=diff, scores=scores) ################################################### ### code chunk number 26: Tutorial.Rnw:348-351 ################################################### sum(scores[nodes(module)]) sum(scores[nodes(module)]>0) sum(scores[nodes(module)]<0) ################################################### ### code chunk number 27: Tutorial.Rnw:372-377 ################################################### library(BioNet) library(DLBCL) library(ALL) data(ALL) data(interactome) ################################################### ### code chunk number 28: Tutorial.Rnw:384-385 ################################################### ALL ################################################### ### code chunk number 29: Tutorial.Rnw:393-394 ################################################### interactome ################################################### ### code chunk number 30: Tutorial.Rnw:404-406 ################################################### mapped.eset = mapByVar(ALL, network=interactome, attr="geneID") mapped.eset[1:5,1:5] ################################################### ### code chunk number 31: Tutorial.Rnw:411-412 ################################################### length(intersect(rownames(mapped.eset), nodes(interactome))) ################################################### ### code chunk number 32: Tutorial.Rnw:421-426 ################################################### network = subNetwork(rownames(mapped.eset), interactome) network network <- largestComp(network) network <- rmSelfLoops(network) network ################################################### ### code chunk number 33: Tutorial.Rnw:437-445 ################################################### library(limma) design= model.matrix(~ -1+ factor(c(substr(unlist(ALL$BT), 0, 1)))) colnames(design)= c("B", "T") contrast.matrix <- makeContrasts(B-T, levels=design) contrast.matrix fit <- lmFit(mapped.eset, design) fit2 <- contrasts.fit(fit, contrast.matrix) fit2 <- eBayes(fit2) ################################################### ### code chunk number 34: Tutorial.Rnw:450-451 ################################################### pval = fit2$p.value[,1] ################################################### ### code chunk number 35: Tutorial.Rnw:462-464 ################################################### fb <- fitBumModel(pval, plot=FALSE) fb ################################################### ### code chunk number 36: Tutorial.Rnw:467-471 ################################################### dev.new(width=13, height=7) par(mfrow=c(1,2)) hist(fb) plot(fb) ################################################### ### code chunk number 37: Tutorial.Rnw:483-484 ################################################### scores <- scoreNodes(network=network, fb=fb, fdr=1e-14) ################################################### ### code chunk number 38: Tutorial.Rnw:490-492 ################################################### writeHeinzEdges(network=network, file="ALL_edges_001", use.score=FALSE) writeHeinzNodes(network=network, file="ALL_nodes_001", node.scores = scores) ################################################### ### code chunk number 39: Tutorial.Rnw:511-513 ################################################### datadir <- file.path(.path.package("BioNet"), "extdata") module <- readHeinzGraph(node.file=file.path(datadir, "ALL_nodes_001.txt.0.hnz"), network=network) ################################################### ### code chunk number 40: Tutorial.Rnw:518-523 ################################################### nodeDataDefaults(module, attr="diff") <- "" nodeData(module, n=nodes(module), attr="diff") <- fit2$coefficients[nodes(module),1] nodeDataDefaults(module, attr="score") <- "" nodeData(module, n=nodes(module), attr="score") <- scores[nodes(module)] nodeData(module)[1] ################################################### ### code chunk number 41: Tutorial.Rnw:529-530 ################################################### saveNetwork(module, file="ALL_module", type="XGMML")