symulacje

##########################################################################################
###                                                    Ocena wariancji estymatora metodą zależnych grup losowych                                                         ###
##########################################################################################

data(chickwts)
hist(chickwts$weight)
x=sample(chickwts$weight,70)
z=c(0,0)
N=10
n=70/N
for (i in 1:N) z[i]=mean(x[((i-1)*n+1):(i*n)])
mean(z)
var1=var(z)/N
z0=mean(chickwts$weight)
var2=sum((z-z0)^2)/(N*(N-1))
print('Odchylenia standardowe - oceny')
sqrt(var1)
sqrt(var2)