[R-br] Exibição de matriz de correlação dada uma condição

ASANTOS alexandresantosbr em yahoo.com.br
Quarta Setembro 2 22:35:44 BRT 2015


Boa noite Helio,

        O problema é que a função some também com variáveis 
significantes mas inversamente correlacionadas, se altero 
ifelse(x>=0.08,NA,x), tenho o problema inverso. Veja as alterações que 
fiz, também criei uma estrutura de correlação um pouco melhor nos dados 
artificiais. Outro problema para colorir estão sendo os NA's, mas acho 
que estamos chegando lá, segue o CRM:

#Variáveis artificiais
dados<-runif(1000)
dados<-(dados*100)+c(1:500, 500:1)
dados<-matrix(dados,nrow=100,ncol=10)
dados<-as.data.frame(dados)
names(dados)<-c("var1", 
"var2","var3","var4","var5","var6","var7","var8","var9","var10")


#Função by Walmes
pn <- function(X){crossprod(!is.na(X))}

cor.prob <- function(X){
      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
      }


correla <- round(cor.prob(dados),4)


#Função by Helio

sig<-0.05## Significância
r<-0.30 ## Correlação

cor.probC <- function(x){

     results <- x
     ifelse(row(x)< col(x),ifelse(x<=sig,NA,x),ifelse(x<=r,NA,x))
}

cr <- cor.probC(correla)
colnames(cr)  <- names(dados)
rownames(cr) <- names(dados)
cr
#

Abraço,

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================










Em 02/09/2015 12:09, Hélio Gallo Rocha escreveu:
> var1<-sort(rnorm(30,1), decreasing = TRUE)
> var2<-rnorm(30,10)
> var3<-sort(rnorm(30,22))
> var4<-sort(rnorm(30,10))
> var5<-sort(rnorm(30,0.5), decreasing = TRUE)
> var6<-sort(rnorm(30,25), decreasing = TRUE)
> var7<-rnorm(30,10)
> dados<-cbind(var1,var2,var3,var4,var5,var6,var7)



Mais detalhes sobre a lista de discussão R-br