<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">A função 'corr' é específica para ser utilizada dentro da função 'boot', então no seu programa troque o boot::corr pelo stats::cor. Uma outra coisa, dentro do loop vc roda a função 'sa' separadamente para as colunas e isso é um problema, pois vc quebra qualquer estrutura de relação entre as variáveis. Em vez disso, modifique a função 'sa' para reamostrar as linhas do seu data frame e coloque 'da' como o input da função no loop. Att.,Rubem --- Em qui, 14/7/11, Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> escreveu: <blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"> De: Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Assunto: Re: [R-br] Algoritmo do boot.ci do livro do Davison, A.C. and Hinkley, D.V.

 (1997)!<br>Para: "R Brasil" <r-br@listas.c3sl.ufpr.br><br>Data: Quinta-feira, 14 de Julho de 2011, 14:44<br><br><div id="yiv759341773"><div style="color:#000;background-color:#fff;font-family:arial, helvetica, sans-serif;font-size:12pt;"><div style="font-family:arial, helvetica, sans-serif;font-size:12pt;"><span>Verdade!! Valeu. Agora uma coisa que estou quebrando a cabeça é a seguinte. Porque que o boot gerado pela função boot está extremamente diferente do boot gerado na 'unha'?. Vejamos um CMR.</span></div><div style="font-family:arial, helvetica, sans-serif;font-size:12pt;"><span><br></span></div><div><span><div>x <- c(0.75,2.55,1.40,0.80,2.60)</div><div>y <- c(0.8249,0.797,0.7762,0.7797,0.7483)</div><div><br></div><div>da <- data.frame(x,y)</div><div>corr(da)</div><div><br></div><div>mean(boot(da,corr,10000)$t)<br></div><div><br></div><div>#Esta função eu fiz para evitar NAs, já que algumas amostragens são todas iguais,

 gerando um desvio padrão zero.</div><div>sa <- function(x){</div><div>  y <- sample(x,replace=T)</div><div> 

 if(y[1:1] == y[2:2] && y[1:1] == y[3:3] && y[1:1] == y[4:4] && y[1:1] == y[5:5]){</div><div>    z <- sample(x,replace=T)}else{</div><div>      y}</div><div>} </div><div><br></div><div>boot_da  <- numeric(10000)</div><div>n <- 10000)</div><div>for(i in 1:n){</div><div>  boot_da[[i]]  <- corr(matrix(c(sa(da[,1]),sa(da[,2])),ncol=2))</div><div>   }  </div><div>boot_da</div><div><br></div><div>mean(boot_da) </div><div><br></div><div>Se alguém puder ajudar de novo, agradeço.</div><div><br></div><div>Abraços!</div><div><br></div><div>Allaman</div><div>(S,f,P)</div></span></div><div style="font-family:arial, helvetica, sans-serif;font-size:12pt;"><span><br></span></div><div style="font-family:arial, helvetica, sans-serif;font-size:12pt;"><span><br></span></div><div style="font-family:arial, helvetica, sans-serif;font-size:12pt;"><span><br></span></div><div

 style="font-family:arial, helvetica, sans-serif;font-size:12pt;"><span><br></span></div><div style="font-family:arial, helvetica, sans-serif;font-size:12pt;"> </div><div style="font-family:arial, helvetica, sans-serif;font-size:12pt;background-color:transparent;color:rgb(115, 115, 115);" align="center"><b>M.Sc Ivan Bezerra Allaman</b> <br>Zootecnista<br><font style="background-color:transparent;" face="comic sans ms" size="2"><font style="background-color:transparent;">Doutorando em Produção Animal/Aquicultura - UFLA</font> <br></font><font style="background-color:transparent;" face="comic sans ms" size="2">email e msn - ivanalaman@yahoo.com.br <br><font size="1">Tel: (35)3826-6608/9900-2924</font></font></div></div></div><br>-----Anexo incorporado-----<br><br><div class="plainMail">_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br"

 href="/mc/compose?to=R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</div></blockquote></td></tr></table>