1 |
J |
|
K |
|
x |
|
est |
|
alpha |
|
grp |
|
dif |
|
nboot |
|
plotit |
|
BA |
|
hoch |
|
... |
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 26 27 28 29 30 31 32 33 34 35 36 37 | ##---- 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 (J, K, x, est = tmean, alpha = 0.05, grp = NA, dif = TRUE,
nboot = NA, plotit = FALSE, BA = FALSE, hoch = FALSE, ...)
{
JK <- J * K
if (is.matrix(x))
x <- listm(x)
if (!is.na(grp[1])) {
yy <- x
for (j in 1:length(grp)) x[[j]] <- yy[[grp[j]]]
}
if (!is.list(x))
stop("Data must be stored in list mode or a matrix.")
for (j in 1:JK) {
xx <- x[[j]]
x[[j]] <- xx[!is.na(xx)]
}
temp <- con2way(J, K)
conA <- temp$conA
conB <- temp$conB
conAB <- temp$conAB
ncon <- max(nrow(conA), nrow(conB), nrow(conAB))
FacA <- rmmcppb(x, con = conA, est = est, plotit = plotit,
dif = dif, grp = grp, nboot = nboot, BA = TRUE, hoch = FALSE,
...)
FacB <- rmmcppb(x, con = conB, est = est, plotit = plotit,
dif = dif, grp = grp, nboot = nboot, BA = TRUE, hoch = FALSE,
...)
FacAB <- rmmcppb(x, con = conAB, est = est, plotit = plotit,
dif = dif, grp = grp, nboot = nboot, BA = TRUE, hoch = FALSE,
...)
list(Factor.A = FacA, Factor.B = FacB, Factor.AB = FacAB)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.