Vectorize {base} | R Documentation |
Vectorize
creates a function wrapper that vectorizes the
action of its argument FUN
.
Vectorize(FUN, vectorize.args = arg.names, SIMPLIFY = TRUE,
USE.NAMES = TRUE)
FUN |
function to apply, found via |
vectorize.args |
a character vector of arguments which should be
vectorized. Defaults to all arguments of |
SIMPLIFY |
logical or character string; attempt to reduce the
result to a vector, matrix or higher dimensional array; see
the |
USE.NAMES |
logical; use names if the first ... argument has names, or if it is a character vector, use that character vector as the names. |
The arguments named in the vectorize.args
argument to
Vectorize
are the arguments passed in the ...
list to
mapply
. Only those that are actually passed will be
vectorized; default values will not. See the examples.
Vectorize
cannot be used with primitive functions as they do
not have a value for formals
.
It also cannot be used with functions that have arguments named
FUN
, vectorize.args
, SIMPLIFY
or
USE.NAMES
, as they will interfere with the Vectorize
arguments. See the combn
example below for a workaround.
A function with the same arguments as FUN
, wrapping a call to
mapply
.
# We use rep.int as rep is primitive
vrep <- Vectorize(rep.int)
vrep(1:4, 4:1)
vrep(times = 1:4, x = 4:1)
vrep <- Vectorize(rep.int, "times")
vrep(times = 1:4, x = 42)
f <- function(x = 1:3, y) c(x, y)
vf <- Vectorize(f, SIMPLIFY = FALSE)
f(1:3, 1:3)
vf(1:3, 1:3)
vf(y = 1:3) # Only vectorizes y, not x
# Nonlinear regression contour plot, based on nls() example
require(graphics)
SS <- function(Vm, K, resp, conc) {
pred <- (Vm * conc)/(K + conc)
sum((resp - pred)^2 / pred)
}
vSS <- Vectorize(SS, c("Vm", "K"))
Treated <- subset(Puromycin, state == "treated")
Vm <- seq(140, 310, length.out = 50)
K <- seq(0, 0.15, length.out = 40)
SSvals <- outer(Vm, K, vSS, Treated$rate, Treated$conc)
contour(Vm, K, SSvals, levels = (1:10)^2, xlab = "Vm", ylab = "K")
# combn() has an argument named FUN
combnV <- Vectorize(function(x, m, FUNV = NULL) combn(x, m, FUN = FUNV),
vectorize.args = c("x", "m"))
combnV(4, 1:4)
combnV(4, 1:4, sum)