<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Obrigado Walmes,<br>
<br>
Deu certo, apenas achei estranho o teste apresentar
valores de p negativos, solução dada por Gabor Grothendieck e Bill
Venables abaixo:<br>
<br>
<br>
require(Hmisc)<br>
<br>
x<-rnorm(100)<br>
y<-rnorm(100)<br>
z<-rnorm(100)<br>
dados<-cbind(x,y,z)<br>
<br>
pn <- function(X){crossprod(!is.na(X))}<br>
<br>
cor.prob <- function(X){<br>
pair.SampSize <- pn(X)<br>
above1 <- row(pair.SampSize) < col(pair.SampSize)<br>
pair.df <- pair.SampSize[above1] - 2<br>
R <- cor(X, use="pair")<br>
above2 <- row(R) < col(R)<br>
r2 <- R[above2]^2<br>
Fstat <- (r2 * pair.df)/(1 - r2)<br>
R[above2] <- 1 - pf(Fstat, 1, pair.df)<br>
R<br>
}<br>
<br>
cor.prob(dados) <br>
<br>
Em 12-12-2011 11:11, Walmes Zeviani escreveu:
<blockquote
cite="mid:CAFU=Ekafma0tY9j7jYr778YSwTzdPnC2G2P1QS4J1daabToxTA@mail.gmail.com"
type="cite"><font face="trebuchet ms,sans-serif">Alexandre,<br>
<br>
Use as funções para correlação disponíveis no pacote Hmisc ou
use a função cor.prob() disponível nessa mensagem (<a
moz-do-not-send="true"
href="http://tolstoy.newcastle.edu.au/R/help/05/07/9122.html">http://tolstoy.newcastle.edu.au/R/help/05/07/9122.html</a>).<br>
<br>
À disposição.<br>
Walmes.</font><br>
<font face="trebuchet ms,sans-serif"><br clear="all">
</font><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques
Zeviani</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório
de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de
Estatística - Universidade Federal do Paraná</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41
3361 3573</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361
3600) 1053 1173</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">e-mail: <a
moz-do-not-send="true" href="mailto:walmes@ufpr.br"
target="_blank">walmes@ufpr.br</a></span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">twitter:
@walmeszeviani</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">homepage: <a
moz-do-not-send="true" href="http://www.leg.ufpr.br/%7Ewalmes"
target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user
number: 531218</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="72">--
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</pre>
</body>
</html>