calcul avec boucle

Bonjour,
Je veux calculer 17 valeurs (prev_juste) à partir 517 valeurs ; chaque valeur de "prev_juste" est calculée à partir de 30 valeurs de ces 517 données.

Voilà un echantillon de :

Code:

          r            p
  [1,] 140305  27290
  [2,] 234094 186669
  [3,]  39960 629890
  [4,] 270511 102235
  [5,] 260099  45712
  [6,] 261448  93736
  [7,] 208707  84767
  [8,] 254184  65987
  [9,]  51059  14399
 [10,] 272400 221521
 [11,] 263225 301509
 [12,] 257781 205960
 [13,] 227997 209948
 [14,] 236170 192819
 [15,]  38553  35596
 [16,] 252677 217541
 [17,] 247943 230868
 [18,] 248403 200379
 [19,] 207028 214038


Je vous donne le script que j'ai utilisé:

Code:

sommeR=0
l=0
k=1
i=0
while(i<length(r)-38)
{
for (k in 1)
 {
  v[k]=p[i+k]-r[i+k]
  sommeR=sommeR+r[i+k]
 }
i=i+30
valAbs=abs(v)
mae=mean(valAbs,na.rm=TRUE)
moyR=sommeR/30
prev_juste=(1-mae/moyR)*100 
sommeR=0
l=l+1
print(prev_juste)
}

J'ai eu un resultat mais il ne correspond pas à ce que recherche (j'ai recours au calcul manuel pour confirmer l'obtention de la premiére valeur) i.e j'obtiens 64.68 et pas 96.65 comme affiché ci-dessous :

R m'affiche:

Code:

[1] 96.65835
[1] 98.00383
[1] 97.85386
[1] 96.77248
[1] 96.86521
[1] 97.93428
[1] 97.76311
[1] 97.22779
[1] 96.68919
[1] 97.76993
[1] 97.86475
[1] 97.32567
[1] 97.46467
[1] 97.66616
[1] 98.07713
[1] 97.8267

Je vous remercie d'avance

Oula,

Je veux bien te donner un coup de main,
mais à quoi correspond exactement la valeur que tu calcules?
Sans un minimum d'explications, dur de voir ce qui cloche dans ta boucle...

Bien cordialement