R/runsim_logreg.R

In LindaNab/simexvsmecor: Comparison of Simex and Mecor for Measurement Error Correction

#' Run simulation study PART II (logreg)
#'
#' @param nrep number of replications
#' @param scen_no scenarios numbers to be used
#' @param output_file_loc directory where output of simulation will be saved
#' @return a file named scen_no_* with the output will be saved in
#' output_file_loc, which constitutes of the uncorrected effect, the corrected
#' effect by means of mecor and the corrected effect by means of simex, the
#' estimated R-squared of the outcome model and the seed no.
#' @export
run_sim_logreg <- function(nrep,
                           scen_no,
                           output_file_loc = "./") {
  output_file <- paste0(output_file_loc, "logreg_scen_no", scen_no, ".Rds")
  for (i in 1:nrep) {
    seed_no <- sample(1:1e6, 1)
    cat("\f")
    print(paste0("Scen #", scen_no, " Run #", i, " with seed: ", seed_no))
    data <- generate_data_logreg_scen_no(seed_no,
                                         scen_no)
    effects <- get_est_effects_logreg(data)
    r_squared <- get_r_squared_logreg(data)
    crude_effect <- get_crude_logreg(data)
    output <- c(effects,
                r_squared = r_squared,
                crude_effect = crude_effect,
                seed_no = seed_no)
    save_output(output,
                output_file)
  }
}