Question

假设我有以下数据框，其中每个观察代表2D空间中特定点的给定变量：

data = data.frame(col1 = c(1,2,3,4,5,6,7,8), col2 = c(2,3,'NA',5,6,7,8,9), col3 = c(3,4,5,6,7,8,9,10), col4 = c(2,3,4,1,2,6,7,8),
                  col5 = c(2,3,'NA','NA',6,7,8,9), col6 = c(1,2,3,5,6,7,8,9), col7 = c(1,2,3,4,6,7,'NA','NA'), col8 = c(1,2,3,4,5,6,'NA','NA'))

> print(data)
  col1 col2 col3 col4 col5 col6 col7 col8
1    1    2    3    2    2    1    1    1
2    2    3    4    3    3    2    2    2
3    3   NA    5    4   NA    3    3    3
4    4    5    6    1   NA    5    4    4
5    5    6    7    2    6    6    6    5
6    6    7    8    6    7    7    7    6
7    7    8    9    7    8    8   NA   NA
8    8    9   10    8    9    9   NA   NA

数据帧为8x8（8列8行），总共36个。

我需要通过计算每个2x2观测组的平均值来将其缩小到4x4数据帧。

2x2 NAs组应该返回NA值，而如果在2x2组内则存在＆lt;平均值必须排除4个NAs，例如平均值为2,2,2，NA = 2。

这是我想要的输出：

newcol1  newcol2  newcol3  newcol4
   2        3        2       1.50
   4        4        4       3.50
   6       5.75     6.50      6
   8       8.50     8.50      NA

我认为我可以通过for循环来解决这个问题，而这正是我尝试过的，但没有成功：

a = 1
b = 2

for (i in 1:15) {
      test[[i]] = mean(c(data[a,a], data[a,b], data[b,a], data[b,b]))
      test[[i]] = mean(c(data[a+i,a+i], data[a+i,b+i], data[b+i,a+i], data[b+i,b+i]))
}

我在网上搜索了很多，但我找不到任何类似的问题或解决方案。

有什么建议吗？

是否有任何R包可以进行这种空间分析？

Answer 1

来自关联帖子，使用@flodel answer，将sum替换为mean：

a <- as.matrix(data)

#convert matrix from char to num, ignore warnings, they are due to NAs
class(a) <- "numeric" 

#adapted from @flodel https://stackoverflow.com/a/16884987/680068
res <- tapply(a, list((row(a) + 1L) %/% 2L, (col(a) + 1L) %/% 2L), mean, na.rm = TRUE)

# remove NANs
res[ is.nan(res) ] <- NA
res
#   1    2   3   4
# 1 2 3.00 2.0 1.5
# 2 4 4.00 4.0 3.5
# 3 6 5.75 6.5 6.0
# 4 8 8.50 8.5  NA

Answer 2

这是一个愚蠢的解决方案（不是非常普遍），但适用于有限的例子。我认为这是一种更聪明的方式。

R8toR4 <- function(mat) {
  stopifnot(identical(dim(mat), c(8L, 8L)))

  out <- matrix(nrow = 4, ncol = 4)

  for (i in 1:nrow(out)) {
    for (j in 1:ncol(out)) {
      submatrix <- mat[i * 2 - c(1, 0), j * 2 - c(1, 0)]
      if (all(is.na(submatrix))) {
        out[i, j] <- NA
      } else {
        out[i, j] <- mean(submatrix, na.rm = TRUE)
      }
    }
  }

  out
}


DataMatrix <- as.matrix(data)
R8toR4(DataMatrix)
     [,1] [,2] [,3] [,4]
[1,]    2 3.00  2.0  1.5
[2,]    4 4.00  4.0  3.5
[3,]    6 5.75  6.5  6.0
[4,]    8 8.50  8.5   NA

通过获取2x2空间平均值，将数据帧从8x8缩小到4x4

2 个答案: