# create the standard error function
se <- function(x, ...) sd(x)/sqrt(length(x))