1 |
x |
|
y |
|
est |
|
nboot |
|
pr |
|
... |
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 = NA, est = median, nboot = 100, pr = TRUE, ...)
{
if (!is.na(y[1]))
x <- cbind(x, y)
x <- elimna(x)
if (pr)
print("Taking Bootstrap Samples. Please wait.")
data <- matrix(sample(nrow(x), size = nrow(x) * nboot, replace = TRUE),
nrow = nboot)
xmat <- matrix(x[data, 1], nrow = nboot, ncol = length(x))
ymat <- matrix(x[data, 2], nrow = nboot, ncol = length(x))
bvec <- apply(xmat, 1, FUN = est, ...) - apply(ymat, 1, FUN = est,
...)
se <- sqrt(var(bvec))
se
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.