mbamv2:

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

Usage

mbamv2(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)
        if (q == 1) {
            out <- lsfit(x[md2 < med], Y[md2 < med])
        }
        else {
            out <- lsfit(x[md2 < med, ], Y[md2 < med])
        }
        FIT <- out$coef[1] + x %*% out$coef[-1]
        RES <- Y - FIT
        par(mfrow = c(2, 1))
        plot(FIT, Y)
        points(FIT[md2 < med], Y[md2 < med], pch = 15)
        abline(0, 1)
        identify(FIT, Y)
        plot(FIT, RES)
        points(FIT[md2 < med], RES[md2 < med], pch = 15)
        abline(0, 0)
        identify(FIT, RES)
    }
  }

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