| manly.boot {equivalence} | R Documentation |
~~ A concise (1-5 lines) description of what the function does. ~~
manly.boot(x, i, Epsilon = 1)
x |
~~Describe x here~~ |
i |
~~Describe i here~~ |
Epsilon |
~~Describe Epsilon here~~ |
~~ If necessary, more details than the description above ~~
~Describe the value returned If it is a LIST, use
comp1 |
Description of 'comp1' |
comp2 |
Description of 'comp2' |
...
....
~~further notes~~
~Make other sections like Warning with section{Warning }{....} ~
~~who you are~~
~put references to the literature/web site here ~
~~objects to See Also as help, ~~~
##---- 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)
}