{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "![](https://github.com/bigdata-icict/ETL-Dataiku-DSS/raw/master/tutoriais/pcdas_1.5.png)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Notebook para criação de tabela de indicadores da PNS - S 2019 Pré-natal - Parte 1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Bibliotecas Utilizadas"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Lendo pacotes necessários\n",
    "library(survey)\n",
    "library(ggplot2)\n",
    "library(dplyr)\n",
    "library(tictoc)\n",
    "library(foreign)\n",
    "library(forcats)\n",
    "library(tidyverse)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Carregando microdados da PNS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>293726</li><li>1087</li></ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 293726\n",
       "\\item 1087\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 293726\n",
       "2. 1087\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "[1] 293726   1087"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Carregando banco de dados para R versão 3.5.0 ou superior\n",
    "load(\"<coloque aqui o caminho para o arquivo dos microdados formato RDATA PNS 2019>\")\n",
    "\n",
    "#conferindo as dimensões (número de linhas e colunas)\n",
    "dim(\"<Coloque aqui o nome do arquivo RDATA PNS 2019>\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Definição do peso e filtragem de respondentes do questionário"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "    Min.  1st Qu.   Median     Mean  3rd Qu.     Max. \n",
       " 0.00562  0.26621  0.54401  1.00000  1.12765 61.09981 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Selecionando registros válidos e calculando peso amostral - summary de verificação\n",
    "pns2019.1<-  <Coloque aqui o nome do arquivo RDATA>  %>% filter(V0025A==1) \n",
    "pns2019.1<-pns2019.1 %>% mutate(peso_morador_selec=((V00291*(90846/168426190))))\n",
    "pns2019.1<-pns2019.1 %>% filter(!is.na(peso_morador_selec))\n",
    "summary(pns2019.1$peso_morador_selec)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Criação de variáveis dos indicadores"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Sim</dt><dd>2734</dd><dt>Não</dt><dd>84</dd><dt>NA's</dt><dd>88028</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Sim] 2734\n",
       "\\item[Não] 84\n",
       "\\item[NA's] 88028\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Sim\n",
       ":   2734Não\n",
       ":   84NA's\n",
       ":   88028\n",
       "\n"
      ],
      "text/plain": [
       "  Sim   Não  NA's \n",
       " 2734    84 88028 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Sim</dt><dd>2657</dd><dt>Não</dt><dd>77</dd><dt>NA's</dt><dd>88112</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Sim] 2657\n",
       "\\item[Não] 77\n",
       "\\item[NA's] 88112\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Sim\n",
       ":   2657Não\n",
       ":   77NA's\n",
       ":   88112\n",
       "\n"
      ],
      "text/plain": [
       "  Sim   Não  NA's \n",
       " 2657    77 88112 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Sim</dt><dd>2356</dd><dt>Não</dt><dd>378</dd><dt>NA's</dt><dd>88112</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Sim] 2356\n",
       "\\item[Não] 378\n",
       "\\item[NA's] 88112\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Sim\n",
       ":   2356Não\n",
       ":   378NA's\n",
       ":   88112\n",
       "\n"
      ],
      "text/plain": [
       "  Sim   Não  NA's \n",
       " 2356   378 88112 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Sim</dt><dd>1822</dd><dt>Não</dt><dd>251</dd><dt>NA's</dt><dd>88773</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Sim] 1822\n",
       "\\item[Não] 251\n",
       "\\item[NA's] 88773\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Sim\n",
       ":   1822Não\n",
       ":   251NA's\n",
       ":   88773\n",
       "\n"
      ],
      "text/plain": [
       "  Sim   Não  NA's \n",
       " 1822   251 88773 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Sim</dt><dd>2056</dd><dt>Não</dt><dd>678</dd><dt>NA's</dt><dd>88112</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Sim] 2056\n",
       "\\item[Não] 678\n",
       "\\item[NA's] 88112\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Sim\n",
       ":   2056Não\n",
       ":   678NA's\n",
       ":   88112\n",
       "\n"
      ],
      "text/plain": [
       "  Sim   Não  NA's \n",
       " 2056   678 88112 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Desfechos - Indicadores\n",
    "\n",
    "# 1. Proporção de mulheres que realizaram pré-natal - S001P.\n",
    "pns2019.1$S001P <- NA\n",
    "pns2019.1$S001P[pns2019.1$C006==2 & pns2019.1$C008>=18 & pns2019.1$S068>0] <- 2\n",
    "pns2019.1$S001P[pns2019.1$C006==2 & pns2019.1$C008>=18 & pns2019.1$S068==1] <- 1\n",
    "pns2019.1$S001P<-factor(pns2019.1$S001P, levels=c(1,2), labels=c(\"Sim\",\"Não\"))\n",
    "summary(pns2019.1$S001P)\n",
    "\n",
    "\n",
    "# 2. Proporção de mulheres que realizaram pré-natal e que possuíam caderneta/cartão da gestante - S002P.\n",
    "pns2019.1$S002P <- NA\n",
    "pns2019.1$S002P[pns2019.1$C006==2 & pns2019.1$C008>=18 & pns2019.1$S068==1] <- 2\n",
    "pns2019.1$S002P[pns2019.1$C006==2 & pns2019.1$C008>=18 & pns2019.1$S068==1 & pns2019.1$S075==1] <- 1\n",
    "pns2019.1$S002P<-factor(pns2019.1$S002P, levels=c(1,2), labels=c(\"Sim\",\"Não\"))\n",
    "summary(pns2019.1$S002P)\n",
    "\n",
    "# 3. Proporção de mulheres que iniciaram pré-natal com menos de 13 semanas ou até 3 meses de gestação - S003P.\n",
    "pns2019.1$S003P <- NA\n",
    "pns2019.1$S003P[pns2019.1$C006==2 & pns2019.1$C008>=18 & pns2019.1$S068==1] <- 2\n",
    "pns2019.1$S003P[pns2019.1$C006==2 & pns2019.1$C008>=18 & pns2019.1$S068==1 & (pns2019.1$S06902 <= 3 | pns2019.1$S06901 < 13)] <- 1\n",
    "pns2019.1$S003P<-factor(pns2019.1$S003P, levels=c(1,2), labels=c(\"Sim\",\"Não\"))\n",
    "summary(pns2019.1$S003P)\n",
    "\n",
    "# 4. Proporção de mulheres que tiveram 6 ou mais consultas de pré-natal entre as gestantes com parto a termo ou pós-termo - S004P.\n",
    "pns2019.1$S004P <- NA\n",
    "pns2019.1$S004P[pns2019.1$C006==2 & pns2019.1$C008>=18 & pns2019.1$S068==1 & pns2019.1$S11801>=37] <- 2\n",
    "pns2019.1$S004P[pns2019.1$C006==2 & pns2019.1$C008>=18 & pns2019.1$S068==1 & pns2019.1$S11801>=37 & (pns2019.1$S070>=6 & pns2019.1$S070<8)] <- 1\n",
    "pns2019.1$S004P<-factor(pns2019.1$S004P, levels=c(1,2), labels=c(\"Sim\",\"Não\"))\n",
    "summary(pns2019.1$S004P)\n",
    "\n",
    "# 5. Proporção de mulheres que realizaram a maioria das consultas de pré-natal em estabelecimentos públicos de saúde - S005P.\n",
    "pns2019.1$S005P <- NA\n",
    "pns2019.1$S005P[pns2019.1$C006==2 & pns2019.1$C008>=18 & pns2019.1$S068==1] <- 2\n",
    "pns2019.1$S005P[pns2019.1$C006==2 & pns2019.1$C008>=18 & pns2019.1$S068==1  & (pns2019.1$S071 <=2 | pns2019.1$S071 <=4)] <- 1\n",
    "pns2019.1$S005P<-factor(pns2019.1$S005P, levels=c(1,2), labels=c(\"Sim\",\"Não\"))\n",
    "summary(pns2019.1$S005P)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Definições de abrangências"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Situação urbana ou rural"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>urbano</dt><dd>69873</dd><dt>rural</dt><dd>20973</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[urbano] 69873\n",
       "\\item[rural] 20973\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "urbano\n",
       ":   69873rural\n",
       ":   20973\n",
       "\n"
      ],
      "text/plain": [
       "urbano  rural \n",
       " 69873  20973 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Situação Urbano ou Rural\n",
    "pns2019.1 <- pns2019.1 %>% rename(Sit_Urbano_Rural=V0026)\n",
    "pns2019.1$Sit_Urbano_Rural<-factor(pns2019.1$Sit_Urbano_Rural, levels=c(1,2), labels=c(\"urbano\", \"rural\"))\n",
    "summary(pns2019.1$Sit_Urbano_Rural)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### UF"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Rondônia</dt><dd>2176</dd><dt>Acre</dt><dd>2380</dd><dt>Amazonas</dt><dd>3479</dd><dt>Roraima</dt><dd>2238</dd><dt>Pará</dt><dd>3853</dd><dt>Amapá</dt><dd>1554</dd><dt>Tocantins</dt><dd>1922</dd><dt>Maranhão</dt><dd>5080</dd><dt>Piauí</dt><dd>2740</dd><dt>Ceará</dt><dd>4265</dd><dt>Rio Grande do Norte</dt><dd>2962</dd><dt>Paraíba</dt><dd>3158</dd><dt>Pernambuco</dt><dd>4083</dd><dt>Alagoas</dt><dd>2987</dd><dt>Sergipe</dt><dd>2610</dd><dt>Bahia</dt><dd>3659</dd><dt>Minas Gerais</dt><dd>5209</dd><dt>Espírito Santo</dt><dd>3541</dd><dt>Rio de Janeiro</dt><dd>4966</dd><dt>São Paulo</dt><dd>6114</dd><dt>Paraná</dt><dd>3967</dd><dt>Santa Catarina</dt><dd>3738</dd><dt>Rio Grande do Sul</dt><dd>3767</dd><dt>Mato Grosso do Sul</dt><dd>2863</dd><dt>Mato Grosso</dt><dd>2468</dd><dt>Goiás</dt><dd>2702</dd><dt>Distrito Federal</dt><dd>2365</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Rondônia] 2176\n",
       "\\item[Acre] 2380\n",
       "\\item[Amazonas] 3479\n",
       "\\item[Roraima] 2238\n",
       "\\item[Pará] 3853\n",
       "\\item[Amapá] 1554\n",
       "\\item[Tocantins] 1922\n",
       "\\item[Maranhão] 5080\n",
       "\\item[Piauí] 2740\n",
       "\\item[Ceará] 4265\n",
       "\\item[Rio Grande do Norte] 2962\n",
       "\\item[Paraíba] 3158\n",
       "\\item[Pernambuco] 4083\n",
       "\\item[Alagoas] 2987\n",
       "\\item[Sergipe] 2610\n",
       "\\item[Bahia] 3659\n",
       "\\item[Minas Gerais] 5209\n",
       "\\item[Espírito Santo] 3541\n",
       "\\item[Rio de Janeiro] 4966\n",
       "\\item[São Paulo] 6114\n",
       "\\item[Paraná] 3967\n",
       "\\item[Santa Catarina] 3738\n",
       "\\item[Rio Grande do Sul] 3767\n",
       "\\item[Mato Grosso do Sul] 2863\n",
       "\\item[Mato Grosso] 2468\n",
       "\\item[Goiás] 2702\n",
       "\\item[Distrito Federal] 2365\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Rondônia\n",
       ":   2176Acre\n",
       ":   2380Amazonas\n",
       ":   3479Roraima\n",
       ":   2238Pará\n",
       ":   3853Amapá\n",
       ":   1554Tocantins\n",
       ":   1922Maranhão\n",
       ":   5080Piauí\n",
       ":   2740Ceará\n",
       ":   4265Rio Grande do Norte\n",
       ":   2962Paraíba\n",
       ":   3158Pernambuco\n",
       ":   4083Alagoas\n",
       ":   2987Sergipe\n",
       ":   2610Bahia\n",
       ":   3659Minas Gerais\n",
       ":   5209Espírito Santo\n",
       ":   3541Rio de Janeiro\n",
       ":   4966São Paulo\n",
       ":   6114Paraná\n",
       ":   3967Santa Catarina\n",
       ":   3738Rio Grande do Sul\n",
       ":   3767Mato Grosso do Sul\n",
       ":   2863Mato Grosso\n",
       ":   2468Goiás\n",
       ":   2702Distrito Federal\n",
       ":   2365\n",
       "\n"
      ],
      "text/plain": [
       "           Rondônia                Acre            Amazonas             Roraima \n",
       "               2176                2380                3479                2238 \n",
       "               Pará               Amapá           Tocantins            Maranhão \n",
       "               3853                1554                1922                5080 \n",
       "              Piauí               Ceará Rio Grande do Norte             Paraíba \n",
       "               2740                4265                2962                3158 \n",
       "         Pernambuco             Alagoas             Sergipe               Bahia \n",
       "               4083                2987                2610                3659 \n",
       "       Minas Gerais      Espírito Santo      Rio de Janeiro           São Paulo \n",
       "               5209                3541                4966                6114 \n",
       "             Paraná      Santa Catarina   Rio Grande do Sul  Mato Grosso do Sul \n",
       "               3967                3738                3767                2863 \n",
       "        Mato Grosso               Goiás    Distrito Federal \n",
       "               2468                2702                2365 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Estados - UFs\n",
    "pns2019.1 <- pns2019.1 %>% rename(Unidades_da_Federacao=V0001)\n",
    "pns2019.1$Unidades_da_Federacao<-factor(pns2019.1$Unidades_da_Federacao, 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),\n",
    "                                       label=c(\"Rondônia\",\"Acre\",\"Amazonas\",\"Roraima\",\"Pará\",\"Amapá\",\"Tocantins\",\"Maranhão\",\"Piauí\",\"Ceará\",\n",
    "                                        \"Rio Grande do Norte\",\"Paraíba\",\"Pernambuco\",\"Alagoas\",\"Sergipe\",\"Bahia\",\n",
    "                                        \"Minas Gerais\",\"Espírito Santo\",\"Rio de Janeiro\",\"São Paulo\",\n",
    "                                        \"Paraná\",\"Santa Catarina\",\"Rio Grande do Sul\", \n",
    "                                        \"Mato Grosso do Sul\",\"Mato Grosso\",\"Goiás\",\"Distrito Federal\"))\n",
    "summary(pns2019.1$Unidades_da_Federacao)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Grandes Regiões"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Norte</dt><dd>17602</dd><dt>Nordeste</dt><dd>31544</dd><dt>Sudeste</dt><dd>19830</dd><dt>Sul</dt><dd>11472</dd><dt>Centro-Oeste</dt><dd>10398</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Norte] 17602\n",
       "\\item[Nordeste] 31544\n",
       "\\item[Sudeste] 19830\n",
       "\\item[Sul] 11472\n",
       "\\item[Centro-Oeste] 10398\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Norte\n",
       ":   17602Nordeste\n",
       ":   31544Sudeste\n",
       ":   19830Sul\n",
       ":   11472Centro-Oeste\n",
       ":   10398\n",
       "\n"
      ],
      "text/plain": [
       "       Norte     Nordeste      Sudeste          Sul Centro-Oeste \n",
       "       17602        31544        19830        11472        10398 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Grandes Regiões\n",
    "pns2019.1 <- pns2019.1 %>% \n",
    "  mutate(GrandesRegioes = fct_collapse(Unidades_da_Federacao, \n",
    "                                          `Norte` = c(\"Rondônia\",\"Acre\",\"Amazonas\",\"Roraima\",\"Pará\", \"Amapá\",\"Tocantins\"),\n",
    "                                          `Nordeste` = c(\"Maranhão\", \"Piauí\", \"Ceará\", \"Rio Grande do Norte\", \"Paraíba\",\"Pernambuco\", \"Alagoas\",\"Sergipe\",\"Bahia\"),\n",
    "                                          `Sudeste` = c(\"Minas Gerais\", \"Espírito Santo\",\"Rio de Janeiro\", \"São Paulo\"), \n",
    "                                          `Sul` = c(\"Paraná\", \"Santa Catarina\", \"Rio Grande do Sul\"),\n",
    "                                          `Centro-Oeste`= c(\"Mato Grosso do Sul\",\"Mato Grosso\", \"Goiás\", \"Distrito Federal\")))\n",
    "summary(pns2019.1$GrandesRegioes)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Capital"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Porto Velho</dt><dd>2176</dd><dt>Rio Branco</dt><dd>2380</dd><dt>Manaus</dt><dd>3479</dd><dt>Boa Vista</dt><dd>2238</dd><dt>Belém</dt><dd>3853</dd><dt>Macapá</dt><dd>1554</dd><dt>Palmas</dt><dd>1922</dd><dt>São Luís</dt><dd>5080</dd><dt>Teresina</dt><dd>2740</dd><dt>Fortaleza</dt><dd>4265</dd><dt>Natal</dt><dd>2962</dd><dt>João Pessoa</dt><dd>3158</dd><dt>Recife</dt><dd>4083</dd><dt>Maceió</dt><dd>2987</dd><dt>Aracaju</dt><dd>2610</dd><dt>Salvador</dt><dd>3659</dd><dt>Belo Horizonte</dt><dd>5209</dd><dt>Vitória</dt><dd>3541</dd><dt>Rio de Janeiro</dt><dd>4966</dd><dt>São Paulo</dt><dd>6114</dd><dt>Curitiba</dt><dd>3967</dd><dt>Florianópolis</dt><dd>3738</dd><dt>Porto Alegre</dt><dd>3767</dd><dt>Campo Grande</dt><dd>2863</dd><dt>Cuiabá</dt><dd>2468</dd><dt>Goiânia</dt><dd>2702</dd><dt>Brasília</dt><dd>2365</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Porto Velho] 2176\n",
       "\\item[Rio Branco] 2380\n",
       "\\item[Manaus] 3479\n",
       "\\item[Boa Vista] 2238\n",
       "\\item[Belém] 3853\n",
       "\\item[Macapá] 1554\n",
       "\\item[Palmas] 1922\n",
       "\\item[São Luís] 5080\n",
       "\\item[Teresina] 2740\n",
       "\\item[Fortaleza] 4265\n",
       "\\item[Natal] 2962\n",
       "\\item[João Pessoa] 3158\n",
       "\\item[Recife] 4083\n",
       "\\item[Maceió] 2987\n",
       "\\item[Aracaju] 2610\n",
       "\\item[Salvador] 3659\n",
       "\\item[Belo Horizonte] 5209\n",
       "\\item[Vitória] 3541\n",
       "\\item[Rio de Janeiro] 4966\n",
       "\\item[São Paulo] 6114\n",
       "\\item[Curitiba] 3967\n",
       "\\item[Florianópolis] 3738\n",
       "\\item[Porto Alegre] 3767\n",
       "\\item[Campo Grande] 2863\n",
       "\\item[Cuiabá] 2468\n",
       "\\item[Goiânia] 2702\n",
       "\\item[Brasília] 2365\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Porto Velho\n",
       ":   2176Rio Branco\n",
       ":   2380Manaus\n",
       ":   3479Boa Vista\n",
       ":   2238Belém\n",
       ":   3853Macapá\n",
       ":   1554Palmas\n",
       ":   1922São Luís\n",
       ":   5080Teresina\n",
       ":   2740Fortaleza\n",
       ":   4265Natal\n",
       ":   2962João Pessoa\n",
       ":   3158Recife\n",
       ":   4083Maceió\n",
       ":   2987Aracaju\n",
       ":   2610Salvador\n",
       ":   3659Belo Horizonte\n",
       ":   5209Vitória\n",
       ":   3541Rio de Janeiro\n",
       ":   4966São Paulo\n",
       ":   6114Curitiba\n",
       ":   3967Florianópolis\n",
       ":   3738Porto Alegre\n",
       ":   3767Campo Grande\n",
       ":   2863Cuiabá\n",
       ":   2468Goiânia\n",
       ":   2702Brasília\n",
       ":   2365\n",
       "\n"
      ],
      "text/plain": [
       "   Porto Velho     Rio Branco         Manaus      Boa Vista          Belém \n",
       "          2176           2380           3479           2238           3853 \n",
       "        Macapá         Palmas       São Luís       Teresina      Fortaleza \n",
       "          1554           1922           5080           2740           4265 \n",
       "         Natal    João Pessoa         Recife         Maceió        Aracaju \n",
       "          2962           3158           4083           2987           2610 \n",
       "      Salvador Belo Horizonte        Vitória Rio de Janeiro      São Paulo \n",
       "          3659           5209           3541           4966           6114 \n",
       "      Curitiba  Florianópolis   Porto Alegre   Campo Grande         Cuiabá \n",
       "          3967           3738           3767           2863           2468 \n",
       "       Goiânia       Brasília \n",
       "          2702           2365 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Capital\n",
    "pns2019.1<- pns2019.1 %>% mutate(Capital= fct_collapse(Unidades_da_Federacao,\n",
    "                                        `Porto Velho`= \"Rondônia\", \n",
    "                                        `Boa Vista`= \"Roraima\",              \n",
    "                                        `Rio Branco`= \"Acre\", \n",
    "                                        `Manaus` = \"Amazonas\",\n",
    "                                        `Belém` = \"Pará\" ,\n",
    "                                        `Macapá`= \"Amapá\",\n",
    "                                        `Palmas` = \"Tocantins\",\n",
    "                                        `São Luís` = \"Maranhão\",\n",
    "                                        `Teresina`= \"Piauí\" ,\n",
    "                                        `Fortaleza`= \"Ceará\",\n",
    "                                        `Natal`= \"Rio Grande do Norte\",\n",
    "                                        `João Pessoa`= \"Paraíba\",\n",
    "                                        `Recife`= \"Pernambuco\",\n",
    "                                        `Maceió`= \"Alagoas\",\n",
    "                                        `Aracaju`= \"Sergipe\",\n",
    "                                        `Salvador`= \"Bahia\",\n",
    "                                        `Belo Horizonte`= \"Minas Gerais\",\n",
    "                                        `Vitória`= \"Espírito Santo\",\n",
    "                                        `Rio de Janeiro`= \"Rio de Janeiro\",\n",
    "                                        `São Paulo`= \"São Paulo\",\n",
    "                                        `Curitiba`= \"Paraná\",\n",
    "                                        `Florianópolis`= \"Santa Catarina\",\n",
    "                                        `Porto Alegre`= \"Rio Grande do Sul\",\n",
    "                                        `Campo Grande`=  \"Mato Grosso do Sul\",\n",
    "                                        `Cuiabá`= \"Mato Grosso\",\n",
    "                                        `Goiânia` = \"Goiás\",\n",
    "                                        `Brasília`= \"Distrito Federal\"))\n",
    "summary(pns2019.1$Capital)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Faixa Etária"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>18 a 24 anos</dt><dd>8145</dd><dt>25 a 29 anos</dt><dd>7249</dd><dt>30 a 39 anos</dt><dd>18150</dd><dt>40 anos ou mais</dt><dd>54987</dd><dt>NA's</dt><dd>2315</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[18 a 24 anos] 8145\n",
       "\\item[25 a 29 anos] 7249\n",
       "\\item[30 a 39 anos] 18150\n",
       "\\item[40 anos ou mais] 54987\n",
       "\\item[NA's] 2315\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "18 a 24 anos\n",
       ":   814525 a 29 anos\n",
       ":   724930 a 39 anos\n",
       ":   1815040 anos ou mais\n",
       ":   54987NA's\n",
       ":   2315\n",
       "\n"
      ],
      "text/plain": [
       "   18 a 24 anos    25 a 29 anos    30 a 39 anos 40 anos ou mais            NA's \n",
       "           8145            7249           18150           54987            2315 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Faixas Etárias\n",
    "pns2019.1 <- pns2019.1 %>% mutate(fx_idade_S=cut(C008,\n",
    "  breaks = c(18,25,30,40,120),\n",
    "  labels = c(\"18 a 24 anos\", \"25 a 29 anos\", \"30 a 39 anos\", \"40 anos ou mais\"), \n",
    "  ordered_result = TRUE, right = FALSE))\n",
    "summary(pns2019.1$fx_idade_S)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Raça"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Branca</dt><dd>33133</dd><dt>Preta</dt><dd>10345</dd><dt>Parda</dt><dd>45994</dd><dt>NA's</dt><dd>1374</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Branca] 33133\n",
       "\\item[Preta] 10345\n",
       "\\item[Parda] 45994\n",
       "\\item[NA's] 1374\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Branca\n",
       ":   33133Preta\n",
       ":   10345Parda\n",
       ":   45994NA's\n",
       ":   1374\n",
       "\n"
      ],
      "text/plain": [
       "Branca  Preta  Parda   NA's \n",
       " 33133  10345  45994   1374 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Raça\n",
    "pns2019.1 <- pns2019.1 %>% mutate(Raca= ifelse(C009==1, 1, \n",
    "                        ifelse(C009==2, 2, \n",
    "                            ifelse(C009==4, 3, 9))))\n",
    "\n",
    "pns2019.1$Raca<-factor(pns2019.1$Raca, levels=c(1,2,3),labels=c(\"Branca\", \"Preta\", \"Parda\"))\n",
    "\n",
    "summary(pns2019.1$Raca)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Renda per capita"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Até 1/2 SM</dt><dd>23697</dd><dt>1/2 até 1 SM</dt><dd>26406</dd><dt>1 até 2 SM</dt><dd>22466</dd><dt>2 até 3 SM</dt><dd>7612</dd><dt>Mais de 3 SM</dt><dd>10643</dd><dt>NA's</dt><dd>22</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Até 1/2 SM] 23697\n",
       "\\item[1/2 até 1 SM] 26406\n",
       "\\item[1 até 2 SM] 22466\n",
       "\\item[2 até 3 SM] 7612\n",
       "\\item[Mais de 3 SM] 10643\n",
       "\\item[NA's] 22\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Até 1/2 SM\n",
       ":   236971/2 até 1 SM\n",
       ":   264061 até 2 SM\n",
       ":   224662 até 3 SM\n",
       ":   7612Mais de 3 SM\n",
       ":   10643NA's\n",
       ":   22\n",
       "\n"
      ],
      "text/plain": [
       "  Até 1/2 SM 1/2 até 1 SM   1 até 2 SM   2 até 3 SM Mais de 3 SM         NA's \n",
       "       23697        26406        22466         7612        10643           22 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Rendimento domiciliar per capita\n",
    "pns2019.1 <- pns2019.1 %>% mutate(rend_per_capita = ifelse(VDF004 %in% 1:2, 1, \n",
    "                        ifelse(VDF004%in% 3, 2, \n",
    "                        ifelse(VDF004%in% 4, 3,\n",
    "                        ifelse(VDF004%in% 5, 4, \n",
    "                        ifelse(is.na(VDF004)==TRUE, NA_real_, 5))))))\n",
    "\n",
    "pns2019.1$rend_per_capita<-factor(pns2019.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\",\n",
    "                                                                                   \"2 até 3 SM\",\"Mais de 3 SM\"))\n",
    "summary(pns2019.1$rend_per_capita)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Escolaridade"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".dl-inline {width: auto; margin:0; padding: 0}\n",
       ".dl-inline>dt, .dl-inline>dd {float: none; width: auto; display: inline-block}\n",
       ".dl-inline>dt::after {content: \":\\0020\"; padding-right: .5ex}\n",
       ".dl-inline>dt:not(:first-of-type) {padding-left: .5ex}\n",
       "</style><dl class=dl-inline><dt>Fundamental incompleto ou equivalente</dt><dd>36276</dd><dt>Médio incompleto ou equivalente</dt><dd>13520</dd><dt>Superior incompleto ou equivalente</dt><dd>27433</dd><dt>Superior completo</dt><dd>13617</dd></dl>\n"
      ],
      "text/latex": [
       "\\begin{description*}\n",
       "\\item[Fundamental incompleto ou equivalente] 36276\n",
       "\\item[Médio incompleto ou equivalente] 13520\n",
       "\\item[Superior incompleto ou equivalente] 27433\n",
       "\\item[Superior completo] 13617\n",
       "\\end{description*}\n"
      ],
      "text/markdown": [
       "Fundamental incompleto ou equivalente\n",
       ":   36276Médio incompleto ou equivalente\n",
       ":   13520Superior incompleto ou equivalente\n",
       ":   27433Superior completo\n",
       ":   13617\n",
       "\n"
      ],
      "text/plain": [
       "Fundamental incompleto ou equivalente       Médio incompleto ou equivalente \n",
       "                                36276                                 13520 \n",
       "   Superior incompleto ou equivalente                     Superior completo \n",
       "                                27433                                 13617 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Escolaridade\n",
    "pns2019.1 <- pns2019.1 %>% mutate(gescol = ifelse(VDD004A %in% 1:2, 1, \n",
    "                        ifelse(VDD004A%in% 3:4, 2, \n",
    "                        ifelse(VDD004A%in% 5:6, 3,4\n",
    "                        ))))\n",
    "\n",
    "pns2019.1$gescol<-factor(pns2019.1$gescol, levels=c(1,2,3,4), \n",
    "                                  labels=c(\"Fundamental incompleto ou equivalente\",\"Médio incompleto ou equivalente\",\n",
    "                                           \"Superior incompleto ou equivalente\",\"Superior completo\"))\n",
    "summary(pns2019.1$gescol)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Criando indicadores"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Filtrando base de indicadores"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "     V0024            UPA_PNS      peso_morador_selec  S001P        S002P      \n",
       " 1210010: 1167   140001681:   18   Min.   : 0.00562   Sim : 2734   Sim : 2657  \n",
       " 1410011:  792   140003815:   18   1st Qu.: 0.26621   Não :   84   Não :   77  \n",
       " 2710111:  779   140005777:   18   Median : 0.54401   NA's:88028   NA's:88112  \n",
       " 2410011:  745   140006746:   18   Mean   : 1.00000                            \n",
       " 5010011:  738   140007081:   18   3rd Qu.: 1.12765                            \n",
       " 3210011:  711   140007715:   18   Max.   :61.09981                            \n",
       " (Other):85914   (Other)  :90738                                               \n",
       "  S003P        S004P        S005P            C008             C006      \n",
       " Sim : 2356   Sim : 2073   Sim : 2056   Min.   : 15.00   Min.   :1.000  \n",
       " Não :  378   Não :  661   Não :  678   1st Qu.: 32.00   1st Qu.:1.000  \n",
       " NA's:88112   NA's:88112   NA's:88112   Median : 45.00   Median :2.000  \n",
       "                                        Mean   : 46.39   Mean   :1.529  \n",
       "                                        3rd Qu.: 60.00   3rd Qu.:2.000  \n",
       "                                        Max.   :107.00   Max.   :2.000  \n",
       "                                                                        \n",
       "      C009            S068           S11801           S118      \n",
       " Min.   :1.000   Min.   :1.00    Min.   : 2.00   Min.   :1.0    \n",
       " 1st Qu.:1.000   1st Qu.:1.00    1st Qu.:38.00   1st Qu.:1.0    \n",
       " Median :4.000   Median :1.00    Median :39.00   Median :1.0    \n",
       " Mean   :2.679   Mean   :1.03    Mean   :38.16   Mean   :1.1    \n",
       " 3rd Qu.:4.000   3rd Qu.:1.00    3rd Qu.:40.00   3rd Qu.:1.0    \n",
       " Max.   :9.000   Max.   :2.00    Max.   :43.00   Max.   :2.0    \n",
       "                 NA's   :87936   NA's   :88227   NA's   :87936  \n",
       "     V0031       Sit_Urbano_Rural    Unidades_da_Federacao      GrandesRegioes \n",
       " Min.   :1.000   urbano:69873     São Paulo     : 6114     Norte       :17602  \n",
       " 1st Qu.:1.000   rural :20973     Minas Gerais  : 5209     Nordeste    :31544  \n",
       " Median :2.000                    Maranhão      : 5080     Sudeste     :19830  \n",
       " Mean   :2.605                    Rio de Janeiro: 4966     Sul         :11472  \n",
       " 3rd Qu.:4.000                    Ceará         : 4265     Centro-Oeste:10398  \n",
       " Max.   :4.000                    Pernambuco    : 4083                         \n",
       "                                  (Other)       :61129                         \n",
       "           Capital                fx_idade_S        Raca      \n",
       " São Paulo     : 6114   18 a 24 anos   : 8145   Branca:33133  \n",
       " Belo Horizonte: 5209   25 a 29 anos   : 7249   Preta :10345  \n",
       " São Luís      : 5080   30 a 39 anos   :18150   Parda :45994  \n",
       " Rio de Janeiro: 4966   40 anos ou mais:54987   NA's  : 1374  \n",
       " Fortaleza     : 4265   NA's           : 2315                 \n",
       " Recife        : 4083                                         \n",
       " (Other)       :61129                                         \n",
       "     rend_per_capita                                    gescol     \n",
       " Até 1/2 SM  :23697   Fundamental incompleto ou equivalente:36276  \n",
       " 1/2 até 1 SM:26406   Médio incompleto ou equivalente      :13520  \n",
       " 1 até 2 SM  :22466   Superior incompleto ou equivalente   :27433  \n",
       " 2 até 3 SM  : 7612   Superior completo                    :13617  \n",
       " Mais de 3 SM:10643                                                \n",
       " NA's        :   22                                                \n",
       "                                                                   "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#Selecionando variáveis para cálculo de indicadores no survey\n",
    "pns2019Ssurvey<- pns2019.1 %>% select(\"V0024\",\"UPA_PNS\",\"peso_morador_selec\", \"S001P\", \"S002P\", \"S003P\",\n",
    "       \"S004P\", \"S005P\", \"C008\", \"C006\", \"C009\", \"S068\", \"S11801\", \"S118\",\n",
    "       \"V0031\", \"Sit_Urbano_Rural\", \"Unidades_da_Federacao\",\n",
    "       \"GrandesRegioes\",  \"Capital\", \"fx_idade_S\", \"Raca\",\n",
    "       \"rend_per_capita\", \"gescol\") #\"Sexo\", \n",
    "summary(pns2019Ssurvey)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Exporta tabela filtrada com os dados específicos - Módulo S 2019 - Parte 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Salvando csv para cálculo de indicadores no survey\n",
    "diretorio_saida <- \"<coloque aqui o diretório>\"\n",
    "write.csv(pns2019Ssurvey, file.path(diretorio_saida, \"pns2019Ssurvey.csv\"))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Cria plano amostral complexo"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "desPNS=svydesign(id=~UPA_PNS, strat=~V0024, weight=~peso_morador_selec, nest=TRUE, \n",
    "                 data=pns2019Ssurvey)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "#survey design S001P\n",
    "desPNSS001P=subset(desPNS, C006==2 & C008>=18 & S068>0)\n",
    "desPNSS001P_C=subset(desPNS, C006==2 & C008>=18 & S068>0 & V0031==1)\n",
    "desPNSS001P_R=subset(desPNS, C006==2 & C008>=18 & S068>0 & !is.na(Raca))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [],
   "source": [
    "#survey design S002P, S003P, S005P\n",
    "desPNSS002P=subset(desPNS, C006==2 & C008>=18 & S068==1)\n",
    "desPNSS002P_C=subset(desPNS, C006==2 & C008>=18 & S068==1 & V0031==1)\n",
    "desPNSS002P_R=subset(desPNS, C006==2 & C008>=18 & S068==1 & !is.na(Raca)) "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "#survey design S004P\n",
    "desPNSS004P=subset(desPNS, C006==2 & C008>=18 & (S068==1 & S11801>=37))\n",
    "desPNSS004P_C=subset(desPNS, C006==2 & C008>=18 & (S068==1 & S11801>=37) & V0031==1) \n",
    "desPNSS004P_R=subset(desPNS, C006==2 & C008>=18 & (S068==1 & S11801>=37) & !is.na(Raca))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Criação da tabela de indicadores\n",
    "Essa tabela é responsável por unir os indicadores no formato do painel de indicadores"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [],
   "source": [
    "matrizIndicadores = data.frame()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Definição de variáveis para iteração dos indicadores"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [],
   "source": [
    "ListaIndicadores = c(~S001P, ~S002P, ~S003P, ~S004P, ~S005P)\n",
    "ListaIndicadoresTexto = c(\"S001P\", \"S002P\", \"S003P\", \"S004P\", \"S005P\" )\n",
    "ListaTotais = c('Brasil','Capital')\n",
    "Ano <- \"2019\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "ListaDominiosS001 = c(~Raca,~rend_per_capita,~fx_idade_S,~Sit_Urbano_Rural,\n",
    "                      ~Unidades_da_Federacao,~GrandesRegioes,~Capital,\n",
    "                      ~gescol)\n",
    "ListaDominiosTextoS001= c(\"raça\",\"rend_per_capita\",\"fx_idade_S\",\"urb_rur\",\n",
    "                          \"uf\",\"região\",\"capital\",\"gescol\") "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Preenchendo a tabela de indicadores\n",
    "Essas iterações rodam por indicador, abrangência e por design"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "#Cálculo dos indicadores usando o pacote survey \n",
    "i <- 0\n",
    "#Para cada indicador\n",
    "for( indicador in ListaIndicadores){\n",
    "    i <- i + 1; j <- 1\n",
    "    if (ListaIndicadoresTexto[i]== \"S001P\" | \n",
    "        ListaIndicadoresTexto[i]== \"S002P\" | \n",
    "        ListaIndicadoresTexto[i]== \"S003P\" | \n",
    "        ListaIndicadoresTexto[i]== \"S004P\" | \n",
    "        ListaIndicadoresTexto[i]== \"S005P\"){\n",
    "        ListaDominios<-ListaDominiosS001\n",
    "        ListaDominiosTexto<-ListaDominiosTextoS001\n",
    "        } else {\n",
    "                    ListaDominios<-ListaDominiosS001\n",
    "ListaDominiosTexto<-ListaDominiosTextoS001\n",
    "                }\n",
    "    #Para cada dominio\n",
    "    for (dominio in ListaDominios){\n",
    "#design especifico para capital que é subconjunto do dataframe total\n",
    "         if (ListaDominiosTexto[j]==\"capital\"){\n",
    "#designs especificos por variavel que são subconjuntos do dataset total\n",
    "            if (ListaIndicadoresTexto[i]== \"S002P\" |\n",
    "                ListaIndicadoresTexto[i]== \"S003P\" |\n",
    "                ListaIndicadoresTexto[i]== \"S005P\"){\n",
    "              dataframe_indicador<-svyby( indicador , dominio , desPNSS002P_C , svymean,vartype= c(\"ci\",\"cv\")) \n",
    "              } else if (ListaIndicadoresTexto[i]== \"S001P\") {\n",
    "                  dataframe_indicador<-svyby( indicador , dominio , desPNSS001P_C , svymean,vartype= c(\"ci\",\"cv\"))\n",
    "                 } else if (ListaIndicadoresTexto[i]== \"S004P\") {\n",
    "                    dataframe_indicador<-svyby( indicador , dominio , desPNSS004P_C, svymean,vartype= c(\"ci\",\"cv\"))\n",
    "                        }\n",
    "               } else if  (ListaDominiosTexto[j]==\"raça\"){\n",
    "                    if (ListaIndicadoresTexto[i]== \"S002P\" |\n",
    "                        ListaIndicadoresTexto[i]== \"S003P\" |\n",
    "                        ListaIndicadoresTexto[i]== \"S005P\"){\n",
    "                            dataframe_indicador<-svyby( indicador , dominio , desPNSS002P_R , svymean,vartype= c(\"ci\",\"cv\"))\n",
    "                        } else if\n",
    "(ListaIndicadoresTexto[i]== \"S001P\"){\n",
    "                            dataframe_indicador<-svyby( indicador , dominio , desPNSS001P_R , svymean,vartype= c(\"ci\",\"cv\"))\n",
    "                        } else if\n",
    "(ListaIndicadoresTexto[i]== \"S004P\"){\n",
    "                            dataframe_indicador<-svyby( indicador , dominio , desPNSS004P_R , svymean,vartype= c(\"ci\",\"cv\"))\n",
    "                        }\n",
    "               #design geral para o subconjunto maior que 15 anos\n",
    "                 } else { \n",
    "                        if (ListaIndicadoresTexto[i]== \"S002P\" |\n",
    "                            ListaIndicadoresTexto[i]== \"S003P\" |\n",
    "                            ListaIndicadoresTexto[i]== \"S005P\"){\n",
    "                                dataframe_indicador<-svyby( indicador , dominio , desPNSS002P , svymean,vartype= c(\"ci\",\"cv\"))\n",
    "                            } else if\n",
    "(ListaIndicadoresTexto[i]== \"S001P\"){\n",
    "                            dataframe_indicador<-svyby( indicador , dominio , desPNSS001P , svymean,vartype= c(\"ci\",\"cv\"))\n",
    "                        } else if\n",
    "(ListaIndicadoresTexto[i]== \"S004P\")\n",
    "                     {\n",
    "                                dataframe_indicador<-svyby( indicador , dominio , desPNSS004P , svymean,vartype= c(\"ci\",\"cv\"))\n",
    "                            }\n",
    "                        }\n",
    "#União do dataframe de indicadores no formato do painel disponibilizado para PNS\n",
    "dataframe_indicador<-data.frame(dataframe_indicador)\n",
    "               colnames(dataframe_indicador) <- c(\"abr_nome\",\"Sim\",\"Nao\",\"LowerS\",\"LowerN\",\"UpperS\",\"UpperN\",\"cvS\",\"cvN\")\n",
    "               dataframe_indicador$Indicador <- ListaIndicadoresTexto[i]\n",
    "               dataframe_indicador$abr_tipo <- ListaDominiosTexto[j]\n",
    "               dataframe_indicador$Ano <- Ano\n",
    "               dataframe_indicador <- dataframe_indicador %>% select(\"abr_tipo\",\"abr_nome\",\"Ano\",\"Indicador\",\"Sim\",\"LowerS\",\"UpperS\",\"cvS\")\n",
    "               matrizIndicadores <-rbind(matrizIndicadores,dataframe_indicador)\n",
    "               j <- j + 1\n",
    "    }\n",
    "}"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Criando a tabela pela abrangência total"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "matriz_totais <- data.frame()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Preenchendo a tabela com as abrangencia Brasil e total das capitais"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [],
   "source": [
    "i <- 0\n",
    "#para cada indicador\n",
    "for(indicador in ListaIndicadores){\n",
    "    i <- i + 1\n",
    "    #para os totais Brasil e total das capitais\n",
    "    for(total in ListaTotais){\n",
    "        dataframe_indicador <- data.frame()\n",
    "        dataframe_indicador_S <- data.frame()\n",
    "        #Uso do design que é subconjunto do dataset para cada Capital\n",
    "        if (total == \"Capital\"){\n",
    "                #Indicadores que são subconjunto do dataset total\n",
    "                if (ListaIndicadoresTexto[i]== \"S002P\" |\n",
    "                    ListaIndicadoresTexto[i]== \"S003P\" |\n",
    "                    ListaIndicadoresTexto[i]== \"S005P\"){\n",
    "                            dataframe_indicador <- svymean(indicador,desPNSS002P_C)\n",
    "                        } else if\n",
    "(ListaIndicadoresTexto[i]== \"S001P\"){\n",
    "                            dataframe_indicador <- svymean(indicador,desPNSS001P_C)\n",
    "                        } else if\n",
    "(ListaIndicadoresTexto[i]== \"S004P\"){\n",
    "                            dataframe_indicador <- svymean(indicador,desPNSS004P_C)\n",
    "                    }\n",
    "                } else {\n",
    "                if (ListaIndicadoresTexto[i]== \"S002P\" |\n",
    "                    ListaIndicadoresTexto[i]== \"S003P\" |\n",
    "                    ListaIndicadoresTexto[i]== \"S005P\"){\n",
    "                            dataframe_indicador <- svymean(indicador,desPNSS002P)\n",
    "                        } else if (ListaIndicadoresTexto[i]== \"S001P\"){\n",
    "                            dataframe_indicador <- svymean(indicador,desPNSS001P)\n",
    "                        } else if\n",
    "(ListaIndicadoresTexto[i]== \"S004P\"){\n",
    "                            dataframe_indicador <- svymean(indicador,desPNSS004P)\n",
    "                    }\n",
    "                }\n",
    "        intervalo_confianca <- confint(dataframe_indicador)\n",
    "        coeficiente_variacao <- cv(dataframe_indicador)\n",
    "        dataframe_indicador <- cbind(data.frame(dataframe_indicador),data.frame(intervalo_confianca))\n",
    "        dataframe_indicador <- cbind(data.frame(dataframe_indicador),data.frame(coeficiente_variacao))\n",
    "        \n",
    "        dataframe_indicador <- dataframe_indicador %>% \n",
    "                               select('mean','X2.5..','X97.5..',coeficiente_variacao) \n",
    "        dataframe_indicador_S <- dataframe_indicador %>% \n",
    "                                 slice(1)\n",
    "                \n",
    "        colnames(dataframe_indicador_S) <- c('Sim','LowerS','UpperS', 'cvS')\n",
    "        dataframe_indicador_S$Indicador <- ListaIndicadoresTexto[i]\n",
    "        \n",
    "        dataframe_indicador_S$abr_tipo <- \"total\"\n",
    "        dataframe_indicador_S$abr_nome <- total\n",
    "        dataframe_indicador_S$Ano <- Ano \n",
    "        dataframe_indicador_S <- dataframe_indicador_S %>% \n",
    "                             select(\"abr_tipo\",\"abr_nome\",\"Ano\",\"Indicador\",\"Sim\",\"LowerS\",\"UpperS\",'cvS')\n",
    "        \n",
    "        matriz_totais <-rbind(matriz_totais,dataframe_indicador_S)\n",
    "    }\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 10 × 8</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>abr_tipo</th><th scope=col>abr_nome</th><th scope=col>Ano</th><th scope=col>Indicador</th><th scope=col>Sim</th><th scope=col>LowerS</th><th scope=col>UpperS</th><th scope=col>cvS</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>S001PSim</th><td>total</td><td>Brasil </td><td>2019</td><td>S001P</td><td>0.9802181</td><td>0.9735800</td><td>0.9868562</td><td>0.003455195</td></tr>\n",
       "\t<tr><th scope=row>S001PSim1</th><td>total</td><td>Capital</td><td>2019</td><td>S001P</td><td>0.9806496</td><td>0.9721002</td><td>0.9891990</td><td>0.004448070</td></tr>\n",
       "\t<tr><th scope=row>S002PSim</th><td>total</td><td>Brasil </td><td>2019</td><td>S002P</td><td>0.9624012</td><td>0.9492620</td><td>0.9755404</td><td>0.006965703</td></tr>\n",
       "\t<tr><th scope=row>S002PSim1</th><td>total</td><td>Capital</td><td>2019</td><td>S002P</td><td>0.9341446</td><td>0.9050284</td><td>0.9632609</td><td>0.015902776</td></tr>\n",
       "\t<tr><th scope=row>S003PSim</th><td>total</td><td>Brasil </td><td>2019</td><td>S003P</td><td>0.8814929</td><td>0.8630816</td><td>0.8999041</td><td>0.010656560</td></tr>\n",
       "\t<tr><th scope=row>S003PSim1</th><td>total</td><td>Capital</td><td>2019</td><td>S003P</td><td>0.8898081</td><td>0.8547462</td><td>0.9248701</td><td>0.020104431</td></tr>\n",
       "\t<tr><th scope=row>S004PSim</th><td>total</td><td>Brasil </td><td>2019</td><td>S004P</td><td>0.8939233</td><td>0.8748669</td><td>0.9129797</td><td>0.010876579</td></tr>\n",
       "\t<tr><th scope=row>S004PSim1</th><td>total</td><td>Capital</td><td>2019</td><td>S004P</td><td>0.9086737</td><td>0.8785255</td><td>0.9388218</td><td>0.016927946</td></tr>\n",
       "\t<tr><th scope=row>S005PSim</th><td>total</td><td>Brasil </td><td>2019</td><td>S005P</td><td>0.6712686</td><td>0.6383040</td><td>0.7042333</td><td>0.025055548</td></tr>\n",
       "\t<tr><th scope=row>S005PSim1</th><td>total</td><td>Capital</td><td>2019</td><td>S005P</td><td>0.5467601</td><td>0.4924781</td><td>0.6010420</td><td>0.050653646</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 10 × 8\n",
       "\\begin{tabular}{r|llllllll}\n",
       "  & abr\\_tipo & abr\\_nome & Ano & Indicador & Sim & LowerS & UpperS & cvS\\\\\n",
       "  & <chr> & <chr> & <chr> & <chr> & <dbl> & <dbl> & <dbl> & <dbl>\\\\\n",
       "\\hline\n",
       "\tS001PSim & total & Brasil  & 2019 & S001P & 0.9802181 & 0.9735800 & 0.9868562 & 0.003455195\\\\\n",
       "\tS001PSim1 & total & Capital & 2019 & S001P & 0.9806496 & 0.9721002 & 0.9891990 & 0.004448070\\\\\n",
       "\tS002PSim & total & Brasil  & 2019 & S002P & 0.9624012 & 0.9492620 & 0.9755404 & 0.006965703\\\\\n",
       "\tS002PSim1 & total & Capital & 2019 & S002P & 0.9341446 & 0.9050284 & 0.9632609 & 0.015902776\\\\\n",
       "\tS003PSim & total & Brasil  & 2019 & S003P & 0.8814929 & 0.8630816 & 0.8999041 & 0.010656560\\\\\n",
       "\tS003PSim1 & total & Capital & 2019 & S003P & 0.8898081 & 0.8547462 & 0.9248701 & 0.020104431\\\\\n",
       "\tS004PSim & total & Brasil  & 2019 & S004P & 0.8939233 & 0.8748669 & 0.9129797 & 0.010876579\\\\\n",
       "\tS004PSim1 & total & Capital & 2019 & S004P & 0.9086737 & 0.8785255 & 0.9388218 & 0.016927946\\\\\n",
       "\tS005PSim & total & Brasil  & 2019 & S005P & 0.6712686 & 0.6383040 & 0.7042333 & 0.025055548\\\\\n",
       "\tS005PSim1 & total & Capital & 2019 & S005P & 0.5467601 & 0.4924781 & 0.6010420 & 0.050653646\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 10 × 8\n",
       "\n",
       "| <!--/--> | abr_tipo &lt;chr&gt; | abr_nome &lt;chr&gt; | Ano &lt;chr&gt; | Indicador &lt;chr&gt; | Sim &lt;dbl&gt; | LowerS &lt;dbl&gt; | UpperS &lt;dbl&gt; | cvS &lt;dbl&gt; |\n",
       "|---|---|---|---|---|---|---|---|---|\n",
       "| S001PSim | total | Brasil  | 2019 | S001P | 0.9802181 | 0.9735800 | 0.9868562 | 0.003455195 |\n",
       "| S001PSim1 | total | Capital | 2019 | S001P | 0.9806496 | 0.9721002 | 0.9891990 | 0.004448070 |\n",
       "| S002PSim | total | Brasil  | 2019 | S002P | 0.9624012 | 0.9492620 | 0.9755404 | 0.006965703 |\n",
       "| S002PSim1 | total | Capital | 2019 | S002P | 0.9341446 | 0.9050284 | 0.9632609 | 0.015902776 |\n",
       "| S003PSim | total | Brasil  | 2019 | S003P | 0.8814929 | 0.8630816 | 0.8999041 | 0.010656560 |\n",
       "| S003PSim1 | total | Capital | 2019 | S003P | 0.8898081 | 0.8547462 | 0.9248701 | 0.020104431 |\n",
       "| S004PSim | total | Brasil  | 2019 | S004P | 0.8939233 | 0.8748669 | 0.9129797 | 0.010876579 |\n",
       "| S004PSim1 | total | Capital | 2019 | S004P | 0.9086737 | 0.8785255 | 0.9388218 | 0.016927946 |\n",
       "| S005PSim | total | Brasil  | 2019 | S005P | 0.6712686 | 0.6383040 | 0.7042333 | 0.025055548 |\n",
       "| S005PSim1 | total | Capital | 2019 | S005P | 0.5467601 | 0.4924781 | 0.6010420 | 0.050653646 |\n",
       "\n"
      ],
      "text/plain": [
       "          abr_tipo abr_nome Ano  Indicador Sim       LowerS    UpperS   \n",
       "S001PSim  total    Brasil   2019 S001P     0.9802181 0.9735800 0.9868562\n",
       "S001PSim1 total    Capital  2019 S001P     0.9806496 0.9721002 0.9891990\n",
       "S002PSim  total    Brasil   2019 S002P     0.9624012 0.9492620 0.9755404\n",
       "S002PSim1 total    Capital  2019 S002P     0.9341446 0.9050284 0.9632609\n",
       "S003PSim  total    Brasil   2019 S003P     0.8814929 0.8630816 0.8999041\n",
       "S003PSim1 total    Capital  2019 S003P     0.8898081 0.8547462 0.9248701\n",
       "S004PSim  total    Brasil   2019 S004P     0.8939233 0.8748669 0.9129797\n",
       "S004PSim1 total    Capital  2019 S004P     0.9086737 0.8785255 0.9388218\n",
       "S005PSim  total    Brasil   2019 S005P     0.6712686 0.6383040 0.7042333\n",
       "S005PSim1 total    Capital  2019 S005P     0.5467601 0.4924781 0.6010420\n",
       "          cvS        \n",
       "S001PSim  0.003455195\n",
       "S001PSim1 0.004448070\n",
       "S002PSim  0.006965703\n",
       "S002PSim1 0.015902776\n",
       "S003PSim  0.010656560\n",
       "S003PSim1 0.020104431\n",
       "S004PSim  0.010876579\n",
       "S004PSim1 0.016927946\n",
       "S005PSim  0.025055548\n",
       "S005PSim1 0.050653646"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "matriz_totais"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Unindo tabela de indicadores e de totais"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [],
   "source": [
    "matrizIndicadores<-rbind(matrizIndicadores,matriz_totais)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Visualizando tabela de indicadores"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 395 × 8</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>abr_tipo</th><th scope=col>abr_nome</th><th scope=col>Ano</th><th scope=col>Indicador</th><th scope=col>Sim</th><th scope=col>LowerS</th><th scope=col>UpperS</th><th scope=col>cvS</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;fct&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>Branca</th><td>raça           </td><td>Branca             </td><td>2019</td><td>S001P</td><td>0.9848342</td><td>0.9728009</td><td>0.9968675</td><td>6.234103e-03</td></tr>\n",
       "\t<tr><th scope=row>Preta</th><td>raça           </td><td>Preta              </td><td>2019</td><td>S001P</td><td>0.9706214</td><td>0.9496967</td><td>0.9915462</td><td>1.099923e-02</td></tr>\n",
       "\t<tr><th scope=row>Parda</th><td>raça           </td><td>Parda              </td><td>2019</td><td>S001P</td><td>0.9802893</td><td>0.9718436</td><td>0.9887350</td><td>4.395738e-03</td></tr>\n",
       "\t<tr><th scope=row>Até 1/2 SM</th><td>rend_per_capita</td><td>Até 1/2 SM         </td><td>2019</td><td>S001P</td><td>0.9748696</td><td>0.9634649</td><td>0.9862743</td><td>5.968842e-03</td></tr>\n",
       "\t<tr><th scope=row>1/2 até 1 SM</th><td>rend_per_capita</td><td>1/2 até 1 SM       </td><td>2019</td><td>S001P</td><td>0.9809000</td><td>0.9691323</td><td>0.9926676</td><td>6.120930e-03</td></tr>\n",
       "\t<tr><th scope=row>1 até 2 SM</th><td>rend_per_capita</td><td>1 até 2 SM         </td><td>2019</td><td>S001P</td><td>0.9900064</td><td>0.9798280</td><td>1.0001847</td><td>5.245551e-03</td></tr>\n",
       "\t<tr><th scope=row>2 até 3 SM</th><td>rend_per_capita</td><td>2 até 3 SM         </td><td>2019</td><td>S001P</td><td>0.9874788</td><td>0.9684689</td><td>1.0064888</td><td>9.822109e-03</td></tr>\n",
       "\t<tr><th scope=row>Mais de 3 SM</th><td>rend_per_capita</td><td>Mais de 3 SM       </td><td>2019</td><td>S001P</td><td>0.9850570</td><td>0.9559967</td><td>1.0141173</td><td>1.505188e-02</td></tr>\n",
       "\t<tr><th scope=row>18 a 24 anos</th><td>fx_idade_S     </td><td>18 a 24 anos       </td><td>2019</td><td>S001P</td><td>0.9780003</td><td>0.9636430</td><td>0.9923577</td><td>7.490081e-03</td></tr>\n",
       "\t<tr><th scope=row>25 a 29 anos</th><td>fx_idade_S     </td><td>25 a 29 anos       </td><td>2019</td><td>S001P</td><td>0.9761319</td><td>0.9599660</td><td>0.9922979</td><td>8.449762e-03</td></tr>\n",
       "\t<tr><th scope=row>30 a 39 anos</th><td>fx_idade_S     </td><td>30 a 39 anos       </td><td>2019</td><td>S001P</td><td>0.9872648</td><td>0.9801366</td><td>0.9943931</td><td>3.683855e-03</td></tr>\n",
       "\t<tr><th scope=row>40 anos ou mais</th><td>fx_idade_S     </td><td>40 anos ou mais    </td><td>2019</td><td>S001P</td><td>0.9584802</td><td>0.9286191</td><td>0.9883414</td><td>1.589556e-02</td></tr>\n",
       "\t<tr><th scope=row>urbano</th><td>urb_rur        </td><td>urbano             </td><td>2019</td><td>S001P</td><td>0.9822630</td><td>0.9752768</td><td>0.9892492</td><td>3.628817e-03</td></tr>\n",
       "\t<tr><th scope=row>rural</th><td>urb_rur        </td><td>rural              </td><td>2019</td><td>S001P</td><td>0.9690309</td><td>0.9495219</td><td>0.9885399</td><td>1.027184e-02</td></tr>\n",
       "\t<tr><th scope=row>Rondônia</th><td>uf             </td><td>Rondônia           </td><td>2019</td><td>S001P</td><td>0.9922432</td><td>0.9813513</td><td>1.0031351</td><td>5.600634e-03</td></tr>\n",
       "\t<tr><th scope=row>Acre</th><td>uf             </td><td>Acre               </td><td>2019</td><td>S001P</td><td>0.9400931</td><td>0.8740058</td><td>1.0061804</td><td>3.586734e-02</td></tr>\n",
       "\t<tr><th scope=row>Amazonas</th><td>uf             </td><td>Amazonas           </td><td>2019</td><td>S001P</td><td>0.9471638</td><td>0.9097590</td><td>0.9845686</td><td>2.014902e-02</td></tr>\n",
       "\t<tr><th scope=row>Roraima</th><td>uf             </td><td>Roraima            </td><td>2019</td><td>S001P</td><td>0.9381625</td><td>0.8778433</td><td>0.9984818</td><td>3.280420e-02</td></tr>\n",
       "\t<tr><th scope=row>Pará</th><td>uf             </td><td>Pará               </td><td>2019</td><td>S001P</td><td>0.9753771</td><td>0.9511006</td><td>0.9996535</td><td>1.269886e-02</td></tr>\n",
       "\t<tr><th scope=row>Amapá</th><td>uf             </td><td>Amapá              </td><td>2019</td><td>S001P</td><td>0.9382793</td><td>0.8674707</td><td>1.0090880</td><td>3.850401e-02</td></tr>\n",
       "\t<tr><th scope=row>Tocantins</th><td>uf             </td><td>Tocantins          </td><td>2019</td><td>S001P</td><td>1.0000000</td><td>1.0000000</td><td>1.0000000</td><td>0.000000e+00</td></tr>\n",
       "\t<tr><th scope=row>Maranhão</th><td>uf             </td><td>Maranhão           </td><td>2019</td><td>S001P</td><td>0.9882002</td><td>0.9774606</td><td>0.9989399</td><td>5.544946e-03</td></tr>\n",
       "\t<tr><th scope=row>Piauí</th><td>uf             </td><td>Piauí              </td><td>2019</td><td>S001P</td><td>1.0000000</td><td>1.0000000</td><td>1.0000000</td><td>0.000000e+00</td></tr>\n",
       "\t<tr><th scope=row>Ceará</th><td>uf             </td><td>Ceará              </td><td>2019</td><td>S001P</td><td>0.9908578</td><td>0.9730003</td><td>1.0087154</td><td>9.195250e-03</td></tr>\n",
       "\t<tr><th scope=row>Rio Grande do Norte</th><td>uf             </td><td>Rio Grande do Norte</td><td>2019</td><td>S001P</td><td>0.9599420</td><td>0.9283055</td><td>0.9915785</td><td>1.681495e-02</td></tr>\n",
       "\t<tr><th scope=row>Paraíba</th><td>uf             </td><td>Paraíba            </td><td>2019</td><td>S001P</td><td>0.9259756</td><td>0.8549651</td><td>0.9969862</td><td>3.912687e-02</td></tr>\n",
       "\t<tr><th scope=row>Pernambuco</th><td>uf             </td><td>Pernambuco         </td><td>2019</td><td>S001P</td><td>0.9653698</td><td>0.9289815</td><td>1.0017582</td><td>1.923183e-02</td></tr>\n",
       "\t<tr><th scope=row>Alagoas</th><td>uf             </td><td>Alagoas            </td><td>2019</td><td>S001P</td><td>0.9641940</td><td>0.9305357</td><td>0.9978523</td><td>1.781065e-02</td></tr>\n",
       "\t<tr><th scope=row>Sergipe</th><td>uf             </td><td>Sergipe            </td><td>2019</td><td>S001P</td><td>1.0000000</td><td>1.0000000</td><td>1.0000000</td><td>1.782482e-17</td></tr>\n",
       "\t<tr><th scope=row>Bahia</th><td>uf             </td><td>Bahia              </td><td>2019</td><td>S001P</td><td>0.9956594</td><td>0.9871376</td><td>1.0041812</td><td>4.366894e-03</td></tr>\n",
       "\t<tr><th scope=row>⋮</th><td>⋮</td><td>⋮</td><td>⋮</td><td>⋮</td><td>⋮</td><td>⋮</td><td>⋮</td><td>⋮</td></tr>\n",
       "\t<tr><th scope=row>João Pessoa4</th><td>capital</td><td>João Pessoa                          </td><td>2019</td><td>S005P</td><td>0.4891024</td><td>0.31333773</td><td>0.6648670</td><td>0.183351143</td></tr>\n",
       "\t<tr><th scope=row>Recife4</th><td>capital</td><td>Recife                               </td><td>2019</td><td>S005P</td><td>0.6154032</td><td>0.35946656</td><td>0.8713398</td><td>0.212189839</td></tr>\n",
       "\t<tr><th scope=row>Maceió4</th><td>capital</td><td>Maceió                               </td><td>2019</td><td>S005P</td><td>0.6770275</td><td>0.53995647</td><td>0.8140986</td><td>0.103297871</td></tr>\n",
       "\t<tr><th scope=row>Aracaju4</th><td>capital</td><td>Aracaju                              </td><td>2019</td><td>S005P</td><td>0.5389670</td><td>0.28151970</td><td>0.7964143</td><td>0.243712643</td></tr>\n",
       "\t<tr><th scope=row>Salvador4</th><td>capital</td><td>Salvador                             </td><td>2019</td><td>S005P</td><td>0.5306442</td><td>0.26607586</td><td>0.7952126</td><td>0.254381995</td></tr>\n",
       "\t<tr><th scope=row>Belo Horizonte4</th><td>capital</td><td>Belo Horizonte                       </td><td>2019</td><td>S005P</td><td>0.3741457</td><td>0.20269372</td><td>0.5455976</td><td>0.233804895</td></tr>\n",
       "\t<tr><th scope=row>Vitória4</th><td>capital</td><td>Vitória                              </td><td>2019</td><td>S005P</td><td>0.4852118</td><td>0.20101229</td><td>0.7694113</td><td>0.298843552</td></tr>\n",
       "\t<tr><th scope=row>Rio de Janeiro9</th><td>capital</td><td>Rio de Janeiro                       </td><td>2019</td><td>S005P</td><td>0.5061534</td><td>0.27485324</td><td>0.7374535</td><td>0.233155491</td></tr>\n",
       "\t<tr><th scope=row>São Paulo8</th><td>capital</td><td>São Paulo                            </td><td>2019</td><td>S005P</td><td>0.3925847</td><td>0.21472169</td><td>0.5704477</td><td>0.231155467</td></tr>\n",
       "\t<tr><th scope=row>Curitiba4</th><td>capital</td><td>Curitiba                             </td><td>2019</td><td>S005P</td><td>0.2620324</td><td>0.03101129</td><td>0.4930534</td><td>0.449830155</td></tr>\n",
       "\t<tr><th scope=row>Florianópolis4</th><td>capital</td><td>Florianópolis                        </td><td>2019</td><td>S005P</td><td>0.5773166</td><td>0.31513885</td><td>0.8394944</td><td>0.231704095</td></tr>\n",
       "\t<tr><th scope=row>Porto Alegre4</th><td>capital</td><td>Porto Alegre                         </td><td>2019</td><td>S005P</td><td>0.4038861</td><td>0.17145033</td><td>0.6363218</td><td>0.293626959</td></tr>\n",
       "\t<tr><th scope=row>Campo Grande4</th><td>capital</td><td>Campo Grande                         </td><td>2019</td><td>S005P</td><td>0.6491680</td><td>0.50650967</td><td>0.7918262</td><td>0.112122258</td></tr>\n",
       "\t<tr><th scope=row>Cuiabá4</th><td>capital</td><td>Cuiabá                               </td><td>2019</td><td>S005P</td><td>0.6268727</td><td>0.37278925</td><td>0.8809562</td><td>0.206799253</td></tr>\n",
       "\t<tr><th scope=row>Goiânia4</th><td>capital</td><td>Goiânia                              </td><td>2019</td><td>S005P</td><td>0.5727847</td><td>0.32629493</td><td>0.8192745</td><td>0.219563127</td></tr>\n",
       "\t<tr><th scope=row>Brasília4</th><td>capital</td><td>Brasília                             </td><td>2019</td><td>S005P</td><td>0.6320122</td><td>0.47297955</td><td>0.7910448</td><td>0.128384529</td></tr>\n",
       "\t<tr><th scope=row>Fundamental incompleto ou equivalente4</th><td>gescol </td><td>Fundamental incompleto ou equivalente</td><td>2019</td><td>S005P</td><td>0.9124938</td><td>0.87271915</td><td>0.9522685</td><td>0.022239685</td></tr>\n",
       "\t<tr><th scope=row>Médio incompleto ou equivalente4</th><td>gescol </td><td>Médio incompleto ou equivalente      </td><td>2019</td><td>S005P</td><td>0.9012396</td><td>0.85482924</td><td>0.9476500</td><td>0.026274048</td></tr>\n",
       "\t<tr><th scope=row>Superior incompleto ou equivalente4</th><td>gescol </td><td>Superior incompleto ou equivalente   </td><td>2019</td><td>S005P</td><td>0.6802039</td><td>0.63250129</td><td>0.7279064</td><td>0.035781180</td></tr>\n",
       "\t<tr><th scope=row>Superior completo4</th><td>gescol </td><td>Superior completo                    </td><td>2019</td><td>S005P</td><td>0.2032387</td><td>0.14803143</td><td>0.2584459</td><td>0.138593119</td></tr>\n",
       "\t<tr><th scope=row>S001PSim</th><td>total  </td><td>Brasil                               </td><td>2019</td><td>S001P</td><td>0.9802181</td><td>0.97357997</td><td>0.9868562</td><td>0.003455195</td></tr>\n",
       "\t<tr><th scope=row>S001PSim1</th><td>total  </td><td>Capital                              </td><td>2019</td><td>S001P</td><td>0.9806496</td><td>0.97210024</td><td>0.9891990</td><td>0.004448070</td></tr>\n",
       "\t<tr><th scope=row>S002PSim</th><td>total  </td><td>Brasil                               </td><td>2019</td><td>S002P</td><td>0.9624012</td><td>0.94926198</td><td>0.9755404</td><td>0.006965703</td></tr>\n",
       "\t<tr><th scope=row>S002PSim1</th><td>total  </td><td>Capital                              </td><td>2019</td><td>S002P</td><td>0.9341446</td><td>0.90502841</td><td>0.9632609</td><td>0.015902776</td></tr>\n",
       "\t<tr><th scope=row>S003PSim</th><td>total  </td><td>Brasil                               </td><td>2019</td><td>S003P</td><td>0.8814929</td><td>0.86308159</td><td>0.8999041</td><td>0.010656560</td></tr>\n",
       "\t<tr><th scope=row>S003PSim1</th><td>total  </td><td>Capital                              </td><td>2019</td><td>S003P</td><td>0.8898081</td><td>0.85474618</td><td>0.9248701</td><td>0.020104431</td></tr>\n",
       "\t<tr><th scope=row>S004PSim</th><td>total  </td><td>Brasil                               </td><td>2019</td><td>S004P</td><td>0.8939233</td><td>0.87486687</td><td>0.9129797</td><td>0.010876579</td></tr>\n",
       "\t<tr><th scope=row>S004PSim1</th><td>total  </td><td>Capital                              </td><td>2019</td><td>S004P</td><td>0.9086737</td><td>0.87852554</td><td>0.9388218</td><td>0.016927946</td></tr>\n",
       "\t<tr><th scope=row>S005PSim</th><td>total  </td><td>Brasil                               </td><td>2019</td><td>S005P</td><td>0.6712686</td><td>0.63830397</td><td>0.7042333</td><td>0.025055548</td></tr>\n",
       "\t<tr><th scope=row>S005PSim1</th><td>total  </td><td>Capital                              </td><td>2019</td><td>S005P</td><td>0.5467601</td><td>0.49247810</td><td>0.6010420</td><td>0.050653646</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 395 × 8\n",
       "\\begin{tabular}{r|llllllll}\n",
       "  & abr\\_tipo & abr\\_nome & Ano & Indicador & Sim & LowerS & UpperS & cvS\\\\\n",
       "  & <chr> & <fct> & <chr> & <chr> & <dbl> & <dbl> & <dbl> & <dbl>\\\\\n",
       "\\hline\n",
       "\tBranca & raça            & Branca              & 2019 & S001P & 0.9848342 & 0.9728009 & 0.9968675 & 6.234103e-03\\\\\n",
       "\tPreta & raça            & Preta               & 2019 & S001P & 0.9706214 & 0.9496967 & 0.9915462 & 1.099923e-02\\\\\n",
       "\tParda & raça            & Parda               & 2019 & S001P & 0.9802893 & 0.9718436 & 0.9887350 & 4.395738e-03\\\\\n",
       "\tAté 1/2 SM & rend\\_per\\_capita & Até 1/2 SM          & 2019 & S001P & 0.9748696 & 0.9634649 & 0.9862743 & 5.968842e-03\\\\\n",
       "\t1/2 até 1 SM & rend\\_per\\_capita & 1/2 até 1 SM        & 2019 & S001P & 0.9809000 & 0.9691323 & 0.9926676 & 6.120930e-03\\\\\n",
       "\t1 até 2 SM & rend\\_per\\_capita & 1 até 2 SM          & 2019 & S001P & 0.9900064 & 0.9798280 & 1.0001847 & 5.245551e-03\\\\\n",
       "\t2 até 3 SM & rend\\_per\\_capita & 2 até 3 SM          & 2019 & S001P & 0.9874788 & 0.9684689 & 1.0064888 & 9.822109e-03\\\\\n",
       "\tMais de 3 SM & rend\\_per\\_capita & Mais de 3 SM        & 2019 & S001P & 0.9850570 & 0.9559967 & 1.0141173 & 1.505188e-02\\\\\n",
       "\t18 a 24 anos & fx\\_idade\\_S      & 18 a 24 anos        & 2019 & S001P & 0.9780003 & 0.9636430 & 0.9923577 & 7.490081e-03\\\\\n",
       "\t25 a 29 anos & fx\\_idade\\_S      & 25 a 29 anos        & 2019 & S001P & 0.9761319 & 0.9599660 & 0.9922979 & 8.449762e-03\\\\\n",
       "\t30 a 39 anos & fx\\_idade\\_S      & 30 a 39 anos        & 2019 & S001P & 0.9872648 & 0.9801366 & 0.9943931 & 3.683855e-03\\\\\n",
       "\t40 anos ou mais & fx\\_idade\\_S      & 40 anos ou mais     & 2019 & S001P & 0.9584802 & 0.9286191 & 0.9883414 & 1.589556e-02\\\\\n",
       "\turbano & urb\\_rur         & urbano              & 2019 & S001P & 0.9822630 & 0.9752768 & 0.9892492 & 3.628817e-03\\\\\n",
       "\trural & urb\\_rur         & rural               & 2019 & S001P & 0.9690309 & 0.9495219 & 0.9885399 & 1.027184e-02\\\\\n",
       "\tRondônia & uf              & Rondônia            & 2019 & S001P & 0.9922432 & 0.9813513 & 1.0031351 & 5.600634e-03\\\\\n",
       "\tAcre & uf              & Acre                & 2019 & S001P & 0.9400931 & 0.8740058 & 1.0061804 & 3.586734e-02\\\\\n",
       "\tAmazonas & uf              & Amazonas            & 2019 & S001P & 0.9471638 & 0.9097590 & 0.9845686 & 2.014902e-02\\\\\n",
       "\tRoraima & uf              & Roraima             & 2019 & S001P & 0.9381625 & 0.8778433 & 0.9984818 & 3.280420e-02\\\\\n",
       "\tPará & uf              & Pará                & 2019 & S001P & 0.9753771 & 0.9511006 & 0.9996535 & 1.269886e-02\\\\\n",
       "\tAmapá & uf              & Amapá               & 2019 & S001P & 0.9382793 & 0.8674707 & 1.0090880 & 3.850401e-02\\\\\n",
       "\tTocantins & uf              & Tocantins           & 2019 & S001P & 1.0000000 & 1.0000000 & 1.0000000 & 0.000000e+00\\\\\n",
       "\tMaranhão & uf              & Maranhão            & 2019 & S001P & 0.9882002 & 0.9774606 & 0.9989399 & 5.544946e-03\\\\\n",
       "\tPiauí & uf              & Piauí               & 2019 & S001P & 1.0000000 & 1.0000000 & 1.0000000 & 0.000000e+00\\\\\n",
       "\tCeará & uf              & Ceará               & 2019 & S001P & 0.9908578 & 0.9730003 & 1.0087154 & 9.195250e-03\\\\\n",
       "\tRio Grande do Norte & uf              & Rio Grande do Norte & 2019 & S001P & 0.9599420 & 0.9283055 & 0.9915785 & 1.681495e-02\\\\\n",
       "\tParaíba & uf              & Paraíba             & 2019 & S001P & 0.9259756 & 0.8549651 & 0.9969862 & 3.912687e-02\\\\\n",
       "\tPernambuco & uf              & Pernambuco          & 2019 & S001P & 0.9653698 & 0.9289815 & 1.0017582 & 1.923183e-02\\\\\n",
       "\tAlagoas & uf              & Alagoas             & 2019 & S001P & 0.9641940 & 0.9305357 & 0.9978523 & 1.781065e-02\\\\\n",
       "\tSergipe & uf              & Sergipe             & 2019 & S001P & 1.0000000 & 1.0000000 & 1.0000000 & 1.782482e-17\\\\\n",
       "\tBahia & uf              & Bahia               & 2019 & S001P & 0.9956594 & 0.9871376 & 1.0041812 & 4.366894e-03\\\\\n",
       "\t⋮ & ⋮ & ⋮ & ⋮ & ⋮ & ⋮ & ⋮ & ⋮ & ⋮\\\\\n",
       "\tJoão Pessoa4 & capital & João Pessoa                           & 2019 & S005P & 0.4891024 & 0.31333773 & 0.6648670 & 0.183351143\\\\\n",
       "\tRecife4 & capital & Recife                                & 2019 & S005P & 0.6154032 & 0.35946656 & 0.8713398 & 0.212189839\\\\\n",
       "\tMaceió4 & capital & Maceió                                & 2019 & S005P & 0.6770275 & 0.53995647 & 0.8140986 & 0.103297871\\\\\n",
       "\tAracaju4 & capital & Aracaju                               & 2019 & S005P & 0.5389670 & 0.28151970 & 0.7964143 & 0.243712643\\\\\n",
       "\tSalvador4 & capital & Salvador                              & 2019 & S005P & 0.5306442 & 0.26607586 & 0.7952126 & 0.254381995\\\\\n",
       "\tBelo Horizonte4 & capital & Belo Horizonte                        & 2019 & S005P & 0.3741457 & 0.20269372 & 0.5455976 & 0.233804895\\\\\n",
       "\tVitória4 & capital & Vitória                               & 2019 & S005P & 0.4852118 & 0.20101229 & 0.7694113 & 0.298843552\\\\\n",
       "\tRio de Janeiro9 & capital & Rio de Janeiro                        & 2019 & S005P & 0.5061534 & 0.27485324 & 0.7374535 & 0.233155491\\\\\n",
       "\tSão Paulo8 & capital & São Paulo                             & 2019 & S005P & 0.3925847 & 0.21472169 & 0.5704477 & 0.231155467\\\\\n",
       "\tCuritiba4 & capital & Curitiba                              & 2019 & S005P & 0.2620324 & 0.03101129 & 0.4930534 & 0.449830155\\\\\n",
       "\tFlorianópolis4 & capital & Florianópolis                         & 2019 & S005P & 0.5773166 & 0.31513885 & 0.8394944 & 0.231704095\\\\\n",
       "\tPorto Alegre4 & capital & Porto Alegre                          & 2019 & S005P & 0.4038861 & 0.17145033 & 0.6363218 & 0.293626959\\\\\n",
       "\tCampo Grande4 & capital & Campo Grande                          & 2019 & S005P & 0.6491680 & 0.50650967 & 0.7918262 & 0.112122258\\\\\n",
       "\tCuiabá4 & capital & Cuiabá                                & 2019 & S005P & 0.6268727 & 0.37278925 & 0.8809562 & 0.206799253\\\\\n",
       "\tGoiânia4 & capital & Goiânia                               & 2019 & S005P & 0.5727847 & 0.32629493 & 0.8192745 & 0.219563127\\\\\n",
       "\tBrasília4 & capital & Brasília                              & 2019 & S005P & 0.6320122 & 0.47297955 & 0.7910448 & 0.128384529\\\\\n",
       "\tFundamental incompleto ou equivalente4 & gescol  & Fundamental incompleto ou equivalente & 2019 & S005P & 0.9124938 & 0.87271915 & 0.9522685 & 0.022239685\\\\\n",
       "\tMédio incompleto ou equivalente4 & gescol  & Médio incompleto ou equivalente       & 2019 & S005P & 0.9012396 & 0.85482924 & 0.9476500 & 0.026274048\\\\\n",
       "\tSuperior incompleto ou equivalente4 & gescol  & Superior incompleto ou equivalente    & 2019 & S005P & 0.6802039 & 0.63250129 & 0.7279064 & 0.035781180\\\\\n",
       "\tSuperior completo4 & gescol  & Superior completo                     & 2019 & S005P & 0.2032387 & 0.14803143 & 0.2584459 & 0.138593119\\\\\n",
       "\tS001PSim & total   & Brasil                                & 2019 & S001P & 0.9802181 & 0.97357997 & 0.9868562 & 0.003455195\\\\\n",
       "\tS001PSim1 & total   & Capital                               & 2019 & S001P & 0.9806496 & 0.97210024 & 0.9891990 & 0.004448070\\\\\n",
       "\tS002PSim & total   & Brasil                                & 2019 & S002P & 0.9624012 & 0.94926198 & 0.9755404 & 0.006965703\\\\\n",
       "\tS002PSim1 & total   & Capital                               & 2019 & S002P & 0.9341446 & 0.90502841 & 0.9632609 & 0.015902776\\\\\n",
       "\tS003PSim & total   & Brasil                                & 2019 & S003P & 0.8814929 & 0.86308159 & 0.8999041 & 0.010656560\\\\\n",
       "\tS003PSim1 & total   & Capital                               & 2019 & S003P & 0.8898081 & 0.85474618 & 0.9248701 & 0.020104431\\\\\n",
       "\tS004PSim & total   & Brasil                                & 2019 & S004P & 0.8939233 & 0.87486687 & 0.9129797 & 0.010876579\\\\\n",
       "\tS004PSim1 & total   & Capital                               & 2019 & S004P & 0.9086737 & 0.87852554 & 0.9388218 & 0.016927946\\\\\n",
       "\tS005PSim & total   & Brasil                                & 2019 & S005P & 0.6712686 & 0.63830397 & 0.7042333 & 0.025055548\\\\\n",
       "\tS005PSim1 & total   & Capital                               & 2019 & S005P & 0.5467601 & 0.49247810 & 0.6010420 & 0.050653646\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 395 × 8\n",
       "\n",
       "| <!--/--> | abr_tipo &lt;chr&gt; | abr_nome &lt;fct&gt; | Ano &lt;chr&gt; | Indicador &lt;chr&gt; | Sim &lt;dbl&gt; | LowerS &lt;dbl&gt; | UpperS &lt;dbl&gt; | cvS &lt;dbl&gt; |\n",
       "|---|---|---|---|---|---|---|---|---|\n",
       "| Branca | raça            | Branca              | 2019 | S001P | 0.9848342 | 0.9728009 | 0.9968675 | 6.234103e-03 |\n",
       "| Preta | raça            | Preta               | 2019 | S001P | 0.9706214 | 0.9496967 | 0.9915462 | 1.099923e-02 |\n",
       "| Parda | raça            | Parda               | 2019 | S001P | 0.9802893 | 0.9718436 | 0.9887350 | 4.395738e-03 |\n",
       "| Até 1/2 SM | rend_per_capita | Até 1/2 SM          | 2019 | S001P | 0.9748696 | 0.9634649 | 0.9862743 | 5.968842e-03 |\n",
       "| 1/2 até 1 SM | rend_per_capita | 1/2 até 1 SM        | 2019 | S001P | 0.9809000 | 0.9691323 | 0.9926676 | 6.120930e-03 |\n",
       "| 1 até 2 SM | rend_per_capita | 1 até 2 SM          | 2019 | S001P | 0.9900064 | 0.9798280 | 1.0001847 | 5.245551e-03 |\n",
       "| 2 até 3 SM | rend_per_capita | 2 até 3 SM          | 2019 | S001P | 0.9874788 | 0.9684689 | 1.0064888 | 9.822109e-03 |\n",
       "| Mais de 3 SM | rend_per_capita | Mais de 3 SM        | 2019 | S001P | 0.9850570 | 0.9559967 | 1.0141173 | 1.505188e-02 |\n",
       "| 18 a 24 anos | fx_idade_S      | 18 a 24 anos        | 2019 | S001P | 0.9780003 | 0.9636430 | 0.9923577 | 7.490081e-03 |\n",
       "| 25 a 29 anos | fx_idade_S      | 25 a 29 anos        | 2019 | S001P | 0.9761319 | 0.9599660 | 0.9922979 | 8.449762e-03 |\n",
       "| 30 a 39 anos | fx_idade_S      | 30 a 39 anos        | 2019 | S001P | 0.9872648 | 0.9801366 | 0.9943931 | 3.683855e-03 |\n",
       "| 40 anos ou mais | fx_idade_S      | 40 anos ou mais     | 2019 | S001P | 0.9584802 | 0.9286191 | 0.9883414 | 1.589556e-02 |\n",
       "| urbano | urb_rur         | urbano              | 2019 | S001P | 0.9822630 | 0.9752768 | 0.9892492 | 3.628817e-03 |\n",
       "| rural | urb_rur         | rural               | 2019 | S001P | 0.9690309 | 0.9495219 | 0.9885399 | 1.027184e-02 |\n",
       "| Rondônia | uf              | Rondônia            | 2019 | S001P | 0.9922432 | 0.9813513 | 1.0031351 | 5.600634e-03 |\n",
       "| Acre | uf              | Acre                | 2019 | S001P | 0.9400931 | 0.8740058 | 1.0061804 | 3.586734e-02 |\n",
       "| Amazonas | uf              | Amazonas            | 2019 | S001P | 0.9471638 | 0.9097590 | 0.9845686 | 2.014902e-02 |\n",
       "| Roraima | uf              | Roraima             | 2019 | S001P | 0.9381625 | 0.8778433 | 0.9984818 | 3.280420e-02 |\n",
       "| Pará | uf              | Pará                | 2019 | S001P | 0.9753771 | 0.9511006 | 0.9996535 | 1.269886e-02 |\n",
       "| Amapá | uf              | Amapá               | 2019 | S001P | 0.9382793 | 0.8674707 | 1.0090880 | 3.850401e-02 |\n",
       "| Tocantins | uf              | Tocantins           | 2019 | S001P | 1.0000000 | 1.0000000 | 1.0000000 | 0.000000e+00 |\n",
       "| Maranhão | uf              | Maranhão            | 2019 | S001P | 0.9882002 | 0.9774606 | 0.9989399 | 5.544946e-03 |\n",
       "| Piauí | uf              | Piauí               | 2019 | S001P | 1.0000000 | 1.0000000 | 1.0000000 | 0.000000e+00 |\n",
       "| Ceará | uf              | Ceará               | 2019 | S001P | 0.9908578 | 0.9730003 | 1.0087154 | 9.195250e-03 |\n",
       "| Rio Grande do Norte | uf              | Rio Grande do Norte | 2019 | S001P | 0.9599420 | 0.9283055 | 0.9915785 | 1.681495e-02 |\n",
       "| Paraíba | uf              | Paraíba             | 2019 | S001P | 0.9259756 | 0.8549651 | 0.9969862 | 3.912687e-02 |\n",
       "| Pernambuco | uf              | Pernambuco          | 2019 | S001P | 0.9653698 | 0.9289815 | 1.0017582 | 1.923183e-02 |\n",
       "| Alagoas | uf              | Alagoas             | 2019 | S001P | 0.9641940 | 0.9305357 | 0.9978523 | 1.781065e-02 |\n",
       "| Sergipe | uf              | Sergipe             | 2019 | S001P | 1.0000000 | 1.0000000 | 1.0000000 | 1.782482e-17 |\n",
       "| Bahia | uf              | Bahia               | 2019 | S001P | 0.9956594 | 0.9871376 | 1.0041812 | 4.366894e-03 |\n",
       "| ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ | ⋮ |\n",
       "| João Pessoa4 | capital | João Pessoa                           | 2019 | S005P | 0.4891024 | 0.31333773 | 0.6648670 | 0.183351143 |\n",
       "| Recife4 | capital | Recife                                | 2019 | S005P | 0.6154032 | 0.35946656 | 0.8713398 | 0.212189839 |\n",
       "| Maceió4 | capital | Maceió                                | 2019 | S005P | 0.6770275 | 0.53995647 | 0.8140986 | 0.103297871 |\n",
       "| Aracaju4 | capital | Aracaju                               | 2019 | S005P | 0.5389670 | 0.28151970 | 0.7964143 | 0.243712643 |\n",
       "| Salvador4 | capital | Salvador                              | 2019 | S005P | 0.5306442 | 0.26607586 | 0.7952126 | 0.254381995 |\n",
       "| Belo Horizonte4 | capital | Belo Horizonte                        | 2019 | S005P | 0.3741457 | 0.20269372 | 0.5455976 | 0.233804895 |\n",
       "| Vitória4 | capital | Vitória                               | 2019 | S005P | 0.4852118 | 0.20101229 | 0.7694113 | 0.298843552 |\n",
       "| Rio de Janeiro9 | capital | Rio de Janeiro                        | 2019 | S005P | 0.5061534 | 0.27485324 | 0.7374535 | 0.233155491 |\n",
       "| São Paulo8 | capital | São Paulo                             | 2019 | S005P | 0.3925847 | 0.21472169 | 0.5704477 | 0.231155467 |\n",
       "| Curitiba4 | capital | Curitiba                              | 2019 | S005P | 0.2620324 | 0.03101129 | 0.4930534 | 0.449830155 |\n",
       "| Florianópolis4 | capital | Florianópolis                         | 2019 | S005P | 0.5773166 | 0.31513885 | 0.8394944 | 0.231704095 |\n",
       "| Porto Alegre4 | capital | Porto Alegre                          | 2019 | S005P | 0.4038861 | 0.17145033 | 0.6363218 | 0.293626959 |\n",
       "| Campo Grande4 | capital | Campo Grande                          | 2019 | S005P | 0.6491680 | 0.50650967 | 0.7918262 | 0.112122258 |\n",
       "| Cuiabá4 | capital | Cuiabá                                | 2019 | S005P | 0.6268727 | 0.37278925 | 0.8809562 | 0.206799253 |\n",
       "| Goiânia4 | capital | Goiânia                               | 2019 | S005P | 0.5727847 | 0.32629493 | 0.8192745 | 0.219563127 |\n",
       "| Brasília4 | capital | Brasília                              | 2019 | S005P | 0.6320122 | 0.47297955 | 0.7910448 | 0.128384529 |\n",
       "| Fundamental incompleto ou equivalente4 | gescol  | Fundamental incompleto ou equivalente | 2019 | S005P | 0.9124938 | 0.87271915 | 0.9522685 | 0.022239685 |\n",
       "| Médio incompleto ou equivalente4 | gescol  | Médio incompleto ou equivalente       | 2019 | S005P | 0.9012396 | 0.85482924 | 0.9476500 | 0.026274048 |\n",
       "| Superior incompleto ou equivalente4 | gescol  | Superior incompleto ou equivalente    | 2019 | S005P | 0.6802039 | 0.63250129 | 0.7279064 | 0.035781180 |\n",
       "| Superior completo4 | gescol  | Superior completo                     | 2019 | S005P | 0.2032387 | 0.14803143 | 0.2584459 | 0.138593119 |\n",
       "| S001PSim | total   | Brasil                                | 2019 | S001P | 0.9802181 | 0.97357997 | 0.9868562 | 0.003455195 |\n",
       "| S001PSim1 | total   | Capital                               | 2019 | S001P | 0.9806496 | 0.97210024 | 0.9891990 | 0.004448070 |\n",
       "| S002PSim | total   | Brasil                                | 2019 | S002P | 0.9624012 | 0.94926198 | 0.9755404 | 0.006965703 |\n",
       "| S002PSim1 | total   | Capital                               | 2019 | S002P | 0.9341446 | 0.90502841 | 0.9632609 | 0.015902776 |\n",
       "| S003PSim | total   | Brasil                                | 2019 | S003P | 0.8814929 | 0.86308159 | 0.8999041 | 0.010656560 |\n",
       "| S003PSim1 | total   | Capital                               | 2019 | S003P | 0.8898081 | 0.85474618 | 0.9248701 | 0.020104431 |\n",
       "| S004PSim | total   | Brasil                                | 2019 | S004P | 0.8939233 | 0.87486687 | 0.9129797 | 0.010876579 |\n",
       "| S004PSim1 | total   | Capital                               | 2019 | S004P | 0.9086737 | 0.87852554 | 0.9388218 | 0.016927946 |\n",
       "| S005PSim | total   | Brasil                                | 2019 | S005P | 0.6712686 | 0.63830397 | 0.7042333 | 0.025055548 |\n",
       "| S005PSim1 | total   | Capital                               | 2019 | S005P | 0.5467601 | 0.49247810 | 0.6010420 | 0.050653646 |\n",
       "\n"
      ],
      "text/plain": [
       "                                       abr_tipo       \n",
       "Branca                                 raça           \n",
       "Preta                                  raça           \n",
       "Parda                                  raça           \n",
       "Até 1/2 SM                             rend_per_capita\n",
       "1/2 até 1 SM                           rend_per_capita\n",
       "1 até 2 SM                             rend_per_capita\n",
       "2 até 3 SM                             rend_per_capita\n",
       "Mais de 3 SM                           rend_per_capita\n",
       "18 a 24 anos                           fx_idade_S     \n",
       "25 a 29 anos                           fx_idade_S     \n",
       "30 a 39 anos                           fx_idade_S     \n",
       "40 anos ou mais                        fx_idade_S     \n",
       "urbano                                 urb_rur        \n",
       "rural                                  urb_rur        \n",
       "Rondônia                               uf             \n",
       "Acre                                   uf             \n",
       "Amazonas                               uf             \n",
       "Roraima                                uf             \n",
       "Pará                                   uf             \n",
       "Amapá                                  uf             \n",
       "Tocantins                              uf             \n",
       "Maranhão                               uf             \n",
       "Piauí                                  uf             \n",
       "Ceará                                  uf             \n",
       "Rio Grande do Norte                    uf             \n",
       "Paraíba                                uf             \n",
       "Pernambuco                             uf             \n",
       "Alagoas                                uf             \n",
       "Sergipe                                uf             \n",
       "Bahia                                  uf             \n",
       "⋮                                      ⋮              \n",
       "João Pessoa4                           capital        \n",
       "Recife4                                capital        \n",
       "Maceió4                                capital        \n",
       "Aracaju4                               capital        \n",
       "Salvador4                              capital        \n",
       "Belo Horizonte4                        capital        \n",
       "Vitória4                               capital        \n",
       "Rio de Janeiro9                        capital        \n",
       "São Paulo8                             capital        \n",
       "Curitiba4                              capital        \n",
       "Florianópolis4                         capital        \n",
       "Porto Alegre4                          capital        \n",
       "Campo Grande4                          capital        \n",
       "Cuiabá4                                capital        \n",
       "Goiânia4                               capital        \n",
       "Brasília4                              capital        \n",
       "Fundamental incompleto ou equivalente4 gescol         \n",
       "Médio incompleto ou equivalente4       gescol         \n",
       "Superior incompleto ou equivalente4    gescol         \n",
       "Superior completo4                     gescol         \n",
       "S001PSim                               total          \n",
       "S001PSim1                              total          \n",
       "S002PSim                               total          \n",
       "S002PSim1                              total          \n",
       "S003PSim                               total          \n",
       "S003PSim1                              total          \n",
       "S004PSim                               total          \n",
       "S004PSim1                              total          \n",
       "S005PSim                               total          \n",
       "S005PSim1                              total          \n",
       "                                       abr_nome                             \n",
       "Branca                                 Branca                               \n",
       "Preta                                  Preta                                \n",
       "Parda                                  Parda                                \n",
       "Até 1/2 SM                             Até 1/2 SM                           \n",
       "1/2 até 1 SM                           1/2 até 1 SM                         \n",
       "1 até 2 SM                             1 até 2 SM                           \n",
       "2 até 3 SM                             2 até 3 SM                           \n",
       "Mais de 3 SM                           Mais de 3 SM                         \n",
       "18 a 24 anos                           18 a 24 anos                         \n",
       "25 a 29 anos                           25 a 29 anos                         \n",
       "30 a 39 anos                           30 a 39 anos                         \n",
       "40 anos ou mais                        40 anos ou mais                      \n",
       "urbano                                 urbano                               \n",
       "rural                                  rural                                \n",
       "Rondônia                               Rondônia                             \n",
       "Acre                                   Acre                                 \n",
       "Amazonas                               Amazonas                             \n",
       "Roraima                                Roraima                              \n",
       "Pará                                   Pará                                 \n",
       "Amapá                                  Amapá                                \n",
       "Tocantins                              Tocantins                            \n",
       "Maranhão                               Maranhão                             \n",
       "Piauí                                  Piauí                                \n",
       "Ceará                                  Ceará                                \n",
       "Rio Grande do Norte                    Rio Grande do Norte                  \n",
       "Paraíba                                Paraíba                              \n",
       "Pernambuco                             Pernambuco                           \n",
       "Alagoas                                Alagoas                              \n",
       "Sergipe                                Sergipe                              \n",
       "Bahia                                  Bahia                                \n",
       "⋮                                      ⋮                                    \n",
       "João Pessoa4                           João Pessoa                          \n",
       "Recife4                                Recife                               \n",
       "Maceió4                                Maceió                               \n",
       "Aracaju4                               Aracaju                              \n",
       "Salvador4                              Salvador                             \n",
       "Belo Horizonte4                        Belo Horizonte                       \n",
       "Vitória4                               Vitória                              \n",
       "Rio de Janeiro9                        Rio de Janeiro                       \n",
       "São Paulo8                             São Paulo                            \n",
       "Curitiba4                              Curitiba                             \n",
       "Florianópolis4                         Florianópolis                        \n",
       "Porto Alegre4                          Porto Alegre                         \n",
       "Campo Grande4                          Campo Grande                         \n",
       "Cuiabá4                                Cuiabá                               \n",
       "Goiânia4                               Goiânia                              \n",
       "Brasília4                              Brasília                             \n",
       "Fundamental incompleto ou equivalente4 Fundamental incompleto ou equivalente\n",
       "Médio incompleto ou equivalente4       Médio incompleto ou equivalente      \n",
       "Superior incompleto ou equivalente4    Superior incompleto ou equivalente   \n",
       "Superior completo4                     Superior completo                    \n",
       "S001PSim                               Brasil                               \n",
       "S001PSim1                              Capital                              \n",
       "S002PSim                               Brasil                               \n",
       "S002PSim1                              Capital                              \n",
       "S003PSim                               Brasil                               \n",
       "S003PSim1                              Capital                              \n",
       "S004PSim                               Brasil                               \n",
       "S004PSim1                              Capital                              \n",
       "S005PSim                               Brasil                               \n",
       "S005PSim1                              Capital                              \n",
       "                                       Ano  Indicador Sim       LowerS    \n",
       "Branca                                 2019 S001P     0.9848342 0.9728009 \n",
       "Preta                                  2019 S001P     0.9706214 0.9496967 \n",
       "Parda                                  2019 S001P     0.9802893 0.9718436 \n",
       "Até 1/2 SM                             2019 S001P     0.9748696 0.9634649 \n",
       "1/2 até 1 SM                           2019 S001P     0.9809000 0.9691323 \n",
       "1 até 2 SM                             2019 S001P     0.9900064 0.9798280 \n",
       "2 até 3 SM                             2019 S001P     0.9874788 0.9684689 \n",
       "Mais de 3 SM                           2019 S001P     0.9850570 0.9559967 \n",
       "18 a 24 anos                           2019 S001P     0.9780003 0.9636430 \n",
       "25 a 29 anos                           2019 S001P     0.9761319 0.9599660 \n",
       "30 a 39 anos                           2019 S001P     0.9872648 0.9801366 \n",
       "40 anos ou mais                        2019 S001P     0.9584802 0.9286191 \n",
       "urbano                                 2019 S001P     0.9822630 0.9752768 \n",
       "rural                                  2019 S001P     0.9690309 0.9495219 \n",
       "Rondônia                               2019 S001P     0.9922432 0.9813513 \n",
       "Acre                                   2019 S001P     0.9400931 0.8740058 \n",
       "Amazonas                               2019 S001P     0.9471638 0.9097590 \n",
       "Roraima                                2019 S001P     0.9381625 0.8778433 \n",
       "Pará                                   2019 S001P     0.9753771 0.9511006 \n",
       "Amapá                                  2019 S001P     0.9382793 0.8674707 \n",
       "Tocantins                              2019 S001P     1.0000000 1.0000000 \n",
       "Maranhão                               2019 S001P     0.9882002 0.9774606 \n",
       "Piauí                                  2019 S001P     1.0000000 1.0000000 \n",
       "Ceará                                  2019 S001P     0.9908578 0.9730003 \n",
       "Rio Grande do Norte                    2019 S001P     0.9599420 0.9283055 \n",
       "Paraíba                                2019 S001P     0.9259756 0.8549651 \n",
       "Pernambuco                             2019 S001P     0.9653698 0.9289815 \n",
       "Alagoas                                2019 S001P     0.9641940 0.9305357 \n",
       "Sergipe                                2019 S001P     1.0000000 1.0000000 \n",
       "Bahia                                  2019 S001P     0.9956594 0.9871376 \n",
       "⋮                                      ⋮    ⋮         ⋮         ⋮         \n",
       "João Pessoa4                           2019 S005P     0.4891024 0.31333773\n",
       "Recife4                                2019 S005P     0.6154032 0.35946656\n",
       "Maceió4                                2019 S005P     0.6770275 0.53995647\n",
       "Aracaju4                               2019 S005P     0.5389670 0.28151970\n",
       "Salvador4                              2019 S005P     0.5306442 0.26607586\n",
       "Belo Horizonte4                        2019 S005P     0.3741457 0.20269372\n",
       "Vitória4                               2019 S005P     0.4852118 0.20101229\n",
       "Rio de Janeiro9                        2019 S005P     0.5061534 0.27485324\n",
       "São Paulo8                             2019 S005P     0.3925847 0.21472169\n",
       "Curitiba4                              2019 S005P     0.2620324 0.03101129\n",
       "Florianópolis4                         2019 S005P     0.5773166 0.31513885\n",
       "Porto Alegre4                          2019 S005P     0.4038861 0.17145033\n",
       "Campo Grande4                          2019 S005P     0.6491680 0.50650967\n",
       "Cuiabá4                                2019 S005P     0.6268727 0.37278925\n",
       "Goiânia4                               2019 S005P     0.5727847 0.32629493\n",
       "Brasília4                              2019 S005P     0.6320122 0.47297955\n",
       "Fundamental incompleto ou equivalente4 2019 S005P     0.9124938 0.87271915\n",
       "Médio incompleto ou equivalente4       2019 S005P     0.9012396 0.85482924\n",
       "Superior incompleto ou equivalente4    2019 S005P     0.6802039 0.63250129\n",
       "Superior completo4                     2019 S005P     0.2032387 0.14803143\n",
       "S001PSim                               2019 S001P     0.9802181 0.97357997\n",
       "S001PSim1                              2019 S001P     0.9806496 0.97210024\n",
       "S002PSim                               2019 S002P     0.9624012 0.94926198\n",
       "S002PSim1                              2019 S002P     0.9341446 0.90502841\n",
       "S003PSim                               2019 S003P     0.8814929 0.86308159\n",
       "S003PSim1                              2019 S003P     0.8898081 0.85474618\n",
       "S004PSim                               2019 S004P     0.8939233 0.87486687\n",
       "S004PSim1                              2019 S004P     0.9086737 0.87852554\n",
       "S005PSim                               2019 S005P     0.6712686 0.63830397\n",
       "S005PSim1                              2019 S005P     0.5467601 0.49247810\n",
       "                                       UpperS    cvS         \n",
       "Branca                                 0.9968675 6.234103e-03\n",
       "Preta                                  0.9915462 1.099923e-02\n",
       "Parda                                  0.9887350 4.395738e-03\n",
       "Até 1/2 SM                             0.9862743 5.968842e-03\n",
       "1/2 até 1 SM                           0.9926676 6.120930e-03\n",
       "1 até 2 SM                             1.0001847 5.245551e-03\n",
       "2 até 3 SM                             1.0064888 9.822109e-03\n",
       "Mais de 3 SM                           1.0141173 1.505188e-02\n",
       "18 a 24 anos                           0.9923577 7.490081e-03\n",
       "25 a 29 anos                           0.9922979 8.449762e-03\n",
       "30 a 39 anos                           0.9943931 3.683855e-03\n",
       "40 anos ou mais                        0.9883414 1.589556e-02\n",
       "urbano                                 0.9892492 3.628817e-03\n",
       "rural                                  0.9885399 1.027184e-02\n",
       "Rondônia                               1.0031351 5.600634e-03\n",
       "Acre                                   1.0061804 3.586734e-02\n",
       "Amazonas                               0.9845686 2.014902e-02\n",
       "Roraima                                0.9984818 3.280420e-02\n",
       "Pará                                   0.9996535 1.269886e-02\n",
       "Amapá                                  1.0090880 3.850401e-02\n",
       "Tocantins                              1.0000000 0.000000e+00\n",
       "Maranhão                               0.9989399 5.544946e-03\n",
       "Piauí                                  1.0000000 0.000000e+00\n",
       "Ceará                                  1.0087154 9.195250e-03\n",
       "Rio Grande do Norte                    0.9915785 1.681495e-02\n",
       "Paraíba                                0.9969862 3.912687e-02\n",
       "Pernambuco                             1.0017582 1.923183e-02\n",
       "Alagoas                                0.9978523 1.781065e-02\n",
       "Sergipe                                1.0000000 1.782482e-17\n",
       "Bahia                                  1.0041812 4.366894e-03\n",
       "⋮                                      ⋮         ⋮           \n",
       "João Pessoa4                           0.6648670 0.183351143 \n",
       "Recife4                                0.8713398 0.212189839 \n",
       "Maceió4                                0.8140986 0.103297871 \n",
       "Aracaju4                               0.7964143 0.243712643 \n",
       "Salvador4                              0.7952126 0.254381995 \n",
       "Belo Horizonte4                        0.5455976 0.233804895 \n",
       "Vitória4                               0.7694113 0.298843552 \n",
       "Rio de Janeiro9                        0.7374535 0.233155491 \n",
       "São Paulo8                             0.5704477 0.231155467 \n",
       "Curitiba4                              0.4930534 0.449830155 \n",
       "Florianópolis4                         0.8394944 0.231704095 \n",
       "Porto Alegre4                          0.6363218 0.293626959 \n",
       "Campo Grande4                          0.7918262 0.112122258 \n",
       "Cuiabá4                                0.8809562 0.206799253 \n",
       "Goiânia4                               0.8192745 0.219563127 \n",
       "Brasília4                              0.7910448 0.128384529 \n",
       "Fundamental incompleto ou equivalente4 0.9522685 0.022239685 \n",
       "Médio incompleto ou equivalente4       0.9476500 0.026274048 \n",
       "Superior incompleto ou equivalente4    0.7279064 0.035781180 \n",
       "Superior completo4                     0.2584459 0.138593119 \n",
       "S001PSim                               0.9868562 0.003455195 \n",
       "S001PSim1                              0.9891990 0.004448070 \n",
       "S002PSim                               0.9755404 0.006965703 \n",
       "S002PSim1                              0.9632609 0.015902776 \n",
       "S003PSim                               0.8999041 0.010656560 \n",
       "S003PSim1                              0.9248701 0.020104431 \n",
       "S004PSim                               0.9129797 0.010876579 \n",
       "S004PSim1                              0.9388218 0.016927946 \n",
       "S005PSim                               0.7042333 0.025055548 \n",
       "S005PSim1                              0.6010420 0.050653646 "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "matrizIndicadores"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Exportando tabela de indicadores calculados - Módulo S 2019"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [],
   "source": [
    "diretorio_saida <- \"<coloque aqui o diretório>\"\n",
    "write.table(matrizIndicadores,file=paste0(diretorio_saida,\"Indicadores_2019S_R.csv\"),sep = \";\",dec = \",\",row.names = FALSE)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R",
   "language": "R",
   "name": "ir"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "4.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}