BOTOX

boxcoxSUPER=function(modelo,data,respuesta,a,b,step){
  library(MASS)
  bmodelo=boxcox(modelo,data=data,lambda=seq(a,b,by=step))
  lambda=bmodelo$x[which.max(bmodelo$y)]
  if(lambda==a || lambda==b){
    warning("Probable máximo estacionario, amplíe el intervalo.")
  }
  print(lambda)
  return((data[[respuesta]]^lambda-1)/lambda)
}

#modelo es un objeto generado por lm()
#data es un data.frame
#respuesta es la variable de respuesta del modelo, debe ser puesta como "character"
#a,b son los límites del intervalo en el cual trabajará boxcox [a,b]
#step es el salto para boxcox

#PILAS SI APARECE EL WARNING: AMPLIAR EL INTERVALO!
#EJEMPLO DE USO: gatitos$transf_boxcox=boxcoxSUPER(modelo1,gatitos,"tamaño",-2,3,0.05)