R中是否有一个函数可以轻松构建重复的单位矩阵(不确定这是否是正确的术语)?
这是我当前正在使用的,但是有点笨拙和不直观:
Ngroups <- 3
NperGroup <- 2
Z <- diag(Ngroups)[rep(1:Ngroups, each = NperGroup), ]
> Z
[,1] [,2] [,3]
[1,] 1 0 0
[2,] 1 0 0
[3,] 0 1 0
[4,] 0 1 0
[5,] 0 0 1
[6,] 0 0 1
这是一种方法
diag(3) %x% c(1,1)
[,1] [,2] [,3]
[1,] 1 0 0
[2,] 1 0 0
[3,] 0 1 0
[4,] 0 1 0
[5,] 0 0 1
[6,] 0 0 1
在你的情况下diag(Ngroups) %x% rep(1,NperGroup)
会做到的。查看?kronecker
更多选项
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句