Nothing
p.int <- function (A, M, delta = 50, N = -1, av = "median", p.adjust.method = "none")
{
if (N < 0) {
N <- 100 * length(A)
}
MavP <- double(N)
#### GERNERATING EMPIRICAL DISTRIBUTION
if (av == "mean") {
for (i in 1:N) {
MavP[i] <- mean(sample(M, 2 * delta + 1))
}
}
if (av == "median") {
for (i in 1:N) {
MavP[i] <- median(sample(M, 2 * delta + 1))
}
}
#### STATISTIC FOR ORIGINAL DATA
Mav <- ma.vector(A, M, delta = delta, av = av)
Mav.l <- length(Mav)
pP <- double(length = length(A)) + NA
MavP.l <- length(MavP)
#### DETERMINING P-VALUES
for (i in 1:Mav.l) {
pP[i] <- sum(MavP >= Mav[i], na.rm = TRUE)/MavP.l
}
nP <- double(length = length(M)) + NA
for (i in 1:Mav.l) {
nP[i] <- sum(MavP <= Mav[i], na.rm = TRUE)/MavP.l
}
### ADJUSTMENT
pP[pP == 0] <- 1/N
nP[nP == 0] <- 1/N
pP.adjust <- p.adjust(pP, method = p.adjust.method)
nP.adjust <- p.adjust(nP, method = p.adjust.method)
pP.adjust[is.na(Mav)] <- NA
nP.adjust[is.na(Mav)] <- NA
list(Pp = pP.adjust, Pn = nP.adjust)
}
#######################################################################
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.