Alpha lattice design in R software

ALPHA LATTICE DESIGN ANALYSIS in R
first method
change into factors using factor function
str(w)
w$Rep=factor(w$Rep)
w$BLOCK=factor(w$BLOCK)
w$Trt=factor(w$Trt)
model
lats=aov(GY~Rep+Trt+Rep:BLOCK,data=w)
summary(lats)
#mean separation
HSD
library(agricolae)
HSD.test(lats,"Trt",alpha = 0.05,console=TRUE)
#LSD
LSD.test(lats,"Trt",alpha = 0.05,console=TRUE)
#Ducan
duncan.test(lats,"Trt",alpha = 0.05,console=TRUE)
#schefee
scheffe.test(lats,"Trt",alpha = 0.05,console=TRUE)
#######################################################
Second method
REML
PBIB.test(w$BLOCK,w$Trt,w$Rep,w$GY,13,method="REML",
test = "lsd", alpha=0.05,console=TRUE, group=TRUE)
#VC
PBIB.test(w$BLOCK,w$Trt,w$Rep,w$GY,13,method="VC",
test = "lsd", alpha=0.05,console=TRUE, group=TRUE)
#######################################################