Discussion :

[mungo32]

Bonjour,

Je n'arrive pas à tester convenablement le résultat d'une division d'un élément de matrice par un autre. Je m'explique…

Je crée trois matrices

Code :

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

1
2
3
4
 
cjmeasure.in <- matrix(rep(0,4*l),,4)				#create a vector to contain measures
minprojin <- matrix(rep(0,2*l),,2)
maxprojin <- matrix(rep(0,2*l),,2)

Elles sont donc en "double". Je remplis les deux premières colonnes de la matrice "cjmeasure.in" ainsi que la première colonne de "minprojin" et "maxprojin" à l'aide de résultats d'un calcul de programmation linéaire.

Ensuite, je remplis les deuxièmes colonnes de "minprojin" et "maxprojin" à l'aide d'un autre problème de programmation linéaire, et complète mes valeurs de "cjmeasure.in" à l'aide de division faites à partir de "minprojin" et "maxprojin". Typiquement :

Code :

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

1
2
3
 
cjmeasure.in[j,3] <- minprojin[j,1]/minprojin[j,2]
cjmeasure.in[j,4] <- maxprojin[j,1]/maxprojin[j,2]

Je souhaite enfin remplacer les éventuels valeurs égales à 1 dans "cjmeasure.in" par "NaN", de sorte qu'il reste de type "double".

Code :

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

1
2
3
4
5
 
if (cjmeasure.in[j,4] == 1)
{
	cjmeasure.in[j,4] <- NaN
}

Et là, alors que la valeur de la colonne 4 de la ligne j vaut 1 (affiché "1" et non "1.0000000" comme dans les autres colonnes, chose assez étonnante!), elle n'est pas remplacée par NaN. Lorsque je teste

Code :

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

cjmeasure.in[j,4] == 1

, la réponse est "FALSE".

Je ne comprend pas pourquoi. Quelqu'un pourrait-il m'aider/me guider?

Merci d'avance,

François

[mungo32]

Après investigations, je me rends compte que mon problème n'est que la conséquence du fait suivant :

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
 
> cjmeasure.in
           phi       pi super_phi super_pi
[1,] 1.0000000 1.444444  0.000000        0
[2,] 1.0000000 1.000000  0.000000        0
[3,] 0.8833333 1.250000  0.000000        0
[4,] 0.6769231 1.000000  1.023622        0
[5,] 1.0000000 1.485294  0.000000        0
[6,] 0.8461538 1.000000  0.000000        0
[7,] 1.0000000 1.181818  0.000000        0
[8,] 0.7416667 1.083333  0.000000        0
> cjmeasure.in[5,1] == 1
  phi 
FALSE
> cjmeasure.in[2,1] == 1
 phi 
TRUE

Dans ma matrice cjmeasure.in, R pense que la valeur de l'élément [5,1] n'est pas égale à 1.
Je ne comprends pas pourquoi. Tandis qu'il me dit que la valeur de l'élément [2,1] est bien égale à 1, elle.

Est-ce que cela peut guider une éventuelle réponse?

[mungo32]

Après investigations, je me rends compte que mon problème n'est que la conséquence du fait suivant :

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
 
> cjmeasure.in
           phi       pi super_phi super_pi
[1,] 1.0000000 1.444444  0.000000        0
[2,] 1.0000000 1.000000  0.000000        0
[3,] 0.8833333 1.250000  0.000000        0
[4,] 0.6769231 1.000000  1.023622        0
[5,] 1.0000000 1.485294  0.000000        0
[6,] 0.8461538 1.000000  0.000000        0
[7,] 1.0000000 1.181818  0.000000        0
[8,] 0.7416667 1.083333  0.000000        0
> cjmeasure.in[5,1] == 1
  phi 
FALSE
> cjmeasure.in[2,1] == 1
 phi 
TRUE

Dans ma matrice cjmeasure.in, R pense que la valeur de l'élément [5,1] n'est pas égale à 1.
Je ne comprends pas pourquoi. Tandis qu'il me dit que la valeur de l'élément [2,1] est bien égale à 1, elle.

Est-ce que cela peut guider une éventuelle réponse?

Nouvelle idée ce matin au réveil :"Puisque ce n'est pas égal à 1, ça vaut combien?"

Premier test : est-ce que la valeur est strictement supérieure à 1 : TRUE
Second : est-ce que la valeur est strictement inférieure à 1,0000001 : TRUE.
Je vous laisse imaginer la suite: ma valeur "cjmeasure.in[5,1]" tourne aux environs de 1.0000000000000002220446 (alors qu'elle devrait être égale à 1).

J'ai perdu beaucoup de temps, mais au moins maintenant je connais la fonction "round()" ainsi que "options(digits=22)"!

Bonne journée à tous,

François