x = c(3, 4, 12, 9, 8, 2, 6, 15, 14, 17, 1, 7, 12)
N = length(x)
n = 4

(br_uzoraka = choose(N, n))

uzorak = list()
s = 1
for (i in 1:N) {
  for (j in 1:N) {
    for (k in 1:N) {
      for (l in 1:N) {
        if((i < j) & (j < k) & (k < l)) {
          uzorak[[s]] = c(x[i], x[j],
                          x[k], x[l])
          s = s + 1
        }
      }
    }
  }
}
length(uzorak)
# verovatnoca svakog od uzoraka je jednaka 1/br_uzoraka

(xn = sapply(uzorak, mean))
(sn2 = sapply(uzorak, var))
(D_xn_ocena = sapply(uzorak, function(t) {
  (N-n) * var(t) / (N*n)
}))
(medijana = sapply(uzorak, median))

# nepristrasnost
(Exn = sum(xn / br_uzoraka))
mean(x)

(Esn2 = sum(sn2 / br_uzoraka))
var(x)