Среднее в трехмерном массиве каждые 12 шагов для третьего измерения

Question

У меня есть массив a

dim (a) 
[1] 126 59 240

Я хотел бы вычислять среднее значение третьего измерения каждые 12 шагов.Я пробовал с:

rollapply(a,3,mean,by=12 )

, но в результате получается вектор длиной 11 и не более.

Это лучше apply и определенная функция?Мой ожидаемый результат - новый массив с dim 126 59 20

Answer 1

Определенно не просто, но это работает.

a = array(1:(126*59*240),c(126,59,240))          # Dummy data
s = seq(1,240,by=12)                             # Sequence of cutpoints for 3rd dim
aux =lapply(s,function(x){
  apply(a[,,x:(x+11)],c(1,2),mean,na.rm=TRUE)    # Get list with mean per each 12 3rd dim
})

aux2 = do.call(cbind,aux)                        # Bind all elements of list by column
Y = array(aux2,dim=c(dim(aux[[1]]),length(aux))) # Reconvert into array

> dim(Y)
[1] 126  59  20