##To create the graphical layout of the D-optimal row-column design
##obtained using the treatment exchange algorithm for
trt.N <- 10 #Number of treatments
col.N <- 10 #Number of arrays
theta <- 0.2 #theta value
OptdesF <- rbind(1:10, c(2:10, 1)) #D-optimal design (loop design)
graphoptrcd.mae(trt.N = 10, col.N = 10, theta = 0.2, OptdesF, Optcrit = "D", cbVal2 = 0)
Run the code above in your browser using DataLab