1 |
x |
|
y |
|
est |
|
pts |
|
fr |
|
... |
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, est = tmean, pts, fr = 1, ...)
{
if (!is.matrix(x))
stop("Predictors are not stored in a matrix.")
if (!is.matrix(pts))
stop("The third argument, pts, must be a matrix.")
library(MASS)
m <- cov.mve(x)
rmd <- 1
nval <- 1
for (i in 1:nrow(pts)) {
rmd[i] <- est(y[near3d(x, pts[i, ], fr, m)], ...)
nval[i] <- length(y[near3d(x, pts[i, ], fr, m)])
}
list(rmd = rmd, nval = nval)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.