Revisit algorithm for assign_trts() #46
Description
The current algorithm doesn't take into account a hierarchical structure with more than one level. Therefore the randomisation for something like below doesn't balance out on the second level of the hierarchy (i.e. rep) but it is on the first top level (i.e. block). A more general approach needed to account for this.
library(edibble)
eg <- design("Balanced Incomplete Block Design") %>%
set_units(rep = 3,
block = nested_in(rep, 12),
unit = nested_in(block, 4)) %>%
set_trts(trt = 48) %>%
allot_trts(trt ~ unit) %>%
assign_trts("random", seed = 32) %>%
serve_table()
with(eg, table(trt,rep))
#> rep
#> trt rep1 rep2 rep3
#> trt1 3 0 0
#> trt10 0 2 1
#> trt11 0 1 2
#> trt12 2 1 0
#> trt13 0 1 2
#> trt14 1 1 1
#> trt15 2 0 1
#> trt16 0 1 2
#> trt17 1 1 1
#> trt18 1 1 1
#> trt19 1 1 1
#> trt2 2 1 0
#> trt20 1 0 2
#> trt21 0 2 1
#> trt22 0 2 1
#> trt23 2 1 0
#> trt24 1 1 1
#> trt25 1 2 0
#> trt26 0 0 3
#> trt27 1 2 0
#> trt28 1 2 0
#> trt29 0 2 1
#> trt3 1 2 0
#> trt30 1 0 2
#> trt31 2 1 0
#> trt32 0 1 2
#> trt33 2 1 0
#> trt34 1 1 1
#> trt35 2 0 1
#> trt36 1 2 0
#> trt37 0 1 2
#> trt38 0 1 2
#> trt39 2 1 0
#> trt4 2 0 1
#> trt40 0 0 3
#> trt41 1 1 1
#> trt42 0 1 2
#> trt43 1 1 1
#> trt44 1 1 1
#> trt45 1 2 0
#> trt46 1 1 1
#> trt47 2 0 1
#> trt48 1 0 2
#> trt5 1 2 0
#> trt6 1 1 1
#> trt7 1 1 1
#> trt8 2 1 0
#> trt9 1 0 2
apply(with(eg, table(trt, block)), 2, table)
#> block
#> block1 block10 block11 block12 block13 block14 block15 block16 block17
#> 0 44 44 44 44 44 44 44 44 44
#> 1 4 4 4 4 4 4 4 4 4
#> block
#> block18 block19 block2 block20 block21 block22 block23 block24 block25
#> 0 44 44 44 44 44 44 44 44 44
#> 1 4 4 4 4 4 4 4 4 4
#> block
#> block26 block27 block28 block29 block3 block30 block31 block32 block33
#> 0 44 44 44 44 44 44 44 44 44
#> 1 4 4 4 4 4 4 4 4 4
#> block
#> block34 block35 block36 block4 block5 block6 block7 block8 block9
#> 0 44 44 44 44 44 44 44 44 44
#> 1 4 4 4 4 4 4 4 4 4Created on 2022-03-21 by the reprex package (v2.0.1)