tvreg2:

Usage Arguments Examples

Usage

1
tvreg2(X, Y, M = 0)

Arguments

X
Y
M

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
##---- 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 = 0) 
{
    X <- as.matrix(X)
    out <- cov.mcd(X)
    center <- out$center
    cov <- out$cov
    rd2 <- mahalanobis(X, center, cov)
    tem <- (100 - M)/100
    val <- quantile(rd2, tem)
    b <- lsfit(X[rd2 <= val, ], Y[rd2 <= val])$coef
    FIT <- X %*% b[-1] + b[1]
    plot(FIT, Y)
    abline(0, 1)
    identify(FIT, Y)
    list(coef = b)
  }

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