Gente, quero fazer uma correlacção entre n variáveis, no Hmisc tem a função rcorr,vejam o exemplo abaixo: require(HmisC) 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)) No meu banco são 90 variáveis como fzer o teste acima sem ter que declarar uma por uma? Algo do tipo rcorr(cbind(x:v)) Bom dia! Edson Lira Estatístico Manaus-Amazonas
Cada uma é um vetor? Não estão todas no mesmo matriz/data.frame? Você testou? require(Hmisc) str(iris) rcorr(cbind(iris[,1], iris[,2])) rcorr(as.matrix(iris[,1:4])) ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Edson Lira -
Walmes Zeviani