Riassegnate == e! = A isTRUE (all.equal())?

Question

A previous post mi ha spinto a postare questa domanda. Sembrerebbe come una best-practice per riassegnare == a isTRUE(all.equal()) (e != a !isTRUE(all.equal()). Mi chiedo se gli altri fanno in pratica? Ho appena realizzato che io uso == e != di fare uguaglianza numerica per tutta la mia base di codice. La mia prima reazione era che ho bisogno di fare un full-scrub e convertirlo in all.equal. Ma in realtà, ogni volta che uso == e != voglio testare l'uguaglianza (indipendentemente dal tipo di dati) .Infatti, non sono sicuro di cosa verrebbero testate queste operazioni per altro che l'uguaglianza, sono sicuro che mi manca qualche concetto qui. Qualcuno può illuminarmi? L'unico argomento che vedo contro questo approccio è che in alcuni casi due numeri non identici sembreranno identici a causa della tolleranza di all.equal. Ma ci viene detto che due numeri che sono in effetti identici potrebbero non passare identical() a causa di come sono memorizzati in memoria. Quindi, in realtà, qual è il punto di non default di all.equal?

Answer 1

Come ha fatto @joran, si verificheranno problemi in virgola mobile con == e != praticamente in qualsiasi altra lingua. Un aspetto importante di essi in R è la parte di vettorizzazione.

Sarebbe molto meglio definire una nuova funzione almostEqual, fuzzyEqual o simile. È un peccato che non ci sia una tale funzione di base. all.equal non è molto efficiente poiché gestisce tutti i tipi di oggetti e restituisce una stringa descrivendo la differenza quando principalmente si desidera solo TRUE o FALSE.

Ecco un esempio di tale funzione. È vettorizzato come ==.

almostEqual <- function(x, y, tolerance=1e-8) { 
    diff <- abs(x - y) 
    mag <- pmax(abs(x), abs(y)) 
    ifelse(mag > tolerance, diff/mag <= tolerance, diff <= tolerance) 
} 

almostEqual(1, c(1+1e-8, 1+2e-8)) # [1] TRUE FALSE 

... Si tratta di circa 2 volte più veloce di all.equal per i valori scalari, e molto più veloce con i vettori.

x <- 1 
y <- 1+1e-8 
system.time(for(i in 1:1e4) almostEqual(x, y)) # 0.44 seconds 
system.time(for(i in 1:1e4) all.equal(x, y)) # 0.93 seconds