lband.fun:

In musto101/wilcox_R: What the package does (short line)

Usage

lband.fun(x, y, crit)

Arguments

x
y
crit

Examples

##---- 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, y, crit) 
{
    pi <- gamma(0.5)^2
    xr <- rank(x)
    yr <- rank(y)
    temp <- apply(cbind(xr, yr), 1, max)
    n <- length(x)
    fj <- NA
    for (i in 1:n) fj[i] <- sum(temp == i)
    v1 <- NA
    for (j in 1:n) v1[j] <- (j - sum(fj[1:j]))/n
    psi <- rep(0, n)
    for (j in 1:n) {
        if (v1[j] > 0) 
            psi[j] <- crit * exp(0 - crit^2/(2 * v1[j]))/sqrt(2 * 
                pi * v1[j]^3)
    }
    res <- mean(fj * psi)
    res
  }

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