agricolae (version 1.1-5)

design.rcbd: Randomized Complete Block Design

Description

It generates Randomized Complete Block Design. "Random" uses the methods of number generation in R. The seed is by set.seed(seed, kinds).

Usage

design.rcbd(trt, r, number = 1, seed = 0, kinds = "Super-Duper", first=FALSE)

Arguments

trt

Treatments

r

Replications or blocks

number

number of first plot

seed

seed

kinds

method for to randomize

first

TRUE or FALSE - randomize rep 1

Value

• trt
  vector, name of the treatments
• r
  vector, numeric
• number
  Numeric
• seed
  Numeric
• kinds
  character
• first
  Logic

Details

kinds <- c("Wichmann-Hill", "Marsaglia-Multicarry", "Super-Duper", "Mersenne-Twister", "Knuth-TAOCP", "user-supplied", "Knuth-TAOCP-2002", "default" )

References

Introduction to Experimental Statistics. Ching Chun Li. McGraw-Hill Book Company, INC, New. York, 1964

See Also

design.crd, design.lsd, random.ab, fact.nk

Examples

library(agricolae)
# 5 treatments and 6 blocks
trt<-c("A","B","C","D","E")
rcbd <-design.rcbd(trt,6,number=101,986,"Wichmann-Hill") # seed = 986
rcbd # field book
# write in hard disk
# write.table(rcbd,"rcbd.txt", row.names=FALSE, sep="\t")
# file.show("rcbd.txt")
# Plots in field model ZIGZAG
plots <-as.numeric(rcbd[,1])
block<- as.numeric(rcbd[,2])
trt <-as.character(rcbd[,3])
nplots<- length(plots)
nblock<- nlevels(rcbd[,2])
ntrt<- nlevels(rcbd[,3])
indice<-1:nplots
for (i in 1:nplots){
if ( 2*floor(block[i]/2) == block[i])indice[i]<-2*block[i]*ntrt-i-ntrt+1
}
plots<-plots[indice]
trt<-trt[indice]
dim(plots)<-c(ntrt,nblock)
dim(trt)<-c(ntrt,nblock)
print(t(plots))
print(t(trt))