###################################################
### chunk number 1: TableExample
###################################################
library(xtable)

samp1 <- rnorm(10)
samp2 <- rnorm(10)
samp3 <- rnorm(10)
samp4 <- rnorm(10)

samps <- cbind(samp1, samp2, samp3, samp4)
means <- apply(samps,2,mean)
samp.means <- rbind(samps, means)

colnames(samp.means) <- c(paste("Sample", 1:4))
dig <- c(0,3,3,3,3)
xtable(samp.means, digits=dig, caption="Samples of $n$=10 individuals from $N(0,1)$, together
with their respective means.")


###################################################
### chunk number 2: FigureExample
###################################################
ms.fxn <- function(x){
samp <- rnorm(x)
mean <- mean(samp)
mean
}

ms <- rep(NA,1000)
for(i in 1:1000){
ms[i] <- ms.fxn(10)
}

hist(ms)