mbamv:

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

Usage

mbamv(x, y, nsamp = 7)

Arguments

x
y
nsamp

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, nsamp = 7) 
{
    x <- as.matrix(x)
    n <- dim(x)[1]
    q <- dim(x)[2]
    covv <- diag(q)
    centers <- sample(n, nsamp)
    for (i in 1:nsamp) {
        md2 <- mahalanobis(x, center = x[centers[i], ], covv)
        med <- median(md2)
        plot(x, y)
        points(x[md2 < med], y[md2 < med], pch = 15)
        abline(lsfit(x[md2 < med], y[md2 < med]))
        identify(x, y)
    }
  }

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