```r
resamp<-replicate(1000,mean(sample(Lshell,10)))
hist(resamp,breaks=20,xlim=c(90,120),main="Sampling distribution of the means of 10 shells")
```

<img src="Quantitative_and_Spatial_files/figure-html/unnamed-chunk-56-1.png" width="672" />