Looking for a functions that performs Tukey's transformation (as it does in r but for larger than 5000 data sets) to identify the best transformation to normalize data for parametric analyses. This would eliminate the guesswork and save everyone time.  The r reference is: 

library(rcompanion)

T_tuk = 
     transformTukey(yourdatanamehere,
                     plotit=FALSE)

    lambda     W Shapiro.p.value
397   -0.1 0.935         0.08248

if (lambda >  0){TRANS = x ^ lambda} 
if (lambda == 0){TRANS = log(x)} 
if (lambda <  0){TRANS = -1 * x ^ lambda}