Me too. Why for example did they use an sapply for column means when they could have just used colMeans with na.rm=T?