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

ASANTOS alexandresantosbr em yahoo.com.br
Segunda Dezembro 12 11:31:27 BRST 2011


   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

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


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