#Bootstrap distribution of median for Mouse observed control group
#Version 1 of program using for loop
#374 Resampling Slide 41

B<-10
origsample<-c(52,104,146,10,50,31,40,27,46)

brep<-rep(0,B)
for(i in 1:B){                             
  bsample<-sample(origsample,length(origsample),replace=T)
  brep[i]<-median(bsample)
  #readline(cat('bootstrap sample:', bsample, 'bootstrap replicate:', brep[i]))
}

hist(brep,breaks=(0:200),freq=FALSE,col="red",axes=FALSE,
  main="Bootstrap distribution for median",xlab="Median", ylab="Probability")
axis(1,at=origsample,font=8)
axis(2)

cat('median of original sample', median(origsample), '\n')
cat('bootstrap estimate of std dev', sd(brep), '\n')