IBGEtabula         package:IBGEPesq         R Documentation(latin1)

_E_s_t_i_m_a_o _d_e _t_o_t_a_i_s _c_o_m _C_V _p_o_r _m_l_t_i_p_l_o_s _d_o_m_n_i_o_s

_D_e_s_c_r_i_p_t_i_o_n:

     Estima totais e CVs, utilizando o pacote survey, permitindo a
     desagregao em mltiplos domnios.

_U_s_a_g_e:

     IBGEtabula(desenho, colunas, dominios, elimina=NA, dominios.vazios=T,
                digitos.cv = 1, total.geral=F, desvio.padrao=F)

_A_r_g_u_m_e_n_t_s:

 desenho: Objeto de classe 'survey.design' obtido pela funo svydesign
          do pacote 'survey'. Contm as informaes do desenho amostral
          e os dados a serem tabulados.

 colunas: Objeto do tipo vetor de characteres com os nomes das colunas
          dos dados contidos em 'desenho' para as quais se deseja
          estimar totais. As variveis indicadas nesta coluna devem ser
          do tipo inteiro com os valores 0 ou 1.  importante ressaltar
          que estas variveis no podem ser do tipo booleano.

dominios: Objeto do tipo vetor de characteres com os nomes das colunas
          dos dados contidos em 'desenho' que definem os domnios de
          estimao, em ordem de precedncia. Ver 'details'.

 elimina: Objeto do tipo lista, com um elemento para cada domnio
          definido no argumento 'dominios'. Cada elemento da lista  um
          vetor de caracteres, com itens que se deseja eliminar da
          tabela final como, por exemplo, sem declarao, ou NA, caso
          nenhum elemento do domnio deva ser eliminado.Se for
          fornecido somente um elemento este ser utilizado em todos os
          domnios.

dominios.vazios: Objeto do tipo lgico (boolean) que determina se
          cruzamentos de domnios vazios devem ser mantidos na tabela
          ou eliminados.

digitos.cv: Nmero de dgitos depois da vrgula para os valores de CV
          em porcentagem.

total.geral: Objeto do tipo lgico (boolean) que determina se deve ser
          calculada uma linha de total geral.

desvio.padrao: Objeto do tipo lgico (boolean) que determina se na
          sada haver tambm uma tabela de desvios padro.

_D_e_t_a_i_l_s:

     Os dados a serem tabulados so definidos dentro do argumento
     'desenho', que deve conter tambm a estrutura do desenho amostral
     utilizado. Este objeto, de classe 'survey.design', pode ser obtido
     utilizando-se a funo svydesign, disponvel no pacote 'survey'.

     Sero calculados totais e CVs para cada varivel definida em
     'colunas', desagregados pelos domnios definidos em 'dominios'. As
     colunas especificadas por esses dois argumentos devem existir nos
     dados contidos em 'desenho'.

     IMPORTANTE: As colunas de domnio devem ser do tipo fator. Os
     'labels' desses fatores sero usados para nomear as linhas da
     tabela final. Os 'labels' no devem conter o sinal de ponto (.)
     nos nomes.

     Itens indesejados de cada domnio devem ser especificados no
     argumento 'elimina'. Este argumento deve ser do tipo lista, com
     nmero de elementos igual ao nmero de domnios.

_V_a_l_u_e:

     Um objeto do tipo lista com dois (ou trs) elementos do tipo
     matriz: 'totais' (e 'desvio.padrao') e 'cv'. O nmero de colunas
     de cada matriz  igual ao tamanho do argumento 'colunas'. As
     linhas so formadas pelos vrios cruzamentos entre os 'domnios',
     formatadas e identadas de acordo com a hierarquia entre eles.

_A_u_t_h_o_r(_s):

     Eloane Gonalves Ramos e Giuseppe de Abreu Antonaci
     <giuseppe.antonaci@ibge.gov.br>

_S_e_e _A_l_s_o:

     Pacote survey.

_E_x_a_m_p_l_e_s:

     library(survey)
     library(IBGEPesq)

     data(ocupacao2005)
     colnames(dados) <- c("UF", "IDADE", "SEXO", "ANOS_ESTUDO", "COND_ATIVIDADE",
      "COND_OCUPACAO", "PREVIDENCIA", "SINDICATO", "ATIVIDADE", "OCUPACAO", "PESOS",
      "gr.anos", "ESP_DOM", "NUM_MOR", "PROJECAO", "INV_FRA", "STRAT", "PSU", "SUBA")

     dados <- dados[,c("UF", "SEXO", "COND_OCUPACAO", "PESOS", "gr.anos",
                    "PROJECAO", "INV_FRA", "STRAT", "PSU", "SUBA")]
     dados <- transform(dados,
             pestrato = PROJECAO,
             pesos = INV_FRA * SUBA,
             brasil = rep(1,nrow(dados)),
       norte = (trunc(UF/10) %in% 1)*1,
       nordeste = (trunc(UF/10) %in% 2)*1,
       sudeste = (trunc(UF/10) %in% 3)*1,
       sul = (trunc(UF/10) %in% 4)*1,
       centroeste = (trunc(UF/10) %in% 5)*1)
       
     dados <- transform(dados, pestrato=factor(pestrato),
                        gr.anos=factor(gr.anos,labels=c("Menos de 1 ano",
                       "1 a 3 anos", "4 a 7 anos", "8 a 10 anos",
                       "11 anos ou mais", "SD")),
                       SEXO = factor(SEXO, labels=c("Homens", "Mulheres")),
                       COND_OCUPACAO = factor(COND_OCUPACAO, labels=c("Ocupados",
                       "Desocupados", "SD"), exclude=NULL))

     levels(dados$pestrato)<-1:length(levels(dados$pestrato))
     pos.estratos <- data.frame(pestrato=levels(dados$pestrato),
             Freq=sort(unique(dados$PROJECAO)))

     options(survey.lonely.psu="average", survey.ultimate.cluster=T)
     dpnad<-svydesign(ids=~PSU, strata=~STRAT, data=dados, nest=TRUE,
             weights=~pesos)
     dpnad<-postStratify(dpnad, ~pestrato, pos.estratos) #ps - estratificao

     colunas <- c("brasil","norte","nordeste","sudeste","sul","centroeste")
     grupos <- c("COND_OCUPACAO","SEXO","gr.anos")
     elimina <- list("SD", NA, "SD")

     estimativa <- IBGEtabula(dpnad, colunas, grupos, elimina, total.geral=T)

