Balanced randomization for EEG O.4

Randomization O.4 - 20 subjects & 32 subjects

1 Randomization Table for 20 Subjects

# Repetate
set.seed(101)
eeg_rand <- psych::block.random(20, c(Condition1 = 2, Order1 = 2)) 
table_eeg_rand <- 
  eeg_rand %>% 
  as.data.frame(.) %>%  
  mutate(Condition1 = case_when(Condition1 == 1 ~ "ogl" ,
                                Condition1 == 2 ~ "ctrl", 
                                TRUE ~ NA_character_)) %>%
  mutate(Order1 = case_when(Order1 == 1 ~ "Instr-Solo" ,
                            Order1 == 2 ~ "Solo-Instr", 
                            TRUE ~ NA_character_)) %>%
  mutate(Condition2 = if_else(Condition1 == "ogl", "ctrl", "ogl")) %>%
  mutate(Order2 = if_else(Order1 == "Instr-Solo", "Solo-Instr", "Instr-Solo")) %>%
  mutate(ID = 1:nrow(.)) %>%
  mutate(Set = rep(NA, nrow(.)),
         Instructor = rep(NA, nrow(.)))
  
table_eeg_rand %>%                              # excel downloadable  DT table
  select(blocks, ID, everything()) %>%
    DT::datatable(
      extensions = 'Buttons',
      options = list(pageLength = 40,
                     dom = 'Bfrtip', 
                     buttons = c('excel', "csv")))

# kable table
# table_eeg_rand %>%                             
#   select(blocks, ID, everything()) %>%
#   knitr::kable(caption = "Randomization Table", align = rep('c', 6))

1.1 Check cross-tabs for each day

table_eeg_rand %>%
  select(Condition1, Order1) %>%
  ftable(row.vars = 1)

           Order1 Instr-Solo Solo-Instr
Condition1                             
ctrl                       5          5
ogl                        5          5

  # knitr::kable()
table_eeg_rand %>%
  select(Condition2, Order2) %>%
  ftable(row.vars = 1)

           Order2 Instr-Solo Solo-Instr
Condition2                             
ctrl                       5          5
ogl                        5          5

# all conditions are balanced by block
# table_eeg_rand %>%
#   filter(blocks == 1) %>%
#   select(Condition1, Instructor1, Order1) %>%
#   ftable(row.vars = 1)

1.2 Check balance

# pairs.panels(eeg_rand)
table_eeg_rand %>%
  select(-ID, -blocks, -Instructor, -Set) %>%
    pairs.panels()

2 Randomization Table for 30 Subjects

# Repetate
set.seed(101)
eeg_rand32 <- psych::block.random(32, c(Condition1 = 2, Order1 = 2, Set1 = 4)) 
table_eeg_rand32 <- 
  eeg_rand32 %>% 
  as.data.frame(.) %>%  
  mutate(Condition1 = case_when(Condition1 == 1 ~ "ogl" ,
                                Condition1 == 2 ~ "ctrl", 
                                TRUE ~ NA_character_)) %>%
  mutate(Order1 = case_when(Order1 == 1 ~ "Instr-Solo" ,
                            Order1 == 2 ~ "Solo-Instr", 
                            TRUE ~ NA_character_)) %>%
  mutate(Set1 = case_when(Set1 == 1 ~ "1-2" ,
                          Set1 == 2 ~ "2-1",
                          Set1 == 3 ~ "1-2f" ,
                          Set1 == 4 ~ "2-1f",
                          TRUE ~ NA_character_)) %>%
  mutate(Condition2 = if_else(Condition1 == "ogl", "ctrl", "ogl")) %>%
  mutate(Order2 = if_else(Order1 == "Instr-Solo", "Solo-Instr", "Instr-Solo")) %>%
  mutate(Set2 = if_else(Set1 %in% c("1-2", "2-1f"), "2-1", "1-2")) %>%
  mutate(ID = 1:nrow(.)) %>%
  mutate(Instructor = rep(NA, nrow(.)))
  
table_eeg_rand32 %>%                              # excel downloadable  DT table
  select(blocks, ID, everything()) %>%
    DT::datatable(
      extensions = 'Buttons',
      options = list(pageLength = 40,
                     dom = 'Bfrtip', 
                     buttons = c('excel', "csv")))

2.1 Check cross-tabs for each day

table_eeg_rand32 %>%
  select(Condition1, Order1, Set1) %>%
  ftable(row.vars = 1)

           Order1 Instr-Solo               Solo-Instr              
           Set1          1-2 1-2f 2-1 2-1f        1-2 1-2f 2-1 2-1f
Condition1                                                         
ctrl                       2    2   2    2          2    2   2    2
ogl                        2    2   2    2          2    2   2    2

  # knitr::kable()
table_eeg_rand32 %>%
  select(Condition2, Order2, Set1) %>%
  ftable(row.vars = 1)

           Order2 Instr-Solo               Solo-Instr              
           Set1          1-2 1-2f 2-1 2-1f        1-2 1-2f 2-1 2-1f
Condition2                                                         
ctrl                       2    2   2    2          2    2   2    2
ogl                        2    2   2    2          2    2   2    2

2.2 Check balance

# pairs.panels(eeg_rand)
table_eeg_rand32 %>%
  select(-ID, -blocks, -Instructor) %>%
    pairs.panels()

3 Session Info

R version 3.5.2 (2018-12-20)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)

Matrix products: default

locale:
[1] LC_COLLATE=Romanian_Romania.1250  LC_CTYPE=Romanian_Romania.1250    LC_MONETARY=Romanian_Romania.1250 LC_NUMERIC=C                     
[5] LC_TIME=Romanian_Romania.1250    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] bindrcpp_0.2.2    DT_0.5            papaja_0.1.0.9842 psych_1.8.10      forcats_0.3.0     stringr_1.3.1     dplyr_0.7.8      
 [8] purrr_0.2.5       readr_1.3.0       tidyr_0.8.2       tibble_1.4.2      ggplot2_3.1.0     tidyverse_1.2.1   pacman_0.5.0     

loaded via a namespace (and not attached):
 [1] tidyselect_0.2.5 xfun_0.4         haven_2.1.0      lattice_0.20-38  colorspace_1.3-2 generics_0.0.2   htmltools_0.3.6  yaml_2.2.0      
 [9] rlang_0.3.0.1    later_0.7.5      pillar_1.3.1     foreign_0.8-71   glue_1.3.1       withr_2.1.2      modelr_0.1.2     readxl_1.1.0    
[17] bindr_0.1.1      plyr_1.8.4       munsell_0.5.0    gtable_0.2.0     cellranger_1.1.0 rvest_0.3.2      htmlwidgets_1.3  knitr_1.21      
[25] httpuv_1.4.5     crosstalk_1.0.0  parallel_3.5.2   broom_0.5.1      Rcpp_1.0.1       xtable_1.8-3     promises_1.0.1   scales_1.0.0    
[33] backports_1.1.3  jsonlite_1.6     mime_0.6         mnormt_1.5-5     digest_0.6.18    hms_0.4.2        stringi_1.2.4    shiny_1.2.0     
[41] grid_3.5.2       cli_1.0.1        tools_3.5.2      magrittr_1.5     lazyeval_0.2.1   crayon_1.3.4     pkgconfig_2.0.2  xml2_1.2.0      
[49] lubridate_1.7.4  assertthat_0.2.1 httr_1.4.0       rstudioapi_0.8   R6_2.4.0         nlme_3.1-137     compiler_3.5.2  

 

---

A work by Claudiu Papasteri

claudiu.papasteri@gmail.com