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

[Leonard de Assis]

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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111212/95fd7f83/attachment.html>