Well, in R you could try this:

    library(MASS)
    boxcox(y~x)
    plot(1/y^2~x) # since the profile likelihood has a maximum near 2

![enter image description here][1]

But it really depends on what you mean by 'better fit to the data'


  [1]: https://i.sstatic.net/lahlV.png