Hi, A recent issue [1] brought up the idea of adding support for complex data in stats.pearsonr. It seems like this would be a useful feature, especially for folks working in signal processing. It wouldn't be too difficult to implement, see [2]. The current estimation of p-values, however, seems to underestimate the significance. That is, the error rate tends to be lower than the p-value. It would be possible to introduce a permutation test to estimate the p-value of the correlation between complex numbers. Is this something worth pursuing? What do you think of introducing an alternative p-value estimation method? Cheers, Chris [1] https://github.com/scipy/scipy/issues/17518 [2] https://en.wikipedia.org/wiki/Complex_random_variable