Discussion :

[Zaazoun]

Bonjour,

Je debute avec R, je doit filtrer une serie temporelle mensuelle de précipitation mensuelle centré réduite De 1962 à 2005 sur R,Je veux filtrer avec le filtre de butterworth du package mFilter (qu'il faut charger).Je voudrais régler les parramatre de ma fonction bwfilter de sorte que le filtre conserve 99% DE L'amplitude des signaux de periode 4 ans.Que dois je taper pour lui indiquer que je veux une periode de 4 ans soit 48 mois?
Voici ce que j'ai déja fais:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
RR<-read.csv2("RRmt.csv",sep="\t",dec=",")#je lis mon fichier csv
RR<-ts(RR,frequency=12,start=c(1962,1)) #je transforme en serie temporelle
plot(RR)
library(mFilter)#je charge mon package
opar<-par(no.readonly=TRUE)
RR.bw<-bwfilter(RR)
plot(RR.bw)
RR.bw1<-bwfilter(RR,drift=TRUE)
RR.bw2<-bwfilter(RR,freq=12,drift=TRUE)
RR.bw3<-bwfilter(RR,freq=48,drift=TRUE)
par(mfrow=c(2,1),mar=c(3,3,2,1),cex=.8)
plot(RR.bw1$x,main="Butterworth filter of precipitation: Trend,drift=TRUE",col=1,ylab="")
lines(RR.bw2$trend,col=3)
lines(RR.bw3$trend,col=4)

Merci d'avance

[pitipoisson]

Bonjour et bienvenue,

Tu peux simplement adapter la fréquence (1/4=0.25) :

Code :

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

1
2
3
4
5
6
tmp <- rnorm(8)
 
tmp.ts <- ts(tmp, frequency=0.25, start=c(1962,1))
tmp.ts
 
plot(tmp.ts)

[Zaazoun]

bonjours
merci pour la réponse.

Par compte c'est pas la fréquence de la série temporelle que je veux changer, mais plutot savoir a quoi correspond la fréquence de la fonction bwfilter, par rapport a la fréquence de la série temporelle...
mais merci bcp pour la rep

[Zaazoun]

Bonjours,
Pour être encore plus précise le filtre doit conserver 99% de l'amplitude des signaux de période 4 ans et 1% de celle des signaux de période 1 ans.