populacija = 1:100
# obelezje: redni broj jedinke
k = 15
n = 20

# 15 puta vadimo psu bp obima n, smestamo ih u listu
uzorci = list()
for(i in 1:15) {
  uzorci[[i]] = sample(populacija, 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] = length(populacija)*mean(uzorci[[i]])
}
t_ocene

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

# 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))]
t_najbliza

# najreprezentativniji uzorak - onaj na kome je vrednost statistike N*Xn jednaka vrednosti t_najbliza
# which(t_ocene == t_najbolja) - vraca indeks iz vektora t_ocene cija je vrednost jednaka vrednosti t_najbliza
najreprezentativniji = which(t_ocene == t_najbliza)
najreprezentativniji_uzorak = uzorci[[najreprezentativniji]]
najreprezentativniji_uzorak