for col_a, col_b in itertools.combinations(columns, 2):
    correlations[col_a + '__' + col_b] = pearsonr(df.loc[:, (col_a)], df.
        loc[:, (col_b)])
result = DataFrame.from_dict(correlations, orient='index')
result.columns = ['PCC', 'p-value']
Вычисление попарной корреляции между всеми столбцами