R/plotNode.R

In ggdendro: Create Dendrograms and Tree Diagrams Using 'ggplot2'

# modified version of stats:::plotNode
plotNode <-  function (x1, x2, subtree, type, center, leaflab, dLeaf, nodePar, 
            edgePar, horiz = FALSE) 
  {
    ddsegments <- NULL
    ddlabels <- list()
    
    wholetree <- subtree
    depth <- 0L
    llimit <- list()
    KK <- integer()
    kk <- integer()
    repeat {
      inner <- !is.leaf(subtree) && x1 != x2
      yTop <- attr(subtree, "height")
      bx <- plotNodeLimit(x1, x2, subtree, center)
      xTop <- bx$x
      depth <- depth + 1L
      llimit[[depth]] <- bx$limit
      hasP <- !is.null(nPar <- attr(subtree, "nodePar"))
      if (!hasP) 
        nPar <- nodePar
      if (getOption("verbose")) {
        cat(if (inner) 
          "inner node"
          else "leaf", ":")
        if (!is.null(nPar)) {
          cat(" with node pars\n")
          utils::str(nPar)
        }
        cat(if (inner) 
          paste(" height", formatC(yTop), "; "), "(x1,x2)= (", 
          formatC(x1, width = 4), ",", formatC(x2, width = 4), 
          ")", "--> xTop=", formatC(xTop, width = 8), "\n", 
          sep = "")
      }
      Xtract <- function(nam, L, default, indx) 
        rep(if (nam %in%  names(L)) L[[nam]] else default, length.out = indx)[indx]
      asTxt <- function(x) if (is.character(x) || is.expression(x) ||  is.null(x))  x else as.character(x)
      i <- if (inner || hasP)  1 else 2
      if (!is.null(nPar)) {
        pch <- Xtract("pch", nPar, default = 1L:2, i)
        cex <- Xtract("cex", nPar, default = c(1, 1), i)
        col <- Xtract("col", nPar, default = par("col"), 
                      i)
        bg <- Xtract("bg", nPar, default = par("bg"), i)
        points(if (horiz) 
          cbind(yTop, xTop)
          else cbind(xTop, yTop), pch = pch, bg = bg, col = col, 
          cex = cex)
      }
      if (leaflab == "textlike") 
        p.col <- Xtract("p.col", nPar, default = "white",  i)
      lab.col <- Xtract("lab.col", nPar, default = par("col"),  i)
      lab.cex <- Xtract("lab.cex", nPar, default = c(1, 1),  i)
      lab.font <- Xtract("lab.font", nPar, default = par("font"),  i)
      lab.xpd <- Xtract("xpd", nPar, default = c(TRUE, TRUE),  i)
      if (is.leaf(subtree)) {
        if (leaflab == "perpendicular") {
          if (horiz) {
            X <- yTop + dLeaf * lab.cex
            Y <- xTop
            srt <- 0
            adj <- c(0, 0.5)
          }
          else {
            Y <- yTop - dLeaf * lab.cex
            X <- xTop
            srt <- 90
            adj <- 1
          }
          nodeText <- asTxt(attr(subtree, "label"))
          # text(X, Y, nodeText, xpd = lab.xpd, srt = srt, 
          #      adj = adj, cex = lab.cex, col = lab.col, font = lab.font)
          ddlabels$xy <- c(ddlabels$xy, X, 0)
          ddlabels$text <- c(ddlabels$text, nodeText)
        }
      }
      else if (inner) {
        # segmentsHV <- function(x0, y0, x1, y1) {
        #   if (horiz) 
        #     segments(y0, x0, y1, x1, col = col, lty = lty,  lwd = lwd)
        #   else segments(x0, y0, x1, y1, col = col, lty = lty,  lwd = lwd)
        # }
        for (k in seq_along(subtree)) {
          child <- subtree[[k]]
          yBot <- attr(child, "height")
          if (getOption("verbose")) 
            cat("ch.", k, "@ h=", yBot, "; ")
          if (is.null(yBot)) 
            yBot <- 0
          xBot <- if (center) 
            mean(bx$limit[k:(k + 1)])
          else bx$limit[k] + .midDend(child)
          hasE <- !is.null(ePar <- attr(child, "edgePar"))
          if (!hasE) 
            ePar <- edgePar
          i <- if (!is.leaf(child) || hasE)  1 else 2
          col <- Xtract("col", ePar, default = par("col"),  i)
          lty <- Xtract("lty", ePar, default = par("lty"),  i)
          lwd <- Xtract("lwd", ePar, default = par("lwd"),  i)
          if (type == "triangle") {
            # segmentsHV(xTop, yTop, xBot, yBot)
            # segmentsHV(xTop, yTop, xBot, yBot)
            # ddsegments <- rbind(ddsegments, data.frame(x = xTop, y = yTop, xend = xBot, yend = yBot))
            ddsegments <- c(ddsegments, xTop, yTop, xBot, yBot)
          }
          else {
            # segmentsHV(xTop, yTop, xBot, yTop)
            # segmentsHV(xBot, yTop, xBot, yBot)
            # ddsegments <- rbind(ddsegments, data.frame(x = xTop, y = yTop, xend = xBot, yend = yTop))
            # ddsegments <- rbind(ddsegments, data.frame(x = xBot, y = yTop, xend = xBot, yend = yBot))
            ddsegments <- c(ddsegments, xTop, yTop, xBot, yTop)
            ddsegments <- c(ddsegments, xBot, yTop, xBot, yBot)
          }
          vln <- NULL
          # if (is.leaf(child) && leaflab == "textlike") {
          #   nodeText <- asTxt(attr(child, "label"))
          #   if (getOption("verbose"))
          #     cat("-- with \"label\"", format(nodeText))
          #   hln <- 0.6 * strwidth(nodeText, cex = lab.cex)/2
          #   vln <- 1.5 * strheight(nodeText, cex = lab.cex)/2
          #   rect(xBot - hln, yBot, xBot + hln, yBot + 2 *
          #          vln, col = p.col)
          #   text(xBot, yBot + vln, nodeText, xpd = lab.xpd,
          #        cex = lab.cex, col = lab.col, font = lab.font)
          # }
          # if (!is.null(attr(child, "edgetext"))) {
          #   edgeText <- asTxt(attr(child, "edgetext"))
          #   if (getOption("verbose"))
          #     cat("-- with \"edgetext\"", format(edgeText))
          #   if (!is.null(vln)) {
          #     mx <- if (type == "triangle")
          #       (xTop + xBot + ((xTop - xBot)/(yTop - yBot)) *
          #          vln)/2
          #     else xBot
          #     my <- (yTop + yBot + 2 * vln)/2
          #   }
          #   else {
          #     mx <- if (type == "triangle")
          #       (xTop + xBot)/2
          #     else xBot
          #     my <- (yTop + yBot)/2
          #   }
          #   p.col <- Xtract("p.col", ePar, default = "white",  i)
          #   p.border <- Xtract("p.border", ePar, default = par("fg"),  i)
          #   p.lwd <- Xtract("p.lwd", ePar, default = lwd,  i)
          #   p.lty <- Xtract("p.lty", ePar, default = lty,  i)
          #   t.col <- Xtract("t.col", ePar, default = col,  i)
          #   t.cex <- Xtract("t.cex", ePar, default = 1,  i)
          #   t.font <- Xtract("t.font", ePar, default = par("font"),  i)
          #   vlm <- strheight(c(edgeText, "h"), cex = t.cex)/2
          #   hlm <- strwidth(c(edgeText, "m"), cex = t.cex)/2
          #   hl3 <- c(hlm[1L], hlm[1L] + hlm[2L], hlm[1L])
          #   if (horiz) {
          #     polygon(my + c(-hl3, hl3), mx + sum(vlm) *
          #               c(-1L:1L, 1L:-1L), col = p.col, border = p.border,
          #             lty = p.lty, lwd = p.lwd)
          #     text(my, mx, edgeText, cex = t.cex, col = t.col,
          #          font = t.font)
          #   }
          #   else {
          #     polygon(mx + c(-hl3, hl3), my + sum(vlm) *
          #               c(-1L:1L, 1L:-1L), col = p.col, border = p.border,
          #             lty = p.lty, lwd = p.lwd)
          #     text(mx, my, edgeText, cex = t.cex, col = t.col,
          #          font = t.font)
          #   }
          # }
        }
      }
      if (inner && length(subtree)) {
        KK[depth] <- length(subtree)
        if (storage.mode(kk) != storage.mode(KK)) 
          storage.mode(kk) <- storage.mode(KK)
        kk[depth] <- 1L
        x1 <- bx$limit[1L]
        x2 <- bx$limit[2L]
        subtree <- subtree[[1L]]
      }
      else {
        repeat {
          depth <- depth - 1L
          if (!depth || kk[depth] < KK[depth]) 
            break
        }
        if (!depth) 
          break
        length(kk) <- depth
        kk[depth] <- k <- kk[depth] + 1L
        x1 <- llimit[[depth]][k]
        x2 <- llimit[[depth]][k + 1L]
        subtree <- wholetree[[kk]]
      }
    }
    # invisible()
    list(
      segments = ddsegments,
      labels = ddlabels
    )
  }