nodelabels.sp <-
function (text, node, adj = c(0.5, 0.5), frame = "rect", pch = NULL, thermo = NULL, pie = NULL, piecol = NULL, col = "black", bg = "lightblue", border = "black", ...) {
    lastPP <- get("last_plot.phylo", envir = .PlotPhyloEnv)
    if (missing(node)) 
	node <- (lastPP$Ntip + 1):length(lastPP$xx)
    XX <- lastPP$xx[node]
    YY <- lastPP$yy[node]
    BOTHlabels.sp(text, node, XX, YY, adj, frame, pch, thermo, pie, piecol, col, bg, border, ...)
}

BOTHlabels.sp <-
function (text, sel, XX, YY, adj, frame, pch, thermo, pie, piecol, col, bg, border, ...) {
    if (missing(text)) 
	text <- NULL
    if (length(adj) == 1) 
	adj <- c(adj, 0.5)
    if (is.null(text) && is.null(pch) && is.null(thermo) && is.null(pie)) 
	text <- as.character(sel)
    frame <- match.arg(frame, c("rect", "circle", "none"))
    args <- list(...)
    CEX <- if ("cex" %in% names(args)) 
	args$cex
    else par("cex")
    if (frame != "none" && !is.null(text)) {
        if (frame == "rect") {
            width <- strwidth(text, units = "inches", cex = CEX)
            height <- strheight(text, units = "inches", cex = CEX)
            if ("srt" %in% names(args)) {
                args$srt <- args$srt%%360
                if (args$srt == 90 || args$srt == 270) {
					tmp <- width
					width <- height
					height <- tmp
                }
                else if (args$srt != 0) 
				warning("only right angle rotation of frame is supported;\n         try  `frame = \"n\"' instead.\n")
            }
            width <- xinch(width)
            height <- yinch(height)
            xl <- XX - xinch(0.10)
            xr <- XX + xinch(-0.01)
            yb <- YY - yinch(0.02)
            yt <- YY + yinch(0.02)
            rect(xl, yb, xr, yt, col = bg, border = border)
        }
        if (frame == "circle") {
            radii <- 0.8 * apply(cbind(strheight(text, units = "inches", 
												 cex = CEX), strwidth(text, units = "inches", 
																	  cex = CEX)), 1, max)
            symbols(XX, YY, circles = radii, inches = max(radii), add = TRUE, bg = bg)
        }
    }
    if (!is.null(thermo)) {
        parusr <- par("usr")
        width <- CEX * (parusr[2] - parusr[1])/40
        height <- CEX * (parusr[4] - parusr[3])/15
        if (is.vector(thermo)) 
		thermo <- cbind(thermo, 1 - thermo)
        thermo <- height * thermo
        xl <- XX - width/2
        xr <- xl + width
        yb <- YY - height/2
        if (is.null(piecol)) 
		piecol <- rainbow(ncol(thermo))
        rect(xl, yb, xr, yb + thermo[, 1], border = NA, col = piecol[1])
        for (i in 2:ncol(thermo)) rect(xl, yb + rowSums(thermo[, 
														1:(i - 1), drop = FALSE]), xr, yb + rowSums(thermo[, 
																									1:i]), border = NA, col = piecol[i])
        rect(xl, yb, xr, yb + height, border = "black")
        segments(xl, YY, xl - width/5, YY)
        segments(xr, YY, xr + width/5, YY)
    }
    if (!is.null(pie)) {
        if (is.vector(pie)) 
		pie <- cbind(pie, 1 - pie)
        xrad <- CEX * diff(par("usr")[1:2])/50
        xrad <- rep(xrad, length(sel))
        for (i in 1:length(sel)) floating.pie.asp(XX[i], YY[i], 
												  pie[i, ], radius = xrad[i], col = piecol)
    }
    if (!is.null(text)) 
	text(XX, YY, text, adj = adj, col = col, ...)
    if (!is.null(pch)) 
	points(XX + adj[1] - 0.5, YY + adj[2] - 0.5, pch = pch, 
		   col = col, bg = bg, ...)
}