Question

Is there some function in R of the form like:

crossApply(v1, v2, func)

that has the same functionality as:

ret = c()
i = 1
for (e1 in v1) {
  for (e2 in v2) {
    ret[i] <- func(e1,e2)
    i <- i + 1
  }
}
return(ret)

Thanks in advance.

Answer 1

I think you may be looking for outer which isn’t exactly what your code does, but it’s close. Specifically, outer will return the matrix (i.e. the outer product) or each combination of the elements of its first two arguments.

You’d probably want to store the result and then extract the lower triangle as a vector. Something like this maybe:

rs <- outer(1:4,-(5:7),"+")
rs[lower.tri(rs,diag = TRUE)]
[1] -4 -3 -2 -1 -4 -3 -2 -4 -3