\begin{Schunk}
\begin{Sinput}
  par(mfcol=c(2,2)) 
 hist(dat$GRIP, prob = TRUE, xlab = "Grip strength", main = "Histogram: Grip Strength")
 curve(dnorm(x, m = sumdat["mean","GRIP"], s = sumdat["sd", "GRIP"]), from =  range(dat$GRIP)[1], to = range(dat$GRIP)[2], add = TRUE) 
 hist(dat$ARM, prob = TRUE, xlab = "Arm strength", main = "Histogram: Arm Strength") 
 curve(dnorm(x, m = sumdat["mean","ARM"], s = sumdat["sd", "ARM"]), from =  range(dat$ARM)[1], to = range(dat$ARM)[2], add = TRUE) 
 hist(dat$RATINGS, prob = TRUE, xlab = "Supervisor Rating", main = "Histogram: Rating") 
 curve(dnorm(x, m = sumdat["mean","RATINGS"], s = sumdat["sd", "RATINGS"]), from = range(dat$RATINGS)[1], to = range(dat$RATINGS)[2], add = TRUE) 
 hist(dat$SIMS, prob = TRUE, xlab="Sims Wrench Test", main="Histogram: Strength Simulation") 
 curve(dnorm(x, m = sumdat["mean","SIMS"], s = sumdat["sd", "SIMS"]), from =  range(dat$SIMS)[1], to = range(dat$SIMS)[2], add = TRUE)
\end{Sinput}
\end{Schunk}
