r - 是否可以在都是矩阵的 A 和 B 上应用外部，以及采用 2 个向量并返回一个标量的函数？

Question

我希望我可以使用类似的东西：

 outer(A,B,myfun)

实现类似：

a<-matrix(nrow=nrow(A),ncol=nrow(B))
for(i in 1:nrow(A))
    for (j in 1:nrow(B))
    {
        a[i,j]<-myfun(A[i,],B[j,])
    }

有没有更好的办法？

Answer 1

无法给出一般性答案，因为详细信息myfun可能会影响答案。

ares <- expand.grid(1:nrow(A), 1:nrow(B))
ares$res <- myfun(A[ares[,1], ), B[res[,2] ) 
# but may need mapply("myfun", A[ares[,1], ), B[res[,2] ) on which Vectorize is based
# or do.call(my.fun, ....)
a <- matrix(ares$res, nrow(A), nrow(B) )