| ptte.boot {equivalence} | R Documentation |
~~ A concise (1-5 lines) description of what the function does. ~~
ptte.boot(x, i, Epsilon = 0.5)
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=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)
}