1 |
x |
|
beta |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | ##---- 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, beta = 0.2)
{
pbvar = 0
x = elimna(x)
w <- abs(x - median(x))
w <- sort(w)
m <- floor((1 - beta) * length(x) + 0.5)
omega <- w[m]
if (omega > 0) {
y <- (x - median(x))/omega
z <- ifelse(y > 1, 1, y)
z <- ifelse(z < (-1), -1, z)
pbvar <- length(x) * omega^2 * sum(z^2)/(length(x[abs(y) <
1]))^2
}
pbvar
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.