perform_correlations() manages all the correlation calculates and
returns a (Z x num of factors) np.array with all the correlations.
Currently only dichotomous metadata is calculated