Discussion :

[JonathanJohns]

Bonjour,

Petit problème paraissant simple mais qui pose des soucis. Je souhaite calculer une moyenne glissante sur cette base de données.



Pas de soucis dans le calcul de la moyenne glissante mais par contre j'aimerai en fait faire une troisième colonne dans ce tableau de données avec les résultats de ma moyenne glissante. Le problème est que les tailles de mes vecteurs ne sont pas les mêmes car les deux premiers jours n'ont pas de moyenne glissante.

J'utilise le code suivant pour calculer ma moyenne glissante :

> roll_mean(dataset[,2], 3, weights = NULL, by = 1L, fill = numeric(0), partial = FALSE, align = c("center", "left", "right"), normalize = TRUE, na.rm = FALSE)

J'obtiens mon vecteur :

[1] 60.17209 60.19096 60.47054 60.41606 60.26796 60.05892 60.03073 60.01870 60.01562 60.23962
[11] 60.22036 60.27191 60.21991 60.31625 60.18935

Au final, j'essaie que mon tableau soit composée de deux NA pour les deux premières lignes de la moyenne glissante ou de créer un nouveau tableau avec la colonne "Date" sans les deux premières lignes et la deuxième colonne avec les valeurs de ma moyenne glissante.

Merci d'avance pour votre aide car je rame un peu là dessus. :/

[mgdondon]

Bonjour,

Voici une proposition :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
> df <- data.frame(V1=c(1,1,1,1,1))
> V2 <- c(2,2,2)
> df$V2 <- c(NA,NA,V2)
> df
  V1 V2
1  1 NA
2  1 NA
3  1  2
4  1  2
5  1  2

Cordialement,

[JonathanJohns]

Bonjour mgdondon,

Merci pour votre réponse. Une question supplémentaire pour savoir s'il est possible de généraliser le code, il m'arrive de calculer ma moyenne glissante sur 3 semaines, 4 semaines, 5 semaines etc.. Pour éviter d'avoir à modifier à chaque fois le nombre de NA dans la variable V2, savez vous s'il est possible de généraliser le code ?

J'entends par là que lorsque je calcule une moyenne sur une période de X temps, je dois créer un vecteur V2 avec (X-1) NA. Par exemple, pour une période de 3 semaines, il faut un vecteur V2 avec deux NA au début. Pour une période de 5 semaines, il faut un vecteur V2 avec 4 NA au début etc...

[mgdondon]

Bonjour Jonathan,

Vous pouvez utiliser la fonction rep() pour générer les données manquantes :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
> df <- data.frame(V1=c(1,1,1,1,1))
> V2 <- c(2,2,2)
> decal <- 2
> df$V2 <- c(rep(NA,decal),V2)
> df
  V1 V2
1  1 NA
2  1 NA
3  1  2
4  1  2
5  1  2

Cordialement,

[JonathanJohns]

Merci de votre réponse, le code fonctionne parfaitement.

Cordialement,

Jonathan