sd <- function(x, na.rm=FALSE) sqrt(var(x, na.rm=na.rm))
