Discussion :

[mantus]

Salut tout le monde;

J'essaie de compiler un petit programme sur R relatif à une série chronologique, seulement étant débutant, un message d'erreur contrecarre mes avancées.

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
19
20
21
22
 
teta<-0
sigma<-0.05
e<-0
x<-0
obs<-0
k<-1
 
while(arret==0) 
{
x[0]=0
x[-1]=0
e[k]=rnorm(1,teta,sigma)
x[k]=-5*x[k-1]-6*x[k-2]+e[k]
if (k<100) {arret=0} else 	{arret=1}
obs[k]<-k
k<-k+1
}
mat<-cbind(obs,x,e)
mat
plot(obs,x)
lines(obs,L1,type='l')

Le message d'erreur suivant est affiché:

Code :

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

1
2
 
Erreur dans xy.coords(x, y) : 'x' and 'y' lengths differ

Le but est de calculer une série de x[k] où k in 1:100. La difficulté principale de
x[k]=-5*x[k-1]-6*x[k-2]+e[k]
est que cette équation prend des valeurs en (k-1) et (k-2) et cela, à partir d'une variable e[k] simulée par rnorm.

Finalement, le but est de représenter graphiquement les obs[k] en absisse, soient les k unité par unité, et les x[k] calculés par l'équation en ordonnée...

Je bloque et ne doute pas que cela reste un problème élémentaire.
Je vous remercie grandement des solutions apportées!

Mantus

[tomaprice]

Ceci devrait résoudre ton problème

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
e=rep(0,102)
x=rep(0,102)
obs=vector()
obs[1]=1
obs[2]=2
 
 
for(i in 3:102){
 
e[i]=rnorm(1,teta,sigma)
x[i]=-5*x[i-1]-6*x[i-2]+e[i]
obs[i]=i
 
}
 
mat<-cbind(obs,x,e)
mat
plot(obs,x)

[mantus]

Merci, cela fonctionne

Bien à vous