忽略语句中列的NA值

问题描述 投票:0回答:1

直到现在,我一直在使用中等大小的数据集进行职业调查(总计约200 mb),如果您要查看它,请使用以下数据:https://drive.google.com/drive/folders/1Od8zlOE3U3DO0YRGnBadFz804OUDnuQZ?usp=sharing

我有以下代码:

hogares<-read.csv("/home/servicio/Escritorio/TR_VIVIENDA01.CSV")
personas<-read.csv("/home/servicio/Escritorio/TR_PERSONA01.CSV")
datos<-merge(hogares,personas)

library(dplyr)


base<-tibble(ID_VIV=datos$ID_VIV, ID_PERSONA=datos$ID_PERSONA, EDAD=datos$EDAD, CONACT=datos$CONACT) 
base$maxage <- ave(base$EDAD, base$ID_VIV, FUN=max)
base$Condición_I<-case_when(base$CONACT==32 & base$EDAD>=60 ~ 1,
                            base$CONACT>=10 & base$EDAD>=60 & base$CONACT<=16 ~ 2,
                            base$CONACT==20 & base$EDAD>=60 | base$CONACT==31 & base$EDAD>=60 | (base$CONACT>=33 & base$CONACT<=35 & base$EDAD>=60) ~ 3)
base <- subset(base, maxage >= 60) 
base<- base %>%  group_by(ID_VIV) %>% mutate(Condición_V = if(n_distinct(Condición_I) > 1) 4 else Condición_I)
base$ID_VIV<-as.character(base$ID_VIV)           
base$ID_PERSONA<-as.character(base$ID_PERSONA)
base

最后得到:

# A tibble: 38,307 x 7
# Groups:   ID_VIV [10,499]
   ID_VIV      ID_PERSONA     EDAD CONACT maxage Condición_I Condición_V
   <chr>       <chr>         <int>  <int>  <int>       <dbl>       <dbl>
 1 10010000007 1001000000701    69     32     69           1           1
 2 10010000008 1001000000803    83     33     83           3           4
 3 10010000008 1001000000802    47     33     83          NA           4
 4 10010000008 1001000000801    47     10     83          NA           4
 5 10010000012 1001000001204     4     NA     60          NA           4
 6 10010000012 1001000001203     2     NA     60          NA           4
 7 10010000012 1001000001201    60     10     60           2           4
 8 10010000012 1001000001202    21     10     60          NA           4
 9 10010000014 1001000001401    67     32     67           1           4
10 10010000014 1001000001402    64     33     67           3           4

Condición_I列的值是每个人(行)劳动条件的代码,其中一些人共享房屋(这就是他们共享ID_VIV的原因),我只关心60岁或以上的人,所有不适用的人都是60岁以上的人,但我不在乎他们的状况(但我需要保留他们),我需要在列[[Condición_V的列中显示以下条件下的另一个值:

Condición_I == 1 ~ 1 Condición_I == 2 ~ 2 Condición_I == 3 ~ 3 Any combination of Condición_I ~ 4

这意味着,如果一个房子中的所有60个[和 __ yo个体都具有Condición_I== 1,那么Condición_V将为1

在存在x.e的情况下,代码3都是正确的。一个人C_I == 1,另一个人C_I == 3在同一屋子里,那么Condición_V将为4我希望得到这样的结果:小报:38,307 x 7# Groups: ID_VIV [10,499] ID_VIV ID_PERSONA EDAD CONACT maxage Condición_I Condición_V <chr> <chr> <int> <int> <int> <dbl> <dbl> 1 10010000007 1001000000701 69 32 69 1 1 2 10010000008 1001000000803 83 33 83 3 3 3 10010000008 1001000000802 47 33 83 NA 3 4 10010000008 1001000000801 47 10 83 NA 3 5 10010000012 1001000001204 4 NA 60 NA 2 6 10010000012 1001000001203 2 NA 60 NA 2 7 10010000012 1001000001201 60 10 60 2 2 8 10010000012 1001000001202 21 10 60 NA 2 9 10010000014 1001000001401 67 32 67 1 4 10 10010000014 1001000001402 64 33 67 3 4

我知道我的错误是在:
`#base<- base %>%  group_by(ID_VIV) %>% mutate(Condición_V = if(n_distinct(Condición_I) > 1) 4 else` Condición_I)

是否有一种方法可以使用该行代码来忽略NA值,或者它是我的最佳选择,所以我不必以我尝试的方式来做,因此,我们将不胜感激任何其他方式或帮助!

r if-statement dplyr tibble
1个回答
0
投票
library(dplyr) base %>% group_by(ID_VIV) %>% mutate(Condición_V = if(n_distinct(na.omit(Condición_I)) > 1) 4 else na.omit(Condición_I)[1]) # A tibble: 10 x 7 # Groups: ID_VIV [4] # ID_VIV ID_PERSONA EDAD CONACT maxage Condición_I Condición_V # <chr> <chr> <int> <int> <int> <int> <dbl> # 1 10010000007 1001000000701 69 32 69 1 1 # 2 10010000008 1001000000803 83 33 83 3 3 # 3 10010000008 1001000000802 47 33 83 NA 3 # 4 10010000008 1001000000801 47 10 83 NA 3 # 5 10010000012 1001000001204 4 NA 60 NA 2 # 6 10010000012 1001000001203 2 NA 60 NA 2 # 7 10010000012 1001000001201 60 10 60 2 2 # 8 10010000012 1001000001202 21 10 60 NA 2 # 9 10010000014 1001000001401 67 32 67 1 4 #10 10010000014 1001000001402 64 33 67 3 4

数据

base <- structure(list(ID_VIV = c("10010000007", "10010000008", "10010000008", 
"10010000008", "10010000012", "10010000012", "10010000012", "10010000012", 
"10010000014", "10010000014"), ID_PERSONA = c("1001000000701", 
"1001000000803", "1001000000802", "1001000000801", "1001000001204", 
"1001000001203", "1001000001201", "1001000001202", "1001000001401", 
"1001000001402"), EDAD = c(69L, 83L, 47L, 47L, 4L, 2L, 60L, 21L, 
67L, 64L), CONACT = c(32L, 33L, 33L, 10L, NA, NA, 10L, 10L, 32L, 
33L), maxage = c(69L, 83L, 83L, 83L, 60L, 60L, 60L, 60L, 67L, 
67L), Condición_I = c(1L, 3L, NA, NA, NA, NA, 2L, NA, 1L, 3L
)), row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", 
"10"), class = "data.frame")
© www.soinside.com 2019 - 2024. All rights reserved.