populacija = 1:4
# Marko 1, Petar 2, Jovan 3, Sasa 6

# obelezje - broj jabuka koje poseduju
vr_obelezja_na_populaciji = c(7,2,2,6)

# Marko i Petar 1/3, Marko i Jovan 1/2, Petar i Jovan 1/6
p_uzorka = function(i,j) {
  if(i==1 && j==2) {
    return(1/3)
  }
  if(i==1 && j==3) {
    return(1/2)
  }
  if(i==2 && j==3) {
    return(1/6)
  }
  else {
    return(0)
  }
}

# Svi uzorci bez ponavljanja obima 2
N = length(populacija)
uzorci = list()
vr_obelezja_na_uzorku = list()
s = 1
for(i in 1:N) {
  for(j in 1:N) {
    if(i<j) { 
      uzorci[[s]] = c(i,j)
      vr_obelezja_na_uzorku[[s]] = c(vr_obelezja_na_populaciji[i], vr_obelezja_na_populaciji[j])
      s = s+1
    }
  }
}
length(uzorci)

# vrednosti koje uzima slucajna velicina t_ocena = N*Xn
t_ocene = c()
for(i in 1:length(uzorci)) {
  t_ocene[i] = N*mean(vr_obelezja_na_uzorku[[i]])
}

t_ocene

# verovatnoce sa kojima uzima te vrednosti
p_t_ocene = c()
for(i in 1:length(uzorci)) {
  p_t_ocene[i] = p_uzorka(uzorci[[i]][1],uzorci[[i]][2])
}
p_t_ocene

E_t_ocene = sum(t_ocene*p_t_ocene)
E_t_ocene

# Suma obelezja na populaciji (ukupan broj jabuka koje poseduju njih cetvorica)
sum(vr_obelezja_na_populaciji)

# Nije nepristrasna