R/thplot2.R

In rgr: Applied Geochemistry EDA

thplot2 <-
function (x, xname = deparse(substitute(x)), ifzero = 0.01, xlow = NA, 
    xhih = NA, yhih = NA, rsd = 5, ptile = 95, main = "", ifalt = FALSE, 
    ...) 
{
    n <- length(x)
    ndup <- n/2
    x1 <- numeric(ndup)
    x2 <- numeric(ndup)
    if (ifalt) {
        for (i in 1:ndup) {
            j <- 2 * (i - 1) + 1
            x1[i] <- x[j]
            x2[i] <- x[j + 1]
        }
    }
    else {
        for (i in 1:ndup) {
            x1[i] <- x[i]
            x2[i] <- x[ndup + i]
        }
    }
    thplot1(x1, x2, xname = xname, ifzero = ifzero, xlow = xlow, 
        xhih = xhih, yhih = yhih, rsd = rsd, ptile = ptile, main = main, 
        ...)
    invisible()
}