for col_a, col_b in itertools.combinations(columns, 2):
    correlations[col_a + '__' + col_b] = pearsonr(df.loc[:, (col_a)], df.
        loc[:, (col_b)])
Calculating pairwise correlation among all columns