1 | deav(alpha = 0.01, k = 5)
|
alpha |
|
k |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | ##---- 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 (alpha = 0.01, k = 5)
{
z <- -log(2 * alpha)
num <- 2 - (2 + 2 * z + z^2) * exp(-z)
den <- (1 - exp(-z)) * (1 - 2 * alpha)
val1 <- num/den
num <- 2 * alpha * z^2
den <- (1 - 2 * alpha)^2
ntmav <- val1 + num/den
zj <- k * log(2)
alphaj <- 0.5 * exp(-zj)
alphaj <- ceiling(100 * alphaj)/100
zj <- -log(2 * alphaj)
num <- 2 - (2 + 2 * zj + zj^2) * exp(-zj)
den <- (1 - exp(-zj)) * (1 - 2 * alphaj)
val1 <- num/den
num <- 2 * alphaj * zj^2
den <- (1 - 2 * alphaj)^2
natmav <- val1 + num/den
return(ntmav, natmav)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.