Discussion :

[Florence Magnin]

Bonjour,

J'essaye de faire un graphe très simple en faisant une moyenne mobile sur une série de données de températures journalières en fonction du temps. J'ai des difficultés sur plusieurs points. Déjà, la fonction à utiliser:

j'ai chargé le package "gregmisc" pour utiliser la fonction "running" qui me paraît plutôt appropriée. Seulement, elle change la longueur de mon vecteur qui ne correspond plus à celle du temps. Je ne peux donc pas faire un graphe de l'un en fonction de l'autre.

Quelqu'un aurait-il quelques conseils pour que je puisse arriver à mes fins sans faire de bricolage insensé?

En vous remerciant.

[Sengar]

Bonjour,

y étant le vecteur des temperatures
x étant le vecteur des dates

Code :

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

1
2
t<-tapply(y,x,mean)
plot(x=as.Date(names(t)),y=t)

Ca marche?

[Florence Magnin]

Merci de la réponse

J'ai essayé, mais j'obtiens un message d'erreur:

Code :

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

"Error in INDEX[[i]] : subscript out of bounds

J'avais déjà essayé en fasiant:

Code :

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

1
2
 a <- running(y,width=31,align="center")
 b <- running(x,width=31,align="center")

Mais j'avoue que j'ai du mal à concevoir d'appliquer une moyenne mobile sur des dates, c'est surtout qui me pose problème.

Quelqu'un pourrait-il m'éclairer et combler quelques lacunes?

Merci beaucoup.

[Sengar]

Comment se présentent tes données initiales?
(length et summary des variables, et str du dataframe s'il y en a un)

Que veux tu faire concrètement?

Tu veux tracer la moyenne des temperatures en fonction de chaque date dont tu disposes?

Ou tu veux tracer la moyenne des temperatures sur une fenetre glissante d'une durée quelconque? (du genre à la date t, tu traces la temperature moyenne pour toutes les dates entre t et t -30; à t+1 la temperature moyenne entre t+1 et t-29, etc...)

[Florence Magnin]

Bonjour,

tout d'abord merci de l'intérêt porté à ma question.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
str: 
'data.frame':   1130 obs. of  16 variables:
 $ dates: POSIXlt, format: "2009-12-15" "2009-12-16" "2009-12-17" ...
 $ 0.3  : num  -7.88 -9.81 -11.53 -13.16 -14.46 ...
 $ 0.5  : num  -7.41 -8.61 -9.91 -11.29 -12.37 ...
 $ 0.7  : num  -7.22 -7.93 -8.9 -9.99 -10.93 ...
 $ 0.9  : num  -7.09 -7.5 -8.21 -9.04 -9.84 ...
 $ 1.1  : num  -7.01 -7.28 -7.78 -8.4 -9.09 ...
 $ 1.4  : num  -6.96 -7.13 -7.44 -7.86 -8.38 ...
 $ 1.7  : num  -6.92 -7.06 -7.27 -7.59 -8 -8.48 -9.14 -9.72 -9.92 -9.96 ...
 $ 2    : num  -6.98 -7.14 -7.32 -7.62 -7.96 -8.44 -9.17 -9.68 -9.8 -9.85 ...
 $ 2.5  : num  -6.85 -7.08 -7.31 -7.67 -7.96 -8.57 -9.28 -9.22 -9.09 -8.93 ...
 $ 3    : num  -5.86 -5.97 -6.08 -6.24 -6.41 -6.63 -7.03 -7.38 -7.47 -7.55 ...
 $ 4    : num  -4.87 -4.94 -4.99 -5.04 -5.1 -5.18 -5.28 -5.42 -5.56 -5.69 ...
 $ 5    : num  -4.28 -4.33 -4.37 -4.4 -4.44 -4.47 -4.52 -4.57 -4.62 -4.68 ...
 $ 7    : num  -3.88 -3.92 -3.94 -3.94 -3.95 -3.97 -3.99 -4.02 -4.02 -4.05 ...
 $ 9    : num  -4 -4.02 -4.02 -4.01 -4.01 -4.01 -4.02 -4.03 -4.03 -4.03 ...
 $ 10   : num  -4.09 -4.11 -4.11 -4.11 -4.11 -4.11 -4.12 -4.12 -4.12 -4.12 ...

ET je cherche juste à tracer la température de la 2ème colonne en fonction de la 1ère (date) pour chaque date dont je dispose, mais en faisant une moyenne mobile des températures sur une quinzaine de jours simplement pour lisser les trop grandes amplitudes.

Cela paraît-il plus clair?

Merci beaucoup.

FM