Dear R users, I have a function I created to generate randomized block designs given below. Once I generate the design, I would like to plot it in a rectangular form to mimic the excel style where the blocks are represented with the treatment values placed exactly where their appropriate X and Y coordinates are.
## A function to generate a RCB design rcbd<-function(b,g,rb,cb,r,c) { # b =number of blocks # g = a vector of treatments # rb = number of rows per blocks # cb =number of columns per block # r = total rows # c = total columns library(foreach) genotypes<-times(b) %do% sample(g,length(g)) block<-rep(1:b,each=length(g)) genotypes<-factor(genotypes) block<-factor(block) ### generate the base design k<-c/cb # number of blocks on the x-axis x<-rep(rep(1:r,each=cb),k) # X-coordinate l<-cb p<-r/rb m<-l+1 d<-l*b/p y<-c(rep(1:l,r),rep(m:d,r)) # Y-coordinate data.frame(x,y,block,genotypes) } set.seed(100) rcbd(b=4,g=1:4,rb=2,cb=2,r=4,c=4) x y block genotypes 1 1 1 1 2 2 1 2 1 1 3 2 1 1 4 4 2 2 1 3 5 3 1 2 2 6 3 2 2 4 7 4 1 2 3 8 4 2 2 1 9 1 3 3 3 10 1 4 3 1 11 2 3 3 2 12 2 4 3 4 13 3 3 4 2 14 3 4 4 4 15 4 3 4 3 16 4 4 4 1 How can I produce a diagram like this one below or a better one for any run of my function? 2 4 2 3 1 3 4 1 3 2 2 3 1 4 4 1 Regards, Laz [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.