box.cox <- function(y, lambda=0) {
  if(lambda==0) log(y) else (y^lambda-1)/lambda }