colSd <- function(x, na.rm = TRUE) {
if (na.rm) {
n <- colSums(!is.na(x))
} else {
n <- nrow(x)
}
colVar <- colMeans(x*x, na.rm = na.rm) - (colMeans(x, na.rm = na.rm))^2
return(sqrt(colVar * n/(n-1)))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.