library(agricolae)

dados <- read.table("vermiwash.txt", header = TRUE)
head(dados)
View(dados)

DOSE <- factor(dados$TREAT)
levels(DOSE)

CUT <- factor(dados$CUT)
levels(CUT)

#com a variável 1
model_var1<-aov(VAR1~DOSE*CUT, data=dados)
summary(model_var1)
cv.model(model_var1)


#
#var1
tukey_var1 <- HSD.test(model_var1, c("DOSE", "CUT"), main="VAR1 ~ DOSE*CUT", console=TRUE)
head(tukey_var1)

#
x11() 
bar.group(tukey_var1$groups,main="Teste de Tukey", ylim=c(0,50), xlab="Tratamentos (Dosagem)") 

#OR

#var1
compara.tukey.agri <- HSD.test(model_var1, "DOSE" ,group=T) 
compara.tukey.agri 

x11() 
bar.group(compara.tukey.agri$groups,main="Teste de Tukey", ylim=c(0,50), xlab="Tratamentos (DOSE)")