1 |
x |
|
y |
|
est |
|
na.rm |
|
plotit |
|
xlab |
|
ylab |
|
... |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 | ##---- 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 = NULL, est = median, na.rm = TRUE, plotit = FALSE,
xlab = "", ylab = "", ...)
{
if (is.null(y)) {
if (ncol(x) != 2)
stop("x should be an n-by-2 matrix")
y = x[, 2]
x = x[, 1]
}
if (na.rm) {
m = elimna(cbind(x, y))
x = m[, 1]
y = m[, 2]
}
temp = as.vector(outer(x, y, FUN = "-"))
val <- est(temp, na.rm = TRUE, ...)
if (plotit)
akerd(temp, xlab = xlab, ylab = ylab)
val
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.