essp:

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

Usage

essp(x, Y, M = 50)

Arguments

x
Y
M

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, Y, M = 50) 
{
    x <- as.matrix(x)
    tval <- M/100
    out <- cov.mcd(x)
    center <- out$center
    cov <- out$cov
    rd2 <- mahalanobis(x, center, cov)
    val <- quantile(rd2, (1 - tval))
    bhat <- lsfit(x[rd2 <= val, ], Y[rd2 <= val])$coef
    ESP <- x %*% bhat[-1]
    plot(ESP, Y)
    identify(ESP, Y)
    return(bhat[-1])
  }

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