#########################################################################
#######################Planejamento do DIC###############################
#########################################################################
require(dae)
r=4 ##numero de repetições
t=3 ##numero de tratamentos
n=r*t
DIC.parc=list(Repeticao=r, Parcela=t)
DIC.nest=list(Parcela="Repeticao")
tratamento=factor(rep(c(1:t), times=r))
DIC=fac.layout(unrandomized = DIC.parc,nested.factors=DIC.nest,randomized = tratamento , seed = 105)
DIC

#########################################################################
############################Análise do DIC###############################
#########################################################################

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

######Análise de Variância#
model=aov(VR~trat)
anova(model)

####Verificar as pressuposições do modelo
####Independencia
require(bstats)
dw.test(model)


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

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

#Teste de Kolmogorov-Smirnov
m=mean(e)
dp=sd(e)
ks.test(e,"pnorm",m,dp)


######Homogeneidade de Variância#########
#teste de Bartlett 
bartlett.test(VR~trat)

#Teste de Levene
require(car)
leveneTest(model)


require(ExpDes)
####Anova com Teste de Tukey####################
crd(trat, VR, quali = TRUE, mcomp = "tukey", sigT = 0.05, sigF = 0.05)

####Anova com Teste de SNK####################
crd(trat, VR, quali = TRUE, mcomp = "snk", sigT = 0.05, sigF = 0.05)

####Anova com Teste de Duncan####################
crd(trat, VR,  quali = TRUE, mcomp = "ducan", sigT = 0.05, sigF = 0.05)

####Anova com Teste de Scott-Knott####################
crd(trat, VR, quali = TRUE, mcomp = "sk", sigT = 0.05, sigF = 0.05)

require(ExpDes)

##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("coelhos.xls",1)
attach(dados)


####Anova com Regressão####################
####Anova com Regressão####################
crd(sub, peso, quali =FALSE,sigF = 0.05)