我有这种矩阵。
我真的很抱歉,但我没有可重现的例子。
表1:
[,1][,2][,3][,4][,5][,6][,7][,8][,9][,10]
[1,] 3 NA NA NA NA NA NA NA NA NA
[2,] 4 2 NA NA NA NA NA NA NA NA
[3,] 4 1 7 NA NA NA NA NA NA NA
[4,] 4 1 2 3 NA NA NA NA NA NA
[5,] 5 2 0 0 5 NA NA NA NA NA
[6,] 2 0 3 3 5 9 NA NA NA NA
[7,] 6 2 0 0 3 4 2 NA NA NA
[8,] 12 6 4 3 1 0 2 6 NA NA
[9,] 16 11 7 6 5 3 4 0 3 NA
[10,] 19 15 13 9 7 6 6 3 3 5
我想创建另一个这样的: 从对角线重建列。
表2:
[,1][,2][,3][,4][,5][,6][,7][,8][,9]
[1,] 3 4 4 5 2 6 12 16 19
[2,] 2 1 1 0 2 6 11 15 NA
[3,] 7 2 0 0 4 7 13 NA NA
[4,] 3 0 3 3 6 9 NA NA NA
[5,] 5 5 3 5 7 NA NA NA NA
[6,] 9 4 0 6 NA NA NA NA NA
[7,] 2 2 4 NA NA NA NA NA NA
[8,] 6 0 3 NA NA NA NA NA NA
[9,] 3 3 NA NA NA NA NA NA NA
[10,] 5 NA NA NA NA NA NA NA NA
Table 2[,1]= Table 1[1,1];[2,2];[3;3]
Table 2[,2]= Table 1 [2,1];[3,2];[4,3]
Table 2[,3]= Table 1 [3,1];[4,2];[5,3]
我尝试了这段代码,但没有成功。
Table2=matrix(NA, ncol=10, nrow=10)
for(i in 0:9)
{
Table2[i+1]=Table1[i+1,i+1]
}
下一步,将计算一个colMean。 所有这些都是用于评估sarima的交叉验证。 这只是一个例子,在实际数据库中,我有超过100 col和100行/
非常感谢
答案 0 :(得分:1)
这个怎么样?
t(apply(t(yr_mat), 1, function(x) c(na.omit(x), rep(NA, sum(is.na(x))))))
采用矩阵的转置可以实现这一目标:
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
V1 3 4 4 4 5 2 6 12 16 19
V2 NA 2 1 1 2 0 2 6 11 15
V3 NA NA 7 2 0 3 0 4 7 13
V4 NA NA NA 3 0 3 0 3 6 9
V5 NA NA NA NA 5 5 3 1 5 7
V6 NA NA NA NA NA 9 4 0 3 6
V7 NA NA NA NA NA NA 2 2 4 6
V8 NA NA NA NA NA NA NA 6 0 3
V9 NA NA NA NA NA NA NA NA 3 3
V10 NA NA NA NA NA NA NA NA NA 5
所以你需要的是滑过'所有的行,所以这些值一直是尾随的NA。该函数通过截断每一行中的所有NA
来执行此操作,然后在良好值之后添加尽可能多的NA。
最后,你必须再次转置这一切,以便它按照你想要的方式进行。
输出结果为:
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
V1 3 4 4 4 5 2 6 12 16 19
V2 2 1 1 2 0 2 6 11 15 NA
V3 7 2 0 3 0 4 7 13 NA NA
V4 3 0 3 0 3 6 9 NA NA NA
V5 5 5 3 1 5 7 NA NA NA NA
V6 9 4 0 3 6 NA NA NA NA NA
V7 2 2 4 6 NA NA NA NA NA NA
V8 6 0 3 NA NA NA NA NA NA NA
V9 3 3 NA NA NA NA NA NA NA NA
V10 5 NA NA NA NA NA NA NA NA NA
答案 1 :(得分:0)
另一种选择是使用' m1'的row
和col
创建数字索引(' indx'),然后{{1}通过' indx'矩阵(' m1'),转换'列表'到'矩阵'来自split
的{{1}}。我们将输出作为字符串,可以转换为'数字'。
stri_list2matrix
stringi