[R-br] Algoritmo do boot.ci do livro do Davison, A.C. and Hinkley, D.V. (1997)!

[Ivan Bezerra Allaman]

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.

x <- c(0.75,2.55,1.40,0.80,2.60)
y <- c(0.8249,0.797,0.7762,0.7797,0.7483)

da <- data.frame(x,y)
corr(da)

mean(boot(da,corr,10000)$t)


#Esta função eu fiz para evitar NAs, já que algumas amostragens são todas iguais, gerando um desvio padrão zero.
sa <- function(x){
  y <- sample(x,replace=T)
  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]){
    z <- sample(x,replace=T)}else{
      y}
} 

boot_da  <- numeric(10000)
n <- 10000)
for(i in 1:n){
  boot_da[[i]]  <- corr(matrix(c(sa(da[,1]),sa(da[,2])),ncol=2))
   }  
boot_da

mean(boot_da) 

Se alguém puder ajudar de novo, agradeço.

Abraços!

Allaman
(S,f,P)




 
M.Sc Ivan Bezerra Allaman 
Zootecnista
Doutorando em Produção Animal/Aquicultura - UFLA 
email e msn - ivanalaman em yahoo.com.br 
Tel: (35)3826-6608/9900-2924
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110714/c0915e34/attachment.html>