r - 从矩阵的对角线计算平均值 + 增量

Question

我有这种矩阵。

非常抱歉，我没有可重现的示例。

表 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 列和 100 行/

非常感谢

Answer 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