rplotsm:

Usage Arguments Examples

Usage

1
rplotsm(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, ...)

Arguments

x
y
est
fr
plotit
pyhat
nboot
atr
nmin
outfun
eout
xlab
ylab
scat
SEED
expand
scale
STAND
varfun
pr
ticktype
theta
phi
...

Examples

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

musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.