Discussion :

[Florence Magnin]

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.

[A. D.]

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 :

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

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 :

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

t<-(length(Tini))*10800

Bonne continuation


Cordialement,

A.D.

[Florence Magnin]

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.