Wraps currently BoxCoxTrans.jl and provides confidence intervals (see doc).
using BoxCoxTransformations
using Distributions
x = rand(Gamma(2,2), 10000) .+ 1;
bc = fit(BoxCoxTransformation, x)
# lambda parameter
bc.lambda
# confidence intervall for lambda
confint(bc)
# transformed data
transform(bc)
See BoxCoxTransformationsMakie.jl for BoxCox-plots.