Considerar matricesd yr condim(d) = J x D ydim(r) = J x R. Deje que la diversión (a, b) sea una función que toma dos vectores de la misma longitud y devuelve algún número.
Quiero tratar las columnas ded yr respectivamente como mis unidades de interés y aplicarouter a ellos.

l siguiente código logra esto mediante la creación de listas de las columnas ded yr y luego usando ambosouter ysapply:

d.cols <- split(d, col(d))
r.cols <- split(r, col(r))
outer(d.cols, r.cols,
      function(x,y) {
           sapply(seq_along(x),
                 function(i) {
                     Fun(x[[i]], y[[i]]) })} )

El código hace lo que quiero y es relativamente eficiente, pero es torpe y poco claro. ¿Hay una mejor manera de lograr lo que estoy tratando de alcanzar?