Hmisc
package's rcorr()
function.
myrho <- function (x, y) { rhosumm <- rcorr(x, y) rho <- round(rhosumm$r[1,2], 3) pval <- round(rhosumm$P[1,2], 4) c(Rho = rho, p.value = pval) }
# Example y <- c(4, 1, 0, 1, 4) z <- c(1, 2, 3, 4, NA) myrho(y, z)