17  Latin square designs

๐Ÿ“– Status of the book

Hi there! This book is a work-in-progress. You may like to come back later when itโ€™s closer to a complete state. If you would like to raise issues or leave feedback, please feel to do this:

lsd <- takeout(menu_lsd(t = 6))
examine_recipe(lsd)
design("Latin Square Design") %>%
  set_units(row = 6,
            col = 6,
            unit = crossed_by(row, col)) %>%
  set_trts(trt = 6) %>%
  allot_trts(trt ~ unit) %>%
  assign_trts("random", seed = 836) %>%
  serve_table()
autoplot(lsd)

17.1 Graeco-Latin square design

graeco <- takeout(menu_graeco(t = 9))
examine_recipe(graeco)
design("Graeco-Latin Square Design") %>%
  set_units(row = 9,
            col = 9,
            unit = crossed_by(row, col)) %>%
  set_trts(trt1 = 9,
           trt2 = 9) %>%
  allot_trts(trt1 ~ unit,
             trt2 ~ unit) %>%
  assign_trts("random", seed = 528) %>%
  serve_table()
autoplot(graeco)

17.2 Hyper-Graeco-Latin square design

hyper_graeco <- takeout(menu_hyper_graeco(t = 6))
examine_recipe(hyper_graeco)
design("Hyper-Graeco-Latin Square Design") %>%
  set_units(block1 = 6,
            block2 = 6,
            block3 = 6,
            block4 = 6,
            unit = crossed_by(block1, block2, block3, block4)) %>%
  set_trts(trt = 6) %>%
  allot_trts(trt ~ unit) %>%
  assign_trts("random", seed = 109) %>%
  serve_table()
autoplot(hyper_graeco)
Warning: Too manu units so 296 units dropped from the plot. If you want to see
all, use `nnode_max = Inf` or use `page = 2` to see the next set.

17.3 Youden square design

youden <- takeout(menu_youden(nc = 4, t = 5))
examine_recipe(youden)
design("Youden Square Design") %>%
  set_units(row = 5,
            col = 4,
            unit = crossed_by(row, col)) %>%
  set_trts(trt = 5) %>%
  allot_trts(trt ~ unit) %>%
  assign_trts("random", seed = 126) %>%
  serve_table()
autoplot(youden)