# Kod 3.10. Losowanie prób bootstrapowych i ocena wariancji średniej z próby
n=20
dane = rnorm(n,mean=10,sd=1)
N_boot=1000
boot_samp=replicate(N_boot,sample(dane,n, TRUE),simplify=FALSE)
x_boot=sapply(boot_samp, mean, simplify=TRUE)
hist(x_boot,prob=TRUE,xlim=c(6,14),breaks=15,main='',ylab='gęstość')
curve(dnorm(x,10,1),add=TRUE)
curve(dnorm(x,10,1/sqrt(n)),add=TRUE,lwd=2)
