R/plotQC.r

In FREddyPro: Post-Processing EddyPro Full Output File

plotQC <-
function (data, var = "co2_flux", qc_var = NULL, xvar = "DOY", 
    lines = FALSE, legendSide = NULL, ylab = var, xlab = "DOY", 
    col = c(1, 2, 3, 4, 5, 6, 7, 8, 9), pch = c(16, 17, 18, 19, 
        20, 21, 22, 23, 24), ...) 
{
    if (is.null(qc_var)) {
        qc = paste("qc_", var, sep = "")
    }
    else {
        qc = qc_var
    }
    plot(data[[var]] ~ data[[xvar]], col = ifelse(data[[qc]] == 
        0, col[1], ifelse(data[[qc]] == 1, col[2], ifelse(data[[qc]] == 
        2, col[3], ifelse(data[[qc]] == 3, col[4], ifelse(data[[qc]] == 
        4, col[5], ifelse(data[[qc]] == 5, col[6], ifelse(data[[qc]] == 
        6, col[7], ifelse(data[[qc]] == 7, col[8], ifelse(data[[qc]] == 
        8, col[9], ifelse(data[[qc]] == 9, col[10], 1)))))))))), 
        pch = ifelse(data[[qc]] == 0, pch[1], ifelse(data[[qc]] == 
            1, pch[2], ifelse(data[[qc]] == 2, pch[3], ifelse(data[[qc]] == 
            3, pch[4], ifelse(data[[qc]] == 4, pch[5], ifelse(data[[qc]] == 
            5, pch[6], ifelse(data[[qc]] == 6, pch[7], ifelse(data[[qc]] == 
            7, pch[8], ifelse(data[[qc]] == 8, pch[9], ifelse(data[[qc]] == 
            9, pch[10], 1)))))))))), ylab = ylab, xlab = xlab, 
        ...)
    if (lines) {
        lines(data[[var]] ~ data[[xvar]], ...)
    }
    if (!is.null(legendSide)) {
        l = list(`0` = "QC=0", `1` = "QC=1", `2` = "QC=2", `3` = "QC=3", 
            `4` = "QC=4", `5` = "QC=5", `6` = "QC=6", `7` = "QC=7", 
            `8` = "QC=8", `9` = "QC=9")
        m = as.vector(sort(unique(data[[qc]])))
        txt = c(do.call("cbind", l[as.character(m)]))
        legend(legendSide, txt, col = c(col[1], col[2], col[3], 
            col[4], col[5], col[6], col[7], col[8], col[9]), 
            pch = c(pch[1], pch[2], pch[3], pch[4], pch[5], pch[6], 
                pch[7], pch[8], pch[9]), horiz = TRUE, bty = "n")
    }
}