R 在数组索引上复制矩阵时出错

这和R的回收规则有关如果在初始化后查看test_array，您会看到以下输出:

test_array <- array(NA, c(4,3,2))
test_matrix <- t(matrix(1:2, nrow = 2, ncol = 4))

test_array
# , , 1
# 
#      [,1] [,2] [,3]
# [1,]   NA   NA   NA
# [2,]   NA   NA   NA
# [3,]   NA   NA   NA
# [4,]   NA   NA   NA
# 
# , , 2
# 
#      [,1] [,2] [,3]
# [1,]   NA   NA   NA
# [2,]   NA   NA   NA
# [3,]   NA   NA   NA
# [4,]   NA   NA   NA

你的test_matrix看起来像这样:

test_matrix
#      [,1] [,2]
# [1,]    1    2
# [2,]    1    2
# [3,]    1    2
# [4,]    1    2

现在，如果你想把它存储到test_array[, 1:3, ](顺便说一句，这与test_array相同)，R复制矩阵，使其符合尺寸.也就是说，每个奇数列将是一个用1填充的列，每个偶数列将用2填充(如果这有助于理解，可以将第二个矩阵的第一列视为总的第四列).

也就是说，如果你要求test_array[, 1, ]，你会得到每个矩阵的第一列.

那么，为什么它工作，如果你使用单一的索引？这是b/c R如果可能的话会悄悄地删除尺寸，除非你指定drop = FALSE:

dim(test_array[, 1, ])
# [1] 4 2
dim(test_array[, 1:3, ])
# [1] 4 3 2
dim(test_array[, 1, , drop = FALSE])
# [1] 4 1 2

总的来说，在赋值的情况下，如果使用多个列，R需要进行一些循环来拟合维度，而对于单个列，维度匹配，不需要循环.

为了解决这个问题，你需要复制你的测试矩阵来考虑最后一个维度:

r <- rep(test_matrix, each = 3)
dim(r) <- dim(test_array)

test_array[, 1:3, ] <- r
test_array
# , , 1
# 
#      [,1] [,2] [,3]
# [1,]    1    1    1
# [2,]    1    1    1
# [3,]    1    1    1
# [4,]    1    1    1
# 
# , , 2
# 
#      [,1] [,2] [,3]
# [1,]    2    2    2
# [2,]    2    2    2
# [3,]    2    2    2
# [4,]    2    2    2