R/plotFit.R

In mccallm/ternarynet: Ternary Network Estimation

plotFit <- function(ternaryFit, type = "interactive", ...) {
  grph <- graphObjMin(ternaryFit)
  deg <- apply(grph, 2, function(x) sum(x > 0))
  el <- cbind(as.vector(grph)[as.vector(grph) > 0], rep(1:ncol(grph), deg))
  el2 <- cbind(geneNames(ternaryFit)[el[, 1]], geneNames(ternaryFit)[el[, 2]])
  if (type == "static") plot(graph.edgelist(el2), ...)
  if (type == "interactive") tkplot(graph.edgelist(el2), ...)
}