Bootstrapping With Missing Data (SRSWR)

 library(sampling)

library(MASS)

#data=read.csv(file.choose())

y=data[,1]; x=data[,2]*data[,4];

DF=data.frame(y,x)

attach(DF)

myr=NULL;zm=c();

N=1000; H=5000; n=30; r=10

for(i in 1:H){

S1=DF[sample(1:N,   n,  replace=T),]

Sr=S1[sample(1:n, r, replace=T),]

myr[i]=mean(Sr$y);

zm[i]=myr[i]

}

msem=(sum((zm-my)^(2)))/H

msem/mse1*100;