首页 > 解决方案 > 创建矩阵,其中所有对角线在反对角方向上的值相同

问题描述

我确定这种类型的矩阵有一个名称,但不确定。在R我想快速转换一个变量x

x = 1:10
x  #[1]  1  2  3  4  5  6  7  8  9 10

成一个矩阵,其中所有对角线(不是主对角线)在一个方向上相同,即反对角方向:

    x
1   1       2       3       4       5       6       7       8       9      10
2   2       3       4       5       6       7       8       9      10      NA
3   3       4       5       6       7       8       9      10      NA      NA
4   4       5       6       7       8       9      10      NA      NA      NA
5   5       6       7       8       9      10      NA      NA      NA      NA
6   6       7       8       9      10      NA      NA      NA      NA      NA
7   7       8       9      10      NA      NA      NA      NA      NA      NA
8   8       9      10      NA      NA      NA      NA      NA      NA      NA
9   9      10      NA      NA      NA      NA      NA      NA      NA      NA
10 10      NA      NA      NA      NA      NA      NA      NA      NA      NA

谢谢

标签: rmatrix

解决方案


我们可以shift使用data.table

library(data.table)
do.call(cbind, shift(x, 0:9, type = 'lead'))
#      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
# [1,]    1    2    3    4    5    6    7    8    9    10
# [2,]    2    3    4    5    6    7    8    9   10    NA
# [3,]    3    4    5    6    7    8    9   10   NA    NA
# [4,]    4    5    6    7    8    9   10   NA   NA    NA
# [5,]    5    6    7    8    9   10   NA   NA   NA    NA
# [6,]    6    7    8    9   10   NA   NA   NA   NA    NA
# [7,]    7    8    9   10   NA   NA   NA   NA   NA    NA
# [8,]    8    9   10   NA   NA   NA   NA   NA   NA    NA
# [9,]    9   10   NA   NA   NA   NA   NA   NA   NA    NA
#[10,]   10   NA   NA   NA   NA   NA   NA   NA   NA    NA

base R,我们可以使用embed

out <- embed(c(x, x), 10)
replace(out, lower.tri(out), NA)[, 10:1]
#      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
# [1,]    1    2    3    4    5    6    7    8    9    10
# [2,]    2    3    4    5    6    7    8    9   10    NA
# [3,]    3    4    5    6    7    8    9   10   NA    NA
# [4,]    4    5    6    7    8    9   10   NA   NA    NA
# [5,]    5    6    7    8    9   10   NA   NA   NA    NA
# [6,]    6    7    8    9   10   NA   NA   NA   NA    NA
# [7,]    7    8    9   10   NA   NA   NA   NA   NA    NA
# [8,]    8    9   10   NA   NA   NA   NA   NA   NA    NA
# [9,]    9   10   NA   NA   NA   NA   NA   NA   NA    NA
#[10,]   10   NA   NA   NA   NA   NA   NA   NA   NA    NA
#[11,]   NA   NA   NA   NA   NA   NA   NA   NA   NA    NA

推荐阅读