# Dados
  setwd("E:\\")
  dados=read.table("dados.csv",sep=";",dec=",",header=TRUE)
  y=dados$y
  x=dados$x

# Pressupostos do modelo de efeitos fixos
  par(mfrow=c(2,2)) ; plot(aov(y ~ x))
  bartlett.test(y ~ x)
  shapiro.test(y)
  shapiro.test( residuals(aov(y~x))  )

# Transformação Box-Cox
  library(MASS)
  boxcox(y ~ x,lambda = seq(-5, 2, 1/10))

  lambda=-1.395
  ynew=((y^lambda)-1)/lambda
  par(mfrow=c(2,2)) ; plot(aov(ynew ~ x))
  bartlett.test(ynew~x)