Dúvida no rotina rcorr no pacote Hmisc
Como faço para colocar os nomes das variáveis na rotina abaixo. Faço a rotina para qa correlação:require(Hmisc) rcorr(cbind(corre$prevalência,corre$precipitação,corre$temp_max, + corre$temp_min) + ,type=c("pearson")) Esta saída é gerada. Como faço para colocar os nomes da variável. [,1] [,2] [,3] [,4] [1,] 1.00 -0.04 -0.15 -0.28 [2,] -0.04 1.00 -0.80 -0.41 [3,] -0.15 -0.80 1.00 0.67 [4,] -0.28 -0.41 0.67 1.00 n= 63 P [,1] [,2] [,3] [,4] [1,] 0.7841 0.2516 0.0249 [2,] 0.7841 0.0000 0.0010 [3,] 0.2516 0.0000 0.0000 [4,] 0.0249 0.0010 0.0000 Edson Lira Estatístico Manaus-Amazonas
Edson, Acho que depende da forma como você fornece os dados. Eu fiz assim x <- matrix(rnorm(1000), ncol=5) x <- as.data.frame(x) str(x) rcorr(as.matrix(x)) E tive o nome nas colunas. O cbind() que você usa não preserva os nomes, veja os dois casos abaixo rcorr(cbind(x$V1,x$V2)) rcorr(cbind(V1=x$V1, V2=x$V2)) Além disso você pode mudar o nome das colunas/linhas de uma matriz com col/rownames(). À disposição. Walmes. Em 24 de março de 2011 11:40, Edson Lira <edinhoestat@yahoo.com.br>escreveu:
Como faço para colocar os nomes das variáveis na rotina abaixo.
Faço a rotina para qa correlação:
require(Hmisc)
rcorr(cbind(corre$prevalência,corre$precipitação,corre$temp_max,
+ corre$temp_min)
+ ,type=c("pearson"))
Esta saída é gerada. Como faço para colocar os nomes da variável.
[,1] [,2] [,3] [,4]
[1,] 1.00 -0.04 -0.15 -0.28
[2,] -0.04 1.00 -0.80 -0.41
[3,] -0.15 -0.80 1.00 0.67
[4,] -0.28 -0.41 0.67 1.00
n= 63
P
[,1] [,2] [,3] [,4]
[1,] 0.7841 0.2516 0.0249
[2,] 0.7841 0.0000 0.0010
[3,] 0.2516 0.0000 0.0000
[4,] 0.0249 0.0010 0.0000
Edson Lira Estatístico Manaus-Amazonas
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-- ========================================================================== 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 ==========================================================================
Valeu Walmes, a rotina rcorr(as.matrix(x)) funcionou. Muito obrigado. Edson Lira Estatístico Manaus-Amazonas --- Em qui, 24/3/11, Walmes Zeviani <walmeszeviani@gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani@gmail.com> Assunto: Re: [R-br] Dúvida no rotina rcorr no pacote Hmisc Para: r-br@listas.c3sl.ufpr.br Data: Quinta-feira, 24 de Março de 2011, 11:54 Edson, Acho que depende da forma como você fornece os dados. Eu fiz assim x <- matrix(rnorm(1000), ncol=5) x <- as.data.frame(x) str(x) rcorr(as.matrix(x)) E tive o nome nas colunas. O cbind() que você usa não preserva os nomes, veja os dois casos abaixo rcorr(cbind(x$V1,x$V2)) rcorr(cbind(V1=x$V1, V2=x$V2)) Além disso você pode mudar o nome das colunas/linhas de uma matriz com col/rownames(). À disposição. Walmes. Em 24 de março de 2011 11:40, Edson Lira <edinhoestat@yahoo.com.br> escreveu: Como faço para colocar os nomes das variáveis na rotina abaixo. Faço a rotina para qa correlação: require(Hmisc) rcorr(cbind(corre$prevalência,corre$precipitação,corre$temp_max, + corre$temp_min) + ,type=c("pearson")) Esta saída é gerada. Como faço para colocar os nomes da variável. [,1] [,2] [,3] [,4] [1,] 1.00 -0.04 -0.15 -0.28 [2,] -0.04 1.00 -0.80 -0.41 [3,] -0.15 -0.80 1.00 0.67 [4,] -0.28 -0.41 0.67 1.00 n= 63 P [,1] [,2] [,3] [,4] [1,] 0.7841 0.2516 0.0249 [2,] 0.7841 0.0000 0.0010 [3,] 0.2516 0.0000 0.0000 [4,] 0.0249 0.0010 0.0000 Edson Lira Estatístico Manaus-Amazonas _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -- ========================================================================== 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 ========================================================================== -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
participantes (2)
-
Edson Lira -
Walmes Zeviani