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