symulacje

##########################################################################################
###                                                                                     Realizacja procedury bootstrap                                                                   ###
##########################################################################################

data(cars)
attach(cars)
N=1000
z=dist/speed
n=length(z); wynik=matrix(0,N,1)
for (i in 1:N)

{wynik[i,1]=mean(sample(z,n,replace=TRUE))}

print('Ocena wartości przeciętnej ilorazu')
mean(wynik)
print('ocena wariancji estymatora')
var(wynik[,1])