R применить функцию к последовательным значениям двух векторов

Question

Я пытаюсь применить функцию к последовательным значениям двух векторов. Проблема может быть упрощена до этого:

x = c(2,4,6,8,10)
y = c(1,2,3,4,5)

for (i in 1:(length(x)-1)) {
    a = mean(x[i:(i+1)])
    b = mean(y[i:(i+1)])
    print(a/b)
}

Мне интересно, что было бы более эффективным способом сделать то же самое, используя функцию из семейства apply? Спасибо!

Answer 1

Ну, может быть, констатирую очевидное, но:

(x[-1] + x[-length(x)]) / (y[-1] + y[-length(x)])
# [1] 2 2 2 2

(это в 15 раз быстрее, чем решение COLMEANS из теста Onyambu)

Answer 2

library(zoo)
rollmean(x,2) / rollmean(y,2)

Answer 3

Вы можете использовать базу R:

a = sapply(1:(length(x)-1),function(x)x:(x+1))
colMeans(structure(x[a],.Dim=dim(a)))/colMeans(structure(y[a],.Dim=dim(a)))
[1] 2 2 2 2

или даже

 tapply(x[a], col(a), mean)/tapply(y[a], col(a), mean)

Кажется, что эти базовые функции работают быстрее, чем свертывание, так как свертывание выполнит свертывание дважды:

microbenchmark::microbenchmark(
   COLMEANS={a=sapply(1:(length(x)-1),function(x)x:(x+1))
   colMeans(structure(x[a],.Dim=dim(a)))/
     colMeans(structure(y[a],.Dim=dim(a)))},
   ROLLAPPLY=rollapply(x,2,mean)/rollapply(y,2,mean),
   TAPPLY={a=sapply(1:(length(x)-1),function(x)x:(x+1))
     tapply(x[a], col(a), mean)/tapply(y[a], col(a), mean)
   },
   ROLLMEANS=rollmean(x,2) / rollmean(y,2)
 )
 Unit: microseconds
      expr      min        lq      mean   median        uq      max neval
  COLMEANS   76.517  109.0040  187.7223  138.499  178.0405 4598.685   100
 ROLLAPPLY 1752.185 1954.8040 2144.2619 2028.543 2211.9260 6244.430   100
    TAPPLY  398.827  519.3725  665.5016  604.224  682.4505 5304.859   100
ROLLMEANS 1366.610 1619.6715 1815.3349 1731.0260 1957.7965 2615.240   100

Answer 4

Вы можете использовать rollapply из библиотеки зоопарка:

  library(zoo)
  rollapply(x,2,mean)/rollapply(y,2,mean)