Vou direto ao ponto.

library(Hmisc)
rcorr(as.matrix(corr), type="pearson") 
z<-rcorr(as.matrix(corr), type="pearson") 
print(z$P, digits = 5)
print(z$r, digits=4)

Claro, vocês já notaram que eu tenho um data frame chamado “corr” que li como uma matriz e calculei a correlação para ela. Os dois últimos comandos nos dão os p-valores das correlações (P) e as ditas cujas (r) com um número de dígitos mais civilizado (ou nem tanto).

É, estou fazendo estas brincadeiras por aqui.

Mais dicas? Bem, vejam isto.