R/randomsample.R

In vincnardelli/crowdpostsampler: Spatial Crowdsourcing Sampling simulatior

#' Title
#'
#' @param model
#'
#' @return
#' @export
#'
#' @examples
randomsample <- function(model){
  if(!is(model,"crowd")) stop("Argument not valid")
  y_list = model$y
  nn = model$nn
  NN = model$NN
  N = sum(NN)
  n = sum(nn)
  sr_size<-n/4

    y_sample <- list()
    for(i in 1:4){
      y_sample[[i]] <- sample(y_list[[i]], sr_size)
    }
    mean <- sapply(y_sample, mean)
    Htest<- mean(mean)
    var_sr<-var(unlist(y_sample))
    var_Htest<- var_sr/(sqrt(n))

  out <- list(method = "Simple Random Sample",
              Htest = Htest,
              var_Htest = var_Htest,
              y_sample = unlist(y_sample),
              model = model
  )
  class(out) <- "crowdsample"
  return(out)
}