Discussion :

[pamram]

Bonjour, j'essaie de faire tourner un programme qui contient une erreur, au niveau de la ligne somme = somme + alpha^(j-1-k)*Xi[k] et en particulier du terme Xi[k]. Je ne sais pas si c'est parce qu'il s'agit d'une série temporelle (je suis une débutante en R), si quelqu'un pouvait m'aider à la résoudre ce serait super !! Merci.

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
23
24
25
26
27
#Lissage Simple
i=1
j=2
k=1
l=1
varls=0
phi1=0.2
phi2=-0.48
alpha=0.001
pas=0.001
somme=0
Xichapls=matrix(nrow=100,ncol=1)
Erreur=matrix(nrow=100,ncol=1)
Erreurmin=0
alphaopt=0
for (i in 1:100){
	Xi <- arima.sim(100,model=list(ar=c(phi1,phi2),ma=c(0,0)))
	#optimisation du alpha
	##### 1. alpha=0.001
	for (j in 2:100){
		for(k in 1:j-1){
		somme = somme + alpha^(j-1-k)*Xi[k]  
		cat(somme,"\n")
		}
		Xichapls[j]=(1-alpha)*somme
	}
}

[hoccha]

Hello,

Je n'ai pas été dans les détails de votre code mais je vois déjà un problème à la ligne 21:

Code :

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

for(k in 1:j-1){

devrait être

Code :

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

for(k in 1:(j-1)){

Et alors votre code ne bloquera plus .... mais je ne sais pas si le résultat est celui voulu.

Pour info, c'est une erreur courante en R. Il ne faut pas oublier que l’opérateur ":" est réalisé avant l'opérateur "-". Dans votre programme, la variable "k" prenait la valeur "0" et donc "Xi[k]" avait une valeur "numeric(0)".