2017-08-01 21:44:47 2017-08-01 21:44:47 R version 3.4.0 (2017-04-21) -- "You Stupid Darkness" 2017-08-01 21:44:47 Copyright (C) 2017 The R Foundation for Statistical Computing 2017-08-01 21:44:47 Platform: x86_64-w64-mingw32/x64 (64-bit) 2017-08-01 21:44:47 2017-08-01 21:44:47 R is free software and comes with ABSOLUTELY NO WARRANTY. 2017-08-01 21:44:47 You are welcome to redistribute it under certain conditions. 2017-08-01 21:44:47 Type 'license()' or 'licence()' for distribution details. 2017-08-01 21:44:47 2017-08-01 21:44:47 Natural language support but running in an English locale 2017-08-01 21:44:47 2017-08-01 21:44:47 R is a collaborative project with many contributors. 2017-08-01 21:44:47 Type 'contributors()' for more information and 2017-08-01 21:44:47 'citation()' on how to cite R or R packages in publications. 2017-08-01 21:44:47 2017-08-01 21:44:47 Type 'demo()' for some demos, 'help()' for on-line help, or 2017-08-01 21:44:47 'help.start()' for an HTML browser interface to help. 2017-08-01 21:44:47 Type 'q()' to quit R. 2017-08-01 21:44:47 2017-08-01 21:44:47 > setwd( Sys.getenv( "RWD" ) ) 2017-08-01 21:44:47 > 2017-08-01 21:44:47 > machine_specific_replacements <- 2017-08-01 21:44:47 + list( 2017-08-01 21:44:47 + 2017-08-01 21:44:47 + # replace the folder path on macnix 2017-08-01 21:44:47 + c( 'path.expand( \"~\" ) , \"PNS\"' , paste0( '"' , getwd() , '"' ) ) , 2017-08-01 21:44:47 + 2017-08-01 21:44:47 + # change other things in the script to be run 2017-08-01 21:44:47 + c( "hello" , "howdy" ) 2017-08-01 21:44:47 + 2017-08-01 21:44:47 + ) 2017-08-01 21:44:47 > 2017-08-01 21:44:47 > source( lodown::syntaxtractor( "pns" , replacements = machine_specific_replacements , setup_test = "setup" ) , echo = TRUE ) 2017-08-01 21:44:47 2017-08-01 21:45:50 2017-08-01 21:45:50 > library(lodown) 2017-08-01 21:45:50 2017-08-01 21:45:50 > lodown("pns", output_dir = file.path("e:/eanthony/r-wd/pns-1501649007")) 2017-08-01 21:45:50 building catalog for pns 2017-08-01 21:45:50 2017-08-01 21:45:56 locally downloading pns 2017-08-01 21:45:56 2017-08-01 21:45:56 'ftp://ftp.ibge.gov.br/PNS/2013/microdados/pns_2013_microdados_2017_03_23.zip' 2017-08-01 21:45:56 cached in 2017-08-01 21:45:56 'E:/eanthony/r-tmp/pns/9be0063137d7f8b42165dc4acdd89694.Rcache' 2017-08-01 21:45:56 copying to 2017-08-01 21:45:56 'E:\eanthony\r-tmp\pns\RtmpUt1m8v\file411c68b2338' 2017-08-01 21:45:56 2017-08-01 21:58:16 pns catalog entry 1 of 1 stored in 'e:/eanthony/r-wd/pns-1501649007/' 2017-08-01 21:58:16 2017-08-01 21:58:16 pns local download completed 2017-08-01 21:58:16 2017-08-01 21:58:16 Warning messages: 2017-08-01 21:58:16 1: In readLines(this_con) : 2017-08-01 21:58:16 incomplete final line found on 'E:/eanthony/r-tmp/pns/RtmpUt1m8v/unzips/Dicionarios_e_input/input_DOMPNS2013.sas' 2017-08-01 21:58:16 2: In readLines(this_con) : 2017-08-01 21:58:16 incomplete final line found on 'E:/eanthony/r-tmp/pns/RtmpUt1m8v/unzips/Dicionarios_e_input/input_PESPNS2013.sas' 2017-08-01 21:58:16 > 2017-08-01 21:58:17 2017-08-01 21:58:17 R version 3.4.0 (2017-04-21) -- "You Stupid Darkness" 2017-08-01 21:58:17 Copyright (C) 2017 The R Foundation for Statistical Computing 2017-08-01 21:58:17 Platform: x86_64-w64-mingw32/x64 (64-bit) 2017-08-01 21:58:17 2017-08-01 21:58:17 R is free software and comes with ABSOLUTELY NO WARRANTY. 2017-08-01 21:58:17 You are welcome to redistribute it under certain conditions. 2017-08-01 21:58:17 Type 'license()' or 'licence()' for distribution details. 2017-08-01 21:58:17 2017-08-01 21:58:17 Natural language support but running in an English locale 2017-08-01 21:58:17 2017-08-01 21:58:17 R is a collaborative project with many contributors. 2017-08-01 21:58:17 Type 'contributors()' for more information and 2017-08-01 21:58:17 'citation()' on how to cite R or R packages in publications. 2017-08-01 21:58:17 2017-08-01 21:58:17 Type 'demo()' for some demos, 'help()' for on-line help, or 2017-08-01 21:58:17 'help.start()' for an HTML browser interface to help. 2017-08-01 21:58:17 Type 'q()' to quit R. 2017-08-01 21:58:17 2017-08-01 21:58:17 > setwd( Sys.getenv( "RWD" ) ) 2017-08-01 21:58:17 > 2017-08-01 21:58:17 > machine_specific_replacements <- 2017-08-01 21:58:17 + list( 2017-08-01 21:58:17 + 2017-08-01 21:58:17 + # replace the folder path on macnix 2017-08-01 21:58:17 + c( 'path.expand( \"~\" ) , \"PNS\"' , paste0( '"' , getwd() , '"' ) ) , 2017-08-01 21:58:17 + 2017-08-01 21:58:17 + # change other things in the script to be run 2017-08-01 21:58:17 + c( "hello" , "howdy" ) 2017-08-01 21:58:17 + 2017-08-01 21:58:17 + ) 2017-08-01 21:58:17 > 2017-08-01 21:58:17 > source( lodown::syntaxtractor( "pns" , replacements = machine_specific_replacements , setup_test = "test" ) , echo = TRUE ) 2017-08-01 21:58:18 2017-08-01 21:58:18 > library(lodown) 2017-08-01 21:58:18 2017-08-01 21:58:18 > pns_cat <- get_catalog("pns", output_dir = file.path("e:/eanthony/r-wd/pns-1501649007")) 2017-08-01 21:58:18 building catalog for pns 2017-08-01 21:58:18 2017-08-01 21:58:24 2017-08-01 21:58:24 > stopifnot(nrow(pns_cat) > 0) 2017-08-01 21:58:24 2017-08-01 21:58:24 > options(survey.lonely.psu = "adjust") 2017-08-01 21:58:24 2017-08-01 21:58:24 > library(survey) 2017-08-01 21:58:24 Loading required package: grid 2017-08-01 21:58:24 Loading required package: Matrix 2017-08-01 21:58:26 Loading required package: survival 2017-08-01 21:58:26 2017-08-01 21:58:26 Attaching package: 'survey' 2017-08-01 21:58:26 2017-08-01 21:58:26 The following object is masked from 'package:graphics': 2017-08-01 21:58:26 2017-08-01 21:58:26 dotchart 2017-08-01 21:58:26 2017-08-01 21:58:26 2017-08-01 21:58:26 > pns_design <- readRDS(file.path("e:/eanthony/r-wd/pns-1501649007", 2017-08-01 21:58:26 + "2013 long questionnaire survey design.rds")) 2017-08-01 22:00:09 2017-08-01 22:00:09 > pns_design <- update(pns_design, one = 1, health_insurance = as.numeric(i001 == 2017-08-01 22:00:09 + 1)) 2017-08-01 22:00:09 2017-08-01 22:00:09 > sum(weights(pns_design, "sampling") != 0) 2017-08-01 22:00:09 [1] 60202 2017-08-01 22:00:09 2017-08-01 22:00:09 > svyby(~one, ~uf, pns_design, unwtd.count) 2017-08-01 22:00:10 uf counts se 2017-08-01 22:00:10 Rondonia Rondonia 1694 0 2017-08-01 22:00:10 Acre Acre 1814 0 2017-08-01 22:00:10 Amazonas Amazonas 2586 0 2017-08-01 22:00:10 Roraima Roraima 1591 0 2017-08-01 22:00:10 Para Para 2004 0 2017-08-01 22:00:10 Amapa Amapa 1332 0 2017-08-01 22:00:10 Tocantins Tocantins 1515 0 2017-08-01 22:00:10 Maranhao Maranhao 1774 0 2017-08-01 22:00:10 Piaui Piaui 1804 0 2017-08-01 22:00:10 Ceara Ceara 2560 0 2017-08-01 22:00:10 Rio Grande do Norte Rio Grande do Norte 1691 0 2017-08-01 22:00:10 Paraiba Paraiba 1943 0 2017-08-01 22:00:10 Pernambuco Pernambuco 2591 0 2017-08-01 22:00:10 Alagoas Alagoas 1748 0 2017-08-01 22:00:10 Sergipe Sergipe 1553 0 2017-08-01 22:00:10 Bahia Bahia 2641 0 2017-08-01 22:00:10 Minas Gerais Minas Gerais 3779 0 2017-08-01 22:00:10 Espirito Santo Espirito Santo 1724 0 2017-08-01 22:00:10 Rio de Janeiro Rio de Janeiro 3486 0 2017-08-01 22:00:10 Sao Paulo Sao Paulo 5305 0 2017-08-01 22:00:10 Parana Parana 3012 0 2017-08-01 22:00:10 Santa Catarina Santa Catarina 1623 0 2017-08-01 22:00:10 Rio Grande do Sul Rio Grande do Sul 2913 0 2017-08-01 22:00:10 Mato Grosso do Sul Mato Grosso do Sul 1809 0 2017-08-01 22:00:10 Mato Grosso Mato Grosso 1476 0 2017-08-01 22:00:10 Goias Goias 2423 0 2017-08-01 22:00:10 Distrito Federal Distrito Federal 1811 0 2017-08-01 22:00:10 2017-08-01 22:00:10 > svytotal(~one, pns_design) 2017-08-01 22:00:12 total SE 2017-08-01 22:00:12 one 146308458 0 2017-08-01 22:00:12 2017-08-01 22:00:12 > svyby(~one, ~uf, pns_design, svytotal) 2017-08-01 22:00:59 uf one se 2017-08-01 22:00:59 Rondonia Rondonia 1189854.3 1.345359e-11 2017-08-01 22:00:59 Acre Acre 479543.9 3.418715e-12 2017-08-01 22:00:59 Amazonas Amazonas 2282655.4 1.890779e-11 2017-08-01 22:00:59 Roraima Roraima 276092.4 2.813633e-12 2017-08-01 22:00:59 Para Para 5216618.5 4.012946e-11 2017-08-01 22:00:59 Amapa Amapa 453577.2 4.166983e-12 2017-08-01 22:00:59 Tocantins Tocantins 987626.7 1.320594e-11 2017-08-01 22:00:59 Maranhao Maranhao 4375124.3 5.527966e-11 2017-08-01 22:00:59 Piaui Piaui 2213132.9 2.648799e-11 2017-08-01 22:00:59 Ceara Ceara 6223082.2 4.228230e-11 2017-08-01 22:00:59 Rio Grande do Norte Rio Grande do Norte 2415675.1 1.981420e-11 2017-08-01 22:00:59 Paraiba Paraiba 2777844.2 3.401668e-11 2017-08-01 22:00:59 Pernambuco Pernambuco 6485947.2 4.879774e-11 2017-08-01 22:00:59 Alagoas Alagoas 2256410.4 2.014613e-11 2017-08-01 22:00:59 Sergipe Sergipe 1531995.8 1.476020e-11 2017-08-01 22:00:59 Bahia Bahia 10668362.4 1.756511e-10 2017-08-01 22:00:59 Minas Gerais Minas Gerais 15342139.8 1.916436e-10 2017-08-01 22:00:59 Espirito Santo Espirito Santo 2795688.0 4.469709e-11 2017-08-01 22:00:59 Rio de Janeiro Rio de Janeiro 12747633.7 8.981232e-11 2017-08-01 22:00:59 Sao Paulo Sao Paulo 33189220.7 4.392020e-10 2017-08-01 22:00:59 Parana Parana 8039355.6 6.829760e-11 2017-08-01 22:00:59 Santa Catarina Santa Catarina 5051440.6 1.015801e-10 2017-08-01 22:00:59 Rio Grande do Sul Rio Grande do Sul 8533867.7 6.045171e-11 2017-08-01 22:00:59 Mato Grosso do Sul Mato Grosso do Sul 1770609.6 3.830371e-11 2017-08-01 22:00:59 Mato Grosso Mato Grosso 2270730.4 2.276871e-11 2017-08-01 22:00:59 Goias Goias 4678041.3 5.301118e-11 2017-08-01 22:00:59 Distrito Federal Distrito Federal 2056187.8 2.249002e-11 2017-08-01 22:00:59 2017-08-01 22:00:59 > svymean(~w00101, pns_design, na.rm = TRUE) 2017-08-01 22:01:01 mean SE 2017-08-01 22:01:01 w00101 71.45 0.1152 2017-08-01 22:01:01 2017-08-01 22:01:01 > svyby(~w00101, ~uf, pns_design, svymean, na.rm = TRUE) 2017-08-01 22:01:44 uf w00101 se 2017-08-01 22:01:44 Rondonia Rondonia 70.57676 0.4500115 2017-08-01 22:01:44 Acre Acre 68.39216 0.4139088 2017-08-01 22:01:44 Amazonas Amazonas 69.56325 0.4035524 2017-08-01 22:01:44 Roraima Roraima 69.30896 0.4834325 2017-08-01 22:01:44 Para Para 67.16428 0.4063297 2017-08-01 22:01:44 Amapa Amapa 69.06946 0.4444721 2017-08-01 22:01:44 Tocantins Tocantins 68.09428 0.5208864 2017-08-01 22:01:44 Maranhao Maranhao 64.33809 0.4498913 2017-08-01 22:01:44 Piaui Piaui 65.32673 0.4205205 2017-08-01 22:01:44 Ceara Ceara 67.38869 0.3611066 2017-08-01 22:01:44 Rio Grande do Norte Rio Grande do Norte 70.20013 0.4876779 2017-08-01 22:01:44 Paraiba Paraiba 68.28660 0.4297494 2017-08-01 22:01:44 Pernambuco Pernambuco 70.33818 0.3672367 2017-08-01 22:01:44 Alagoas Alagoas 69.79392 0.4233045 2017-08-01 22:01:44 Sergipe Sergipe 70.47851 0.4411245 2017-08-01 22:01:44 Bahia Bahia 68.93628 0.5974439 2017-08-01 22:01:44 Minas Gerais Minas Gerais 70.58504 0.4722967 2017-08-01 22:01:44 Espirito Santo Espirito Santo 70.95230 0.4223929 2017-08-01 22:01:44 Rio de Janeiro Rio de Janeiro 73.91826 0.3141405 2017-08-01 22:01:44 Sao Paulo Sao Paulo 73.65539 0.3192958 2017-08-01 22:01:44 Parana Parana 72.54831 0.4019706 2017-08-01 22:01:44 Santa Catarina Santa Catarina 74.85520 0.5407879 2017-08-01 22:01:44 Rio Grande do Sul Rio Grande do Sul 75.01492 0.4049787 2017-08-01 22:01:44 Mato Grosso do Sul Mato Grosso do Sul 73.63023 0.4439434 2017-08-01 22:01:44 Mato Grosso Mato Grosso 73.36957 0.5221943 2017-08-01 22:01:44 Goias Goias 71.58753 0.4368148 2017-08-01 22:01:44 Distrito Federal Distrito Federal 71.82446 0.3726693 2017-08-01 22:01:44 2017-08-01 22:01:44 > svymean(~c006, pns_design) 2017-08-01 22:01:46 mean SE 2017-08-01 22:01:46 c006masculino 0.47104 0 2017-08-01 22:01:46 c006feminino 0.52896 0 2017-08-01 22:01:46 2017-08-01 22:01:46 > svyby(~c006, ~uf, pns_design, svymean) 2017-08-01 22:02:32 uf c006masculino c006feminino 2017-08-01 22:02:32 Rondonia Rondonia 0.4902478 0.5097522 2017-08-01 22:02:32 Acre Acre 0.4770861 0.5229139 2017-08-01 22:02:32 Amazonas Amazonas 0.4888971 0.5111029 2017-08-01 22:02:32 Roraima Roraima 0.4942450 0.5057550 2017-08-01 22:02:32 Para Para 0.4874239 0.5125761 2017-08-01 22:02:32 Amapa Amapa 0.4769455 0.5230545 2017-08-01 22:02:32 Tocantins Tocantins 0.4833958 0.5166042 2017-08-01 22:02:32 Maranhao Maranhao 0.4786299 0.5213701 2017-08-01 22:02:32 Piaui Piaui 0.4798349 0.5201651 2017-08-01 22:02:32 Ceara Ceara 0.4740662 0.5259338 2017-08-01 22:02:32 Rio Grande do Norte Rio Grande do Norte 0.4648216 0.5351784 2017-08-01 22:02:32 Paraiba Paraiba 0.4659286 0.5340714 2017-08-01 22:02:32 Pernambuco Pernambuco 0.4644740 0.5355260 2017-08-01 22:02:32 Alagoas Alagoas 0.4612037 0.5387963 2017-08-01 22:02:32 Sergipe Sergipe 0.4790062 0.5209938 2017-08-01 22:02:32 Bahia Bahia 0.4658573 0.5341427 2017-08-01 22:02:32 Minas Gerais Minas Gerais 0.4730581 0.5269419 2017-08-01 22:02:32 Espirito Santo Espirito Santo 0.4769511 0.5230489 2017-08-01 22:02:32 Rio de Janeiro Rio de Janeiro 0.4521380 0.5478620 2017-08-01 22:02:32 Sao Paulo Sao Paulo 0.4693414 0.5306586 2017-08-01 22:02:32 Parana Parana 0.4736853 0.5263147 2017-08-01 22:02:32 Santa Catarina Santa Catarina 0.4867316 0.5132684 2017-08-01 22:02:32 Rio Grande do Sul Rio Grande do Sul 0.4682191 0.5317809 2017-08-01 22:02:32 Mato Grosso do Sul Mato Grosso do Sul 0.4750854 0.5249146 2017-08-01 22:02:32 Mato Grosso Mato Grosso 0.4918424 0.5081576 2017-08-01 22:02:32 Goias Goias 0.4815873 0.5184127 2017-08-01 22:02:32 Distrito Federal Distrito Federal 0.4522139 0.5477861 2017-08-01 22:02:32 se.c006masculino se.c006feminino 2017-08-01 22:02:32 Rondonia 3.729872e-18 5.551524e-18 2017-08-01 22:02:32 Acre 7.046073e-18 7.046073e-18 2017-08-01 22:02:32 Amazonas 5.589483e-18 6.766625e-18 2017-08-01 22:02:32 Roraima 5.764340e-18 4.505724e-18 2017-08-01 22:02:32 Para 4.643261e-18 5.051309e-18 2017-08-01 22:02:32 Amapa 3.691331e-18 3.691331e-18 2017-08-01 22:02:32 Tocantins 6.828329e-18 6.828329e-18 2017-08-01 22:02:32 Maranhao 6.678935e-18 9.279058e-18 2017-08-01 22:02:32 Piaui 9.024090e-18 8.388295e-18 2017-08-01 22:02:32 Ceara 4.649624e-18 4.649624e-18 2017-08-01 22:02:32 Rio Grande do Norte 6.709916e-18 6.314833e-18 2017-08-01 22:02:32 Paraiba 4.964085e-18 4.964085e-18 2017-08-01 22:02:32 Pernambuco 4.725737e-18 4.725737e-18 2017-08-01 22:02:32 Alagoas 4.709530e-18 5.058202e-18 2017-08-01 22:02:32 Sergipe 6.653677e-18 6.048867e-18 2017-08-01 22:02:32 Bahia 3.891155e-18 5.590533e-18 2017-08-01 22:02:32 Minas Gerais 8.666421e-18 7.800171e-18 2017-08-01 22:02:32 Espirito Santo 5.631134e-18 5.614455e-18 2017-08-01 22:02:32 Rio de Janeiro 5.751597e-18 5.123351e-18 2017-08-01 22:02:32 Sao Paulo 8.962785e-18 8.962785e-18 2017-08-01 22:02:32 Parana 4.130164e-18 4.130164e-18 2017-08-01 22:02:32 Santa Catarina 3.858517e-18 5.212865e-18 2017-08-01 22:02:32 Rio Grande do Sul 4.457178e-18 4.275378e-18 2017-08-01 22:02:32 Mato Grosso do Sul 5.668093e-18 5.668093e-18 2017-08-01 22:02:32 Mato Grosso 1.074540e-17 9.280739e-18 2017-08-01 22:02:32 Goias 3.331133e-18 3.795933e-18 2017-08-01 22:02:32 Distrito Federal 1.003518e-17 9.038829e-18 2017-08-01 22:02:32 2017-08-01 22:02:32 > svytotal(~w00101, pns_design, na.rm = TRUE) 2017-08-01 22:02:34 total SE 2017-08-01 22:02:34 w00101 9.906e+09 25570730 2017-08-01 22:02:34 2017-08-01 22:02:34 > svyby(~w00101, ~uf, pns_design, svytotal, na.rm = TRUE) 2017-08-01 22:03:22 uf w00101 se 2017-08-01 22:03:22 Rondonia Rondonia 81207966 824349.2 2017-08-01 22:03:22 Acre Acre 30954162 319900.5 2017-08-01 22:03:22 Amazonas Amazonas 152089547 1187636.0 2017-08-01 22:03:22 Roraima Roraima 18243829 173749.3 2017-08-01 22:03:22 Para Para 334165128 3572022.5 2017-08-01 22:03:22 Amapa Amapa 29853686 333018.1 2017-08-01 22:03:22 Tocantins Tocantins 64457101 674705.4 2017-08-01 22:03:22 Maranhao Maranhao 272013397 2553807.7 2017-08-01 22:03:22 Piaui Piaui 139196851 1249866.3 2017-08-01 22:03:22 Ceara Ceara 395862742 3601098.4 2017-08-01 22:03:22 Rio Grande do Norte Rio Grande do Norte 162348422 1446338.9 2017-08-01 22:03:22 Paraiba Paraiba 180068058 1545640.8 2017-08-01 22:03:22 Pernambuco Pernambuco 429428903 3538312.2 2017-08-01 22:03:22 Alagoas Alagoas 149906741 1421022.9 2017-08-01 22:03:22 Sergipe Sergipe 104682696 766202.8 2017-08-01 22:03:22 Bahia Bahia 709685792 6563655.0 2017-08-01 22:03:22 Minas Gerais Minas Gerais 997655053 11632534.1 2017-08-01 22:03:22 Espirito Santo Espirito Santo 185559789 2244836.6 2017-08-01 22:03:22 Rio de Janeiro Rio de Janeiro 910188301 5356163.2 2017-08-01 22:03:22 Sao Paulo Sao Paulo 2293690500 17315373.8 2017-08-01 22:03:22 Parana Parana 552312293 5412468.3 2017-08-01 22:03:22 Santa Catarina Santa Catarina 354901662 4964588.6 2017-08-01 22:03:22 Rio Grande do Sul Rio Grande do Sul 614708304 4261031.2 2017-08-01 22:03:22 Mato Grosso do Sul Mato Grosso do Sul 125349738 1095731.9 2017-08-01 22:03:22 Mato Grosso Mato Grosso 160483906 1509042.9 2017-08-01 22:03:22 Goias Goias 319011025 2766915.3 2017-08-01 22:03:22 Distrito Federal Distrito Federal 138013197 1573500.4 2017-08-01 22:03:22 2017-08-01 22:03:22 > svytotal(~c006, pns_design) 2017-08-01 22:03:23 total SE 2017-08-01 22:03:23 c006masculino 68916470 0 2017-08-01 22:03:23 c006feminino 77391988 0 2017-08-01 22:03:23 2017-08-01 22:03:23 > svyby(~c006, ~uf, pns_design, svytotal) 2017-08-01 22:04:08 uf c006masculino c006feminino 2017-08-01 22:04:08 Rondonia Rondonia 583323.4 606530.8 2017-08-01 22:04:08 Acre Acre 228783.7 250760.2 2017-08-01 22:04:08 Amazonas Amazonas 1115983.6 1166671.8 2017-08-01 22:04:08 Roraima Roraima 136457.3 139635.1 2017-08-01 22:04:08 Para Para 2542704.7 2673913.7 2017-08-01 22:04:08 Amapa Amapa 216331.6 237245.6 2017-08-01 22:04:08 Tocantins Tocantins 477414.6 510212.1 2017-08-01 22:04:08 Maranhao Maranhao 2094065.4 2281058.8 2017-08-01 22:04:08 Piaui Piaui 1061938.3 1151194.6 2017-08-01 22:04:08 Ceara Ceara 2950153.1 3272929.1 2017-08-01 22:04:08 Rio Grande do Norte Rio Grande do Norte 1122857.9 1292817.2 2017-08-01 22:04:08 Paraiba Paraiba 1294277.1 1483567.0 2017-08-01 22:04:08 Pernambuco Pernambuco 3012553.8 3473393.4 2017-08-01 22:04:08 Alagoas Alagoas 1040664.8 1215745.6 2017-08-01 22:04:08 Sergipe Sergipe 733835.5 798160.3 2017-08-01 22:04:08 Bahia Bahia 4969934.8 5698427.6 2017-08-01 22:04:08 Minas Gerais Minas Gerais 7257722.9 8084416.9 2017-08-01 22:04:08 Espirito Santo Espirito Santo 1333406.4 1462281.6 2017-08-01 22:04:08 Rio de Janeiro Rio de Janeiro 5763689.7 6983944.0 2017-08-01 22:04:08 Sao Paulo Sao Paulo 15577076.9 17612143.8 2017-08-01 22:04:08 Parana Parana 3808124.4 4231231.2 2017-08-01 22:04:08 Santa Catarina Santa Catarina 2458695.9 2592744.7 2017-08-01 22:04:08 Rio Grande do Sul Rio Grande do Sul 3995719.8 4538147.9 2017-08-01 22:04:08 Mato Grosso do Sul Mato Grosso do Sul 841190.8 929418.9 2017-08-01 22:04:08 Mato Grosso Mato Grosso 1116841.5 1153888.9 2017-08-01 22:04:08 Goias Goias 2252885.1 2425156.2 2017-08-01 22:04:08 Distrito Federal Distrito Federal 929836.8 1126351.0 2017-08-01 22:04:08 se.c006masculino se.c006feminino 2017-08-01 22:04:08 Rondonia 1.055509e-11 8.162454e-12 2017-08-01 22:04:08 Acre 2.778211e-12 1.809925e-12 2017-08-01 22:04:08 Amazonas 1.414805e-11 1.316256e-11 2017-08-01 22:04:08 Roraima 1.981371e-12 2.011455e-12 2017-08-01 22:04:08 Para 2.775186e-11 2.351348e-11 2017-08-01 22:04:08 Amapa 1.868308e-12 3.917436e-12 2017-08-01 22:04:08 Tocantins 9.297487e-12 9.307504e-12 2017-08-01 22:04:08 Maranhao 3.346834e-11 5.114452e-11 2017-08-01 22:04:08 Piaui 2.157676e-11 1.136971e-11 2017-08-01 22:04:08 Ceara 3.061303e-11 3.168642e-11 2017-08-01 22:04:08 Rio Grande do Norte 1.056384e-11 1.740655e-11 2017-08-01 22:04:08 Paraiba 1.683168e-11 3.120145e-11 2017-08-01 22:04:08 Pernambuco 3.920445e-11 3.111754e-11 2017-08-01 22:04:08 Alagoas 1.341281e-11 1.893330e-11 2017-08-01 22:04:08 Sergipe 8.512099e-12 1.272918e-11 2017-08-01 22:04:08 Bahia 1.603945e-10 7.901851e-11 2017-08-01 22:04:08 Minas Gerais 9.357109e-11 1.675816e-10 2017-08-01 22:04:08 Espirito Santo 3.566086e-11 3.271353e-11 2017-08-01 22:04:08 Rio de Janeiro 6.024668e-11 6.726535e-11 2017-08-01 22:04:08 Sao Paulo 3.989328e-10 1.837736e-10 2017-08-01 22:04:08 Parana 4.797330e-11 4.871226e-11 2017-08-01 22:04:08 Santa Catarina 3.861720e-11 8.879706e-11 2017-08-01 22:04:08 Rio Grande do Sul 5.272979e-11 3.273625e-11 2017-08-01 22:04:08 Mato Grosso do Sul 2.008009e-11 2.759123e-11 2017-08-01 22:04:08 Mato Grosso 1.899481e-11 1.655943e-11 2017-08-01 22:04:08 Goias 2.954581e-11 4.184105e-11 2017-08-01 22:04:08 Distrito Federal 7.654015e-12 2.162701e-11 2017-08-01 22:04:08 2017-08-01 22:04:08 > svyquantile(~w00101, pns_design, 0.5, na.rm = TRUE) 2017-08-01 22:04:08 0.5 2017-08-01 22:04:08 w00101 69.7 2017-08-01 22:04:08 2017-08-01 22:04:08 > svyby(~w00101, ~uf, pns_design, svyquantile, 0.5, 2017-08-01 22:04:08 + ci = TRUE, keep.var = TRUE, na.rm = TRUE) 2017-08-01 22:05:03 uf w00101 se 2017-08-01 22:05:03 Rondonia Rondonia 68.90000 0.5612348 2017-08-01 22:05:03 Acre Acre 66.00000 0.4847028 2017-08-01 22:05:03 Amazonas Amazonas 67.30000 0.4081708 2017-08-01 22:05:03 Roraima Roraima 67.70000 0.6291724 2017-08-01 22:05:03 Para Para 66.00000 0.5612348 2017-08-01 22:05:03 Amapa Amapa 67.01431 0.5102135 2017-08-01 22:05:03 Tocantins Tocantins 67.00000 0.6632775 2017-08-01 22:05:03 Maranhao Maranhao 63.00000 0.7468389 2017-08-01 22:05:03 Piaui Piaui 64.10000 0.5357241 2017-08-01 22:05:03 Ceara Ceara 66.10000 0.4462180 2017-08-01 22:05:03 Rio Grande do Norte Rio Grande do Norte 68.60000 0.5357241 2017-08-01 22:05:03 Paraiba Paraiba 66.80000 0.4870273 2017-08-01 22:05:03 Pernambuco Pernambuco 68.30000 0.4591921 2017-08-01 22:05:03 Alagoas Alagoas 67.90000 0.5314039 2017-08-01 22:05:03 Sergipe Sergipe 68.62337 0.4166224 2017-08-01 22:05:03 Bahia Bahia 66.90000 0.7398095 2017-08-01 22:05:03 Minas Gerais Minas Gerais 69.00000 0.5867455 2017-08-01 22:05:03 Espirito Santo Espirito Santo 69.70000 0.6377668 2017-08-01 22:05:03 Rio de Janeiro Rio de Janeiro 72.19962 0.2551067 2017-08-01 22:05:03 Sao Paulo Sao Paulo 71.80000 0.4336814 2017-08-01 22:05:03 Parana Parana 70.90000 0.5102135 2017-08-01 22:05:03 Santa Catarina Santa Catarina 74.00000 0.7908309 2017-08-01 22:05:03 Rio Grande do Sul Rio Grande do Sul 73.10000 0.5612348 2017-08-01 22:05:03 Mato Grosso do Sul Mato Grosso do Sul 71.14348 0.4336814 2017-08-01 22:05:03 Mato Grosso Mato Grosso 71.70000 0.6803694 2017-08-01 22:05:03 Goias Goias 69.00000 0.5102135 2017-08-01 22:05:03 Distrito Federal Distrito Federal 70.10224 0.4404234 2017-08-01 22:05:03 2017-08-01 22:05:03 > svyratio(numerator = ~w00203, denominator = ~w00101, 2017-08-01 22:05:03 + pns_design, na.rm = TRUE) 2017-08-01 22:05:07 Ratio estimator: svyratio.survey.design2(numerator = ~w00203, denominator = ~w00101, 2017-08-01 22:05:07 pns_design, na.rm = TRUE) 2017-08-01 22:05:07 Ratios= 2017-08-01 22:05:07 w00101 2017-08-01 22:05:07 w00203 2.298467 2017-08-01 22:05:07 SEs= 2017-08-01 22:05:07 w00101 2017-08-01 22:05:07 w00203 0.003472622 2017-08-01 22:05:07 2017-08-01 22:05:07 > sub_pns_design <- subset(pns_design, atfi04 == 1) 2017-08-01 22:05:07 2017-08-01 22:05:07 > svymean(~w00101, sub_pns_design, na.rm = TRUE) 2017-08-01 22:05:09 mean SE 2017-08-01 22:05:09 w00101 70.147 0.1826 2017-08-01 22:05:09 2017-08-01 22:05:09 > this_result <- svymean(~w00101, pns_design, na.rm = TRUE) 2017-08-01 22:05:10 2017-08-01 22:05:10 > coef(this_result) 2017-08-01 22:05:10 w00101 2017-08-01 22:05:10 71.4502 2017-08-01 22:05:10 2017-08-01 22:05:10 > SE(this_result) 2017-08-01 22:05:10 w00101 2017-08-01 22:05:10 w00101 0.1151966 2017-08-01 22:05:10 2017-08-01 22:05:10 > confint(this_result) 2017-08-01 22:05:10 2.5 % 97.5 % 2017-08-01 22:05:10 w00101 71.22442 71.67598 2017-08-01 22:05:10 2017-08-01 22:05:10 > cv(this_result) 2017-08-01 22:05:10 w00101 2017-08-01 22:05:10 w00101 0.001612264 2017-08-01 22:05:10 2017-08-01 22:05:10 > grouped_result <- svyby(~w00101, ~uf, pns_design, 2017-08-01 22:05:10 + svymean, na.rm = TRUE) 2017-08-01 22:05:59 2017-08-01 22:05:59 > coef(grouped_result) 2017-08-01 22:05:59 Rondonia Acre Amazonas Roraima 2017-08-01 22:05:59 70.57676 68.39216 69.56325 69.30896 2017-08-01 22:05:59 Para Amapa Tocantins Maranhao 2017-08-01 22:05:59 67.16428 69.06946 68.09428 64.33809 2017-08-01 22:05:59 Piaui Ceara Rio Grande do Norte Paraiba 2017-08-01 22:05:59 65.32673 67.38869 70.20013 68.28660 2017-08-01 22:05:59 Pernambuco Alagoas Sergipe Bahia 2017-08-01 22:05:59 70.33818 69.79392 70.47851 68.93628 2017-08-01 22:05:59 Minas Gerais Espirito Santo Rio de Janeiro Sao Paulo 2017-08-01 22:05:59 70.58504 70.95230 73.91826 73.65539 2017-08-01 22:05:59 Parana Santa Catarina Rio Grande do Sul Mato Grosso do Sul 2017-08-01 22:05:59 72.54831 74.85520 75.01492 73.63023 2017-08-01 22:05:59 Mato Grosso Goias Distrito Federal 2017-08-01 22:05:59 73.36957 71.58753 71.82446 2017-08-01 22:05:59 2017-08-01 22:05:59 > SE(grouped_result) 2017-08-01 22:05:59 [1] 0.4500115 0.4139088 0.4035524 0.4834325 0.4063297 0.4444721 0.5208864 2017-08-01 22:05:59 [8] 0.4498913 0.4205205 0.3611066 0.4876779 0.4297494 0.3672367 0.4233045 2017-08-01 22:05:59 [15] 0.4411245 0.5974439 0.4722967 0.4223929 0.3141405 0.3192958 0.4019706 2017-08-01 22:05:59 [22] 0.5407879 0.4049787 0.4439434 0.5221943 0.4368148 0.3726693 2017-08-01 22:05:59 2017-08-01 22:05:59 > confint(grouped_result) 2017-08-01 22:05:59 2.5 % 97.5 % 2017-08-01 22:05:59 Rondonia 69.69475 71.45877 2017-08-01 22:05:59 Acre 67.58091 69.20341 2017-08-01 22:05:59 Amazonas 68.77230 70.35420 2017-08-01 22:05:59 Roraima 68.36145 70.25647 2017-08-01 22:05:59 Para 66.36789 67.96067 2017-08-01 22:05:59 Amapa 68.19831 69.94061 2017-08-01 22:05:59 Tocantins 67.07336 69.11520 2017-08-01 22:05:59 Maranhao 63.45632 65.21986 2017-08-01 22:05:59 Piaui 64.50253 66.15094 2017-08-01 22:05:59 Ceara 66.68094 68.09645 2017-08-01 22:05:59 Rio Grande do Norte 69.24430 71.15596 2017-08-01 22:05:59 Paraiba 67.44431 69.12889 2017-08-01 22:05:59 Pernambuco 69.61841 71.05795 2017-08-01 22:05:59 Alagoas 68.96426 70.62358 2017-08-01 22:05:59 Sergipe 69.61392 71.34310 2017-08-01 22:05:59 Bahia 67.76531 70.10724 2017-08-01 22:05:59 Minas Gerais 69.65935 71.51072 2017-08-01 22:05:59 Espirito Santo 70.12442 71.78017 2017-08-01 22:05:59 Rio de Janeiro 73.30256 74.53397 2017-08-01 22:05:59 Sao Paulo 73.02958 74.28120 2017-08-01 22:05:59 Parana 71.76046 73.33616 2017-08-01 22:05:59 Santa Catarina 73.79527 75.91512 2017-08-01 22:05:59 Rio Grande do Sul 74.22117 75.80866 2017-08-01 22:05:59 Mato Grosso do Sul 72.76012 74.50035 2017-08-01 22:05:59 Mato Grosso 72.34608 74.39305 2017-08-01 22:05:59 Goias 70.73139 72.44367 2017-08-01 22:05:59 Distrito Federal 71.09404 72.55487 2017-08-01 22:05:59 2017-08-01 22:05:59 > cv(grouped_result) 2017-08-01 22:05:59 Rondonia Acre Amazonas Roraima 2017-08-01 22:05:59 0.006376200 0.006051992 0.005801230 0.006975036 2017-08-01 22:05:59 Para Amapa Tocantins Maranhao 2017-08-01 22:05:59 0.006049789 0.006435146 0.007649488 0.006992611 2017-08-01 22:05:59 Piaui Ceara Rio Grande do Norte Paraiba 2017-08-01 22:05:59 0.006437188 0.005358564 0.006946966 0.006293319 2017-08-01 22:05:59 Pernambuco Alagoas Sergipe Bahia 2017-08-01 22:05:59 0.005221014 0.006065063 0.006258993 0.008666612 2017-08-01 22:05:59 Minas Gerais Espirito Santo Rio de Janeiro Sao Paulo 2017-08-01 22:05:59 0.006691174 0.005953195 0.004249836 0.004334996 2017-08-01 22:05:59 Parana Santa Catarina Rio Grande do Sul Mato Grosso do Sul 2017-08-01 22:05:59 0.005540730 0.007224453 0.005398642 0.006029363 2017-08-01 22:05:59 Mato Grosso Goias Distrito Federal 2017-08-01 22:05:59 0.007117315 0.006101828 0.005188613 2017-08-01 22:05:59 2017-08-01 22:05:59 > degf(pns_design) 2017-08-01 22:05:59 [1] 5480 2017-08-01 22:05:59 2017-08-01 22:05:59 > svyvar(~w00101, pns_design, na.rm = TRUE) 2017-08-01 22:06:03 variance SE 2017-08-01 22:06:03 w00101 243.89 3.2704 2017-08-01 22:06:03 2017-08-01 22:06:03 > svymean(~w00101, pns_design, na.rm = TRUE, deff = TRUE) 2017-08-01 22:06:09 mean SE DEff 2017-08-01 22:06:09 w00101 71.4502 0.1152 3.1181 2017-08-01 22:06:09 2017-08-01 22:06:09 > svymean(~w00101, pns_design, na.rm = TRUE, deff = "replace") 2017-08-01 22:06:14 mean SE DEff 2017-08-01 22:06:14 w00101 71.4502 0.1152 3.1168 2017-08-01 22:06:14 2017-08-01 22:06:14 > svyciprop(~health_insurance, pns_design, method = "likelihood", 2017-08-01 22:06:14 + na.rm = TRUE) 2017-08-01 22:06:21 2.5% 97.5% 2017-08-01 22:06:21 health_insurance 0.303 0.294 0.31 2017-08-01 22:06:21 2017-08-01 22:06:21 > svyttest(w00101 ~ health_insurance, pns_design) 2017-08-01 22:06:23 2017-08-01 22:06:23 Design-based t-test 2017-08-01 22:06:23 2017-08-01 22:06:23 data: w00101 ~ health_insurance 2017-08-01 22:06:23 t = 14.395, df = 5479, p-value < 2.2e-16 2017-08-01 22:06:23 alternative hypothesis: true difference in mean is not equal to 0 2017-08-01 22:06:23 sample estimates: 2017-08-01 22:06:23 difference in mean 2017-08-01 22:06:23 3.668906 2017-08-01 22:06:23 2017-08-01 22:06:23 2017-08-01 22:06:23 > svychisq(~health_insurance + c006, pns_design) 2017-08-01 22:06:25 2017-08-01 22:06:25 Pearson's X^2: Rao & Scott adjustment 2017-08-01 22:06:25 2017-08-01 22:06:25 data: svychisq(~health_insurance + c006, pns_design) 2017-08-01 22:06:25 F = 6.2074, ndf = 1, ddf = 5480, p-value = 0.01275 2017-08-01 22:06:25 2017-08-01 22:06:25 2017-08-01 22:06:25 > glm_result <- svyglm(w00101 ~ health_insurance + c006, 2017-08-01 22:06:25 + pns_design) 2017-08-01 22:06:27 2017-08-01 22:06:27 > summary(glm_result) 2017-08-01 22:06:30 2017-08-01 22:06:30 Call: 2017-08-01 22:06:30 svyglm(formula = w00101 ~ health_insurance + c006, pns_design) 2017-08-01 22:06:30 2017-08-01 22:06:30 Survey design: 2017-08-01 22:06:30 update(pns_design, one = 1, health_insurance = as.numeric(i001 == 2017-08-01 22:06:30 1)) 2017-08-01 22:06:30 2017-08-01 22:06:30 Coefficients: 2017-08-01 22:06:30 Estimate Std. Error t value Pr(>|t|) 2017-08-01 22:06:30 (Intercept) 75.1053 0.1857 404.44 <2e-16 *** 2017-08-01 22:06:30 health_insurance 3.8433 0.2430 15.81 <2e-16 *** 2017-08-01 22:06:30 c006feminino -9.1428 0.2226 -41.08 <2e-16 *** 2017-08-01 22:06:30 --- 2017-08-01 22:06:30 Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 2017-08-01 22:06:30 2017-08-01 22:06:30 (Dispersion parameter for gaussian family taken to be 208.6816) 2017-08-01 22:06:30 2017-08-01 22:06:30 Number of Fisher Scoring iterations: 2 2017-08-01 22:06:30 2017-08-01 22:06:30 2017-08-01 22:06:30 > library(srvyr) 2017-08-01 22:06:30 2017-08-01 22:06:30 Attaching package: 'srvyr' 2017-08-01 22:06:30 2017-08-01 22:06:30 The following object is masked from 'package:stats': 2017-08-01 22:06:30 2017-08-01 22:06:30 filter 2017-08-01 22:06:30 2017-08-01 22:06:30 2017-08-01 22:06:30 > pns_srvyr_design <- as_survey(pns_design) 2017-08-01 22:06:30 2017-08-01 22:06:30 > pns_srvyr_design %>% summarize(mean = survey_mean(w00101, 2017-08-01 22:06:30 + na.rm = TRUE)) 2017-08-01 22:06:32 # A tibble: 1 x 2 2017-08-01 22:06:32 mean mean_se 2017-08-01 22:06:32 2017-08-01 22:06:32 1 71.4502 0.1151966 2017-08-01 22:06:32 2017-08-01 22:06:32 > pns_srvyr_design %>% group_by(uf) %>% summarize(mean = survey_mean(w00101, 2017-08-01 22:06:32 + na.rm = TRUE)) 2017-08-01 22:07:18 # A tibble: 27 x 3 2017-08-01 22:07:18 uf mean mean_se 2017-08-01 22:07:18 2017-08-01 22:07:18 1 Rondonia 70.57676 0.4500115 2017-08-01 22:07:18 2 Acre 68.39216 0.4139088 2017-08-01 22:07:18 3 Amazonas 69.56325 0.4035524 2017-08-01 22:07:18 4 Roraima 69.30896 0.4834325 2017-08-01 22:07:18 5 Para 67.16428 0.4063297 2017-08-01 22:07:18 6 Amapa 69.06946 0.4444721 2017-08-01 22:07:18 7 Tocantins 68.09428 0.5208864 2017-08-01 22:07:18 8 Maranhao 64.33809 0.4498913 2017-08-01 22:07:18 9 Piaui 65.32673 0.4205205 2017-08-01 22:07:18 10 Ceara 67.38869 0.3611066 2017-08-01 22:07:18 # ... with 17 more rows 2017-08-01 22:07:18 There were 27 warnings (use warnings() to see them) 2017-08-01 22:07:18 >