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

[Hélio Gallo Rocha]

Alexandre. Estou num celular,  mas acho que incluindo o caracter &
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:[hidden email]
<http:///user/SendEmail.jtp?type=node&node=4665035&i=0>
         [hidden email]
<http:///user/SendEmail.jtp?type=node&node=4665035&i=1>
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)

_______________________________________________
R-br mailing list
[hidden email] <http:///user/SendEmail.jtp?type=node&node=4665035&i=2>
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.


------------------------------
If you reply to this email, your message will be added to the discussion
below:
http://r-br.2285057.n4.nabble.com/R-br-Exibicao-de-matriz-de-correlacao-dada-uma-condicao-tp4664911p4665035.html
To unsubscribe from R-br, click here
<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4>
.
NAML
<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150903/62ecc309/attachment.html>