1SN <- function(x){
2    20*log10(mean(x)/sd(x))
3}