[R-br] Extrair p valor de uma matriz de correlação [Resolvido]

[Walmes Zeviani]

Será que você não tá fazendo a leitura errada do resultado? Os p-valores
(elementos acima da diagonal) são todos positivos, veja

> # pairwise sample size
> # Gabor G - 11/23/2004 R-help List
>
> pn <- function(X){ crossprod(!is.na(X)) }
> cor.prob <- function(X){
+   ## Correlations Below Main Diagonal
+   ## Significance Tests with Pairwise Deletion
+   ## Above Main Diagonal
+   ## Believe part of this came from Bill Venables
+   pair.SampSize <- pn(X)
+   above1 <- row(pair.SampSize) < col(pair.SampSize)
+   pair.df <- pair.SampSize[above1] - 2
+   R <- cor(X, use="pair")
+   above2 <- row(R) < col(R)
+   r2 <- R[above2]^2
+   Fstat <- (r2*pair.df)/(1-r2)
+   R[above2] <- 1-pf(Fstat, 1, pair.df)
+   R
+ }
>
> mydata <- matrix(rnorm(1000), ncol=10)
> round(cor.prob(mydata),2)
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6] [,7]  [,8] [,9] [,10]
 [1,]  1.00  0.07  0.92  0.71  0.29  0.52 0.05  0.80 0.08  0.77
 [2,]  0.18  1.00  0.01  0.62  0.45  0.53 0.63  0.05 0.89  0.41
 [3,]  0.01 -0.25  1.00  0.63  0.36  0.73 0.94  0.09 0.28  0.67
 [4,] -0.04  0.05 -0.05  1.00  0.11  0.76 0.92  0.26 0.63  0.75
 [5,] -0.11 -0.08  0.09  0.16  1.00  0.13 0.04  0.20 0.78  0.94
 [6,]  0.07  0.06  0.04 -0.03 -0.15  1.00 0.53  0.02 0.14  0.85
 [7,]  0.20  0.05 -0.01  0.01 -0.20 -0.06 1.00  0.92 0.79  0.09
 [8,]  0.03  0.19 -0.17  0.11 -0.13 -0.23 0.01  1.00 0.29  0.02
 [9,]  0.17 -0.01 -0.11  0.05  0.03  0.15 0.03 -0.11 1.00  0.12
[10,]  0.03  0.08  0.04 -0.03  0.01 -0.02 0.17 -0.24 0.16  1.00
>

As correlações estão abaixo da diagonal.

À disposição.
Walmes.

==========================================================================
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 em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111212/1cd2c922/attachment.html>