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) }