# script_rboot # Autor # 01/05/2020 # Obter as réplicas bootstrap rm(list = ls()) # função media com argumento vetor media <- function(vetor){ # comprimento do vetor nv <- length(vetor) # soma sm <- 0 for(i in 1:nv){ sm <- sm+vetor[i] } # média out_media <- sm/nv # saída da função return(out_media) } # função com o argumento dados e B rboot <- function(dados,B){ # vetor para guardar as réplicas bootstrap rboot <- rep(0,B) # réplicas bootstrap for(j in 1:B){ # amostra bootstrap aboot <- sample(dados,length(dados),replace=TRUE) rboot[j] <- media(aboot) } # saída da função return(rboot) } # amostra original DAP <- c(34.15,36.89,35.41,34.01,39.49,34.49,32.61, 34.33,34.16,33.56) # Semente set.seed(12345) # réplicas boostrap rboot(DAP,3)