manly.boot {equivalence}R Documentation

~~function to do ... ~~

Description

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

Usage

manly.boot(x, i, Epsilon = 1)

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=1)
{
    reps <- 50
    x <- x[i]
    mean <- mean(x)
    n <- length(x)

    pmanly1.boot <- function(d, i)
    {  randomize <- 1 - 2*rbinom(n, 1, 0.5)
       d1.boot <- (d[i]+Epsilon)*randomize
       return(mean(d1.boot) >= mean+Epsilon) }
    pmanly2.boot <- function(d, i)
    {  randomize <- 1 - 2*rbinom(n, 1, 0.5)
       d2.boot <- (d[i]-Epsilon)*randomize
       return(mean(d2.boot) <= mean-Epsilon) }

    pvalue <- mean(boot(x, pmanly1.boot, reps)$t) +
              mean(boot(x, pmanly2.boot, reps)$t)

    result <- (pvalue < 0.05)
  }

[Package equivalence version 0.4.1 Index]