我想构建一个维数为T x T的矩阵。在第一行和最后一列中我想要全零。此外,从第二行直到行T以及从第一列直到列T-1,我想要单位矩阵。在T = 4的情况下,它看起来应该是这样的:
1. column 2. column 3. column 4. column
1. row: 0 0 0 0
2. row: 1 0 0 0
3. row: 0 1 0 0
4. row: 0 0 1 0
我希望这是有道理的,
感谢。
答案 0 :(得分:2)
查看? diag
,rbind
和?cbind
:
n <- 4
rbind(rep(0, n), cbind(diag(1, n-1), rep(0, n-1)))
[,1] [,2] [,3] [,4]
[1,] 0 0 0 0
[2,] 1 0 0 0
[3,] 0 1 0 0
[4,] 0 0 1 0
(请注意,T
通常用作TRUE
的快捷方式,因此您应该将其作为变量名称使用,否则您将遇到一些问题......)