
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