1 | regYci.sum(x, y, regfun = tsreg, pts = x, nboot = 100, xout = FALSE, outfun = out, SEED = TRUE, alpha = 0.05, crit = NULL, null.value = 0, plotPV = FALSE, ADJ = FALSE, MC = FALSE, scale = FALSE, span = 0.75, xlab = "X", xlab1 = "X1", xlab2 = "X2", ylab = "p-values", theta = 50, phi = 25, pch = "*", ...)
|
x |
|
y |
|
regfun |
|
pts |
|
nboot |
|
xout |
|
outfun |
|
SEED |
|
alpha |
|
crit |
|
null.value |
|
plotPV |
|
ADJ |
|
MC |
|
scale |
|
span |
|
xlab |
|
xlab1 |
|
xlab2 |
|
ylab |
|
theta |
|
phi |
|
pch |
|
... |
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 (x, y, regfun = tsreg, pts = x, nboot = 100, xout = FALSE,
outfun = out, SEED = TRUE, alpha = 0.05, crit = NULL, null.value = 0,
plotPV = FALSE, ADJ = FALSE, MC = FALSE, scale = FALSE, span = 0.75,
xlab = "X", xlab1 = "X1", xlab2 = "X2", ylab = "p-values",
theta = 50, phi = 25, pch = "*", ...)
{
xy = elimna(cbind(x, y))
x <- as.matrix(x)
p = ncol(x)
if (p > 1)
stop("This function is designed for one independent variable only")
p1 = p + 1
x <- xy[, 1:p]
y <- xy[, p1]
x <- as.matrix(x)
if (xout) {
m <- cbind(x, y)
flag <- outfun(x, plotit = FALSE, ...)$keep
m <- m[flag, ]
x <- m[, 1:p]
y <- m[, p1]
x = as.matrix(x)
}
res = regYci(x, y, regfun = regfun, nboot = nboot, null.value = null.value,
alpha = alpha, crit = crit, ADJ = ADJ, MC = MC, plotPV = plotPV,
...)
xord = order(pts)
outp = cbind(pts[xord], res[xord, ])
dimnames(outp) = list(NULL, c("X", "Pred. Y", "Lower.ci",
"Upper.ci", "p.value"))
outp
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.