<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif">Alexandre, bom dia!</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Se ainda lhe for Ãºtil, segue uma sugestão...</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style=""><div class="gmail_default"><font face="monospace, monospace">### <code r></font></div><div class="gmail_default"><font face="monospace, monospace">var1 <- rnorm(200,5,0.25)</font></div><div class="gmail_default"><font face="monospace, monospace">var2 <- 1:500</font></div><div class="gmail_default"><font face="monospace, monospace">var3 <- rnorm(100,5,0.25)</font></div><div class="gmail_default"><font face="monospace, monospace">var4 <- 500:1</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">cor.prob <- function(X){ <br></font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â  pair.SampSize <- crossprod(!<a href="http://is.na">is.na</a>(X)) </font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â  above1 <- row(pair.SampSize) < col(pair.SampSize) </font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â  pair.df <- pair.SampSize[above1] - 2 </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â  R <- cor(X, use="pair") </font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â  above2 <- row(R) < col(R) </font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â  r2 <- R[above2]^2 </font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â  Fstat <- (r2 * pair.df)/(1 - r2) </font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â  R[above2] <- 1 - pf(Fstat, 1, pair.df) </font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â  R </font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â  } </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">tam Â  <- 100 # amostra <br></font></div><div class="gmail_default"><font face="monospace, monospace">n_sim <- Â 99 # simulações </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">RES <- NULL <br></font></div><div class="gmail_default"><font face="monospace, monospace">for(i in 1:n_sim){</font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â dados Â  <- sapply(c("var1", "var2", "var3", "var4"), function(x) sample(get(x), 100))</font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â correla <- round(cor.prob(dados),4) </font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â res Â  Â  <- cbind(i, correla) </font></div><div class="gmail_default"><font face="monospace, monospace">  Â  Â RES Â  Â  <- rbind(RES, res) </font></div><div class="gmail_default"><font face="monospace, monospace">} </font></div><div class="gmail_default"><font face="monospace, monospace"># </font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">RES[c(1:5,392:396),]</font></div><div class="gmail_default"><font face="monospace, monospace"># Â  Â  Â  i Â  Â var1 Â  Â var2 Â  Â var3 Â  var4</font></div><div class="gmail_default"><font face="monospace, monospace"># var1 Â 1 Â 1.0000 Â 0.8025 Â 0.8451 0.9178</font></div><div class="gmail_default"><font face="monospace, monospace"># var2 Â 1 -0.0253 Â 1.0000 Â 0.6609 0.0106</font></div><div class="gmail_default"><font face="monospace, monospace"># var3 Â 1 -0.0198 Â 0.0444 Â 1.0000 0.4850</font></div><div class="gmail_default"><font face="monospace, monospace"># var4 Â 1 -0.0104 -0.2545 -0.0706 1.0000</font></div><div class="gmail_default"><font face="monospace, monospace"># var1 Â 2 Â 1.0000 Â 0.5200 Â 0.1704 0.8461</font></div><div class="gmail_default"><font face="monospace, monospace"># var4 98 Â 0.1309 Â 0.0332 -0.1361 1.0000</font></div><div class="gmail_default"><font face="monospace, monospace"># var1 99 Â 1.0000 Â 0.1433 Â 0.2047 0.3260</font></div><div class="gmail_default"><font face="monospace, monospace"># var2 99 -0.1474 Â 1.0000 Â 0.5391 0.9524</font></div><div class="gmail_default"><font face="monospace, monospace"># var3 99 -0.1279 -0.0621 Â 1.0000 0.2885</font></div><div class="gmail_default"><font face="monospace, monospace"># var4 99 -0.0992 -0.0060 Â 0.1072 1.0000</font></div><div class="gmail_default"><font face="monospace, monospace">### </code></font></div><div class="gmail_default"><font face="verdana, sans-serif"><br></font></div><div style="font-family:verdana,sans-serif"><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font face="arial, helvetica, sans-serif" style="font-size:small"><div style="font-family:'trebuchet ms',sans-serif;display:inline">​</div></font><font face="arial, helvetica, sans-serif" style="font-size:12.8px"><span style="font-size:small">================================================</span><br>Éder Comunello</font><div style="font-size:12.8px"><span style="font-family:arial,helvetica,sans-serif">Agronomist (UEM), MSc in Environ. Sciences (UEM)</span><font face="arial, helvetica, sans-serif"><br></font></div><div style="font-size:12.8px"><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">DSc in Agricultural Systems Engineering (USP/Esalq)</span><br></div><div><span style="font-size:small">Brazilian Agricultural Research Corporation (</span><font face="arial, helvetica, sans-serif" style="font-size:small">Embrapa)</font></div><div><div><font face="arial, helvetica, sans-serif" style="font-size:small">Dourados, MS, Brazil </font><font face="arial, helvetica, sans-serif" size="2">|<O>|</font></div><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">================================================</span><br></div></div></div><div style="font-size:small"><font face="arial, helvetica, sans-serif">GEO, </font>-22.2752, -54.8182, 408m<span style="font-family:arial,helvetica,sans-serif"><br></span></div><div style="font-size:12.8px"><span style="font-size:small;font-family:arial,helvetica,sans-serif">UTC-04:00 / DST: UTC-03:00</span></div><div><div><div><br></div><div><br></div></div><div style="font-size:small"><br></div></div></div></div></div></div></div></div></div></div></div></div></div>