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

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