[R-br] Extrair p valor de uma matriz de correlação [Resolvido]

[ASANTOS]

   Walmes e Leonard,

             Comparei cor(dados) com cor.prob(dados) e realmente os p 
valores estão na parte superior e o r na parte inferior, dei bobeira é 
só depois vi isto explicitado na função,

Obrigado,



Em 12-12-2011 13:07, Leonard de Assis escreveu:
> a parte de cima da Matriz são os P, a parte de baixo são as correlações
> []s
> Leonard de Assis
> assis<dot>  leonard<at>  gmail<dot>  com
>
> Em 12/12/2011 11:31, ASANTOS escreveu:
>>   Obrigado Walmes,
>>
>>            Deu certo, apenas achei estranho o teste apresentar 
>> valores de p negativos, solução dada por Gabor Grothendieck e Bill 
>> Venables abaixo:
>>
>>
>> require(Hmisc)
>>
>> x<-rnorm(100)
>> y<-rnorm(100)
>> z<-rnorm(100)
>> dados<-cbind(x,y,z)
>>
>> 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
>> }
>>
>> cor.prob(dados)
>>
>> Em 12-12-2011 11:11, Walmes Zeviani escreveu:
>>> Alexandre,
>>>
>>> Use as funções para correlação disponíveis no pacote Hmisc ou use a 
>>> função cor.prob() disponível nessa mensagem 
>>> (http://tolstoy.newcastle.edu.au/R/help/05/07/9122.html).
>>>
>>> À disposição.
>>> Walmes.
>>>
>>> ==========================================================================
>>> Walmes Marques Zeviani
>>> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 
>>> 49.231759 W)
>>> Departamento de Estatística - Universidade Federal do Paraná
>>> fone: (+55) 41 3361 3573
>>> VoIP: (3361 3600) 1053 1173
>>> e-mail: walmes em ufpr.br <mailto:walmes em ufpr.br>
>>> twitter: @walmeszeviani
>>> homepage: http://www.leg.ufpr.br/~walmes 
>>> <http://www.leg.ufpr.br/%7Ewalmes>
>>> linux user number: 531218
>>> ==========================================================================
>>>
>>>
>>> _______________________________________________
>>> 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.
>>
>>
>> -- 
>> Alexandre DOS SANTOS
>> Engenheiro Florestal, Msc.
>> Laboratório de Entomologia Florestal
>> Departamento de Entomologia
>> Universidade Federal de Lavras
>> Caixa Postal 3037
>> 37200-000 - Lavras - Minas Gerais - Brasil
>> Tel: +55 35 92230304
>>
>>
>> _______________________________________________
>> 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.
>
>
> _______________________________________________
> 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.


-- 
Alexandre DOS SANTOS
Engenheiro Florestal, Msc.
Laboratório de Entomologia Florestal
Departamento de Entomologia
Universidade Federal de Lavras
Caixa Postal 3037
37200-000 - Lavras - Minas Gerais - Brasil
Tel: +55 35 92230304

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