deca = read.table("deca.txt")
View(deca)

N = length(deca$ulica)
n = 200

# jedinke koje vadimo iz populacije
indeks = sample(deca$ulica, n, replace = F)
# vrednosti obelezja na jedinkama iz uzorka
uzorak = deca$br_dece[indeks]
uzorak

t_ocena = N*mean(uzorak)
t_ocena

# P(t u I_poverenja_t_90)=0.9
alpha = 0.1
z = qnorm(1-alpha/2)
sn = sqrt(var(uzorak))
I_poverenja_t_90 = c(t_ocena-z*sn*sqrt(N*(N-n)/n), t_ocena+z*sn*sqrt(N*(N-n)/n))
I_poverenja_t_90

sum(deca$br_dece)
# upada u interval poverenja