将代表扩展到矩阵？

Question

如果您在矩阵上调用rep ，它会重复其元素而不是整个矩阵。 传统的解决方法是调用rep(list(theMatrix),...) 。 我想扩展rep以便它自动执行此操作。

我试图使用

rep.matrix<-function(x,...) rep(list(x),...)

这确实将rep.matrix添加到methods(rep)

> methods(rep)
[1] rep.bibentry*       rep.Date            rep.factor          rep.matrix         
[5] rep.numeric_version rep.POSIXct         rep.POSIXlt         rep.roman*         
see '?methods' for accessing help and source code

但是，在矩阵上调用 rep 似乎并没有发送到rep.matrix 。

> rep(diag(5),3)
 [1] 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0
[42] 0 0 1 0 0 0 0 0 1 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1

尽管直接调用rep.matrix没有错误。

> rep.matrix(diag(5),3)
[[1]]
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    0    0    0    0
[2,]    0    1    0    0    0
[3,]    0    0    1    0    0
[4,]    0    0    0    1    0
[5,]    0    0    0    0    1

[[2]]
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    0    0    0    0
[2,]    0    1    0    0    0
[3,]    0    0    1    0    0
[4,]    0    0    0    1    0
[5,]    0    0    0    0    1

[[3]]
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    0    0    0    0
[2,]    0    1    0    0    0
[3,]    0    0    1    0    0
[4,]    0    0    0    1    0
[5,]    0    0    0    0    1

如果我创建和使用rep.array而不是rep.matrix ，我会得到相同的结果。

我的错误在哪里？ 为什么rep不发送到rep.matrix ？ 我是否以某种方式使用了错误的 object 系统？

Answer 1

这是无法做到的。 关于这个事实，文档并不太公开，但rep是一个内部通用 function，如此处所述。 这意味着它只会分派给is.object返回TRUE的东西。 矩阵不是这样的东西，因此您不能将rep扩展到矩阵并将其分派给它们。