Nothing
ace <- function (x, y, wt = rep(1, nrow(x)), cat = NULL, mon = NULL,
lin = NULL, circ = NULL, delrsq = 0.01)
{
x <- as.matrix(x)
if (delrsq <= 0) {
cat("delrsq must be positive")
return()
}
iy <- ncol(x) + 1
l <- matrix(1, ncol = iy)
if (!is.null(circ)) {
for (i in 1:length(circ)) {
if (circ[i] < 0 || circ[i] > ncol(x)) {
cat("bad circ= specification")
return()
}
if (circ[i] == 0) {
nncol <- iy
}
else {
nncol <- circ[i]
}
if (l[nncol] != 2 & l[nncol] != 1) {
cat("conflicting transformation specifications")
return()
}
l[nncol] <- 2
}
}
if (length(mon)) {
for (i in 1:length(mon)) {
if (mon[i] < 0 || mon[i] > ncol(x)) {
cat("bad mon= specification")
return()
}
if (mon[i] == 0) {
nncol <- iy
}
else {
nncol <- mon[i]
}
if (l[nncol] != 3 && l[nncol] != 1) {
cat("conflicting transformation specifications")
return()
}
l[nncol] <- 3
}
}
if (length(lin)) {
for (i in 1:length(lin)) {
if (lin[i] < 0 || lin[i] > ncol(x)) {
cat("bad lin= specification")
return()
}
if (lin[i] == 0) {
nncol <- iy
}
else {
nncol <- lin[i]
}
if (l[nncol] != 4 && l[nncol] != 1) {
cat("conflicting transformation specifications")
return()
}
l[nncol] <- 4
}
}
if (length(cat)) {
for (i in 1:length(cat)) {
if (cat[i] < 0 || cat[i] > ncol(x)) {
cat("bad cat= specification")
return()
}
if (cat[i] == 0) {
nncol <- iy
}
else {
nncol <- cat[i]
}
if (l[nncol] != 5 && l[nncol] != 1) {
cat("conflicting transformation specifications")
return()
}
l[nncol] <- 5
}
}
tx <- x
ty <- y
m <- matrix(0, nrow = nrow(x), ncol = iy)
z <- matrix(0, nrow = nrow(x), ncol = 12)
z <- as.matrix(z)
ns <- 1
mode(x) <- "double"
mode(y) <- "double"
mode(tx) <- "double"
mode(ty) <- "double"
mode(wt) <- "double"
mode(delrsq) <- "double"
mode(z) <- "double"
junk <- .Fortran("mace", p = as.integer(ncol(x)), n = as.integer(nrow(x)),
x = t(x), y = y, w = as.double(wt), l = as.integer(l),
delrsq = delrsq, ns = as.integer(ns), tx = tx, ty = ty,
rsq = double(1), ierr = integer(1), m = as.integer(m),
z = z, PACKAGE = "acepack")
return(junk)
}
avas <- function (x, y, wt = rep(1, nrow(x)), cat = NULL, mon = NULL,
lin = NULL, circ = NULL, delrsq = 0.01, yspan = 0)
{
x <- as.matrix(x);
if (delrsq <= 0) {
cat("delrsq must be positive")
return()
}
iy <- ncol(x) + 1
l <- matrix(1, ncol = iy)
if (length(circ)) {
for (i in 1:length(circ)) {
if (circ[i] < 0 || circ[i] > ncol(x)) {
cat("bad circ= specification")
return()
}
if (circ[i] == 0) {
nncol <- iy
}
else {
nncol <- circ[i]
}
if (l[nncol] != 2 & l[nncol] != 1) {
cat("conflicting transformation specifications")
return()
}
l[nncol] <- 2
}
}
if (length(mon)) {
for (i in 1:length(mon)) {
if (mon[i] < 0 || mon[i] > ncol(x)) {
cat("bad mon= specification")
return()
}
if (mon[i] == 0) {
nncol <- iy
}
else {
nncol <- mon[i]
}
if (l[nncol] != 3 && l[nncol] != 1) {
cat("conflicting transformation specifications")
return()
}
l[nncol] <- 3
}
}
if (length(lin)) {
for (i in 1:length(lin)) {
if (lin[i] < 0 || lin[i] > ncol(x)) {
cat("bad lin= specification")
return()
}
if (lin[i] == 0) {
nncol <- iy
}
else {
nncol <- lin[i]
}
if (l[nncol] != 4 && l[nncol] != 1) {
cat("conflicting transformation specifications")
return()
}
l[nncol] <- 4
}
}
if (length(cat)) {
for (i in 1:length(cat)) {
if (cat[i] < 0 || cat[i] > ncol(x)) {
cat("bad cat= specification")
return()
}
if (cat[i] == 0) {
nncol <- iy
}
else {
nncol <- cat[i]
}
if (l[nncol] != 5 && l[nncol] != 1) {
cat("conflicting transformation specifications")
return()
}
l[nncol] <- 5
}
}
tx <- x
ty <- y
m <- matrix(0, nrow = nrow(x), ncol = ncol(x) + 2)
z <- matrix(0, nrow = nrow(x), ncol = 17)
z <- as.matrix(z)
iters <- matrix(0, nrow = 100, ncol = 2)
mode(x) <- "double"
mode(y) <- "double"
mode(tx) <- "double"
mode(ty) <- "double"
mode(wt) <- "double"
mode(m) <- "integer"
mode(l) <- "integer"
mode(delrsq) <- "double"
mode(z) <- "double"
mode(yspan) <- "double"
mode(iters) <- "double"
junk <- .Fortran("avas", as.integer(ncol(x)), as.integer(nrow(x)),
x, y, wt, l, delrsq, tx = tx, ty = ty, rsq = double(1),
ierr = integer(1), m, z, yspan = yspan, niter = integer(1),
iters = iters, PACKAGE = "acepack")
junk$iters <- junk$iters[1:junk$niter, ]
return(list(x = t(x), y = y, tx = junk$tx, ty = junk$ty, rsq = junk$rsq,
l=l, m, yspan = junk$yspan, iters = junk$iters, niters = junk$niter))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.