# NOT RUN {
##To create the graphical layout of the D-optimal block design
##obtained using the treatment exchange algorithm for
trt.N <- 10 #Number of treatments
blk.N <- 10 #Number of blocks
alpha <- 0.1 #alpha value
beta <- 0.1 #beta value
OptdesF <- rbind(1:10, c(2:10,1)) #Bayesian A-optimal block design (loop design)
graphoptBbd(trt.N = 10, blk.N = 10, alpha = 0.1, beta = 0.1, OptdesF, Optcrit = "A")
# }
Run the code above in your browser using DataCamp Workspace