将新列添加到 R 中矩阵列表中的每个矩阵

Question

我有一个包含 61 个矩阵的矩阵列表。下面是 61 个矩阵之一的示例。

list_matrix_Tanzania_Mod500
[[60]]
      LU_2000
 [1,]      12
 [2,]      12
 [3,]      12
 [4,]      12
 [5,]      12
 [6,]      12
 [7,]      12
 [8,]      12
 [9,]      12
[10,]      11
[11,]      11
[12,]      12
[13,]      11
[14,]      12
[15,]      12 

我想向矩阵列表的每个矩阵添加 10 个新列。添加这些新列后，我想为每个列指定一个标题名称，其值为NA："ID", "LU_1990", "CHLU_90_00", "LU_2005", "CHLU_00_05", "Tile", "UNIQ_ID", "AREA", "Sour_90_00", "Sour_00_05"

我现在想用几个矩阵填充矩阵列表 list_matrix_Tanzania_Modis：LU_Mod2000, LU_Mod2005 and list_matrix_Tanzania。我试过了，但该列没有很好地更新。例如，当我为 LU_Mod2005 运行 Mapcode 行时，LU_2000 的值会从矩阵中删除。

columnsToTransfer1 <- 'LU_2000'
columnsToTransfer2 <- 'LU_2005'
columnsToTransfer3<- ('ID', 'Tile', 'UNIQ_ID')

    Map(function(x,y,z) {x[,z] <- y[,z];x},
                list_matrix_Tanzania_Modis, LU_Mod2000, columnsToTransfer1)
    Map(function(x,y,z) {x[,z] <- y[,z];x},
        list_matrix_Tanzania_Modis, LU_Mod2005, columnsToTransfer2)
Map(function(x,y,z) {x[,z] <- y[,z];x},
    list_matrix_Tanzania_Modis, list_matrix_Tanzania, columnsToTransfer3)

这个技巧有效，但我想将矩阵列表的名称保留为 list_matrix_Tanzania_Modis

list_matrix_Tanzania_Modis1<- Map(function(x,y,z) {x[,z] <- y[,z];x},
            list_matrix_Tanzania_Modis, LU_Mod2000, columnsToTransfer1)
list_matrix_Tanzania_Modis2<-Map(function(x,y,z) {x[,z] <- y[,z];x},
    list_matrix_Tanzania_Modis1, LU_Mod2005, columnsToTransfer2)

Answer 1

如果我理解正确，一个简单的cbind 在创建一个正确尺寸的空矩阵后应该可以工作：

示例输入：

L <- list(matrix(1:3, ncol = 1, dimnames = list(NULL, "a")),
          matrix(1:5, ncol = 1, dimnames = list(NULL, "b")))
L
# [[1]]
#      a
# [1,] 1
# [2,] 2
# [3,] 3
# 
# [[2]]
#      b
# [1,] 1
# [2,] 2
# [3,] 3
# [4,] 4
# [5,] 5

您要添加的新列：

newcols <- c("ID", "LU_1990", "CHLU_90_00", "LU_2005", "CHLU_00_05",
             "Tile", "UNIQ_ID", "AREA", "Sour_90_00", "Sour_00_05")

使用lapply 循环浏览列表，将相关的空列添加到每个列表项。

myNewList <- lapply(L, function(x) {
  M <- matrix(NA, nrow = nrow(x), ncol = 10, dimnames = list(NULL, newcols))
  cbind(x, M)
})

myNewList
# [[1]]
#      a ID LU_1990 CHLU_90_00 LU_2005 CHLU_00_05 Tile UNIQ_ID AREA Sour_90_00 Sour_00_05
# [1,] 1 NA      NA         NA      NA         NA   NA      NA   NA         NA         NA
# [2,] 2 NA      NA         NA      NA         NA   NA      NA   NA         NA         NA
# [3,] 3 NA      NA         NA      NA         NA   NA      NA   NA         NA         NA
# 
# [[2]]
#      b ID LU_1990 CHLU_90_00 LU_2005 CHLU_00_05 Tile UNIQ_ID AREA Sour_90_00 Sour_00_05
# [1,] 1 NA      NA         NA      NA         NA   NA      NA   NA         NA         NA
# [2,] 2 NA      NA         NA      NA         NA   NA      NA   NA         NA         NA
# [3,] 3 NA      NA         NA      NA         NA   NA      NA   NA         NA         NA
# [4,] 4 NA      NA         NA      NA         NA   NA      NA   NA         NA         NA
# [5,] 5 NA      NA         NA      NA         NA   NA      NA   NA         NA         NA