Consider the following data frame:
dat <- data.frame( ID = c(1:200), var1 = rnorm(200), var2 = rnorm(200), var3 =