Discussion :

[enicnath]

Bonjour
Je souhaite tester l'égalité entre des valeurs, sur deux variables qui possèdent des NA.
Le problème est que je n'obtiens pas du tout ce que je veux. Exemple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
X <- data.frame(Id=c(1,2,3,4,5),V1=c(0,0,1,NA,1),V2=c(NA,0,0,1,1))
X[X$V1 != X$V2,]

Le résultat est le suivant :

Id V1 V2
NA NA NA NA
3 3 1 0
NA.1 NA NA NA

Si j'y met un which(...)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

X[which(X$V1 != X$V2),]

J'obtiens :

Id V1 V2
3 3 1 0

Ce n'est pas non plus ce que je désire. Je voudrais avoir

Id V1 V2
1 1 0 NA
3 3 1 0
4 4 NA 1

Bien sur, on peut aussi ajouter des contrôles sur les NA, de la façon suivante :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

X[which( (is.na(X$V1) & !is.na(X$V2)) | (is.na(X$V2) & !is.na(X$V1)) | (X$V1 != X$V2)),]

Mais n'y a-t-il pas plus simple ?

[xavier-Pierre]

bonjour ,
essaye

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

X[X$V1 != X$V2|is.na(X$V1 != X$V2),]

cdt

[enicnath]

c'est pas mal, mais pas tout à fait suffisant.
Exemple, j'ajoute une observation où les deux variables sont non renseignées

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
X <- data.frame(Id=c(1,2,3,4,5,6),V1=c(0,0,1,NA,NA,1),V2=c(NA,0,0,NA,1,1))
X[X$V1 != X$V2|is.na(X$V1 != X$V2),]

Ta solution produit ceci

Id V1 V2
1 1 0 NA
3 3 1 0
4 4 NA NA
5 5 NA 1

Et on voit que la ligne 4 contient NA et NA, alors que je ne voudrais pas qu'elle apparaisse.

[xavier-Pierre]

mais ça ne fonctionne pas si on a 2 NA

[enicnath]

mais ça ne fonctionne pas si on a 2 NA

Nos réponses se sont croisées.
Oui ça ne fonctionne pas dans le cas de 2 NA
Pas d'autre solution simple ?

[olivier.decourt]

Bonjour.
N'ayant aucune idée simple, j'ai construit une fonction bien bourrine et pas du tout vectorisée (c'est lamentable).

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
X <- data.frame(Id=c(1,2,3,4,5,6),V1=c(0,0,1,NA,1,NA),V2=c(NA,0,0,1,1,NA))
 
 
comparaison <- function(x,y){
 
  if (is.na(x) & is.na(y)) {
    res <- TRUE 
  } else {
    if (is.na(x) | is.na(y)) {
      res <- FALSE
    } else {
      res <- (x == y)
    }
  }
  return (res)
}
 
 
 
purrr::map2_lgl(X$V1, X$V2, comparaison)

Et puis j'ai eu honte et j'ai repris mes conditions en vectorisant.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
comparaison2 <- function (v1, v2) {
  (is.na(v1) & is.na(v2)) |
    (!is.na(v1) & !is.na(v2) & v1 == v2)
}
comparaison2(X$V1, X$V2)

Après je serais surpris que ce genre de fonction n'existe dans aucun package... mais Google ne m'en a pas proposé (j'ai cherché "elementwise vector comparison r function").