R/plotFit.R

In 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), ...)
}