Is this what you were after:
> ## 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)
> ans <- rcbd(b=4,g=1:4,rb=2,cb=2,r=4,c=4)
>
> result <- matrix(nrow = max(ans$x), ncol = max(ans$y))
> result[cbind(ans$y, ans$x)] <- ans$genotypes
> result
[,1] [,2] [,3] [,4]
[1,] 2 4 2 3
[2,] 1 3 4 1
[3,] 3 2 2 3
[4,] 1 4 4 1
>
Jim Holtman
Data Munger Guru
What is the problem that you are trying to solve?
Tell me what you want to do, not how you want to do it.
On Sun, Sep 22, 2013 at 1:07 AM, Laz <lmra...@ufl.edu> wrote:
> 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.
______________________________________________
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.
