24

在今天的答案中引用的@hadley 关于函数的文章的提示下，我决定重新审视一个关于outer函数如何工作（或不工作）的持久性难题。为什么会失败：

outer(0:5, 0:6, sum) # while outer(0:5, 0:6, "+") succeeds

这显示了我认为outer 应该如何处理如下函数sum：

 Outer <- function(x,y,fun) {
   mat <- matrix(NA, length(x), length(y))
   for (i in seq_along(x)) {
            for (j in seq_along(y)) {mat[i,j] <- fun(x[i],y[j])} }
   mat}

>  Outer(0:5, 0:6, `+`)
     [,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,]    0    1    2    3    4    5    6
[2,]    1    2    3    4    5    6    7
[3,]    2    3    4    5    6    7    8
[4,]    3    4    5    6    7    8    9
[5,]    4    5    6    7    8    9   10
[6,]    5    6    7    8    9   10   11

好的，我的索引并没有与该示例完全对齐，但修复起来并不难。问题是为什么像这样的函数sum应该能够接受两个参数并返回适合矩阵元素的（原子）值，而在传递给base::outer函数时却不能这样做？

所以@agstudy 为更紧凑的版本提供了灵感，Outer他的版本更加紧凑​​：

 Outer <- function(x,y,fun) {
       mat <- matrix(mapply(fun, rep(x, length(y)), 
                                 rep(y, each=length(x))),
                     length(x), length(y))

然而，问题仍然存在。术语“矢量化”在这里有些模棱两可，我认为“二元”更正确，因为sin并且cos在该术语的通常意义上是“矢量化”的。outer期望以可以使用非二元函数的方式扩展其论点是否存在基本的逻辑障碍。

这是另一个outer错误，可能与我对这个问题缺乏了解类似：

> Vectorize(sum)
function (..., na.rm = FALSE)  .Primitive("sum")
>  outer(0:5, 0:6, function(x,y) Vectorize(sum)(x,y) )
Error in outer(0:5, 0:6, function(x, y) Vectorize(sum)(x, y)) : 
  dims [product 42] do not match the length of object [1]
4

1 回答 1

32

outer(0:5, 0:6, sum) 不起作用，因为sum没有“矢量化”（在返回与其两个参数长度相同的向量的意义上）。这个例子应该解释差异：

 sum(1:2,2:3)
  8
 1:2 + 2:3
 [1] 3 5

sum您可以使用mapply例如矢量化：

identical(outer(0:5, 0:6, function(x,y)mapply(sum,x,y)),
          outer(0:5, 0:6,'+'))
TRUE

PS：一般在使用之前outer我使用browser在调试模式下创建我的函数：

outer(0:2, 1:3, function(x,y)browser())
Called from: FUN(X, Y, ...)
Browse[1]> x
[1] 0 1 2 0 1 2 0 1 2
Browse[1]> y
[1] 1 1 1 2 2 2 3 3 3
Browse[1]> sum(x,y)
[1] 27          ## this give an error 
Browse[1]> x+y  
[1] 1 2 3 2 3 4 3 4 5 ## this is vectorized
于 2013-08-07T18:09:54.193 回答