symulacje

Porównanie oszacowań metodami bootstrap, DRG I IRG

x=rnorm(1000,100,5)
n=50; G=5
### metoda bootstrap
B=10000
sx=sample(x,n)
zB=c()
for (i in 1:B)
{zB=c(zB,mean(sample(sx,n,replace=TRUE)))}
var(zB)
### metoda DRG
sx=sample(x,n)
zDRG=c()
for (g in 1:G)
{zDRG=c(zDRG,mean(sx[(1+(g-1)*10):(g*10)]))}
var(zDRG)/G
### metoda IRG
zIRG=c()
for (g in 1:G)
{zIRG=c(zIRG,mean(sample(x,10)))}
var(zIRG)/G