It would be good to include an introduction to vectorization, especially with respect to R's random number generator functions. Right now we introduce it tangentially in discussing rnorm() but might be better to have a section of its own.
I have material to use.