1 | bivar(x)
|
x |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | ##---- 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)
{
m <- median(x)
u <- abs((x - m)/(9 * qnorm(0.75) * mad(x)))
av <- ifelse(u < 1, 1, 0)
top <- length(x) * sum(av * (x - m)^2 * (1 - u^2)^4)
bot <- sum(av * (1 - u^2) * (1 - 5 * u^2))
bi <- top/bot^2
bi
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.