[R-br] Exibição de matriz de correlação dada uma condição
Hélio Gallo Rocha
heliogallorocha em gmail.com
Quarta Setembro 2 13:09:20 BRT 2015
Boa tarde Alexandre,
Alexandre
A linha que sugeri é para ser rodada após a função do Walmes, ficando desta
forma:
#Variáveis artificiais
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)
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),2)
cor.probC <- function(x){
results <- x
ifelse(row(x)< col(x),ifelse(x<=0.05,NA,x),ifelse(x<=0.08,NA,x))
}
cr <- cor.probC(correla)
colnames(cr) <- names(dados)
rownames(cr) <- names(dados)
cr
Tentei colorir as células de acordo com as exigencias para a correlação e o
p-valor, mas não tenho conhecimento para isso.
Outra coisa, seria muito bom colocarmos lado a lado, em colunas, os valores
da correlação e o p-valor, facilitaria muito a análise...
Abraço
Em 2 de setembro de 2015 00:25, ASANTOS [via R-br] <
ml-node+s2285057n4665017h95 em n4.nabble.com> escreveu:
> Boa noite Hélio,
>
> Estou mexendo ainda, mas ainda não funcionou, veja o CRM abaixo,
> comparando a sua função e a do Walmes:
>
> #Variáveis artificiais
> 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)
>
> #Função by IFSULDEMINAS
> cor.probC <- function(x){
>
> results <- x
> ifelse(row(x)< col(x),ifelse(x<=0.05,NA,x),ifelse(x<=0.08,NA,x))
> }
> #
>
> #Teste
> cor.probC(dados)
>
>
> require(Hmisc)## Pacote
>
>
> #Função que dispõem p valor na diagonal superior e r valor na inferior
>
> 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
> }
>
> #Resultado da matriz
> cor.prob(dados)
> #
>
> --
> ======================================================================
> 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)[hidden email] <http:///user/SendEmail.jtp?type=node&node=4665017&i=0>
> [hidden email] <http:///user/SendEmail.jtp?type=node&node=4665017&i=1>
> Lattes: http://lattes.cnpq.br/1360403201088680
> ======================================================================
>
>
>
>
>
>
>
> Em 29/08/2015 11:35, Hélio Gallo Rocha escreveu:
>
> Alexandre,
>
> Conseguimos resolver, mas não rodando dentro da função, testa ai
>
> correla <- cor.prob(dados)
> correla
>
>
> corr <- ifelse(row(correla)<
> col(correla),ifelse(correla<=0.05,NA,correla),ifelse(correla<=0.08,NA,correla));corr
>
>
>
> Em 24 de agosto de 2015 16:50, ASANTOS [via R-br] <[hidden email]
> <http:///user/SendEmail.jtp?type=node&node=4665017&i=2>> escreveu:
>
>> Boa tarde Helio,
>>
>> Ainda não estou quebrando a cabeça nas partes above1 <-
>> row(pair.SampSize) < col(pair.SampSize) e above2 <- row(R) < col(R) que é
>> onde vão as condições, mas esse é um caso para Programming master em R. Se
>> funcionar essa seria uma boa função para seleção de variáveis multivariadas
>> quando se trabalha com um p muito grande,
>>
>> Vou continuar tentando,
>>
>> Abraços,
>>
>>
>> --
>> ======================================================================
>> 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: <a moz-do-not-send="true" href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112 (TIM) <a moz-do-not-send="true" href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970 (VIVO)
>> [hidden email] <http:///user/SendEmail.jtp?type=node&node=4664976&i=0>
>> [hidden email] <http:///user/SendEmail.jtp?type=node&node=4664976&i=1>
>> Lattes: http://lattes.cnpq.br/1360403201088680
>> ======================================================================
>>
>>
>>
>> Em 24/08/2015 15:35, Hélio Gallo Rocha escreveu:
>>
>> Alexandre,
>>
>> Conseguiu resolver esta questão?
>>
>> Olhei o site sugerido pelo Walmes, mas ñ consegui inserir na função.
>>
>> Se der resultado, posta na lista, pois a é bem iteressante.
>>
>> Abraço
>>
>> Em 19 de agosto de 2015 21:38, Walmes Zeviani 2 [via R-br] <[hidden
>> email] <http:///user/SendEmail.jtp?type=node&node=4664976&i=2>> escreveu:
>>
>>> Isso pode servir de inspiração.
>>>
>>>
>>> http://stackoverflow.com/questions/9439619/using-r-replace-all-values-in-a-matrix-0-1-with-0
>>>
>>> À disposição.
>>> Walmes.
>>>
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> [hidden email] <http:///user/SendEmail.jtp?type=node&node=4664945&i=0>
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem ( <http://www.leg.ufpr.br/r-br-guia>
>>> 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-tp4664911p4664945.html
>>> To unsubscribe from R-br, click here.
>>> 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>
>>>
>>
>>
>>
>> --
>> Hélio Gallo Rocha
>> IFSULDEMINAS - Câmpus Muzambinho
>>
>>
>> _______________________________________________
>> R-br mailing list[hidden email] <http:///user/SendEmail.jtp?type=node&node=4664976&i=3>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.
>>
>>
>>
>> _______________________________________________
>> R-br mailing list
>> [hidden email] <http:///user/SendEmail.jtp?type=node&node=4664976&i=4>
>> 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-tp4664911p4664976.html
>> To unsubscribe from R-br, click here.
>> 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>
>>
>
>
>
> --
> Hélio Gallo Rocha
> IFSULDEMINAS - Câmpus Muzambinho
>
>
> _______________________________________________
> R-br mailing list[hidden email] <http:///user/SendEmail.jtp?type=node&node=4665017&i=3>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.
>
>
>
> _______________________________________________
> R-br mailing list
> [hidden email] <http:///user/SendEmail.jtp?type=node&node=4665017&i=4>
> 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-tp4664911p4665017.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>
>
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150902/8bfa9506/attachment.html>
Mais detalhes sobre a lista de discussão R-br