我是R的新手，因此我的问题可能很简单，但是尽管如此，我还是花了很多时间试图弄清楚我做错了什么，但毫无用处。在过去一周中，我在此站点上搜索其他问题/答案时发现了很多帮助（谢谢！），但是作为新手，通常很难解释别人的代码。
我正在尝试构建多个数据文件的3维数组，每个文件具有相同的尺寸57x57。
# read in 100 files
Files = lapply(Sys.glob('File*.txt'), read.table, sep='\t', as.is=TRUE)
# convert to dataframes
Files = lapply(Files[1:100], as.data.frame)
# check dimensions of first file (it's the same for all)
dim(Files[[1]])
[1] 57 57
# build empty array
Array = array(dim=c(57,57,100))
# read in the first data frame
Array[,,1] = Files[1]
# read in the second data frame
Array[,,2] = Files[2]
Error in Array[, , 2] = Files[2] : incorrect number of subscripts
# if I check...
Array[,,1] = Files[1]
Error in Array[, , 1] : incorrect number of dimensions
# The same thing happens when I do it in a loop:
x = 0
for(i in 1:100){
Array[,,x+1] = Files[[i]]
x = x + 1
}
Error in Array[, , 1] = Files[[1]] :
incorrect number of subscripts
在进行分配之前，您需要将数据帧转换为矩阵：
l <- list(data.frame(x=1:2, y=3:4), data.frame(x=5:6, y=7:8))
arr <- array(dim=c(2, 2, 2))
arr[,,1] <- as.matrix(l[[1]])
arr[,,2] <- as.matrix(l[[2]])
arr
# , , 1
#
# [,1] [,2]
# [1,] 1 3
# [2,] 2 4
#
# , , 2
#
# [,1] [,2]
# [1,] 5 7
# [2,] 6 8
实际上，您可以将unlist函数应用到要合并的矩阵列表中，从而在一行中构建数组：
arr2 <- array(unlist(lapply(l, as.matrix)), dim=c(dim(l[[1]]), length(l)))
all.equal(arr, arr2)
# [1] TRUE
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句