R/a_error_muestra.R

In pelishk/upax_library: Librería de análisis

#' cálculo de tamaño de muestra
#'
#' Función para calcular el tamaño de una muestra
#'
#' @param n numeric - tamano de muestra
#' @param pob numeric - tamano de la población. Default 1000000
#' @param confianza numeric - nivel de confianza. Default = 95
#' @param margen numeric - porcentaje esperado de la variable más importante del estudio. Default = 50
#' @param deff numeric - efecto de diseno. Default = 1.25
#' @param nr numeric - no respuesta. Degault = 0

#' @return un objeto con el error de la muestra y los parámetros utilizados 

#' 
a_error_muestra<-function(n,pob,confianza = 95,margen = 50,deff = 1.25,nr = 0){
  z      =  qnorm(.5+confianza/200)
  m      =  margen/100
  nr     =  nr/100
  
  a <- (pob-n) * (z^2 * m * (1-m) * deff)
  b <- (1-nr) * n * (pob-1)
  error <- sqrt(a/b) * 100
  
  
  
  
  NC     = c(paste(confianza, "%", sep = ""))
  IC     = c(paste(round(error,2), "%", sep = ""))
  NR     = c(paste(nr*100, "%", sep = ""))
  M      = c(paste(margen, "%", sep = ""))
  METODO = c("Calculos del error de muestreo")
  RES    = data.frame(pob, NC, IC, M, round(n, digits = 0),deff,NR)
  pre    = structure(list(POB = "Poblacion", NC = "Nivel de confianza",
                          IC = "Intervalo de confianza (+/-)",
                          MDE = "Margen de error", TM = "Tamano de muestra",
                          DEFF = "Efecto de diseno", NR = "No respuesta esperada",
                          method = METODO), class = "power.htest")
  names(RES) = c("POB","NC", "IC", "MDE", "TM","DEFF","NR")
  print(pre)
  print(RES)
}