deav:

Usage Arguments Examples

Usage

1
deav(alpha = 0.01, k = 5)

Arguments

alpha
k

Examples

 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)
  }

musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.