12 Dez
2011
12 Dez
'11
13:42
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@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================