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

ASANTOS alexandresantosbr em yahoo.com.br
Quarta Setembro 2 00:31:06 BRT 2015


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 em yahoo.com.br
         alexandre.santos em 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 em n4.nabble.com 
> <mailto:ml-node+s2285057n4664976h24 em 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-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) 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
>     <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 em 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.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150901/86e536ec/attachment.html>


Mais detalhes sobre a lista de discussão R-br