Primer círculo de R infierno. 0,1 = 0,3/3

Duplicar posibles:
Numeric comparison difficulty in R Primer círculo de R infierno. 0,1 = 0,3/3

Hola a todos,

Según "R Inferno" papel. Estoy ahora en el primer círculo del infierno R Aquí es donde los paganos esperan 0.1 == 0.3/3. Paper recomienda utilizar la función all.equal para estos casos; sin embargo, debo verificar las condiciones "> =" o "< =". Con el ejemplo actual de ellos fracasan:

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

¿Hay una función similar a all.equal que comprueba las desigualdades?

Gracias,

Ilya

Fuente

2010-11-02 ilya

Eso no es un error específico de I, pero un error común que es causada por la representación binaria de los números de punto flotante. Ver http://en.wikipedia.org/wiki/Floating_point#Accuracy_problems por ejemplo. – schnaader

Sé por qué ocurre un error. Me gustaría saber si hay una función R que puede verificar las desigualdades o si necesito escribir una con los puntos flotantes en mente. – ilya

sí, necesito hacer mejores búsquedas – ilya

ver estas preguntas:

En términos generales, se puede hacer frente a este mediante la inclusión de un nivel de tolerancia de acuerdo con el segundo enlace anterior.

Fuente

2010-11-02 15:26:02 Shane

Vea el segundo enlace que responde directamente a la pregunta. – Shane

Gracias Shane, en el segundo enlace encontré tu publicación con la creación de un nuevo operador, exactamente lo que quería. ¡¡Gracias!! – ilya

@ilya ¡Me alegra que pudiéramos ayudar! – Shane

Por favor ver el R FAQ entry Why doesn't R think these numbers are equal y las referencias en él.

Fuente

2010-11-02 15:20:25

Esto es fascinante, pero no una respuesta para el OP. – mob

Claro está, o lees otra entrada 7.31. El ejemplo que hay 'sqrt (2) * sqrt (2) == 2' sale como FALSO y es el mismo problema aquí: representación de los números de coma flotante. –

Podría intentar el uso juicioso de zapsmall() que parece funcionar según el comportamiento que está buscando. No sé si esto funciona en todas las situaciones. por ejemplo,

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

Fuente

2010-11-02 15:28:45 kmm

La prueba principal de all.equal es si abs(x-y) < tolerance para algunos valores x y y y algunos pequeños tolerance. pruebas de desigualdad equivalentes comprobarían:

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

Fuente

2010-11-02 15:30:08 mob

Primer círculo de R infierno. 0,1 = 0,3/3

Respuesta

Cuestiones relacionadas