Write an R function moment (k,x, centered=(k>1)) that takes a numeric vector and returns the kth sample moment (either ˆμk or ˆμ k depending on the value of centered). Test your function by comparing its output to the output of mean() and var (). How does your function handle missing data (NA values)?