Red and Black : Search for the Best

Data

class_roll <- read_excel("../data/class_roll0303_deid.xlsx", 
                        range = "C2:I162")
str(class_roll)
class_roll <- class_roll[, -6]
names(class_roll) <- c("dept", "id", "name", "year", "email", "cell_no")
class(class_roll$id) <- "character"
class(class_roll$year) <- "character"
str(class_roll)
source("./RnB_2023_id.R")

Search for Best Configuration

M1 <- 1
M2 <- 1000
Xsum <- sapply(M1:M2, red_and_black)
names(Xsum) <- M1:M2

Summary

Xsum %>%
  summary %>%
  round(2) 

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    0.25    2.86    4.51    5.18    6.79   20.32

Xsum %>%
  sd %>%
  round(2)

## [1] 3.18

Xsum %>%
  `<=`(0.1) %>%
  which %>%
  `[`(Xsum, .) %>%
  round(2)

## named numeric(0)

Xmin <- names(Xsum[which(Xsum == min(Xsum))])
Xmin

## [1] "213" "340" "403"

Plot

Randomization

set.seed(Xmin)
N <- nrow(class_roll) 
class_roll$group <- 
  sample(1:N) %%
  2 %>%
  factor(levels = c(0, 1), labels = c("Red", "Black"))
red_and_black(Xmin)

## [1] 0.2471627

학번

	2012 이전	2013	2014	2015	2016	2017
Red	7	13	10	7	12	31
Black	7	13	10	6	14	30

Comments

R 코드를 설명에 포함하려면 class_roll %$%처럼 backtick(`) 으로 시작과 끝을 묶어줍니다. $는 LaTeX 으로 수식을 나타낼 때 사용하는 특수 문자이기 때문에 R 코드에서 %$% 를 설명하고 싶을 때 그대로 사용하면 수식으로 바뀔 수 있으니까 반드시 ` 으로 묶어 주세요.