Exibição de matriz de correlação dada uma condição
Prezados listeiros, No meu CRM abaixo tenho uma matriz de correlação com o valor de p na diagonal superior o e valor r de Pearson da diagonal inferior, no entanto eu gostaria de exibir uma matriz de correlação apenas com os p valores maiores ou igual 0.05 e r maior ou igual a 0.08 isso é possível? Segue CRM: require(Hmisc)## Pacote #Variáveis artificiais var1<-sort(rnorm(500,1), decreasing = TRUE) var2<-rnorm(500,10) var3<-sort(rnorm(500,22)) var4<-sort(rnorm(500,10)) var5<-sort(rnorm(500,0.5), decreasing = TRUE) var6<-sort(rnorm(500,25), decreasing = TRUE) var7<-rnorm(500,10) dados<-cbind(var1,var2,var3,var4,var5,var6,var7) #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) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
Alexandre, Salvo se eu me perdi na leitura do seu texto e na sua função, basta você decidir o quê você vai colocar nas células da matriz quando os valores não estiverem como você deseja e usar um teste if com cláusula else nas atribuições, ou uma construção similar à que você já faz para determinar se vai escrever na diagonal superior ou inferior com seus vetores above1 e above2. HTH 2015-08-16 23:50 GMT-03:00 ASANTOS <alexandresantosbr@yahoo.com.br>:
Prezados listeiros,
No meu CRM abaixo tenho uma matriz de correlação com o valor de p na diagonal superior o e valor r de Pearson da diagonal inferior, no entanto eu gostaria de exibir uma matriz de correlação apenas com os p valores maiores ou igual 0.05 e r maior ou igual a 0.08 isso é possível?
Segue CRM:
require(Hmisc)## Pacote
#Variáveis artificiais var1<-sort(rnorm(500,1), decreasing = TRUE) var2<-rnorm(500,10) var3<-sort(rnorm(500,22)) var4<-sort(rnorm(500,10)) var5<-sort(rnorm(500,0.5), decreasing = TRUE) var6<-sort(rnorm(500,25), decreasing = TRUE) var7<-rnorm(500,10) dados<-cbind(var1,var2,var3,var4,var5,var6,var7)
#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) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
Isso pode servir de inspiração. http://stackoverflow.com/questions/9439619/using-r-replace-all-values-in-a-m... À disposição. Walmes.
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] < ml-node+s2285057n4664945h63@n4.nabble.com> escreveu:
Isso pode servir de inspiração.
http://stackoverflow.com/questions/9439619/using-r-replace-all-values-in-a-m...
À 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) 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... 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
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: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br 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] <ml-node+s2285057n4664945h63@n4.nabble.com <mailto:ml-node+s2285057n4664945h63@n4.nabble.com>> escreveu:
Isso pode servir de inspiração.
http://stackoverflow.com/questions/9439619/using-r-replace-all-values-in-a-m...
À 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) 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...
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
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
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] < ml-node+s2285057n4664976h24@n4.nabble.com> 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: (+55) 65 8132-8112 (TIM) (+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-m...
À 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) 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... 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... 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
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) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br 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] <ml-node+s2285057n4664976h24@n4.nabble.com <mailto:ml-node+s2285057n4664976h24@n4.nabble.com>> 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:(+55) 65 8132-8112 <tel:%28%2B55%29%2065%208132-8112> (TIM)(+55) 65 9686-6970 <tel:%28%2B55%29%2065%209686-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-m...
À 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) 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...
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...
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
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
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@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-m...
À 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... 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... 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... 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
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@yahoo.com.br alexandre.santos@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)
Vou ver se tem como Em 02/09/2015 22:36, "ASANTOS [via R-br]" < ml-node+s2285057n4665035h53@n4.nabble.com> escreveu:
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... 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>
Alexandre. Estou num celular, mas acho que incluindo o caracter & no campo relativo a correlação, vai funcionar Em 02/09/2015 22:36, "ASANTOS [via R-br]" < ml-node+s2285057n4665035h53@n4.nabble.com> escreveu:
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... 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>
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... 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>
participantes (4)
-
ASANTOS -
Cesar Rabak -
Hélio Gallo Rocha -
Walmes Zeviani