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
	}
}