06_consumo_interno_importado_WIOD.R

#Consumo das Famílias Interno e Importado - a partir da WIOD
library(data.table)


consumobr <- function(ano = 2003, pasta = "sourcedata/wiots/") {
  if (ano > 2014) { fonte <- paste0(pasta,"WIOT2014_October16_ROW.rds")}
  else {fonte <- paste0(pasta,"WIOT",ano,"_October16_ROW.rds")}
  

consumofam2003 <- readRDS(fonte)

consumobrasfam2003 <- consumofam2003%>%select(IndustryCode,Country,BRA57)

consumointernobr2003 <- consumobrasfam2003%>%filter(Country == "BRA")%>%select(-Country)
names(consumointernobr2003)[2] <- "interno"

consumoimportadobr2003 <- consumobrasfam2003%>%
  filter(Country != "BRA")%>%
  group_by(IndustryCode)%>%
  summarize(importado = sum(BRA57))%>%
  filter(IndustryCode %in% consumointernobr2003$IndustryCode)

ciM <- consumointernobr2003%>%left_join(consumoimportadobr2003)

ciM[ciM == 0] <- 0.0000000000000001

ciM[,2:3] <- as.data.frame(prop.table(as.array(as.matrix(ciM[,-1])),1))

ciM$ano <- ano

ciM
}


consumos <- rbindlist(lapply(c(2003,2009,2018),consumobr))


consumos$IndustryCode <- substr(consumos$IndustryCode,1,3)



consumos$traded <- consumos$importado >= 0.1





write_csv2(consumos,"dados/proporcao-consumo-importado-geral-wiod-2003-2009-2014-8.csv")







###ABANDONADO - FEITA A TRADUÇÃO MANUALMENTE E JÁ EMBUTIDA EM GASTOS_SCN
# 
# tradisic <- read_xls("sourcedata/tradutores/CNAE20_Correspondencia_Cnae20xIsic4.xls",skip=3)[-1,-4]
# names(tradisic) <- c("cnae2","desc_cnae","isic4","desc_isic","obs")
# 
# tradisic <- tradisic %>%mutate(cnaeag = ifelse(is.na(as.numeric(cnae2)),cnae2,""),
#   isicsec = ifelse(is.na(as.numeric(isic4)),isic4,""),
#                                )
# 
# 
# for (i in 1:nrow(tradisic)) {
#   c <- length(tradisic)
#   if(tradisic[i,c]==""){
#     tradisic[i,c] <- tradisic[i-1,c]
#   }
#   if(tradisic[i,"isic4"] != tradisic[i,c]) {
#   tradisic[i,"isic4"] <- paste0(tradisic[i,c],tradisic[i,"isic4"])
#   tradisic[i,"cnaeag"] <- paste0(tradisic[i,c],tradisic[i,"cnae2"])
#   }
# }
# 
# tradisicrs <- tradisic[nchar(tradisic$isic4)< 4,]
# 
# 
# 
# 
# consumos <- consumos%>%left_join(tradisicrs,by = c("IndustryCode"= "isic4"))
# 
# 
# unique(consumos$cnae2)
# 
# 
# 
# ###Não suficiente sigamos
# cnaescn <- read_xls("sourcedata/tradutores/cnae_10_ativ147_ativ55.xls", skip=2)
# 
# 
# for (i in nrow(cnaescn):1) {
#   if(is.na(cnaescn[i,]$`Descrição Atividade 147`)) {
#     cnaescn[i,]$`Descrição Atividade 147` <- cnaescn[i+1,]$`Descrição Atividade 147`
#     cnaescn[i,]$`Descrição Atividade 56` <- cnaescn[i+1,]$`Descrição Atividade 56`
#   }
# }
# 
# cnaescna <- cnaescn[!is.na(cnaescn$DENOMINAÇÃO)&(!is.na(cnaescn$SEÇÃO)|!is.na(cnaescn$DIVISÃO)),]
# 
# cnaescna$codigo <- ""
# 
# for (i in 1:nrow(cnaescna)) {
#   if(is.na(cnaescna[i,1])) {
#     cnaescna[i,1] <- cnaescna[i-1,1]
#     cnaescna[i,"codigo"] <- paste0(cnaescna[i,1],cnaescna[i,2])
#   }
#   if(is.na(cnaescna[i,2])) {
#     cnaescna[i,"codigo"] <- cnaescna[i,1]
#   }
#   
# }
# 
# cnaescna <- cnaescna%>%select(codigo,DENOMINAÇÃO,`Descrição Atividade 147`,`Descrição Atividade 56`)