#########################################################################
#######################Planejamento do DQL###############################
#########################################################################
require(agricolae)
t = 5 ##numero de tratamento
trat=seq(1,t,1)
DQL=design.lsd(trat, number = 1, seed = 12345)
Croqui1=matrix(DQL[,4],t,t)
Croqui1

#########################################################################
############################Análise do DQL###############################
#########################################################################

##Mudar o diretório###
setwd("J:/Veterinaria/Conjuntos de Dados")

##Chamar o pacote xlsx para ler o arquivo de dados
require(xlsx)


##Lendo o arquivo de dados
dados=read.xlsx("anestesico.xls",1)
dados
attach(dados)
summary(dados)


######Análise de Variância#####################################3
model=aov(VR~Coluna+Linha+tratamento)
anova(model)

##########Verificar as pressuposições do modelo#####################33

####Independencia
require(bstats)
dw.test(model)

##Obter o erros do modelo#####################
e=resid(model)
hist(e)

####Verificar a normalidade##################
#Teste de Shapiro-Wilk
shapiro.test(e) 

######Homogeneidade de Variância#########
#teste de Bartlett 
bartlett.test(VR~Coluna+Linha+tratamento)

###Teste de aditividade dos efeitos de tratamento e bloco
require(dae)
model1=aov(VR~Coluna+Linha+tratamento+Error(Coluna*Linha))
tukey.1df(model1, dados, error.term = "Coluna:Linha")


require(ExpDes)
####Anova com Teste de Tukey####################
latsd(tratamento, Linha,Coluna,VR, quali = TRUE, mcomp = "tukey", sigT = 0.05, sigF = 0.05)

####Anova com Teste de SNK####################
latsd(tratamento, Linha,Coluna,VR, quali = TRUE, mcomp = "snk", sigT = 0.05, sigF = 0.05)

####Anova com Teste de Duncan####################
latsd(tratamento, Linha,Coluna,VR, quali = TRUE, mcomp = "ducan", sigT = 0.05, sigF = 0.05)

####Anova com Teste de Scott-Knott####################
latsd(tratamento, Linha,Coluna,VR, quali = TRUE, mcomp = "sk", sigT = 0.05, sigF = 0.05)