title: How do I calculate lambda to use scipy.special.boxcox1p function for my entire dataframe of 500 columns?
body: <p>I have a dataframe with total sales of around 500 product categories in each row. So there are 500 columns in my dataframe. I am trying to find the highest correlated category with my another dataframe columns.
  So I will use Pearson correlation method for this.
  But the Total sales for all the categories are highly skewed data, with the skewness level ranging from 10 to 40 for all the category columns. So I want to log transform this sales data using boxcox transformation.
  Since, my sales data has 0 values as well, I want to use boxcox1p function.
  Can somebody help me, how do I calculate lambda for boxcox1p function, since it is a mandatory parameter for this function?
  Also, Is this the correct approach for my problem statement to find highly correlated categories?</p>