Problème dans la gestion des arrondis sur R

Bonjour à tous et à toutes :D

je n'ai pas trouvé de réponse à mes questions sur le forum.. J'espère ne pas faire de doublons en tout cas :oops::oops:

J'aurais une petite question sur R à propos de comment il gère les arrondis, les décimales, etc.

Je suis en train de faire des calculs qui demandent beaucoup de précision... et je vois que R bugue à ce niveau là.

Par exemple si j'écris, pour un certain p et x que j'ai calculé :
(p-x-x)-(p-2*x)

R me donne : -5.551115e-17
(alors que ça devrait être parfaitement 0)

Ou encore, si j'écris (pour d'autres valeurs de x, y et p):
x + y + p - x -y == p

Il me répond : FALSE

(il ne va pas me donner ça pour n'importe quelle valeur de p, x et y, mais dans certains cas, si !)

Bref, je ne sais pas ce qui se passe, mais cette imprécision rend mes calculs inexacts. :calim2:

Je ne sais pas si je suis claire... J'imagine que R fait des arrondis et que c'est ça qui rend mes calculs faux... :(

Est-ce que vous pourriez m'aider ?

Ou alors je passe à un autre logiciel ? Python ? (que je ne connais pas... et j'ai beaucoup de retard dans ma thèse... donc je ne sais pas si c'est une bonne idée...)

Merci d'avance de votre aide ! :D

En effet R ne semble pas au top au niveau de la précision.
En cas de doute tu peux comparer les valeurs obtenues avec celle de la précision de ta machine.
Pour l'obtenir :

Code:

---

1 2 3

.Machine$double.eps # [1] 2.220446e-16

---

La valeur obtenue sera différente en fonction de la machine.
Dans ton cas '-5.551115e-17' est probablement inférieur à la précision de ta machine, du coup tu peux le considérer comme égale à 0 même si il n'affiche pas ce résultat.

Merci beaucou MamounMob !

J'imagine qu'il n'y a pas moyen d'augmenter la précision de R ?

Est-ce vous savez si le logiciel Python (assez facile d'utilisation d'après ce que j'ai compris) est plus précis ?

Merci d'avance !

je ne sais si il est possible de modifier la précision ou si Python est plus précis. Mais peut que le langage C est meilleure pour le coup, il faudrait creuser la question.
Bonne chance

ok, merci :)