x = dat$x
N = length(x)
k = 15 # broj uzoraka
n = 20 # obimi uzoraka

# 15 puta biramo psu bp obima n, smestamo ih u listu
uzorci = list()
for(i in 1:15) {
  uzorci[[i]] = sample(x, n, replace = F)
}
uzorci

# realizovane vrednosti statistike t_ocena = N*Xn na uzorcima
t_ocene = c()
for(i in 1:length(uzorci)) {
  t_ocene[i] = N * mean(uzorci[[i]])
}
t_ocene

# suma obelezja na populaciji
(t = sum(x))

# which.min(abs(t_ocene-t)) - vraca indeks vrednosti iz vektora t_ocene
# gde je apsoulutna vrednost razlike t_ocene - t najmanja
(t_najbliza = t_ocene[which.min(abs(t_ocene - t))])

# which(t_ocene == t_najbliza) - vraca indeks vrednosti iz vektora t_ocene koja je jednaka vrednosti t_najbliza
najreprezentativniji = which(t_ocene == t_najbliza)
(najreprezentativniji_uzorak = uzorci[[najreprezentativniji]])