1 | ratmn(x, k1 = 6, k2 = 6)
|
x |
|
k1 |
|
k2 |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 | ##---- 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, k1 = 6, k2 = 6)
{
madd <- mad(x, constant = 1)
med <- median(x)
LM <- sum(x < (med - k1 * madd))
nmUM <- sum(x > (med + k2 * madd))
n <- length(x)
ll <- ceiling((100 * LM)/n)
hh <- ceiling((100 * (nmUM))/n)
tem <- sort(x)
ln <- floor((ll * n)/100)
un <- floor((n * (100 - hh))/100)
low <- ln + 1
val1 <- tem[low]
val2 <- tem[un]
rtmn <- mean(x[(x >= val1) & (x <= val2)])
trmn
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.