covid19br / nowcasting

Repositório para análises de nowcasting
11 stars 1 forks source link

read.sivep seleiona campo de geocode com NA #20

Closed piklprado closed 4 years ago

piklprado commented 4 years ago

Debugando a read.sivep para município vi que a linhas 55-57 do codigo rea.sive.R onde há

f (escala == "municipio") {
      if(residentes)
          dados <- dados[dados$co_mun_res == as.numeric(geocode), ]

Está selecionando também linhas que têm NA no campo co_mun_res. Não esperava isso, e não entendi pq a indexação está se comportando desta maneira. Mas testei com duas alternativas que selecionaram apenas os geocodes e não os campos com NA. Seria na linha 57:

dados <- dados[dados$co_mun_res == as.numeric(geocode)&!is.na(dados$co_mun_res), ]

ou

dados <- dplyr::filter(dados, co_mun_res == as.numeric(geocode)) %>%as.data.frame()

Estou perdendo algo? Se é isso , seria o mesmo para as demais linha análogas de seleção. Aliás, acho bom verificar todas as linhas de seleção de campos que podem conter NA que pdoe estar acontecendo o mesmo problema.