1 |
xy |
|
tr |
|
alpha |
|
plotit |
|
plot.dif |
|
pts |
|
sm |
|
pr |
|
xout |
|
outfun |
|
LP |
|
npts |
|
p.crit |
|
nreps |
|
SCAT |
|
xlab |
|
ylab |
|
pch1 |
|
pch2 |
|
... |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 | ##---- 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 (xy, tr = 0.2, alpha = 0.05, plotit = FALSE, plot.dif = FALSE,
pts = NA, sm = FALSE, pr = TRUE, xout = FALSE, outfun = out,
LP = TRUE, npts = 25, p.crit = NULL, nreps = 2000, SCAT = TRUE,
xlab = "X", ylab = "Y", pch1 = "*", pch2 = "+", ...)
{
xy1 = elimna(xy[, 1:2])
xy2 = elimna(xy[, 3:4])
x1 = xy1[, 1]
y1 = xy1[, 2]
x2 = xy2[, 1]
y2 = xy2[, 2]
if (ncol(as.matrix(x1)) > 1)
stop("One covariate only is allowed with this function")
res1 = ancova(x1, y1, x2, y2, pr = FALSE, plotit = FALSE)$output
pts = seq(res1[1, 1], res1[5, 1], length.out = npts)
res = ancova(x1, y1, x2, y2, tr = tr, alpha = alpha, plotit = FALSE,
pr = FALSE, pts = pts, skip.crit = TRUE)$output
res.out = min(res[, 9])
res.out
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.