Notebook para criação de tabela de indicadores da PNS - módulo I 2013 Cobertura de Plano de Saúde¶
Bibliotecas Utilizadas¶
In [3]:
#Lendo pacotes necessários
library(survey)
library(ggplot2)
library(dplyr)
library(foreign)
library(forcats)
library(tidyverse)
Carregando microdados da PNS¶
In [4]:
#Carregando banco de dados para R versão 3.5.0 ou superior
load("<coloque aqui o caminho para o arquivo dos microdados formato RDATA PNS 2013>")
#conferindo as dimensões (número de linhas e colunas)
dim("<Coloque aqui o nome do arquivo RDATA PNS 2013>")
Definição do peso e filtragem de respondentes do questionário¶
In [5]:
#Selecionando registros válidos e calculando peso amostral - summary de verificação
pns2013.1<- <Coloque aqui o nome do arquivo RDATA> %>% mutate(peso_moradores=((V00281*205546/199551444)))
pns2013.1<-pns2013.1 %>% filter(!is.na(peso_moradores))
summary(pns2013.1$peso_moradores)
Criação de variáveis dos indicadores¶
In [6]:
#Desfechos - Indicadores
#Posse de plano de saúde médico ou odontológico particular - I001P
pns2013.1 <- pns2013.1 %>% mutate(I001P= ifelse(I001==1, 1,2))
pns2013.1$I001P<-factor(pns2013.1$I001P, levels=c(1,2), labels=c("Sim", "Não"))
summary(pns2013.1$I001P)
##Posse de plano de saúde médico particular - I002P
pns2013.1 <- pns2013.1 %>% mutate(I002P = ifelse(I001==2 ,2,
ifelse(I001==1 & I002==1 & I003==1,2,
ifelse(I001==1 & I002==1 & I003==2,1,
ifelse(I001==1 & I002>1 & I003==1,1,
ifelse(I001==1 & I002>1 & I003==2,1,2))))))
pns2013.1$I002P<-factor(pns2013.1$I002P, levels=c(1,2), labels=c("Sim","Não"))
summary(pns2013.1$I002P)
#Posse de plano de saúde odontológico particular - I003P - Plano amostral distinto
pns2013.1 <- pns2013.1 %>% mutate(I003P= if_else(I003==1,1,2))
pns2013.1$I003P<-factor(pns2013.1$I003P, levels=c(1,2), labels=c("Sim", "Não"))
summary(pns2013.1$I003P)
Definições de abrangências¶
Situação urbana ou rural¶
In [7]:
#Situação Urbano ou Rural
pns2013.1 <- pns2013.1 %>% rename(Situação_Urbano_Rural=V0026)
pns2013.1$Situação_Urbano_Rural<-factor(pns2013.1$Situação_Urbano_Rural, levels=c(1,2), labels=c("Urbano", "Rural"))
summary(pns2013.1$Situação_Urbano_Rural)
Sexo¶
In [8]:
#Sexo
pns2013.1 <- pns2013.1 %>% rename(Sexo=C006)
pns2013.1$Sexo<-factor(pns2013.1$Sexo, levels=c(1,2), labels=c("Masculino", "Feminino"))
summary(pns2013.1$Sexo)
UF¶
In [9]:
#Estados - UFs
pns2013.1 <- pns2013.1 %>% rename(Unidades_da_Federação=V0001)
pns2013.1$Unidades_da_Federação<-factor(pns2013.1$Unidades_da_Federação, levels=c(11,12,13,14,15,16,17,21,22,23,24,25,26,27,28,29,31,32,33,35,41,42,43,50,51,52,53),
label=c("Rondônia","Acre","Amazonas","Roraima","Pará","Amapá","Tocantins","Maranhão","Piauí","Ceará",
"Rio Grande do Norte","Paraíba","Pernambuco","Alagoas","Sergipe","Bahia",
"Minas Gerais","Espírito Santo","Rio de Janeiro","São Paulo",
"Paraná","Santa Catarina","Rio Grande do Sul",
"Mato Grosso do Sul","Mato Grosso","Goiás","Distrito Federal"))
summary(pns2013.1$Unidades_da_Federação)
Grandes Regiões¶
In [10]:
#Grandes Regiões
pns2013.1 <- pns2013.1 %>%
mutate(GrandesRegioes = fct_collapse(Unidades_da_Federação,
`Norte` = c("Rondônia","Acre","Amazonas","Roraima","Pará", "Amapá","Tocantins"),
`Nordeste` = c("Maranhão", "Piauí", "Ceará", "Rio Grande do Norte", "Paraíba","Pernambuco", "Alagoas","Sergipe","Bahia"),
`Sudeste` = c("Minas Gerais", "Espírito Santo","Rio de Janeiro", "São Paulo"),
`Sul` = c("Paraná", "Santa Catarina", "Rio Grande do Sul"),
`Centro-Oeste`= c("Mato Grosso do Sul","Mato Grosso", "Goiás", "Distrito Federal")))
summary(pns2013.1$GrandesRegioes)
Capital¶
In [11]:
#Capital
pns2013.1<- pns2013.1 %>% mutate(Capital= fct_collapse(Unidades_da_Federação,
`Porto Velho`= "Rondônia",
`Boa Vista`= "Roraima",
`Rio Branco`= "Acre",
`Manaus` = "Amazonas",
`Belém` = "Pará" ,
`Macapá`= "Amapá",
`Palmas` = "Tocantins",
`São Luís` = "Maranhão",
`Teresina`= "Piauí" ,
`Fortaleza`= "Ceará",
`Natal`= "Rio Grande do Norte",
`João Pessoa`= "Paraíba",
`Recife`= "Pernambuco",
`Maceió`= "Alagoas",
`Aracaju`= "Sergipe",
`Salvador`= "Bahia",
`Belo Horizonte`= "Minas Gerais",
`Vitória`= "Espírito Santo",
`Rio de Janeiro`= "Rio de Janeiro",
`São Paulo`= "São Paulo",
`Curitiba`= "Paraná",
`Florianópolis`= "Santa Catarina",
`Porto Alegre`= "Rio Grande do Sul",
`Campo Grande`= "Mato Grosso do Sul",
`Cuiabá`= "Mato Grosso",
`Goiânia` = "Goiás",
`Brasília`= "Distrito Federal"))
summary(pns2013.1$Capital)
Faixa Etária¶
In [12]:
#Faixas Etárias
pns2013.1 <- pns2013.1 %>% mutate(faixa_idade = ifelse(C008<=0,1,
ifelse(C008%in% 1:4, 2,
ifelse(C008%in% 5:14, 3,
ifelse(C008%in% 15:24, 4,
ifelse(C008%in% 25:34, 5,
ifelse(C008%in% 35:44, 6,
ifelse(C008%in% 45:59, 7,8))))))))
pns2013.1$faixa_idade<-factor(pns2013.1$faixa_idade, levels=c(1,2,3,4,5,6,7,8), labels=c("Menor de 1 ano","1 a 4 anos","5 a 14 anos","15 a 24 anos","25 a 34 anos","35 a 44 anos","45 a 59 anos","60 anos ou mais"))
summary(pns2013.1$faixa_idade)
Raça¶
In [13]:
#Raça
pns2013.1 <- pns2013.1 %>% mutate(Raca= ifelse(C009==1, 1,
ifelse(C009==2, 2,
ifelse(C009==4, 3, 9))))
pns2013.1$Raca<-factor(pns2013.1$Raca, levels=c(1,2,3),labels=c("Branca", "Preta", "Parda"))
summary(pns2013.1$Raca)
Renda per capita¶
In [14]:
#Rendimento domiciliar per capita - Outro plano amostral
pns2013.1 <- pns2013.1 %>% mutate(VDF004= ifelse(VDF003 %in% 0:170, 1,
ifelse(VDF003%in% 170:339, 2,
ifelse(VDF003%in% 340:678, 3,
ifelse(VDF003%in% 679:1356, 4,
ifelse(VDF003%in% 1357:2034, 5,
ifelse(VDF003%in% 2035:3390, 6,
ifelse(VDF003=="NA","NA",7))))))))
pns2013.1 <- pns2013.1 %>% mutate(rend_per_capita= ifelse(VDF004 %in% 1:2, 1,
ifelse(VDF004==3, 2,
ifelse(VDF004==4, 3,
ifelse(VDF004==5, 4,
ifelse(VDF004=="NA","NA",5))))))
pns2013.1$rend_per_capita<-factor(pns2013.1$rend_per_capita, levels=c(1,2,3,4,5), labels=c("Até 1/2 SM","1/2 até 1 SM","1 até 2 SM",
"2 até 3 SM","Mais de 3 SM"))
summary(pns2013.1$rend_per_capita)
Criando indicadores¶
Filtrando base de indicadores¶
In [15]:
#Selecionando variáveis para cálculo de indicadores no survey
pns2013survey<- pns2013.1 %>% select("V0024","UPA_PNS","peso_moradores", "I001P","I002P","I003P", "V0031","C009","VDF003",
"Situação_Urbano_Rural","Sexo","Unidades_da_Federação", "GrandesRegioes",
"Capital","faixa_idade", "Raca","rend_per_capita")
summary(pns2013survey)
Exporta tabela filtrada com os dados específicos - Módulo I 2013¶
In [17]:
#Salvando csv para cálculo de indicadores no survey
diretorio_saida <- "<coloque aqui o diretório>"
write.csv(pns2013survey, file.path(diretorio_saida, "pns2013Isurvey.csv"))
Cria plano amostral complexo¶
In [18]:
#survey design
desPNSI=svydesign(id=~UPA_PNS, strat=~V0024, weight=~peso_moradores, nest=TRUE, data=pns2013survey)
desPNSIC=subset(desPNSI, V0031==1)
desPNSIR=subset(desPNSI, C009!=9)
desPNSIRE=subset(desPNSI, VDF003!=9)
In [19]:
#survey design - I003P
desPNSI003P=subset(desPNSI, !is.na(I003P))
desPNSI003P_R=subset(desPNSI, !is.na(I003P) & C009!=9)
desPNSI003P_C=subset(desPNSI, V0031==1 & !is.na(I003P))
desPNSI003P_RE=subset(desPNSI, VDF003!=9 & !is.na(I003P))
Criação da tabela de indicadores¶
Essa tabela é responsável por unir os indicadores no formato do painel de indicadores
In [20]:
matrizIndicadores = data.frame()
Definição de variáveis para iteração dos indicadores¶
In [21]:
ListaIndicadores = c(~I001P,~I002P,~I003P)
ListaIndicadoresTexto = c("I001P","I002P","I003P")
ListaDominios = c(~Sexo,~Raca,~rend_per_capita,~faixa_idade,~Situação_Urbano_Rural,~Unidades_da_Federação,~GrandesRegioes,~Capital)
ListaDominiosTexto = c("Sexo","raça","rend_per_capita","fx_idade","urb_rur","uf","região","capital")
ListaTotais = c('Brasil','Capital')
Ano <- "2013"
Preenchendo a tabela de indicadores¶
In [22]:
#Cálculo dos indicadores usando o pacote survey
i <- 0
for( indicador in ListaIndicadores){
i <- i + 1
j <- 1
for (dominio in ListaDominios){
if (ListaDominiosTexto[j]=="capital"){
if(ListaIndicadoresTexto[i] == "I003P"){
dataframe_indicador<-svyby( indicador , dominio , desPNSI003P_C , svymean,vartype= c("ci","cv"))
}else{
dataframe_indicador<-svyby( indicador , dominio , desPNSIC , svymean,vartype= c("ci","cv"))
}
}else if (ListaDominiosTexto[j]=="raça"){
if(ListaIndicadoresTexto[i] == "I003P"){
dataframe_indicador<-svyby( indicador , dominio , desPNSI003P_R , svymean,vartype= c("ci","cv"))
}else{
dataframe_indicador<-svyby( indicador , dominio , desPNSIR , svymean,vartype= c("ci","cv"))
}
}else if (ListaDominiosTexto[j]=="rend_per_capita"){
if(ListaIndicadoresTexto[i] == "I003P"){
dataframe_indicador<-svyby( indicador , dominio , desPNSI003P_RE , svymean,vartype= c("ci","cv"))
}else{
dataframe_indicador<-svyby( indicador , dominio , desPNSIRE , svymean,vartype= c("ci","cv"))
}
}else {
if(ListaIndicadoresTexto[i] == "I003P"){
dataframe_indicador<-svyby( indicador , dominio , desPNSI003P , svymean,vartype= c("ci","cv"))
}else{
dataframe_indicador<-svyby( indicador , dominio , desPNSI , svymean,vartype= c("ci","cv"))
}
}
dataframe_indicador<-data.frame(dataframe_indicador)
colnames(dataframe_indicador) <- c("abr_nome","Sim","Nao","LowerS","LowerN","UpperS","UpperN","cvS","cvN")
dataframe_indicador$Indicador <- ListaIndicadoresTexto[i]
dataframe_indicador$abr_tipo <- ListaDominiosTexto[j]
dataframe_indicador$Ano <- Ano
dataframe_indicador <- dataframe_indicador %>% select("abr_tipo","abr_nome","Ano","Indicador","Sim","LowerS","UpperS","cvS")
matrizIndicadores <-rbind(matrizIndicadores,dataframe_indicador)
j <- j + 1
}
}
Criando a tabela pela abrangência total¶
In [23]:
matriz_totais <- data.frame()
Preenchendo a tabela com as abrangencia Brasil e total das capitais¶
In [24]:
i=0
for(indicador in ListaIndicadores){
i <- i+1
for(total in ListaTotais){
dataframe_indicador <- data.frame()
dataframe_indicador_S <- data.frame()
if (total == "Capital"){
if(ListaIndicadoresTexto[i] == "I003P"){
dataframe_indicador <- svymean(indicador,desPNSI003P_C)
}else{
dataframe_indicador <- svymean(indicador,desPNSIC)
}
} else {
if(ListaIndicadoresTexto[i] == "I003P"){
dataframe_indicador <- svymean(indicador,desPNSI003P)
}else{
dataframe_indicador <- svymean(indicador,desPNSI)
}
}
intervalo_confianca <- confint(dataframe_indicador)
coeficiente_variacao <- cv(dataframe_indicador)
dataframe_indicador <- cbind(data.frame(dataframe_indicador),data.frame(intervalo_confianca))
dataframe_indicador <- cbind(data.frame(dataframe_indicador),data.frame(coeficiente_variacao))
dataframe_indicador <- dataframe_indicador %>%
select('mean','X2.5..','X97.5..',coeficiente_variacao)
dataframe_indicador_S <- dataframe_indicador %>%
slice(1)
colnames(dataframe_indicador_S) <- c('Sim','LowerS','UpperS', 'cvS')
dataframe_indicador_S <- dataframe_indicador_S %>%
select('Sim','LowerS','UpperS','cvS')
dataframe_indicador_S$Indicador <- ListaIndicadoresTexto[i]
dataframe_indicador_S$abr_tipo <- "total"
dataframe_indicador_S$abr_nome <- total
dataframe_indicador_S$Ano <- Ano
dataframe_indicador_S <- dataframe_indicador_S %>%
select("abr_tipo","abr_nome","Ano","Indicador","Sim","LowerS","UpperS",'cvS')
matriz_totais <-rbind(matriz_totais,dataframe_indicador_S)
}
}
Unindo tabela de indicadores e de totais¶
In [25]:
matrizIndicadores<-rbind(matrizIndicadores,matriz_totais)
Visualizando a tabela de indicadores¶
In [26]:
matrizIndicadores
Exportando tabela de indicadores calculados - Módulo I 2013¶
In [28]:
diretorio_saida <- "<coloque aqui o diretório>"
write.table(matrizIndicadores,file=paste0(diretorio_saida, Indicadores_2013I_R.csv",sep = ";",dec = ",",row.names = FALSE)