сравнение фрейма данных - процент несовпадающих случаев - PullRequest
0 голосов
/ 06 марта 2020

У меня есть два человека, которые кодируют одни и те же переменные для психологического теста (более 400 переменных), и мне нужно сравнить наборы данных. Мне нужны два результата:

  1. Мне нужно видеть только указанные c случаи с несовпадением и
  2. В качестве конечного результата мне нужен процент несоответствия для каждой переменной.

вот что я подразумеваю под «процентом несоответствия по переменной»:

A <- tibble( ID = c(1:10),  
         v1 = rep(1),
         v2 = rep(2), 
         v3 = rep(3))


B <- tibble( ID = c(1:10),  
             v1 = c(1,1,1,1,10,1,1,1,1,1),
             v2 = c(30,2,2,2,51,2,2,2,2,40), 
             v3 = c(3,3,3,3,3,3,3,65,3,90))

A;B

# A tibble: 10 x 4
      ID    v1    v2    v3
   <int> <dbl> <dbl> <dbl>
 1     1     1     2     3
 2     2     1     2     3
 3     3     1     2     3
 4     4     1     2     3
 5     5     1     2     3
 6     6     1     2     3
 7     7     1     2     3
 8     8     1     2     3
 9     9     1     2     3
10    10     1     2     3
# A tibble: 10 x 4
      ID    v1    v2    v3
   <int> <dbl> <dbl> <dbl>
 1     1     1    30     3
 2     2     1     2     3
 3     3     1     2     3
 4     4     1     2     3
 5     5    10    51     3
 6     6     1     2     3
 7     7     1     2     3
 8     8     1     2    65
 9     9     1     2     3
10    10     1    40    90

Как сравнить набор данных A с набором данных B, чтобы получить такой результат:

result<- tibble(variables = c("v1", "v2", "v3"),
                 n.mismatch = c(1,3,2),
                 percentage.mismatch = c(0.10, 0.30, 0.20))
result
# A tibble: 3 x 3
  variables n.mismatch percentage.mismatch
  <chr>          <dbl>               <dbl>
1 v1                 1                 0.1
2 v2                 3                 0.3
3 v3                 2                 0.2

1 Ответ

1 голос
/ 06 марта 2020

Мы можем использовать Map для сравнения значений столбца.

as.data.frame(Map(function(x, y) {
   inds <- x != y
   c(n.mismatch = sum(inds), percentage.mismatch = mean(inds))
}, A[-1], B[-1]))

#                    v1  v2  v3
#n.mismatch          1.0 3.0 2.0
#percentage.mismatch 0.1 0.3 0.2

Аналогично, в tidyverse мы можем использовать map2

purrr::map2_df(A[-1], B[-1], ~{
   inds = .x != .y
   tibble::tibble(n.mismatch = sum(inds), percentage.mismatch = mean(inds))
 }, .id = "variables")

#  variables n.mismatch percentage.mismatch
#  <chr>          <int>               <dbl>
#1 v1                 1                 0.1
#2 v2                 3                 0.3
#3 v3                 2                 0.2
...