Problème de test sur valeurs de matrices

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:

---

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:

---

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:

---

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:

---

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

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

Code:

---

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?

Citation:

---

Envoyé par mungo32

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

Code:

---

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