renderSubNetSimple(NodeName, sampleName, GeneName, intome = intome, gisticCopyCalls = NULL, resultObj, fileOut = NULL)
Arguments
| NodeName |
|
| sampleName |
|
| GeneName |
|
| intome |
|
| gisticCopyCalls |
|
| resultObj |
|
| fileOut |
|
Details
Value
References
Note
See also
Examples
#> function (NodeName, sampleName, GeneName, intome = intome, gisticCopyCalls = NULL,
#> resultObj, fileOut = NULL)
#> {
#> library(Rgraphviz)
#> library(BioNet)
#> if (!is.null(gisticCopyCalls)) {
#> cellFrame <- gisticCopyCalls[, c("Gene.Symbol", "NodeName",
#> sampleName)]
#> }
#> mutCopyFrame <- resultObj$mutCopyFrames[[sampleName]]
#> cellResults <- resultObj$cellResults
#> nodenet <- c(NodeName, intersect(as.character(inEdges(NodeName,
#> intome)[[1]]), as.character(mutCopyFrame$NodeName)))
#> brcamut <- intersect(nodenet, as.character(mutCopyFrame[mutCopyFrame$mutations !=
#> 0, "NodeName"]))
#> labeltab <- geneIntTable[geneIntTable$NodeName %in% nodenet,
#> ]
#> labels <- as.character(labeltab$Gene)
#> names(labels) <- labeltab$NodeName
#> nodenetwork <- subNetwork(nodenet, intome)
#> nodeRenderInfo(nodenetwork) <- list(shape = "circle", iheight = 0.5,
#> iwidth = 0.5, fixedsize = FALSE, label = as.list(labels))
#> RPPAshape <- list()
#> RPPAshape[[NodeName]] <- "diamond"
#> nodeRenderInfo(nodenetwork) <- list(shape = RPPAshape)
#> if (!is.null(gisticCopyCalls)) {
#> nodecopyframe <- cellFrame[as.character(cellFrame$NodeName) %in%
#> nodenet, ]
#> rownames(nodecopyframe) <- nodecopyframe$NodeName
#> colnames(nodecopyframe) <- c("GeneName", "NodeName",
#> "copyStatus")
#> }
#> else {
#> nodecopyframe <- mutCopyFrame[, c("GeneName", "NodeName",
#> "copyStatus")]
#> }
#> nodecolors <- list()
#> for (nd in rownames(nodecopyframe)) {
#> ndcol <- "white"
#> if (nodecopyframe[nd, "copyStatus"] > 0) {
#> ndcol <- "pink"
#> }
#> if (nodecopyframe[nd, "copyStatus"] < 0) {
#> ndcol <- "lightgreen"
#> }
#> nodecolors[[nd]] <- ndcol
#> }
#> for (nd in brcamut) {
#> nodecolors[[nd]] <- "lightblue"
#> }
#> cellRes <- cellResults[[sampleName]]
#> NodDegree <- cellRes[cellRes$NodeName == NodeName, "degree"]
#> NeighborMuts <- cellRes[cellRes$NodeName == NodeName, "neighborVec"]
#> pval <- cellRes[cellRes$NodeName == NodeName, "pvalue"]
#> nodenetwork <- layoutGraph(nodenetwork)
#> if (NeighborMuts > 0 & length(nodes(nodenetwork)) > 1) {
#> nodeRenderInfo(nodenetwork) <- list(fill = nodecolors)
#> if (is.null(fileOut)) {
#> fileOut <- paste(sampleName, "-", GeneName, ".svg",
#> sep = "")
#> }
#> svg(height = 7, width = 7, filename = fileOut)
#> renderGraph(nodenetwork)
#> plotTitle <- paste(NodeName, " (", sampleName, " ", ", p=",
#> pval, ", n=", NeighborMuts, ", d=", NodDegree, ")",
#> sep = "")
#> title(plotTitle)
#> dev.off()
#> }
#> }
#> <environment: 0x10bde39c0>