formuler une boucle

Bonjour,

j'essaye de calculer un ensemble de valeurs à partir d'une colonne de nombres (nommée "Tini") extraite d'un data.frame nommé D:

Code:

---

1 2 3 4

Tini <- D[[2]] str(Tini) num [1:7585] -3.37 -1.5 -3.96 -7.7 -9 -9.16 -9.2 -9.14 -6.88 -7.07

---

Pour mon calcul j'utilise une boucle mais j'ai du mal à la formuler:

Code:

---

1 2 3 4 5 6 7 8 9 10 11

t <- 10800 Tc <- c()   for (n in Tini[1:n-1]){   T <- (Tini[n+1]-Tini[n])*erfc(z/2*sqrt(k*t))   Tc <- c(Tc,T)   t <- t+10800 }

---

Mais j'obtiens des messages d'erreur me disant que "n" est introuvable.
J'ai essayé en le nommant "i" mais j'obtiens le même message.

Ainsi, l'idée de ma boucle est de calculer T pour chacun des "n" éléments de Tini, à partir du premier élément de Tini et jusqu'à la fin (moins le dernier donc : "n-1");

t est le temps (3h donc 10800 s);

Je le définis une première fois avant la boucle, puis une deuxième fois en fin
de boucle pour que l'incrémentation de temps se fasse à chaque nouveau calcul de T.

Tc est un vecteur qui contiendra l'ensemble des valeurs calculées dans cette boucle.

Je pense avoir exposé l'idée de mon calcul, mais je peine à l'écrire correctement.
Quelqu'un pourrait-il apporter quelques corrections à cette formulation?
Si besoin de plus d'infos pas de problème.
En vous remerciant d'avance.

Bonjour,

Je ne suis pas bien sûre d'avoir compris ce que vous cherchez à faire, cependant je pense que le problème est dû à la manière dont vous écrivez le critère pour la boucle.

Est-ce que ceci ne fonctionnerait pas mieux?

Code:

---

1 2 3

for( n in 1:(length(Tini) - 1) ) { ... }

---

Autre petite remarque, que représente z et k?

Enfin, je ne vois pas bien à quoi sert le calcul/l'incrémentation de t au sein de la boucle? Pourquoi ne pas tout simplement faire :

Code:

---

t<-(length(Tini))*10800

---

Bonne continuation :)


Cordialement,

A.D.

Super merci.
En effet j'avais du mal à écrire le critère de la boucle et cela semble fonctionner désormais.
Et puis oui pour l'incrément je n'y avais pas pensé.
Merci.