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

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