N = 25
x = c(15.07, 8.60, 10.10, 11.22, 7.63, 10.01, 10.00, 4.74, 13.05,
      11.80, 8.12, 9.49, 11.52, 9.22, 9.27, 5.64, 11.66, 10.37,
      10.82, 5.42, 14.95, 10.46, 8.84, 16.09, 9.86)
y = c(58.23, 44.73, 39.92, 59.28, 41.90, 47.14, 55.89, 27.64, 62.04,
      45.88, 42.11, 43.42, 62.14, 57.01, 47.02, 42.27, 54.60, 54.51,
      49.42, 29.31, 56.18, 53.88, 58.25, 69.47, 47.77)

# formiramo sve uzorke bez ponavljanja obima 10
uzorci_10_ind = combn(1:N, 10, simplify = FALSE)
uzorci_10 = lapply(uzorci_10_ind, function(t) {
  rbind(x[t], y[t])
})

# vrednosti koje uzima slucajna velicina R_ocena = xn/yn na uzorcima
# posto je u pitanju prosto slucajno uzorkovanje, svaki od uzoraka ima 
# jednaku verovatnocu (1/broj uzoraka) 
R_ocena = c()
for(i in 1:length(uzorci_10)) {
  R_ocena[i] = mean(uzorci_10[[i]][1, ] / uzorci_10[[i]][2, ])
}

(ER_ocena = sum(R_ocena / length(uzorci_10)))

# kolicnik obelezja populacije
(R = mean(x) / mean(y))