# Group-divisible NPBIB on v = m*n = 8 treatments (m = 4 rows, n = 2),
# with q = 2 sub-blocks per block; merge the first two rows into v2 = 2 controls
gd <- list(c(1,3), c(2,4), c(5,7), c(6,8), c(1,5), c(2,6),
c(3,7), c(4,8), c(1,7), c(2,8), c(3,5), c(4,6))
d <- construct_method3(gd, m = 4, n = 2, q = 2, merge_rows = c(1, 2), v2 = 2)
d
Run the code above in your browser using DataLab