# script_epjack
# Autor
# 01/05/2020
# Determinar a estimativa jackknife do erro padrão
rm(list = ls())

# Semente
set.seed(12345)

# 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 vetor
epjack <- function(vetor){
  # comprimento do vetor
  n <- length(vetor)
  # vetor para guardar as réplicas jackknife
  rjack <- rep(0,n) 
  # réplicas jackknife
  for(i in 1:n){
    # amostra jackknife
    ajack <- vetor[-i]
    rjack[i] <- media(ajack)
  }
  # estimativa jackknife
  ejack <- media(rjack)
  # soma da estimativa jackknife do erro padrão
  d2 <- 0
  for(j in 1:n){
    d2 <- d2+((rjack[j]-ejack)^2)
  }
  # estimativa jackknife do erro padrão
  out <- sqrt(((n-1)/n)*d2)
  # saída da função
  return(out)
}

# amostra original
DAP <- c(34.15,36.89,35.41,34.01,39.49,34.49,32.61,
         34.33,34.16,33.56)
# estimativa jackknife do erro padrão 
epjack(DAP)