我是使用R的新手,因此我的问题可能是一个简单的问题,但尽管如此,我花了很多时间试图弄清楚我做错了什么,但无济于事。在过去的一周里,我在这个网站上发现了很多帮助,通过其他问题/答案搜索(谢谢!),但作为一个新手,通常很难解释别人的代码。
我正在尝试构建多个数据文件的三维数组,每个文件具有相同的尺寸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