Note, uses denominator n-1

gm_sd(x, na.rm = FALSE)