Rodrigo, tem uma rotina no pacote Hmisc que talvez possa te ajudar.

require(Hmisc)
?rcorr

x <- c(-2, -1, 0, 1, 2)
y <- c(4,   1, 0, 1, 4)
z <- c(1,   2, 3, 4, NA)
v <- c(1,   2, 3, 4, 5)
rcorr(cbind(x,y,z,v))

[  ]'s.
Edson Lira
Estatístico
Manaus-Amazonas


Em Segunda-feira, 24 de Março de 2014 23:04, Mauro Sznelwar <sznelwar@uol.com.br> escreveu:
Este código não funciona.
 
O seu loop for tá sobre escrevendo resultados em
cor.mult<- cor.test(x, y[,i])
e portanto somente o ultimo está  sendo guardado

(note tb que sua matrz de dados esá reciclando valores)

Uma sugestão:

x <- 1:10
y <- matrix(rnorm(200), nrow=10)
apply(y, 2, function(y, x=x) cor.test(x,y)$p.value



_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.