ptte.boot {equivalence}R Documentation

~~function to do ... ~~

Description

~~ A concise (1-5 lines) description of what the function does. ~~

Usage

ptte.boot(x, i, Epsilon = 0.5)

Arguments

x ~~Describe x here~~
i ~~Describe i here~~
Epsilon ~~Describe Epsilon here~~

Details

~~ If necessary, more details than the description above ~~

Value

~Describe the value returned If it is a LIST, use

comp1 Description of 'comp1'
comp2 Description of 'comp2'

...

Warning

....

Note

~~further notes~~

~Make other sections like Warning with section{Warning }{....} ~

Author(s)

~~who you are~~

References

~put references to the literature/web site here ~

See Also

~~objects to See Also as help, ~~~

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, i, Epsilon=0.5)
{
    x <- x[i]
    mean <- mean(x, na.rm = TRUE)
    std <- sd(x, na.rm = TRUE)
    n = sum(!is.na(x))

    df1 <- 1
    df2 <- n - 1
    cutoff <- c(0)
    ncp <- Epsilon^2 * n

    findquant <- function(alpha, q, df1, df2, ncp) alpha - pf(q,
        df1, df2, ncp)
    cutoff <- sqrt(uniroot(findquant, c(0, 1e+23), alpha = 0.05,
        df1 = df1, df2 = min(df2, 4000), ncp = ncp)$root)
    tstat <- mean/(std/sqrt(n))
    result <- (abs(tstat) < cutoff)
  }

[Package equivalence version 0.4.1 Index]