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 : Sélectionner tout - Visualiser dans une fenêtre à part
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 : Sélectionner tout - Visualiser dans une fenêtre à part
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.