3 votes

Moyenne mobile par groupe R data.table

Je veux calculer une moyenne mobile annuelle par groupe en commençant par la première ligne du groupe et en terminant par la dernière ligne. Exemple ci-dessous...

Group <- c(rep("a",5), rep("b",5))
Sales <- c(2,4,3,3,5,9,7,8,10,11)
Result <- c(2,3,3,3,3.4,9,8,8,8.5,9)
df <- data.frame(Group, Sales, Result)

La colonne Résultat correspond à ce que j'attends de la moyenne mobile.

11voto

nachti Points 549

Utilisation de cumsum :

dt <- as.data.table(df)
dt[, res := cumsum(Sales)/(1:.N), by = Group]
dt
    Group Sales Result res
 1:     a     2    2.0 2.0
 2:     a     4    3.0 3.0
 3:     a     3    3.0 3.0
 4:     a     3    3.0 3.0
 5:     a     5    3.4 3.4
 6:     b     9    9.0 9.0
 7:     b     7    8.0 8.0
 8:     b     8    8.0 8.0
 9:     b    10    8.5 8.5
10:     b    11    9.0 9.0

ou avec rollapplyr de la zoo paquet :

dt[, res := rollapplyr(Sales, 1:.N, mean), by = Group]

ou avec la base R :

ave(df$Sales, df$Group, FUN = function(x) cumsum(x) / seq_along(x))

2voto

akrun Points 148302

Nous pouvons utiliser dplyr con zoo . Le site %>% relie la lhs avec la rhs et il est très facile de le comprendre et de l'exécuter.

library(dplyr)
library(zoo)
df %>%
   group_by(Group) %>% 
   mutate(Sales = rollapplyr(Sales, row_number(), mean))  
#    Group Sales Result
#    <fctr> <dbl>  <dbl>
#1       a   2.0    2.0
#2       a   3.0    3.0
#3       a   3.0    3.0
#4       a   3.0    3.0
#5       a   3.4    3.4
#6       b   9.0    9.0
#7       b   8.0    8.0
#8       b   8.0    8.0
#9       b   8.5    8.5
#10      b   9.0    9.0

Prograide.com

Prograide est une communauté de développeurs qui cherche à élargir la connaissance de la programmation au-delà de l'anglais.
Pour cela nous avons les plus grands doutes résolus en français et vous pouvez aussi poser vos propres questions ou résoudre celles des autres.

Powered by:

X