R/probability.R

In neckl004/NecklenJacobTools: Jacob Necklen's Tools

#' Calculate weighted mean, variance and standard deviation
           #'
           #' @param d data matrix
           #'
           #' @return list
           #'
           #' @export
           #'
           #' @examples
           #' probs=runif(10)
           #' probability(list(x = runif(10), p=(probs/sum(probs))))
           #'
           probability = function(d)
           {
              mu = sum(d$x * d$p)
              vari = sum(((d$x-mu)^2)*d$p)
              sd = sqrt(vari)
              list(mu = mu,vari = vari,sd = sd)
           }