R/svyEDER.R

In jcms2665/svyINEGI: Paquete para hacer las precisiones estadisticas con datos de la ENOE

#' @title Precisiones estadísticas, EDER 2011 
#' @name errEDER
#' @description Pruebas de significancia estadística de los indicadores de la Encuesta Demográfica Retrospectiva (EDER) 2011
#' @param u Nombre de la variable. 
#' @param data Base de datos.
#' @return Matriz de resultados.
#' @author Julio C.
#' @note Cualquier sugerencia es bienvenida
#' @details Para cada categoria, la funcion genera: la estimacion del indicador (Estim); error estandar (ErrStd); coeficiente de variación (CV); limite inferior (LInf(90%));limite superior (LSup(90%)); evaluacion significancia (Signif?)
#' @export

errEDER<-
  function(u, data = parent.frame()){
    
    if(!any(names(data) == u[1]))
      stop("[!] La variable ",u," no está en la base")
    
    if(!any(names(data) == "FOLIO"))
      stop("[!] Se requiere la variable FOLIO")
    
    if(!is.numeric(data[,u]))
      stop("[!] La variable ",u," debe ser numérica")
    
    mer<-merge(CAR_EST,data[,c("FOLIO",u)],by=c("FOLIO"),all.x=TRUE)    
    mer$UNO<-1
    
    options(survey.lonely.psu="adjust")
    
    mydesign <- svydesign(id = ~UPM,data = mer,weight = ~FAC,strata = ~EST_DIS,nest=TRUE)
    sy1<-round(svyby(~UNO,~mer[,u], mydesign, svytotal), digits = 0)
    cv1<-round(data.frame(cv(sy1)*100), digits = 3)
    li1<-round(data.frame(confint(sy1, level = 0.90)), digits = 0)
    
    cc <- matrix(NA,nrow=nrow(sy1),ncol=1)
    for(j in 1:nrow(sy1)){
      if(cv1[j,1]<15){cc[j,1]<-"Y"}else{cc[j,1]<-"N"}
    }
    r<-data.frame(sy1,cv1,li1,cc)
    names(r) <- c("Categ","Estim","ErrStd","CV","LInf(90%)","LSup(90%)","Signif?")
    return(r)
    
  }