1 |
x |
|
q |
|
alpha |
|
op |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | ##---- 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, q = 0.5, alpha = 0.05, op = 3)
{
if (sum(duplicated(x) > 0))
stop("Duplicate values detected; use hdpb")
n <- length(x)
xsort <- sort(x)
iq <- floor(q * n + 0.5)
qest <- xsort[iq]
se <- qse(x, q, op = op)
crit <- qnorm(1 - alpha/2)
ci.low <- qest - crit * se
ci.up <- qest + crit * se
list(ci.low = ci.low, ci.up = ci.up, q.est = qest)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.