# *******************************************************************
# 1. Definição da trilha de dados
# *******************************************************************
setwd("C:\\dados\\")


# *******************************************************************
# 2. Leitura dos dados
# *******************************************************************
X<-read.table("c:\\dados\\Plano1.dat",header=F)
# Listar dados
X


# *******************************************************************
# 3. Anova em Blocos ao acaso
# *******************************************************************
Trat<-as.factor(X[,1])
Blocos<-as.factor(X[,2])
#Número de fatores (tratamento e bloco)
col_inicial<-3
#rotina ANOVA todas as variáveis
for(i in col_inicial:ncol(X))
{
Modelo<-aov(X[,i]~Blocos+Trat)
print(summary(Modelo))
}


# *******************************************************************
# 4a. Teste de Normalidade  (Shapiro-Wilk normality test)
# *******************************************************************
variavel<-3
Trat<-as.factor(X[,1])
Blocos<-as.factor(X[,2])
Modelo<-aov(X[,variavel]~Blocos+Trat)
summary (Modelo)
resid<-residuals(Modelo)
shapiro.test (resid)


# *******************************************************************
# 4b. Teste de Normalidade  (Shapiro-Wilk normality test)
# *******************************************************************
variavel<-3
Trat<-as.factor(X[,1])
Blocos<-as.factor(X[,2])
Modelo<-aov(X[,variavel]~Blocos+Trat)
anova<-anova(Modelo)
anova
resid<-residuals(Modelo)
shapiro.test (resid)