R/a_tamano_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 nivel.confianza numeric - nivel de confianza. Default = 95
#' @param margen numeric - porcentaje esperado de la variable más importante del estudio. Default = 50
#' @param i.c numeric - intervalo de confianza. Default = 3.5
#' @param poblacion numeric - tamano de la población. Default 1000000
#' @param deff numeric - efecto de diseno. Default = 1.25
#' @param nr numeric - no respuesta. Degault = 0

#' @return un objeto con el tamaño de la muestra y los parámetros utilizados 
#' 
#' @details para obtener una muestra "normal" hay que fijar el efecto de diseño en 1
#' 
a_tamano_muestra<-function(nivel.confianza=95,margen=50,i.c=3.5,poblacion=1000000,deff=1.25,nr=0){
  z      =  qnorm(.5+nivel.confianza/200)
  m      =  margen/100
  int    =  i.c/100
  nr     =  nr/100
  tm     =  (z^2*(m)*(1-m)*deff)/(int^2*(1-nr))
  p.tm   =  tm/(1 + ((tm-1)/poblacion))
  NC     = c(paste(nivel.confianza, "%", sep = ""))
  IC     = c(paste(i.c, "%", sep = ""))
  NR     = c(paste(nr*100, "%", sep = ""))
  M      = c(paste(margen, "%", sep = ""))
  METODO = c("Calculos del tamnano de muestra")
  RES    = data.frame(poblacion, NC, IC, M, round(p.tm, 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)
}