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.