我正在研究R中的一些数据,这些数据包括由三个空间维度和时间维度组成的四维数组:x,y,z,t.对于我的一些分析,我想获得一组空间坐标x,y,z的时间维度中的所有数据.到目前为止,我已经使用哪个函数来获得感兴趣的空间位置的索引.但是当我去获取与空间位置相对应的时间维度中的所有相关数据时,我找不到优雅的R解决方案并且已经使用了repmat,一个移植的MATLAB函数.
a4d <- array(rnorm(10000), rep(10,4)) #x, y, z, t
#arbitrary set of 3d spatial indices x, y, z (here, using high values at first timepoint)
indices <- which(a4d[,,,1] > 2, arr.ind=TRUE)
str(indices)
# int [1:20, 1:3] 10 2 6 5 8 2 6 8 2 10 ...
# - attr(*, "dimnames")=List of 2
# ..$ : NULL
# ..$ : chr [1:3] "dim1" "dim2" "dim3"
#Now, I would like to use these indices to get data x, y, z for …