plot_robpca:

In musto101/wilcox_R: What the package does (short line)

Usage

plot_robpca(robpca.obj, classic = 0, labod = 3, labsd = 3)

Arguments

robpca.obj
classic
labod
labsd

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (robpca.obj, classic = 0, labod = 3, labsd = 3) 
{
    diagnosticplot <- !any(robpca.obj$od <= as.vector(1e-06, 
        mode(robpca.obj$od)))
    if (diagnosticplot == T) {
        xmax <- max(max(robpca.obj$sd), robpca.obj$cutoff$sd)
        ymax <- max(max(robpca.obj$od), robpca.obj$cutoff$od)
        plot(robpca.obj$sd, robpca.obj$od, xlab = "Score distance", 
            ylab = "Orthogonal distance", xlim = c(0, xmax), 
            ylim = c(0, ymax), type = "p")
        abline(v = robpca.obj$cutoff$sd)
        abline(h = robpca.obj$cutoff$od)
        givelabel(robpca.obj, labod, labsd)
    }
    else {
        ymax <- max(max(robpca.obj$sd), robpca.obj$cutoff$sd)
        plot(robpca.obj$sd, xlab = "Index", ylab = "Score distance", 
            ylim = c(0, ymax), type = "p")
        abline(h = robpca.obj$cutoff$sd)
        givelabel(robpca.obj, labod = 0, labsd, indexplot = 1)
    }
    title("ROBPCA")
    if (classic == 1) {
        diagnosticplot <- !any(robpca.obj$classic$od <= as.vector(1e-06, 
            mode(robpca.obj$classic$od)))
        if (diagnosticplot == T) {
            xmax <- max(max(robpca.obj$classic$sd), robpca.obj$classic$cutoff$sd)
            ymax <- max(max(robpca.obj$classic$od), robpca.obj$classic$cutoff$od)
            plot(robpca.obj$classic$sd, robpca.obj$classic$od, 
                xlab = "Score distance", ylab = "Orthogonal distance", 
                xlim = c(0, xmax), ylim = c(0, ymax), type = "p")
            abline(v = robpca.obj$classic$cutoff$sd)
            abline(h = robpca.obj$classic$cutoff$od)
            givelabel(robpca.obj$classic, labod, labsd)
        }
        else {
            ymax <- max(max(robpca.obj$classic$sd), robpca.obj$classic$cutoff$sd)
            plot(robpca.obj$classic$sd, xlab = "Index", ylab = "Score distance", 
                ylim = c(0, ymax), type = "p")
            abline(h = robpca.obj$classic$cutoff$sd)
            givelabel(robpca.obj$classic, labod = 0, labsd, indexplot = 1)
        }
        title("CPCA")
    }
    invisible(robpca.obj)
  }

musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.