sem <- function(x, na.rm=F){
  sem<-sd(x, na.rm = na.rm)/sqrt(length(x))
  return(sem)
}