1 |
x |
|
y |
|
est |
|
fr |
|
plotit |
|
pyhat |
|
nboot |
|
atr |
|
nmin |
|
outfun |
|
eout |
|
xlab |
|
ylab |
|
scat |
|
SEED |
|
expand |
|
scale |
|
STAND |
|
varfun |
|
pr |
|
ticktype |
|
theta |
|
phi |
|
... |
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 38 39 40 41 42 43 44 45 46 | ##---- 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, y, est = tmean, fr = 1, plotit = TRUE, pyhat = FALSE,
nboot = 40, atr = 0, nmin = 0, outfun = outpro, eout = FALSE,
xlab = "X", ylab = "Y", scat = TRUE, SEED = TRUE, expand = 0.5,
scale = FALSE, STAND = TRUE, varfun = pbvar, pr = TRUE, ticktype = "simple",
theta = 50, phi = 25, ...)
{
x <- as.matrix(x)
xx <- cbind(x, y)
xx <- elimna(xx)
x <- xx[, 1:ncol(x)]
x <- as.matrix(x)
y <- xx[, ncol(x) + 1]
if (ncol(x) == 1) {
val <- runmbo(x, y, est = est, scat = scat, fr = fr,
plotit = plotit, pyhat = TRUE, STAND = STAND, xlab = xlab,
ylab = ylab, eout = eout, nboot = nboot, outfun = outfun,
SEED = SEED, atr = atr, ...)
}
if (ncol(x) > 1) {
if (ncol(x) == 2 && !scale) {
if (pr) {
print("scale=F is specified.")
print("If there is dependence, use scale=T")
}
}
if (ncol(x) > 2)
plotit <- F
val <- run3bo(x, y, est = est, fr = fr, nmin = nmin,
plotit = plotit, pyhat = TRUE, phi = phi, theta = theta,
xlab = xlab, ylab = ylab, ticktype = ticktype, STAND = STAND,
SEED = SEED, expand = expand, scale = scale, nboot = nboot,
...)
val <- val$output
}
E.power <- varfun(val[!is.na(val)])/varfun(y)
if (!pyhat)
val <- NULL
E.power = as.numeric(E.power)
list(Strength.Assoc = sqrt(E.power), Explanatory.Power = E.power,
yhat = val)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.