Premier cercle de R enfer. 0,1 = 0,3/3

Possible en double:
Numeric comparison difficulty in R Premier cercle de R enfer. 0,1 = 0,3/3

Bonjour à tous,

Selon "R Inferno" papier. Je suis maintenant dans le premier cercle de R enfer. C'est là que les païens s'attendent à 0,1 == 0,3/3. Le papier recommande d'utiliser la fonction all.equal pour de tels cas, cependant je dois vérifier les conditions "> =" ou "< =". Avec l'exemple actuel sur d'entre eux échouent:

> .1 >= .3/3 
[1] TRUE 
> .1 <= .3/3 
[1] FALSE

Y at-il une fonction similaire à tous.qualité qui vérifie les inégalités?

Merci,

Ilya

Source

2010-11-02 ilya

Ce n'est pas une erreur spécifique à R, mais une erreur commune provoquée par la représentation binaire des nombres à virgule flottante. Voir http://en.wikipedia.org/wiki/Floating_point#Accuracy_problems par exemple. – schnaader

Je sais pourquoi une erreur se produit. Je voudrais savoir s'il existe une fonction R capable de vérifier les inégalités ou si je dois en écrire une avec des points flottants en tête. – ilya

oui, je dois faire de meilleures recherches – ilya

Voir ces questions:

D'une manière générale, vous pouvez traiter en incluant un niveau de tolérance selon le second lien ci-dessus.

Source

2010-11-02 15:26:02 Shane

Voir le deuxième lien qui répond directement à la question. – Shane

Merci Shane, en deuxième lien, j'ai trouvé votre message avec la création d'un nouvel opérateur - exactement ce que je voulais. Je vous remercie!! – ilya

@ilya Heureux nous pourrions aider! – Shane

S'il vous plaît voir le R FAQ entry Why doesn't R think these numbers are equal et les références qui s'y trouvent.

Source

2010-11-02 15:20:25

C'est fascinant, mais pas une réponse à l'OP. – mob

Sûr est, ou vous lisez une autre entrée 7.31. L'exemple 'sqrt (2) * sqrt (2) == 2' apparaît comme FALSE et est le même problème ici: représentation de nombres à virgule flottante. –

Vous pourriez essayer l'utilisation judicieuse de zapsmall() qui semble donner le comportement que vous recherchez. Je ne sais pas si cela fonctionne dans toutes les situations. par exemple,

.1 >= zapsmall(.3/3) 
[1] TRUE 
> .1 <= zapsmall(.3/3) 
[1] TRUE

Source

2010-11-02 15:28:45 kmm

Le test principal de all.equal est de savoir si abs(x-y) < tolerance pour certaines valeurs x et y et une petite tolerance. Les tests d'inégalité équivalents vérifient:

x <= y:   x-y < tolerance 
x < y:   x-y < -tolerance 
x >= y:   x-y > -tolerance 
x > y:   x-y > tolerance

Source

2010-11-02 15:30:08 mob

Répondre

Questions connexes