Generalized/power mean with R

powermean.r

# Power mean
## Examples:
## p = 1, arithmetic mean
## p = 2, quadratic mean (RMS)
pmean <- function(x, p) {
  library(TTR)
  return(tail(TTR::SMA(x^p, length(x))^(1/p), n=1))
}