wddplot:

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

Usage

wddplot(x)

Arguments

x

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 (x) 
{
    n <- dim(x)[1]
    wt <- 0 * (1:n)
    p <- dim(x)[2]
    center <- apply(x, 2, mean)
    cov <- var(x)
    md2 <- mahalanobis(x, center, cov)
    out <- cov.mcd(x)
    center <- out$center
    cov <- out$cov
    rd2 <- mahalanobis(x, center, cov)
    md <- sqrt(md2)
    rd <- sqrt(rd2)
    const <- sqrt(qchisq(0.5, p))/median(rd)
    rd <- const * rd
    wt[rd < sqrt(qchisq(0.975, p))] <- 1
    MD <- md[wt > 0]
    RD <- rd[wt > 0]
    plot(MD, RD)
  }

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