R/weighted_mean.R

In tobiasschoch/robsurvey: Robust Survey Statistics Estimation

# weighted mean (Hajek estimator)
weighted_mean <- function(x, w, na.rm = FALSE)
{
    dat <- .check_data_weights(x, w, na.rm)
    if (is.null(dat))
        NA
    else
        sum(dat$x * dat$w) / sum(dat$w)
}
# weighted total
weighted_total <- function(x, w, na.rm = FALSE)
{
    dat <- .check_data_weights(x, w, na.rm)
    if (is.null(dat))
        NA
    else
        sum(dat$w * dat$x)
}