1 |
x |
|
y |
|
est |
|
... |
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, est = tukmed, ...)
{
x <- as.matrix(x)
y <- as.matrix(y)
FLAG <- F
if (ncol(x) != ncol(y))
stop("x and y should have the same number of columns")
if (ncol(x) == 1 && ncol(y) == 1)
FLAG <- T
if (FLAG)
val <- loc2dif(x, y, est = est, ...)
if (!FLAG) {
J <- (ncol(x)^2 - ncol(x))/2
mat <- matrix(NA, ncol = ncol(x), nrow = nrow(x) * nrow(y))
for (j in 1:ncol(x)) mat[, j] <- as.vector(outer(x[,
j], y[, j], FUN = "-"))
val <- est(mat, ...)
}
val
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.