prihodi = read.table("baze\\prihodi.txt")
podaci = prihodi$V1
# varijacioni niz
podaci = sort(podaci) 
n = length(podaci)  
 
# ocena metodom zamene
Tn = sum(podaci*((2*(1:n) - 1)/n - 1))/sum(podaci)

# ocene dobijene na osnovu uzoraka bez i-te opservacije, i = 1,2,...,n 
T_i = rep(0, n)
for (i in 1:n) {
  y = podaci[-i]
  T_i[i] = sum(y*((2*(1:(n-1)) - 1)/(n - 1) - 1))/sum(y)
}
Tn_dot = mean(T_i)

# Dzeknajf ocena
(Tn_jack = n*Tn - (n - 1)*Tn_dot)