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

##---- Should be DIRECTLY executable !! ---- ##-- ==> Define data, use random, ##-- or do help(data=index) for the standard data sets. ## The function is currently defined as 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() } }
#> 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>