data-raw/t_res.R

In meghapsimatrix/SimHelpers: Helper Functions for Simulation Studies

## code to prepare `DATASET` dataset goes here

library(tidyverse)
library(broom)
library(kableExtra)

set.seed(20191228)

# function to create normally distributed data for each group to run t test

generate_dat <- function(n = 50, effect_x){

  dat <- tibble(group_1 = rnorm(n, effect_x, 1),
                group_2 = rnorm(n, 0, 1))

  return(dat)

}


# function to calculate t-test, pulls out estimate of the mean difference, p val and ci

estimate_t <- function(sim_dat, effect_x){

  res <- tidy(t.test(sim_dat$group_1, sim_dat$group_2)) %>%
    select(est = estimate, p_val = p.value, lower_bound = conf.low, upper_bound = conf.high) %>%
    mutate(true_param = effect_x)

  return(res)

}


# generating 1000 iterations

t_res <-
  map_dfr(1:1000, ~ {
    dat <- generate_dat(effect_x = .5)
    estimate_t(dat, effect_x = .5)
  })

usethis::use_data(t_res, overwrite = TRUE)