Cor.test() takes vectors x and y as arguments, but I have an entire matrix of data that I want to test, pairwise. Cor() t
"The accepted solution (corr.test function in the psych package) works, but is extremely slow for large matrices."
If you use ci=FALSE, then the speed is much faster.
By default, confidence intervals are found. However, this leads to a slight slowdown of speed. So, for just the rs, ts and ps, set ci=FALSE.