g_mean<- function(x) {
  n = sum(!is.na(x))
  res <- prod(x, na.rm = T)^(1/n)
  res
}